From johnrahill at college.harvard.edu Sat Nov 4 14:35:43 2017 From: johnrahill at college.harvard.edu (John Rahill) Date: Sat, 4 Nov 2017 17:35:43 -0400 Subject: [CIG-SEISMO] specfem question Message-ID: Hello, My name is John Rahill and I am a senior at Harvard University. I am using SPECFEM2D/3D for my senior thesis project to model high frequency wave sin elastic media in the near subsurface. I have been running various simulations with different parameters thus far, but I am having trouble executing the ./bin/xspecfem2D command. Currently in the bin folder there is only xmeshfem2D and xmeshfem2D.dSYM. I am wondering whether I set up the program fully / correctly, or if I am just trying to execute the command incorrectly or from the wrong directory. Thank you very much for the help, John -------------- next part -------------- An HTML attachment was scrubbed... URL: From jas11 at Princeton.EDU Sun Nov 5 12:31:53 2017 From: jas11 at Princeton.EDU (James A. Smith) Date: Sun, 5 Nov 2017 20:31:53 +0000 Subject: [CIG-SEISMO] specfem question In-Reply-To: References: Message-ID: Hi John, What happens if you type 'make xspecfem2D' in your specfem2d directory? This should build the xspecfem2D executable in your bin folder. James ________________________________ From: CIG-SEISMO [cig-seismo-bounces at geodynamics.org] on behalf of John Rahill [johnrahill at college.harvard.edu] Sent: Saturday, November 04, 2017 4:35 PM To: cig-seismo at geodynamics.org Subject: [CIG-SEISMO] specfem question Hello, My name is John Rahill and I am a senior at Harvard University. I am using SPECFEM2D/3D for my senior thesis project to model high frequency wave sin elastic media in the near subsurface. I have been running various simulations with different parameters thus far, but I am having trouble executing the ./bin/xspecfem2D command. Currently in the bin folder there is only xmeshfem2D and xmeshfem2D.dSYM. I am wondering whether I set up the program fully / correctly, or if I am just trying to execute the command incorrectly or from the wrong directory. Thank you very much for the help, John -------------- next part -------------- An HTML attachment was scrubbed... URL: From giuseppe.digiulio at ingv.it Mon Nov 6 04:26:32 2017 From: giuseppe.digiulio at ingv.it (Giuseppe Di Giulio) Date: Mon, 6 Nov 2017 13:26:32 +0100 Subject: [CIG-SEISMO] specfem2d receivers on surface with an external mesh Message-ID: ​​​​​ profilo_geo.msh.gz ​​​Dear SpecFem developers I am using Specfem2D 7.0 with an external mesh (built by gmsh) aimed at investigating the topographic effect. My will is to fix the position of the receivers on the topographyc surface. My understanding is that the parameter "enreg_surf_same_vertical" in the Par_file cannot be "true" with an external mesh. I try to set multiple receiver lines in the Par file where (xdeb,zdeb) are the coordinates of my receivers. The receiver lines (nreceiversets parameter) are fixed to the same number of receivers, and nrec is equal to 1 for each receiver line (and for this I do not care of xfin and zfin; see the end of message). Unfortunately it does not work, the error message after the run of xmeshfem2D is "At line 326 of file ../shared/read_value_parameters.f90 Fortran runtime error: Bad real number in item 1 of list input" Any suggestion to place receivers on topographic surface with an external 2D mesh is welcome. The Par_file and my mesh are available as attach Thanks Giuseppe # first receiver line nrec = 1 # number of receivers xdeb = 10. # first receiver x in meters zdeb = 880. # first receiver z in meters xfin = 2500. # last receiver x in meters (ignored if onlyone receiver) zfin = 870.000427 # last receiver z in meters (ignored if onlyone receiver) enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface # second receiver line nrec = 1 # number of receivers xdeb = 1000. # first receiver x in meters zdeb = 1000. # first receiver z in meters xfin = 2500. # last receiver x in meters (ignored if onlyone receiver) zfin = 870.000427 # last receiver z in meters (ignored if onlyone receiver) enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface -------------- next part -------------- An HTML attachment was scrubbed... URL: From giuseppe.digiulio at ingv.it Mon Nov 6 04:29:05 2017 From: giuseppe.digiulio at ingv.it (Giuseppe Di Giulio) Date: Mon, 6 Nov 2017 13:29:05 +0100 Subject: [CIG-SEISMO] specfem2d receivers on surface with an external mesh In-Reply-To: References: Message-ID: Sorry I probably forgot the Par_file in my previous message Giuseppe 2017-11-06 13:26 GMT+01:00 Giuseppe Di Giulio : > ​​​​​ > profilo_geo.msh.gz > > ​​​Dear SpecFem developers > > I am using Specfem2D 7.0 with an external mesh (built by gmsh) aimed at > investigating the topographic effect. > > My will is to fix the position of the receivers on the topographyc surface. > > My understanding is that the parameter "enreg_surf_same_vertical" in the > Par_file cannot be "true" with an external mesh. > > I try to set multiple receiver lines in the Par file where (xdeb,zdeb) are > the coordinates of my receivers. The receiver lines (nreceiversets > parameter) are fixed to the same number of receivers, and nrec is equal to > 1 for each receiver line (and for this I do not care of xfin and zfin; see > the end of message). > > Unfortunately it does not work, the error message after the run of > xmeshfem2D is > > "At line 326 of file ../shared/read_value_parameters.f90 > Fortran runtime error: Bad real number in item 1 of list input" > > Any suggestion to place receivers on topographic surface with an external > 2D mesh is welcome. > > The Par_file and my mesh are available as attach > > Thanks > > Giuseppe > > > # first receiver line > nrec = 1 # number of receivers > xdeb = 10. # first receiver x in > meters > zdeb = 880. # first receiver z in > meters > xfin = 2500. # last receiver x in > meters (ignored if onlyone receiver) > zfin = 870.000427 # last receiver z > in meters (ignored if onlyone receiver) > enreg_surf_same_vertical = .false. # receivers inside the > medium or at the surface > > # second receiver line > nrec = 1 # number of receivers > xdeb = 1000. # first receiver x in > meters > zdeb = 1000. # first receiver z in > meters > xfin = 2500. # last receiver x in > meters (ignored if onlyone receiver) > zfin = 870.000427 # last receiver z > in meters (ignored if onlyone receiver) > enreg_surf_same_vertical = .false. # receivers inside the > medium or at the surface > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Par_file Type: application/octet-stream Size: 13404 bytes Desc: not available URL: From komatitsch at lma.cnrs-mrs.fr Tue Nov 7 06:59:03 2017 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Tue, 7 Nov 2017 15:59:03 +0100 Subject: [CIG-SEISMO] specfem2d receivers on surface with an external mesh In-Reply-To: References: Message-ID: <27ed9d4c-ea07-389a-d5ba-f708d3009da0@lma.cnrs-mrs.fr> Dear Giuseppe, dear all, You are right, option "enreg_surf_same_vertical" in the Par_file currently does not work for external meshes because we have no easy way of knowing where the topography is in an external mesh (we would need to have a list of all element edges that are along the upper surface; could be done, but not done yet; not too difficult). In the meantime the easiest thing to do is probably to zoom on your topography in Gmesh (zooming a lot) and then click on the point where you want your receiver to be located, and write down the coordinates of that point (and then enter that manually in the Par_file). Best regards, Dimitri. On 11/06/2017 01:26 PM, Giuseppe Di Giulio wrote: > ​​​​​ > profilo_geo.msh.gz > > ​​​Dear SpecFem developers > > I am using Specfem2D 7.0 with an external mesh (built by gmsh) aimed at > investigating the topographic effect. > > My will is to fix the position of the receivers on the topographyc surface. > > My understanding is that the parameter "enreg_surf_same_vertical" in the > Par_file cannot be "true" with an external mesh. > > I try to set multiple receiver lines in the Par file where (xdeb,zdeb) > are the coordinates of my receivers. The receiver lines (nreceiversets > parameter) are fixed to the same number of receivers, and  nrec is equal > to 1 for each receiver line (and for this I do not care of xfin and > zfin; see the end of message). > > Unfortunately it does not work, the error message after the run of > xmeshfem2D is > > "At line 326 of file ../shared/read_value_parameters.f90 > Fortran runtime error: Bad real number in item 1 of list input" > > Any suggestion to place receivers on topographic surface with an > external  2D mesh is welcome. > > The Par_file and my mesh are available as attach > > Thanks > > Giuseppe > > > # first receiver line > nrec                            = 1             # number of receivers > xdeb                            = 10.           # first receiver x in meters > zdeb                            = 880.          # first receiver z in meters > xfin                            = 2500.          # last receiver x in > meters (ignored if onlyone receiver) > zfin                            = 870.000427          # last receiver z > in meters (ignored if onlyone receiver) > enreg_surf_same_vertical        = .false.         # receivers inside the > medium or at the surface > > # second receiver line > nrec                            = 1             # number of receivers > xdeb                            = 1000.           # first receiver x in > meters > zdeb                            = 1000.          # first receiver z in > meters > xfin                            = 2500.          # last receiver x in > meters (ignored if onlyone receiver) > zfin                            = 870.000427          # last receiver z > in meters (ignored if onlyone receiver) > enreg_surf_same_vertical        = .false.         # receivers inside the > medium or at the surface > > > > > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > -- Dimitri Komatitsch, CNRS Research Director (DR CNRS) Laboratory of Mechanics and Acoustics, Marseille, France http://komatitsch.free.fr From komatitsch at lma.cnrs-mrs.fr Tue Nov 7 07:33:33 2017 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Tue, 7 Nov 2017 16:33:33 +0100 Subject: [CIG-SEISMO] specfem2d receivers on surface with an external mesh / PS In-Reply-To: <27ed9d4c-ea07-389a-d5ba-f708d3009da0@lma.cnrs-mrs.fr> References: <27ed9d4c-ea07-389a-d5ba-f708d3009da0@lma.cnrs-mrs.fr> Message-ID: Dear Giuseppe, Using a single receiver per set, as in your: # first receiver line nrec = 1 # number of receivers xdeb = 10. # first receiver x in meters zdeb = 880. # first receiver z in meters xfin = 2500. # last receiver x in meters (ignored if onlyone receiver) zfin = 870.000427 # last receiver z in meters (ignored if onlyone receiver) enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface # second receiver line nrec = 1 # number of receivers xdeb = 1000. # first receiver x in meters zdeb = 1000. # first receiver z in meters xfin = 2500. # last receiver x in meters (ignored if onlyone receiver) zfin = 870.000427 # last receiver z in meters (ignored if onlyone receiver) enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface works fine here though. Are you using a recent version of the code? Cheers, Dimitri. On 11/07/2017 03:59 PM, Dimitri Komatitsch wrote: > > Dear Giuseppe, dear all, > > You are right, option "enreg_surf_same_vertical" in the Par_file > currently does not work for external meshes because we have no easy way > of knowing where the topography is in an external mesh (we would need to > have a list of all element edges that are along the upper surface; could > be done, but not done yet; not too difficult). > > In the meantime the easiest thing to do is probably to zoom on your > topography in Gmesh (zooming a lot) and then click on the point where > you want your receiver to be located, and write down the coordinates of > that point (and then enter that manually in the Par_file). > > Best regards, > Dimitri. > > On 11/06/2017 01:26 PM, Giuseppe Di Giulio wrote: >> ​​​​​ >> profilo_geo.msh.gz >> >> >> ​​​Dear SpecFem developers >> >> I am using Specfem2D 7.0 with an external mesh (built by gmsh) aimed >> at investigating the topographic effect. >> >> My will is to fix the position of the receivers on the topographyc >> surface. >> >> My understanding is that the parameter "enreg_surf_same_vertical" in >> the Par_file cannot be "true" with an external mesh. >> >> I try to set multiple receiver lines in the Par file where (xdeb,zdeb) >> are the coordinates of my receivers. The receiver lines (nreceiversets >> parameter) are fixed to the same number of receivers, and  nrec is >> equal to 1 for each receiver line (and for this I do not care of xfin >> and zfin; see the end of message). >> >> Unfortunately it does not work, the error message after the run of >> xmeshfem2D is >> >> "At line 326 of file ../shared/read_value_parameters.f90 >> Fortran runtime error: Bad real number in item 1 of list input" >> >> Any suggestion to place receivers on topographic surface with an >> external  2D mesh is welcome. >> >> The Par_file and my mesh are available as attach >> >> Thanks >> >> Giuseppe >> >> >> # first receiver line >> nrec                            = 1             # number of receivers >> xdeb                            = 10.           # first receiver x in >> meters >> zdeb                            = 880.          # first receiver z in >> meters >> xfin                            = 2500.          # last receiver x in >> meters (ignored if onlyone receiver) >> zfin                            = 870.000427          # last receiver >> z in meters (ignored if onlyone receiver) >> enreg_surf_same_vertical        = .false.         # receivers inside >> the medium or at the surface >> >> # second receiver line >> nrec                            = 1             # number of receivers >> xdeb                            = 1000.           # first receiver x >> in meters >> zdeb                            = 1000.          # first receiver z in >> meters >> xfin                            = 2500.          # last receiver x in >> meters (ignored if onlyone receiver) >> zfin                            = 870.000427          # last receiver >> z in meters (ignored if onlyone receiver) >> enreg_surf_same_vertical        = .false.         # receivers inside >> the medium or at the surface >> >> >> >> >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >> > -- Dimitri Komatitsch, CNRS Research Director (DR CNRS) Laboratory of Mechanics and Acoustics, Marseille, France http://komatitsch.free.fr From giuseppe.digiulio at ingv.it Tue Nov 7 09:36:20 2017 From: giuseppe.digiulio at ingv.it (Giuseppe Di Giulio) Date: Tue, 7 Nov 2017 18:36:20 +0100 Subject: [CIG-SEISMO] specfem2d receivers on surface with an external mesh / PS In-Reply-To: References: <27ed9d4c-ea07-389a-d5ba-f708d3009da0@lma.cnrs-mrs.fr> Message-ID: Hi Dimitri I am using the release 7.0. I will try with the most recent one. Thanks Giuseppe 2017-11-07 16:33 GMT+01:00 Dimitri Komatitsch : > > Dear Giuseppe, > > Using a single receiver per set, as in your: > > # first receiver line > nrec = 1 # number of receivers > xdeb = 10. # first receiver x in > meters > zdeb = 880. # first receiver z in > meters > xfin = 2500. # last receiver x in > meters (ignored if onlyone receiver) > zfin = 870.000427 # last receiver z in > meters (ignored if onlyone receiver) > enreg_surf_same_vertical = .false. # receivers inside the > medium or at the surface > > # second receiver line > nrec = 1 # number of receivers > xdeb = 1000. # first receiver x in > meters > zdeb = 1000. # first receiver z in > meters > xfin = 2500. # last receiver x in > meters (ignored if onlyone receiver) > zfin = 870.000427 # last receiver z in > meters (ignored if onlyone receiver) > enreg_surf_same_vertical = .false. # receivers inside the > medium or at the surface > > works fine here though. Are you using a recent version of the code? > > Cheers, > Dimitri. > > On 11/07/2017 03:59 PM, Dimitri Komatitsch wrote: > >> >> Dear Giuseppe, dear all, >> >> You are right, option "enreg_surf_same_vertical" in the Par_file >> currently does not work for external meshes because we have no easy way of >> knowing where the topography is in an external mesh (we would need to have >> a list of all element edges that are along the upper surface; could be >> done, but not done yet; not too difficult). >> >> In the meantime the easiest thing to do is probably to zoom on your >> topography in Gmesh (zooming a lot) and then click on the point where you >> want your receiver to be located, and write down the coordinates of that >> point (and then enter that manually in the Par_file). >> >> Best regards, >> Dimitri. >> >> On 11/06/2017 01:26 PM, Giuseppe Di Giulio wrote: >> >>> ​​​​​ >>> profilo_geo.msh.gz >> gv.it/file/d/19zHUn9T6deq0q1KeVzeJs0jqAw3uFkaR/view?usp=drive_web> >>> ​​​Dear SpecFem developers >>> >>> I am using Specfem2D 7.0 with an external mesh (built by gmsh) aimed at >>> investigating the topographic effect. >>> >>> My will is to fix the position of the receivers on the topographyc >>> surface. >>> >>> My understanding is that the parameter "enreg_surf_same_vertical" in the >>> Par_file cannot be "true" with an external mesh. >>> >>> I try to set multiple receiver lines in the Par file where (xdeb,zdeb) >>> are the coordinates of my receivers. The receiver lines (nreceiversets >>> parameter) are fixed to the same number of receivers, and nrec is equal to >>> 1 for each receiver line (and for this I do not care of xfin and zfin; see >>> the end of message). >>> >>> Unfortunately it does not work, the error message after the run of >>> xmeshfem2D is >>> >>> "At line 326 of file ../shared/read_value_parameters.f90 >>> Fortran runtime error: Bad real number in item 1 of list input" >>> >>> Any suggestion to place receivers on topographic surface with an >>> external 2D mesh is welcome. >>> >>> The Par_file and my mesh are available as attach >>> >>> Thanks >>> >>> Giuseppe >>> >>> >>> # first receiver line >>> nrec = 1 # number of receivers >>> xdeb = 10. # first receiver x in >>> meters >>> zdeb = 880. # first receiver z in >>> meters >>> xfin = 2500. # last receiver x in >>> meters (ignored if onlyone receiver) >>> zfin = 870.000427 # last receiver z >>> in meters (ignored if onlyone receiver) >>> enreg_surf_same_vertical = .false. # receivers inside >>> the medium or at the surface >>> >>> # second receiver line >>> nrec = 1 # number of receivers >>> xdeb = 1000. # first receiver x in >>> meters >>> zdeb = 1000. # first receiver z in >>> meters >>> xfin = 2500. # last receiver x in >>> meters (ignored if onlyone receiver) >>> zfin = 870.000427 # last receiver z >>> in meters (ignored if onlyone receiver) >>> enreg_surf_same_vertical = .false. # receivers inside >>> the medium or at the surface >>> >>> >>> >>> >>> >>> _______________________________________________ >>> CIG-SEISMO mailing list >>> CIG-SEISMO at geodynamics.org >>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >>> >>> >> > -- > Dimitri Komatitsch, CNRS Research Director (DR CNRS) > Laboratory of Mechanics and Acoustics, Marseille, France > http://komatitsch.free.fr > -------------- next part -------------- An HTML attachment was scrubbed... URL: From emanuel.kaestle at upmc.fr Mon Nov 13 07:39:31 2017 From: emanuel.kaestle at upmc.fr (Emanuel Kaestle) Date: Mon, 13 Nov 2017 16:39:31 +0100 Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) Message-ID: <1f76121a-b588-b5ff-bbb9-15b0e5a087b5@upmc.fr> Hi there, I am having some touble compiling SPECFEM3D. Execution of the configure script returns me the following error /## --- ## ## MPI ## ## --- ## checking mpi.h usability... no checking mpi.h presence... no checking for mpi.h... no configure: error: MPI header not found; try setting MPI_INC./ I start it with the command /./configure FC=gfortran MPIFC=mpif90/ Changing it to /./configure FC=gfortran MPIFC=mpif90 MPI_INC='/usr/include/mpich/mpi.h'/ doesn't change anything. It returns the same error. The mpi.h header should be present, at least I can find it with the /find/ command: //usr/src/linux-headers-4.8.0-53/include/linux/mpi.h// ///usr/src/linux-headers-4.10.0-28/include/linux/mpi.h// ///usr/include/mpich/mpi.h/ The small script from the manual (program main / include 'mpif.h' / ....) runs without problems if I compile it with mpif90 (not with gfortran). What is the correct way to set MPI_INC? Thanks for your help! Emanuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From rama.kishan.v.malladi at intel.com Mon Nov 13 07:50:09 2017 From: rama.kishan.v.malladi at intel.com (Malladi, Rama Kishan V) Date: Mon, 13 Nov 2017 15:50:09 +0000 Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) In-Reply-To: <1f76121a-b588-b5ff-bbb9-15b0e5a087b5@upmc.fr> References: <1f76121a-b588-b5ff-bbb9-15b0e5a087b5@upmc.fr> Message-ID: Typically MPI_INC is the path to the location where the header “mpi.h” is available. So, try setting MPI_INC to /usr/include/mpich From: CIG-SEISMO [mailto:cig-seismo-bounces at geodynamics.org] On Behalf Of Emanuel Kaestle Sent: Monday, November 13, 2017 9:10 PM To: cig-seismo at geodynamics.org Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) Hi there, I am having some touble compiling SPECFEM3D. Execution of the configure script returns me the following error ## --- ## ## MPI ## ## --- ## checking mpi.h usability... no checking mpi.h presence... no checking for mpi.h... no configure: error: MPI header not found; try setting MPI_INC. I start it with the command ./configure FC=gfortran MPIFC=mpif90 Changing it to ./configure FC=gfortran MPIFC=mpif90 MPI_INC='/usr/include/mpich/mpi.h' doesn't change anything. It returns the same error. The mpi.h header should be present, at least I can find it with the find command: /usr/src/linux-headers-4.8.0-53/include/linux/mpi.h /usr/src/linux-headers-4.10.0-28/include/linux/mpi.h /usr/include/mpich/mpi.h The small script from the manual (program main / include 'mpif.h' / ....) runs without problems if I compile it with mpif90 (not with gfortran). What is the correct way to set MPI_INC? Thanks for your help! Emanuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From emanuel.kaestle at upmc.fr Mon Nov 13 08:21:00 2017 From: emanuel.kaestle at upmc.fr (Emanuel Kaestle) Date: Mon, 13 Nov 2017 17:21:00 +0100 Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) In-Reply-To: References: <1f76121a-b588-b5ff-bbb9-15b0e5a087b5@upmc.fr> Message-ID: <1095a4b3-f347-e03e-0fa4-d6685400949f@upmc.fr> Thanks for your quick reply! But why does the following command not do the trick (independent whether I use inverted commas or not)? /./////configure FC=gfortran MPIFC=mpif90 MPI_INC='/usr/include/mpich/mpi.h//'////// /The configure script recognizes that MPI_INC is set. Before, I got the error: /gfortran:error: unrecognized command line option '-showme:incdir/s' Now I just get the error /MPI header not found; .../ Thanks again for any suggestions, Emanuel / / On 11/13/2017 04:50 PM, Malladi, Rama Kishan V wrote: > > /Typically//MPI_INC is the path to the location where the header > “mpi.h” is available. So, try setting MPI_INC to /usr/include/mpich/ > > *From:*CIG-SEISMO [mailto:cig-seismo-bounces at geodynamics.org] *On > Behalf Of *Emanuel Kaestle > *Sent:* Monday, November 13, 2017 9:10 PM > *To:* cig-seismo at geodynamics.org > *Subject:* [CIG-SEISMO] compilation error (configure: error: MPI > header not found; try setting MPI_INC) > > Hi there, > > I am having some touble compiling SPECFEM3D. Execution of the > configure script returns me the following error > > /## --- ## > ## MPI ## > ## --- ## > checking mpi.h usability... no > checking mpi.h presence... no > checking for mpi.h... no > configure: error: MPI header not found; try setting MPI_INC./ > > I start it with the command /./configure FC=gfortran MPIFC=mpif90/ > > Changing it to /./configure FC=gfortran MPIFC=mpif90 > MPI_INC='/usr/include/mpich/mpi.h'/ doesn't change anything. It > returns the same error. > > The mpi.h header should be present, at least I can find it with the > /find/ command: > > //usr/src/linux-headers-4.8.0-53/include/linux/mpi.h > /usr/src/linux-headers-4.10.0-28/include/linux/mpi.h > /usr/include/mpich/mpi.h/ > > The small script from the manual (program main / include 'mpif.h' / > ....) runs without problems if I compile it with mpif90 (not with > gfortran). > > What is the correct way to set MPI_INC? > > Thanks for your help! > > Emanuel > > > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo -------------- next part -------------- An HTML attachment was scrubbed... URL: From emanuel.kaestle at upmc.fr Mon Nov 13 08:49:14 2017 From: emanuel.kaestle at upmc.fr (Emanuel Kaestle) Date: Mon, 13 Nov 2017 17:49:14 +0100 Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) In-Reply-To: <1095a4b3-f347-e03e-0fa4-d6685400949f@upmc.fr> References: <1f76121a-b588-b5ff-bbb9-15b0e5a087b5@upmc.fr> <1095a4b3-f347-e03e-0fa4-d6685400949f@upmc.fr> Message-ID: <809f72ab-6745-c274-609e-5908196711ab@upmc.fr> Now it works, I was not reading your answer carefully enough, Rama Kishan. Thank you! On 11/13/2017 05:21 PM, Emanuel Kaestle wrote: > > Thanks for your quick reply! > > But why does the following command not do the trick (independent > whether I use inverted commas or not)? > /./////configure FC=gfortran MPIFC=mpif90 > MPI_INC='/usr/include/mpich/mpi.h//'////// > > /The configure script recognizes that MPI_INC is set. Before, I got > the error: > /gfortran:error: unrecognized command line option '-showme:incdir/s' > > Now I just get the error /MPI header not found; .../ > > Thanks again for any suggestions, > > Emanuel > > / > / > > > On 11/13/2017 04:50 PM, Malladi, Rama Kishan V wrote: >> >> /Typically//MPI_INC is the path to the location where the header >> “mpi.h” is available. So, try setting MPI_INC to /usr/include/mpich/ >> >> *From:*CIG-SEISMO [mailto:cig-seismo-bounces at geodynamics.org] *On >> Behalf Of *Emanuel Kaestle >> *Sent:* Monday, November 13, 2017 9:10 PM >> *To:* cig-seismo at geodynamics.org >> *Subject:* [CIG-SEISMO] compilation error (configure: error: MPI >> header not found; try setting MPI_INC) >> >> Hi there, >> >> I am having some touble compiling SPECFEM3D. Execution of the >> configure script returns me the following error >> >> /## --- ## >> ## MPI ## >> ## --- ## >> checking mpi.h usability... no >> checking mpi.h presence... no >> checking for mpi.h... no >> configure: error: MPI header not found; try setting MPI_INC./ >> >> I start it with the command /./configure FC=gfortran MPIFC=mpif90/ >> >> Changing it to /./configure FC=gfortran MPIFC=mpif90 >> MPI_INC='/usr/include/mpich/mpi.h'/ doesn't change anything. It >> returns the same error. >> >> The mpi.h header should be present, at least I can find it with the >> /find/ command: >> >> //usr/src/linux-headers-4.8.0-53/include/linux/mpi.h >> /usr/src/linux-headers-4.10.0-28/include/linux/mpi.h >> /usr/include/mpich/mpi.h/ >> >> The small script from the manual (program main / include 'mpif.h' / >> ....) runs without problems if I compile it with mpif90 (not with >> gfortran). >> >> What is the correct way to set MPI_INC? >> >> Thanks for your help! >> >> Emanuel >> >> >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > > > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo -------------- next part -------------- An HTML attachment was scrubbed... URL: From J.Wookey at bristol.ac.uk Mon Nov 13 08:27:26 2017 From: J.Wookey at bristol.ac.uk (James Wookey) Date: Mon, 13 Nov 2017 16:27:26 +0000 Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) In-Reply-To: <1095a4b3-f347-e03e-0fa4-d6685400949f@upmc.fr> References: <1f76121a-b588-b5ff-bbb9-15b0e5a087b5@upmc.fr> <1095a4b3-f347-e03e-0fa4-d6685400949f@upmc.fr> Message-ID: <95A14355-FCEE-4CC0-9C07-8036777C4D2D@bristol.ac.uk> It looks like you are pointing at the file specifically, not the directory which contains it. I think MPI_INC should be the directory containing mpi.h, i.e.: ./configure FC=gfortran MPIFC=mpif90 MPI_INC=/usr/include/mpich > On 13 Nov 2017, at 16:21, Emanuel Kaestle wrote: > > Thanks for your quick reply! > > But why does the following command not do the trick (independent whether I use inverted commas or not)? > ./configure FC=gfortran MPIFC=mpif90 MPI_INC='/usr/include/mpich/mpi.h' > > The configure script recognizes that MPI_INC is set. Before, I got the error: > gfortran:error: unrecognized command line option '-showme:incdirs' > > Now I just get the error MPI header not found; ... > Thanks again for any suggestions, > > Emanuel > > > > On 11/13/2017 04:50 PM, Malladi, Rama Kishan V wrote: >> Typically MPI_INC is the path to the location where the header “mpi.h” is available. So, try setting MPI_INC to /usr/include/mpich >> >> >> >> From: CIG-SEISMO [mailto:cig-seismo-bounces at geodynamics.org] On Behalf Of Emanuel Kaestle >> Sent: Monday, November 13, 2017 9:10 PM >> To: cig-seismo at geodynamics.org >> Subject: [CIG-SEISMO] compilation error (configure: error: MPI header not found; try setting MPI_INC) >> >> >> >> Hi there, >> >> I am having some touble compiling SPECFEM3D. Execution of the configure script returns me the following error >> >> ## --- ## >> ## MPI ## >> ## --- ## >> checking mpi.h usability... no >> checking mpi.h presence... no >> checking for mpi.h... no >> configure: error: MPI header not found; try setting MPI_INC. >> >> I start it with the command ./configure FC=gfortran MPIFC=mpif90 >> >> Changing it to ./configure FC=gfortran MPIFC=mpif90 MPI_INC='/usr/include/mpich/mpi.h' doesn't change anything. It returns the same error. >> >> The mpi.h header should be present, at least I can find it with the find command: >> >> /usr/src/linux-headers-4.8.0-53/include/linux/mpi.h >> /usr/src/linux-headers-4.10.0-28/include/linux/mpi.h >> /usr/include/mpich/mpi.h >> >> The small script from the manual (program main / include 'mpif.h' / ....) runs without problems if I compile it with mpif90 (not with gfortran). >> >> What is the correct way to set MPI_INC? >> >> Thanks for your help! >> >> Emanuel >> >> >> >> _______________________________________________ >> CIG-SEISMO mailing list >> >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo From tjesser at ucdavis.edu Mon Nov 13 17:21:16 2017 From: tjesser at ucdavis.edu (Tyler Esser) Date: Mon, 13 Nov 2017 17:21:16 -0800 Subject: [CIG-SEISMO] SW4 2.0 Released Message-ID: We are pleased to announce the release of SW4 2.0. SW4 (Seismic Waves, 4th order) implements substantial capabilities for 3-D seismic modeling, with a free surface condition on the top boundary, absorbing super-grid conditions on the far-field boundaries, and an arbitrary number of point force and/or point moment tensor source terms. SW4 is available at https://geodynamics.org/cig/software/sw4/ and the source repository is available at https://github.com/geodynamics/sw4. Version 2.0 of SW4 implements mesh refinement with hanging nodes. Mesh refinement is currently supported in the Cartesian portion of the mesh, but can be used together with realistic topography and heterogeneous isotropic visco-elastic material models. This release also marks the completion of the migration of SW4Lite into its own repository at https://github.com/geodynamics/sw4lite. Thanks to N. Anders Petersson and Bjorn Sjogreen! -------------- next part -------------- An HTML attachment was scrubbed... URL: From jfontiela at uevora.pt Wed Nov 15 03:49:42 2017 From: jfontiela at uevora.pt (=?utf-8?Q?Jo=C3=A3o_Fontiela?=) Date: Wed, 15 Nov 2017 11:49:42 +0000 Subject: [CIG-SEISMO] Unsuccessuful SW4 test Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47@uevora.pt> Dear all, We installed SW4 (latest release) successfully nevertheless when we execute test_sw4.py (we have installed also python3) we got the following error: [root at terra pytest]# ./test_sw4.py Running all tests for level 0 ... Starting test # 1 in directory: meshrefine with input file: refine-el-1.in ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 32512 aborting test run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in >& refine-el-1.out DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) test_sw4 was unsuccessful Make.inc configuration is: FC = gfortran CXX = mpicxx EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran Can anyone help understand what is going on? Yours, João Fontiela ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela -------------- next part -------------- An HTML attachment was scrubbed... URL: From petersson1 at llnl.gov Wed Nov 15 16:33:16 2017 From: petersson1 at llnl.gov (Petersson, Anders) Date: Thu, 16 Nov 2017 00:33:16 +0000 Subject: [CIG-SEISMO] =?utf-8?q?Unsuccessuful_SW4_test_=28Jo=C3=A3o_Fontie?= =?utf-8?b?bGEp?= Message-ID: Hi Joao, Sorry to hear about your difficulties testing SW4. You have some strange characters in the output from the test script. The first line is: Starting test # 1 in directory: meshrefine with input file: refine-el-1.in The hyperlink should not be there (). Can you can run the test manually from the pytest directory? > cd pytest > mpirun –np 1 ../optimize/sw4 reference/meshrefine/refine-el-1.in Good luck, Anders On 11/15/17, 12:05 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org" wrote: Send CIG-SEISMO mailing list submissions to cig-seismo at geodynamics.org To subscribe or unsubscribe via the World Wide Web, visit http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo or, via email, send a message with subject or body 'help' to cig-seismo-request at geodynamics.org You can reach the person managing the list at cig-seismo-owner at geodynamics.org When replying, please edit your Subject line so it is more specific than "Re: Contents of CIG-SEISMO digest..." Today's Topics: 1. Unsuccessuful SW4 test (João Fontiela) ---------------------------------------------------------------------- Message: 1 Date: Wed, 15 Nov 2017 11:49:42 +0000 From: João Fontiela To: cig-seismo at geodynamics.org Subject: [CIG-SEISMO] Unsuccessuful SW4 test Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47 at uevora.pt> Content-Type: text/plain; charset="utf-8" Dear all, We installed SW4 (latest release) successfully nevertheless when we execute test_sw4.py (we have installed also python3) we got the following error: [root at terra pytest]# ./test_sw4.py Running all tests for level 0 ... Starting test # 1 in directory: meshrefine with input file: refine-el-1.in ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 32512 aborting test run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in >& refine-el-1.out DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) test_sw4 was unsuccessful Make.inc configuration is: FC = gfortran CXX = mpicxx EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran Can anyone help understand what is going on? Yours, João Fontiela ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ Subject: Digest Footer _______________________________________________ CIG-SEISMO mailing list CIG-SEISMO at geodynamics.org http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo ------------------------------ End of CIG-SEISMO Digest, Vol 118, Issue 8 ****************************************** From tjesser at ucdavis.edu Wed Nov 15 17:33:26 2017 From: tjesser at ucdavis.edu (Tyler Esser) Date: Wed, 15 Nov 2017 17:33:26 -0800 Subject: [CIG-SEISMO] =?utf-8?q?Unsuccessuful_SW4_test_=28Jo=C3=A3o_Fontie?= =?utf-8?b?bGEp?= In-Reply-To: References: Message-ID: Anders, the hyperlink text was probably added by the mailing list. I'd guess that the filename, refine-el-1.in, was automatically treated as a URL. I expect my email will be similarly mangled. João, I think I was able to replicate your issue. It looks like the error was caused by the way the Python script calls mpirun. In pytest/test_sw4.py, lines 200 - 211 setup the mpirun command and handle the resulting error. In particular, line 206: status = os.system(run_cmd) > executes a string created on line 201: run_cmd = sw4_mpi_run + ' ' + sw4_input_file + ' >& ' + sw4_stdout_file > I think the culprit is the output redirect: ' >& ', as apparently that syntax is only supported in csh and tcsh shells. ( https://unix.stackexchange.com/a/176224) It would be helpful to have more information on your system environment so that we can figure out whether you are using one of those system shells, or if your problem comes from something else. However, the Python documentation recommends against using os.system() to execute commands, instead recommending the subprocess python module: ( https://docs.python.org/3/library/os.html#os.system). The subprocess module should allow output redirection in a system agnostic way, which I think will prevent your issue. Anders, if you agree with this assessment, I can write a patch using the subprocess module and submit it to the Github repo. Tyler Esser P.S. Here's the error I got when running on my system: tyler at Latitude:~/Downloads/sw4-v2.0/pytest$ ./test_sw4.py -v > pytest_dir = /home/tyler/Downloads/sw4-v2.0/pytest > sw4_base_dir = /home/tyler/Downloads/sw4-v2.0 > optimize_dir = /home/tyler/Downloads/sw4-v2.0/optimize > reference_dir = /home/tyler/Downloads/sw4-v2.0/pytest/reference os.uname= posix.uname_result(sysname='Linux', nodename='Latitude', > release='4.4.0-98-generic', version='#121-Ubuntu SMP Tue Oct 10 14:24:03 > UTC 2017', machine='x86_64') > node_name= Latitude > sys_name= Linux > mpirun_cmd = mpirun -np 1 > Running all tests for level 0 ... > Starting test # 1 in directory: meshrefine with input file: refine-el-1.in > sh: 1: Syntax error: Bad fd number > ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 512 > aborting test > run_cmd= mpirun -np 1 /home/tyler/Downloads/sw4-v2.0/optimize/sw4 > /home/tyler/Downloads/sw4-v2.0/pytest/reference/meshrefine/refine-el-1.in > >& refine-el-1.out > DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) > test_sw4 was unsuccessful On Wed, Nov 15, 2017 at 4:33 PM, Petersson, Anders wrote: > Hi Joao, > Sorry to hear about your difficulties testing SW4. > You have some strange characters in the output from the test script. The > first line is: > Starting test # 1 in directory: meshrefine with input file: refine-el-1.in > > The hyperlink should not be there (). > > Can you can run the test manually from the pytest directory? > > cd pytest > > mpirun –np 1 ../optimize/sw4 reference/meshrefine/refine-el-1.in > > Good luck, > Anders > > On 11/15/17, 12:05 PM, "CIG-SEISMO on behalf of cig-seismo-request@ > geodynamics.org" cig-seismo-request at geodynamics.org> wrote: > > Send CIG-SEISMO mailing list submissions to > cig-seismo at geodynamics.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > or, via email, send a message with subject or body 'help' to > cig-seismo-request at geodynamics.org > > You can reach the person managing the list at > cig-seismo-owner at geodynamics.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of CIG-SEISMO digest..." > > > Today's Topics: > > 1. Unsuccessuful SW4 test (João Fontiela) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Nov 2017 11:49:42 +0000 > From: João Fontiela > To: cig-seismo at geodynamics.org > Subject: [CIG-SEISMO] Unsuccessuful SW4 test > Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47 at uevora.pt> > Content-Type: text/plain; charset="utf-8" > > Dear all, > > > We installed SW4 (latest release) successfully nevertheless when we > execute test_sw4.py (we have installed also python3) we got the following > error: > > [root at terra pytest]# ./test_sw4.py > Running all tests for level 0 ... > Starting test # 1 in directory: meshrefine with input file: > refine-el-1.in > ERROR: Test refine-el-1.in : sw4 returned > non-zero exit status= 32512 aborting test > run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 > /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in < > http://refine-el-1.in/> >& refine-el-1.out > DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d > OPTION) > test_sw4 was unsuccessful > > Make.inc configuration is: > FC = gfortran > CXX = mpicxx > EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran > > Can anyone help understand what is going on? > > Yours, > João Fontiela > > ********************* > ********************* > João Fontiela > PhD Student > skype: fontiela > email: fontiela at gmail.com > mobile phone +351 917 862 830 > > Profissional contacts: > email: jfontiela at uevora.pt > Tel: (+351) 266 745 354 ext: 5395 or 5413 > Fax: (+351) 266 745 394 > > Profissional address: > Universidade de Évora > Instituto de Ciências da Terra > Colégio Luís Verney > Rua Romão Ramalho, 59 > 7000-671 Évora - Portugal > > ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: attachments/20171115/901b5718/attachment-0001.html> > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > > ------------------------------ > > End of CIG-SEISMO Digest, Vol 118, Issue 8 > ****************************************** > > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo -------------- next part -------------- An HTML attachment was scrubbed... URL: From petersson1 at llnl.gov Mon Nov 20 13:43:41 2017 From: petersson1 at llnl.gov (Petersson, Anders) Date: Mon, 20 Nov 2017 21:43:41 +0000 Subject: [CIG-SEISMO] =?utf-8?q?Unsuccessuful_SW4_test_=28Jo=C3=A3o_Fontie?= =?utf-8?b?bGEp?= Message-ID: Hi Joao, We have improved the test_sw4.py script and included it in a new release of SW4, with version 2.01. Note that only the test script was changed, and the source code is identical to version 2.0. The new version can be downloaded from the CIG website, or directly from github. It would be good to know whether the new script works for you. Thanks, Anders On 11/15/17, 12:00 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org" wrote: Send CIG-SEISMO mailing list submissions to cig-seismo at geodynamics.org To subscribe or unsubscribe via the World Wide Web, visit http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo or, via email, send a message with subject or body 'help' to cig-seismo-request at geodynamics.org You can reach the person managing the list at cig-seismo-owner at geodynamics.org When replying, please edit your Subject line so it is more specific than "Re: Contents of CIG-SEISMO digest..." Today's Topics: 1. Unsuccessuful SW4 test (João Fontiela) ---------------------------------------------------------------------- Message: 1 Date: Wed, 15 Nov 2017 11:49:42 +0000 From: João Fontiela To: cig-seismo at geodynamics.org Subject: [CIG-SEISMO] Unsuccessuful SW4 test Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47 at uevora.pt> Content-Type: text/plain; charset="utf-8" Dear all, We installed SW4 (latest release) successfully nevertheless when we execute test_sw4.py (we have installed also python3) we got the following error: [root at terra pytest]# ./test_sw4.py Running all tests for level 0 ... Starting test # 1 in directory: meshrefine with input file: refine-el-1.in ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 32512 aborting test run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in >& refine-el-1.out DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) test_sw4 was unsuccessful Make.inc configuration is: FC = gfortran CXX = mpicxx EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran Can anyone help understand what is going on? Yours, João Fontiela ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ Subject: Digest Footer _______________________________________________ CIG-SEISMO mailing list CIG-SEISMO at geodynamics.org http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo ------------------------------ End of CIG-SEISMO Digest, Vol 118, Issue 8 ****************************************** From jfontiela at uevora.pt Mon Nov 20 15:17:30 2017 From: jfontiela at uevora.pt (=?utf-8?Q?Jo=C3=A3o_Fontiela?=) Date: Mon, 20 Nov 2017 23:17:30 +0000 Subject: [CIG-SEISMO] =?utf-8?q?Unsuccessuful_SW4_test_=28Jo=C3=A3o_Fontie?= =?utf-8?b?bGEp?= In-Reply-To: References: Message-ID: <73966844-7F32-4904-AB80-6A9A788F4D0D@uevora.pt> Hi Anders, Thank you very much for your availability to fix this problem I will test the new release of test_sw4.py and let you know if it works. Thank you, João Fontiela ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > No dia 20/11/2017, às 21:43, Petersson, Anders > escreveu: > > Hi Joao, > We have improved the test_sw4.py script and included it in a new release of SW4, with version 2.01. Note that only the test script was changed, and the source code is identical to version 2.0. The new version can be downloaded from the CIG website, or directly from github. > It would be good to know whether the new script works for you. > Thanks, > Anders > > > On 11/15/17, 12:00 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org " on behalf of cig-seismo-request at geodynamics.org > wrote: > > Send CIG-SEISMO mailing list submissions to > cig-seismo at geodynamics.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > or, via email, send a message with subject or body 'help' to > cig-seismo-request at geodynamics.org > > You can reach the person managing the list at > cig-seismo-owner at geodynamics.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of CIG-SEISMO digest..." > > > Today's Topics: > > 1. Unsuccessuful SW4 test (João Fontiela) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Nov 2017 11:49:42 +0000 > From: João Fontiela > To: cig-seismo at geodynamics.org > Subject: [CIG-SEISMO] Unsuccessuful SW4 test > Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47 at uevora.pt> > Content-Type: text/plain; charset="utf-8" > > Dear all, > > > We installed SW4 (latest release) successfully nevertheless when we execute test_sw4.py (we have installed also python3) we got the following error: > > [root at terra pytest]# ./test_sw4.py > Running all tests for level 0 ... > Starting test # 1 in directory: meshrefine with input file: refine-el-1.in > ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 32512 aborting test > run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in >& refine-el-1.out > DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) > test_sw4 was unsuccessful > > Make.inc configuration is: > FC = gfortran > CXX = mpicxx > EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran > > Can anyone help understand what is going on? > > Yours, > João Fontiela > > ********************* > ********************* > João Fontiela > PhD Student > skype: fontiela > email: fontiela at gmail.com > mobile phone +351 917 862 830 > > Profissional contacts: > email: jfontiela at uevora.pt > Tel: (+351) 266 745 354 ext: 5395 or 5413 > Fax: (+351) 266 745 394 > > Profissional address: > Universidade de Évora > Instituto de Ciências da Terra > Colégio Luís Verney > Rua Romão Ramalho, 59 > 7000-671 Évora - Portugal > > ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > > ------------------------------ > > End of CIG-SEISMO Digest, Vol 118, Issue 8 > ****************************************** > > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela -------------- next part -------------- An HTML attachment was scrubbed... URL: From jfontiela at uevora.pt Tue Nov 21 06:11:33 2017 From: jfontiela at uevora.pt (=?utf-8?Q?Jo=C3=A3o_Fontiela?=) Date: Tue, 21 Nov 2017 14:11:33 +0000 Subject: [CIG-SEISMO] =?utf-8?q?Unsuccessuful_SW4_test_=28Jo=C3=A3o_Fontie?= =?utf-8?b?bGEp?= In-Reply-To: <73966844-7F32-4904-AB80-6A9A788F4D0D@uevora.pt> References: <73966844-7F32-4904-AB80-6A9A788F4D0D@uevora.pt> Message-ID: <717C625F-DBF8-48A8-853B-C95B3A3470EB@uevora.pt> Hi Anders, This time I was able run SW4 tests nevertheless the first test failed and the remaining ones passed. The test log: Running all tests for level 0 ... Starting test # 1 in directory: meshrefine with input file: refine-el-1.in INFO: compare_one_line: col= 0 test= 8.7590587e-06 base= 8.7588796e-06 err= 2.0447820746341383e-05 ERROR: compare_one_line: err= 2.0447820746341383e-05 > tolerance= 1e-05 INFO: compare_one_line: col= 1 test= 1.4151418e-06 base= 1.4151566e-06 err= 1.0458206533397305e-05 ERROR: compare_one_line: err= 1.0458206533397305e-05 > tolerance= 1e-05 Test # 1 Input file: refine-el-1.in FAILED Starting test # 2 in directory: meshrefine with input file: refine-att-1.in Test # 2 Input file: refine-att-1.in PASSED Starting test # 3 in directory: meshrefine with input file: refine-att-2nd-1.in Test # 3 Input file: refine-att-2nd-1.in PASSED Starting test # 4 in directory: attenuation with input file: tw-att-1.in Test # 4 Input file: tw-att-1.in PASSED Starting test # 5 in directory: attenuation with input file: tw-att-2.in Test # 5 Input file: tw-att-2.in PASSED Starting test # 6 in directory: attenuation with input file: tw-topo-att-1.in Test # 6 Input file: tw-topo-att-1.in PASSED Starting test # 7 in directory: pointsource with input file: pointsource-sg-1.in Test # 7 Input file: pointsource-sg-1.in PASSED Starting test # 8 in directory: twilight with input file: flat-twi-1.in Test # 8 Input file: flat-twi-1.in PASSED Starting test # 9 in directory: twilight with input file: flat-twi-2.in Test # 9 Input file: flat-twi-2.in PASSED Starting test # 10 in directory: twilight with input file: gauss-twi-1.in Test # 10 Input file: gauss-twi-1.in PASSED Starting test # 11 in directory: twilight with input file: gauss-twi-2.in Test # 11 Input file: gauss-twi-2.in PASSED Starting test # 12 in directory: lamb with input file: lamb-1.in Test # 12 Input file: lamb-1.in PASSED Out of 12 tests, 11 passed and 1 failed. I thing the problem is related with input file and not with SW4. My point is if I can use SW4 without concerns? Yours, João Fontiela ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > No dia 20/11/2017, às 23:17, João Fontiela escreveu: > > Hi Anders, > > > Thank you very much for your availability to fix this problem > I will test the new release of test_sw4.py and let you know if it works. > > Thank you, > João Fontiela > > ********************* > ********************* > João Fontiela > PhD Student > skype: fontiela > email: fontiela at gmail.com > mobile phone +351 917 862 830 > > Profissional contacts: > email: jfontiela at uevora.pt > Tel: (+351) 266 745 354 ext: 5395 or 5413 > Fax: (+351) 266 745 394 > > Profissional address: > Universidade de Évora > Instituto de Ciências da Terra > Colégio Luís Verney > Rua Romão Ramalho, 59 > 7000-671 Évora - Portugal > > ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > >> No dia 20/11/2017, às 21:43, Petersson, Anders > escreveu: >> >> Hi Joao, >> We have improved the test_sw4.py script and included it in a new release of SW4, with version 2.01. Note that only the test script was changed, and the source code is identical to version 2.0. The new version can be downloaded from the CIG website, or directly from github. >> It would be good to know whether the new script works for you. >> Thanks, >> Anders >> >> >> On 11/15/17, 12:00 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org " on behalf of cig-seismo-request at geodynamics.org > wrote: >> >> Send CIG-SEISMO mailing list submissions to >> cig-seismo at geodynamics.org >> >> To subscribe or unsubscribe via the World Wide Web, visit >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >> or, via email, send a message with subject or body 'help' to >> cig-seismo-request at geodynamics.org >> >> You can reach the person managing the list at >> cig-seismo-owner at geodynamics.org >> >> When replying, please edit your Subject line so it is more specific >> than "Re: Contents of CIG-SEISMO digest..." >> >> >> Today's Topics: >> >> 1. Unsuccessuful SW4 test (João Fontiela) >> >> >> ---------------------------------------------------------------------- >> >> Message: 1 >> Date: Wed, 15 Nov 2017 11:49:42 +0000 >> From: João Fontiela >> To: cig-seismo at geodynamics.org >> Subject: [CIG-SEISMO] Unsuccessuful SW4 test >> Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47 at uevora.pt> >> Content-Type: text/plain; charset="utf-8" >> >> Dear all, >> >> >> We installed SW4 (latest release) successfully nevertheless when we execute test_sw4.py (we have installed also python3) we got the following error: >> >> [root at terra pytest]# ./test_sw4.py >> Running all tests for level 0 ... >> Starting test # 1 in directory: meshrefine with input file: refine-el-1.in >> ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 32512 aborting test >> run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in >& refine-el-1.out >> DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) >> test_sw4 was unsuccessful >> >> Make.inc configuration is: >> FC = gfortran >> CXX = mpicxx >> EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran >> >> Can anyone help understand what is going on? >> >> Yours, >> João Fontiela >> >> ********************* >> ********************* >> João Fontiela >> PhD Student >> skype: fontiela >> email: fontiela at gmail.com >> mobile phone +351 917 862 830 >> >> Profissional contacts: >> email: jfontiela at uevora.pt >> Tel: (+351) 266 745 354 ext: 5395 or 5413 >> Fax: (+351) 266 745 394 >> >> Profissional address: >> Universidade de Évora >> Instituto de Ciências da Terra >> Colégio Luís Verney >> Rua Romão Ramalho, 59 >> 7000-671 Évora - Portugal >> >> ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela >> >> -------------- next part -------------- >> An HTML attachment was scrubbed... >> URL: >> >> ------------------------------ >> >> Subject: Digest Footer >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >> >> ------------------------------ >> >> End of CIG-SEISMO Digest, Vol 118, Issue 8 >> ****************************************** >> >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > > > > ********************* > ********************* > João Fontiela > PhD Student > skype: fontiela > email: fontiela at gmail.com > mobile phone +351 917 862 830 > > Profissional contacts: > email: jfontiela at uevora.pt > Tel: (+351) 266 745 354 ext: 5395 or 5413 > Fax: (+351) 266 745 394 > > Profissional address: > Universidade de Évora > Instituto de Ciências da Terra > Colégio Luís Verney > Rua Romão Ramalho, 59 > 7000-671 Évora - Portugal > > ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo -------------- next part -------------- An HTML attachment was scrubbed... URL: From petersson1 at llnl.gov Tue Nov 21 09:13:26 2017 From: petersson1 at llnl.gov (Petersson, Anders) Date: Tue, 21 Nov 2017 17:13:26 +0000 Subject: [CIG-SEISMO] =?utf-8?q?Unsuccessuful_SW4_test_=28Jo=C3=A3o_Fontie?= =?utf-8?b?bGEp?= Message-ID: <4B0A684C-884D-47B2-B58E-38441881F72F@llnl.gov> Hi Joao, First, I am glad to hear that the test script worked for you. The errors reported by the test script for the first test are very small and I would not worry about it. Your installation is fine. The reason that the test failed is due to round-off errors, related to the linear system that is solved within the mesh refinement algorithm. This should not be noticeable in practical seismic simulations. Good luck with your SW4 project! /Anders On 11/21/17, 7:07 AM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org" wrote: Send CIG-SEISMO mailing list submissions to cig-seismo at geodynamics.org To subscribe or unsubscribe via the World Wide Web, visit http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo or, via email, send a message with subject or body 'help' to cig-seismo-request at geodynamics.org You can reach the person managing the list at cig-seismo-owner at geodynamics.org When replying, please edit your Subject line so it is more specific than "Re: Contents of CIG-SEISMO digest..." Today's Topics: 1. Re: (João Fontiela) ---------------------------------------------------------------------- Message: 1 Date: Tue, 21 Nov 2017 14:11:33 +0000 From: João Fontiela To: cig-seismo at geodynamics.org Subject: Re: [CIG-SEISMO] Unsuccessuful SW4 test (João Fontiela) Message-ID: <717C625F-DBF8-48A8-853B-C95B3A3470EB at uevora.pt> Content-Type: text/plain; charset="utf-8" Hi Anders, This time I was able run SW4 tests nevertheless the first test failed and the remaining ones passed. The test log: Running all tests for level 0 ... Starting test # 1 in directory: meshrefine with input file: refine-el-1.in INFO: compare_one_line: col= 0 test= 8.7590587e-06 base= 8.7588796e-06 err= 2.0447820746341383e-05 ERROR: compare_one_line: err= 2.0447820746341383e-05 > tolerance= 1e-05 INFO: compare_one_line: col= 1 test= 1.4151418e-06 base= 1.4151566e-06 err= 1.0458206533397305e-05 ERROR: compare_one_line: err= 1.0458206533397305e-05 > tolerance= 1e-05 Test # 1 Input file: refine-el-1.in FAILED Starting test # 2 in directory: meshrefine with input file: refine-att-1.in Test # 2 Input file: refine-att-1.in PASSED Starting test # 3 in directory: meshrefine with input file: refine-att-2nd-1.in Test # 3 Input file: refine-att-2nd-1.in PASSED Starting test # 4 in directory: attenuation with input file: tw-att-1.in Test # 4 Input file: tw-att-1.in PASSED Starting test # 5 in directory: attenuation with input file: tw-att-2.in Test # 5 Input file: tw-att-2.in PASSED Starting test # 6 in directory: attenuation with input file: tw-topo-att-1.in Test # 6 Input file: tw-topo-att-1.in PASSED Starting test # 7 in directory: pointsource with input file: pointsource-sg-1.in Test # 7 Input file: pointsource-sg-1.in PASSED Starting test # 8 in directory: twilight with input file: flat-twi-1.in Test # 8 Input file: flat-twi-1.in PASSED Starting test # 9 in directory: twilight with input file: flat-twi-2.in Test # 9 Input file: flat-twi-2.in PASSED Starting test # 10 in directory: twilight with input file: gauss-twi-1.in Test # 10 Input file: gauss-twi-1.in PASSED Starting test # 11 in directory: twilight with input file: gauss-twi-2.in Test # 11 Input file: gauss-twi-2.in PASSED Starting test # 12 in directory: lamb with input file: lamb-1.in Test # 12 Input file: lamb-1.in PASSED Out of 12 tests, 11 passed and 1 failed. I thing the problem is related with input file and not with SW4. My point is if I can use SW4 without concerns? Yours, João Fontiela ********************* ********************* João Fontiela PhD Student skype: fontiela email: fontiela at gmail.com mobile phone +351 917 862 830 Profissional contacts: email: jfontiela at uevora.pt Tel: (+351) 266 745 354 ext: 5395 or 5413 Fax: (+351) 266 745 394 Profissional address: Universidade de Évora Instituto de Ciências da Terra Colégio Luís Verney Rua Romão Ramalho, 59 7000-671 Évora - Portugal ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > No dia 20/11/2017, às 23:17, João Fontiela escreveu: > > Hi Anders, > > > Thank you very much for your availability to fix this problem > I will test the new release of test_sw4.py and let you know if it works. > > Thank you, > João Fontiela > > ********************* > ********************* > João Fontiela > PhD Student > skype: fontiela > email: fontiela at gmail.com > mobile phone +351 917 862 830 > > Profissional contacts: > email: jfontiela at uevora.pt > Tel: (+351) 266 745 354 ext: 5395 or 5413 > Fax: (+351) 266 745 394 > > Profissional address: > Universidade de Évora > Instituto de Ciências da Terra > Colégio Luís Verney > Rua Romão Ramalho, 59 > 7000-671 Évora - Portugal > > ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > >> No dia 20/11/2017, às 21:43, Petersson, Anders > escreveu: >> >> Hi Joao, >> We have improved the test_sw4.py script and included it in a new release of SW4, with version 2.01. Note that only the test script was changed, and the source code is identical to version 2.0. The new version can be downloaded from the CIG website, or directly from github. >> It would be good to know whether the new script works for you. >> Thanks, >> Anders >> >> >> On 11/15/17, 12:00 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org " on behalf of cig-seismo-request at geodynamics.org > wrote: >> >> Send CIG-SEISMO mailing list submissions to >> cig-seismo at geodynamics.org >> >> To subscribe or unsubscribe via the World Wide Web, visit >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >> or, via email, send a message with subject or body 'help' to >> cig-seismo-request at geodynamics.org >> >> You can reach the person managing the list at >> cig-seismo-owner at geodynamics.org >> >> When replying, please edit your Subject line so it is more specific >> than "Re: Contents of CIG-SEISMO digest..." >> >> >> Today's Topics: >> >> 1. Unsuccessuful SW4 test (João Fontiela) >> >> >> ---------------------------------------------------------------------- >> >> Message: 1 >> Date: Wed, 15 Nov 2017 11:49:42 +0000 >> From: João Fontiela >> To: cig-seismo at geodynamics.org >> Subject: [CIG-SEISMO] Unsuccessuful SW4 test >> Message-ID: <2E1502C9-B63E-4C1F-9B8C-88A3EC744B47 at uevora.pt> >> Content-Type: text/plain; charset="utf-8" >> >> Dear all, >> >> >> We installed SW4 (latest release) successfully nevertheless when we execute test_sw4.py (we have installed also python3) we got the following error: >> >> [root at terra pytest]# ./test_sw4.py >> Running all tests for level 0 ... >> Starting test # 1 in directory: meshrefine with input file: refine-el-1.in >> ERROR: Test refine-el-1.in : sw4 returned non-zero exit status= 32512 aborting test >> run_cmd= mpirun -np 1 /share/apps/sw4-master/optimize/sw4 /share/apps/sw4-master/pytest/reference/meshrefine/refine-el-1.in >& refine-el-1.out >> DID YOU USE THE CORRECT SW4 EXECUTABLE? (SPECIFY DIRECTORY WITH -d OPTION) >> test_sw4 was unsuccessful >> >> Make.inc configuration is: >> FC = gfortran >> CXX = mpicxx >> EXTRA_LINK_FLAGS =-L/usr/lib64/atlas -llapack -lblas -lgfortran >> >> Can anyone help understand what is going on? >> >> Yours, >> João Fontiela >> >> ********************* >> ********************* >> João Fontiela >> PhD Student >> skype: fontiela >> email: fontiela at gmail.com >> mobile phone +351 917 862 830 >> >> Profissional contacts: >> email: jfontiela at uevora.pt >> Tel: (+351) 266 745 354 ext: 5395 or 5413 >> Fax: (+351) 266 745 394 >> >> Profissional address: >> Universidade de Évora >> Instituto de Ciências da Terra >> Colégio Luís Verney >> Rua Romão Ramalho, 59 >> 7000-671 Évora - Portugal >> >> ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela >> >> -------------- next part -------------- >> An HTML attachment was scrubbed... >> URL: >> >> ------------------------------ >> >> Subject: Digest Footer >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >> >> ------------------------------ >> >> End of CIG-SEISMO Digest, Vol 118, Issue 8 >> ****************************************** >> >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > > > > ********************* > ********************* > João Fontiela > PhD Student > skype: fontiela > email: fontiela at gmail.com > mobile phone +351 917 862 830 > > Profissional contacts: > email: jfontiela at uevora.pt > Tel: (+351) 266 745 354 ext: 5395 or 5413 > Fax: (+351) 266 745 394 > > Profissional address: > Universidade de Évora > Instituto de Ciências da Terra > Colégio Luís Verney > Rua Romão Ramalho, 59 > 7000-671 Évora - Portugal > > ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ Subject: Digest Footer _______________________________________________ CIG-SEISMO mailing list CIG-SEISMO at geodynamics.org http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo ------------------------------ End of CIG-SEISMO Digest, Vol 118, Issue 11 ******************************************* From yjgao at gfz-potsdam.de Fri Nov 24 09:53:11 2017 From: yjgao at gfz-potsdam.de (yjgao) Date: Fri, 24 Nov 2017 18:53:11 +0100 Subject: [CIG-SEISMO] =?utf-8?q?=28no_subject=29?= Message-ID: <2AD94A28-FEAC-45E1-A7D2-330057086080@gfz-potsdam.de> An HTML attachment was scrubbed... URL: From mzhao at cea-igp.ac.cn Wed Nov 29 01:34:37 2017 From: mzhao at cea-igp.ac.cn (zm) Date: Wed, 29 Nov 2017 17:34:37 +0800 Subject: [CIG-SEISMO] some problems in gravity calculation of SPECFEM3D Message-ID: Hi everyone, I tried to use specfem3d to calculate gravity integration of a very simple model(see attachment 2,the central block has different density).However,it seems the gravity field is quite rough and uneven,look at attachment 1,it seems related to the mesh,but it is a simple and regular mesh,what could it be?Anyone has any ideas,thanks! Best Ming Zhap -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: abnormal_gravity_seismic_snapshot.png Type: image/png Size: 212903 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: gravity_seismic_abonomal_in_middle.jpg Type: image/jpeg Size: 295571 bytes Desc: not available URL: From komatitsch at lma.cnrs-mrs.fr Wed Nov 29 04:43:53 2017 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Wed, 29 Nov 2017 13:43:53 +0100 Subject: [CIG-SEISMO] some problems in gravity calculation of SPECFEM3D In-Reply-To: References: Message-ID: Hi Ming Zhap, Thank you for your message! That is strange, because we used it successfully for instance in http://komatitsch.free.fr/preprints/GJI_Martin_gravimetry_2017.pdf and in http://komatitsch.free.fr/preprints/GJI_Roland_Martin_2013.pdf . Could you configure the code with --enable-debug and try again? Also, note that the recording stations must be above the surface, since the integrals contain a 1/r term the formulation used is not suitable for calculation at the surface, it is fine only above the surface. Thank you, Best regards, Dimitri. On 11/29/2017 10:34 AM, zm wrote: > Hi everyone, > > I tried to use specfem3d to calculate gravity integration of a very > simple model(see attachment 2,the central block has different > density).However,it seems the gravity field is quite rough and > uneven,look at attachment 1,it seems related to the mesh,but it is a > simple and regular mesh,what could it be?Anyone has any ideas,thanks! > > Best > > Ming Zhap > > > > _______________________________________________ > CIG-SEISMO mailing list > CIG-SEISMO at geodynamics.org > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo > -- Dimitri Komatitsch, CNRS Research Director (DR CNRS) Laboratory of Mechanics and Acoustics, Marseille, France http://komatitsch.free.fr