From tsuboi at jamstec.go.jp Mon Jul 2 00:41:53 2018 From: tsuboi at jamstec.go.jp (Seiji Tsuboi) Date: Mon, 2 Jul 2018 16:41:53 +0900 Subject: [CIG-SEISMO] specfem3d_globe-7.0.0 overflow Message-ID: Dear Sir, I have used specfem3d_globe-7.0.0 on intel cluster and have overflow problems. If I set NPROC=16 and NEX=768 (1536 cores), I could calculate seismograms correctly. However if I increase NEX, for example NPROC=17, NEX=816 (1734 cores), xspecfem3D stops with overflow in the very beginning of computation. ATTENUATION is set as true and it seems that the program overflows in "multiply_accel_elastic" routine. I tried both ifort and gfortran but the situation is the same. I would be grateful if someone would let me know how I should solve this issue. Best regards, Seiji Tsuboi -------------- next part -------------- An HTML attachment was scrubbed... URL: From komatitsch at lma.cnrs-mrs.fr Mon Jul 2 14:02:35 2018 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Mon, 2 Jul 2018 23:02:35 +0200 Subject: [CIG-SEISMO] specfem3d_globe-7.0.0 overflow In-Reply-To: References: Message-ID: <07cd1545-4ad2-787b-9d87-e8fdf75822fc@lma.cnrs-mrs.fr> Dear Tsuboi-san, It is probably the time step computed automatically by the code that is slightly above the CFL stability limit in that case. Thus, you probably just need to reduce the value of variable DT a bit (say 10% or so). Best regards, Dimitri. On 07/02/2018 09:41 AM, Seiji Tsuboi wrote: > Dear Sir, > > I have used specfem3d_globe-7.0.0 on intel cluster and have overflow > problems. > If I set NPROC=16 and NEX=768 (1536 cores), I could calculate > seismograms correctly. However if I increase NEX, for example NPROC=17, > NEX=816 (1734 cores), xspecfem3D stops with overflow in the very > beginning of computation. ATTENUATION is set as true and it seems that > the program overflows in "multiply_accel_elastic" routine. I tried both > ifort and gfortran but the situation is the same. I would be grateful if > someone would let me know how I should solve this issue. > > Best regards, > Seiji Tsuboi > > > _______________________________________________ > 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 tsuboi at jamstec.go.jp Mon Jul 2 21:41:39 2018 From: tsuboi at jamstec.go.jp (Seiji Tsuboi) Date: Tue, 3 Jul 2018 13:41:39 +0900 Subject: [CIG-SEISMO] specfem3d_globe-7.0.0 overflow Message-ID: Dear Komatitsch-san Thank you very much for your very quick reply to my inquiry. Since I could not find where I should reduce DT in the program, I have turned both USE_LDDRK and INCREASE_CFL_FOR_LDDRK on. Then computation proceeds for NPROC=17 and NEX=816 case without overflow but the time per iteration becomes almost 6 times longer. It may be necessary to compute for more than one day. Is this what we expect? Should I turn USE_LDDRK false and reduce DT instead? Best regards, Seiji Tsuboi 2018-07-03 6:02 GMT+09:00 Dimitri Komatitsch : > > Dear Tsuboi-san, > > It is probably the time step computed automatically by the code that is > slightly above the CFL stability limit in that case. Thus, you probably > just need to reduce the value of variable DT a bit (say 10% or so). > > Best regards, > Dimitri. > > On 07/02/2018 09:41 AM, Seiji Tsuboi wrote: > >> Dear Sir, >> >> I have used specfem3d_globe-7.0.0 on intel cluster and have overflow >> problems. >> If I set NPROC=16 and NEX=768 (1536 cores), I could calculate seismograms >> correctly. However if I increase NEX, for example NPROC=17, NEX=816 (1734 >> cores), xspecfem3D stops with overflow in the very beginning of >> computation. ATTENUATION is set as true and it seems that the program >> overflows in "multiply_accel_elastic" routine. I tried both ifort and >> gfortran but the situation is the same. I would be grateful if someone >> would let me know how I should solve this issue. >> >> Best regards, >> Seiji Tsuboi >> >> >> _______________________________________________ >> 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 alejandro.duran at sed.ethz.ch Fri Jul 6 00:23:16 2018 From: alejandro.duran at sed.ethz.ch (Duran Neme Javier Alejandro) Date: Fri, 6 Jul 2018 07:23:16 +0000 Subject: [CIG-SEISMO] Question: Acoustic simulation equivalent to a curl-free Elastic one (both homogeneous media)? SPECFEM2D v. 6.1 Message-ID: <79888D86BA6AC7488EDCF1DEAAA176B932CA1960@MBX115.d.ethz.ch> Dear all I am comparing the results obtained in an acoustic simulation with a curl-free elastic one (both homogeneous media). I would want to know if both numerical experiments are equivalent on SPECFEM2D (I’m using version 6.1)? In the acoustic case (vs=0m/s), I sent a plane wave (generated from several pressure point sources aligned horizontally) and I measure the seismograms at a distant point. In the elastic case I sent also a plane P-wave, generated from several explosive sources aligned horizontally (An explosive source should generate only P-waves in a homogeneous elastic medium). For the sources I used moment tensors in this case (Mxx=Mzz=1, Mxz=0) and measure the seismograms at the same point as before. I observed that Curl(u)~0 (which is expected if there are almost no S-waves) and Div(u) is non zero. In the acoustic case I can only use Pressure sources and my receiver is an Hydrophone. In the elastic one my sources are moment tensors and my receiver is a geophone , thus I can measure displacement u, Curl(u) and Div(u). In both scenarios I used Ricker time function with the same amplification factor and frequency. My question is: • Are both experiments equivalent (acoustic <--> elastic with curl(u)=0?) even though sources and receivers are different? If they are equivalent, how can I compare both synthetic data (Pressure and Div(u) ) to make sure they are the same? I tried Pressure=-Bulk_Modulus*Div(u), but they are several orders of magnitude off. • In the acoustic case, which are the units of the source? [Pascals_2D]?=[N/m]? I appreciate your reply. I’m using this information in order to compute the scattering cross section of a localised velocity anomaly in both scenarios. If both are equivalent, I expected to obtain the same cross section, but so far it is not the case. With kind regards, Alejandro -------------- next part -------------- An HTML attachment was scrubbed... URL: From komatitsch at lma.cnrs-mrs.fr Mon Jul 9 14:31:10 2018 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Mon, 9 Jul 2018 23:31:10 +0200 Subject: [CIG-SEISMO] specfem3d_globe-7.0.0 overflow In-Reply-To: References: Message-ID: <274842a7-fec4-74ff-4e65-a383e1f97368@lma.cnrs-mrs.fr> Dear Tsuboi-san, Thank you very much for your answer. To change the value of the time step DT automatically computed by the code, you can edit file src/shared/get_timestep_and_layers.f90 and put something like: DT = DT * 0.90 at the end of that routine. This should fix the problem. Best regards, Dimitri Komatitsch. On 07/03/2018 06:41 AM, Seiji Tsuboi wrote: > Dear Komatitsch-san > > Thank you very much for your very quick reply to my inquiry. > Since I could not find where I should reduce DT in the program, I have > turned both USE_LDDRK and INCREASE_CFL_FOR_LDDRK on. Then computation > proceeds for NPROC=17 and NEX=816 case without overflow but the time per > iteration becomes almost 6 times longer. It may be necessary to compute > for more than one day. Is this what we expect? Should I turn USE_LDDRK > false and reduce DT instead? > > Best regards, > Seiji Tsuboi > > 2018-07-03 6:02 GMT+09:00 Dimitri Komatitsch >: > > > Dear Tsuboi-san, > > It is probably the time step computed automatically by the code that > is slightly above the CFL stability limit in that case. Thus, you > probably just need to reduce the value of variable DT a bit (say 10% > or so). > > Best regards, > Dimitri. > > On 07/02/2018 09:41 AM, Seiji Tsuboi wrote: > > Dear Sir, > > I have used specfem3d_globe-7.0.0 on intel cluster and have > overflow problems. > If I set NPROC=16 and NEX=768 (1536 cores), I could calculate > seismograms correctly. However if I increase NEX, for example > NPROC=17, NEX=816 (1734 cores), xspecfem3D stops with overflow > in the very beginning of computation. ATTENUATION is set as true > and it seems that the program overflows in > "multiply_accel_elastic" routine. I tried both ifort and > gfortran but the situation is the same. I would be grateful if > someone would let me know how I should solve this issue. > > Best regards, > Seiji Tsuboi > > > _______________________________________________ > 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 > > -- Dimitri Komatitsch, CNRS Research Director (DR CNRS) Laboratory of Mechanics and Acoustics, Marseille, France http://komatitsch.free.fr From gwwang80 at 163.com Wed Jul 11 18:26:39 2018 From: gwwang80 at 163.com (=?GBK?B?zfW5887E?=) Date: Thu, 12 Jul 2018 09:26:39 +0800 (CST) Subject: [CIG-SEISMO] Question about SPECFEM3D install Message-ID: <56a7453c.2ebb.1648c18e65c.Coremail.gwwang80@163.com> Hi sir: I meet some problems When I compile the software SPECFEM 3D Cartesian. There is a syntax error in the configure file when i run "./configure",but I can not find the error.Could you give me some advices?My cluster's operating system is linux(Cent os) and CPU is intel. I have no GPU. Thanks for your help! The configure config.log and configure.ac files were sent as attachments. the command is "./configure FC=ifort CC=icc MPIFC=mpiifort --with-mpi" The error message is as follows -------------------------------------------------------------- checking for gcc... icc checking whether we are using the GNU C compiler... yes checking whether icc accepts -g... yes checking for icc option to accept ISO C89... none needed checking for dummy main to link with Fortran libraries... none checking for Fortran name-mangling scheme... lower case, underscore, no extra underscore checking for Fortran flag to compile .f90 files... none checking for Fortran flag needed to allow free-form source... none checking for Fortran flag to compile preprocessed .F files... none checking how to define symbols for preprocessed Fortran... -D checking for Fortran flag to compile preprocessed .F90 files... none checking for gcc... (cached) icc checking whether we are using the GNU C compiler... (cached) yes checking whether icc accepts -g... (cached) yes checking for icc option to accept ISO C89... (cached) none needed ## ------ ## ## SCOTCH ## ## ------ ## checking for scotchfarchinit in -lscotch... yes checking whether Scotch is usable... ./configure: line 6181: syntax error near unexpected token `{' ./configure: line 6181: ` ACX_PTHREAD({ $as_echo "$as_me:$LINENO: result: pthread found" >&5' --------------------------------------------------------------------------- -- Wang Guiwen -- shanxi normal university , China -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: configure Type: application/octet-stream Size: 268395 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: configure.ac Type: application/octet-stream Size: 13987 bytes Desc: not available URL: -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: config.log URL: From enkst000 at gmail.com Mon Jul 16 15:20:39 2018 From: enkst000 at gmail.com (ed cas) Date: Mon, 16 Jul 2018 18:20:39 -0400 Subject: [CIG-SEISMO] Model error in mesh Message-ID: Dear developers and users, I have been trying to produce a readable mesh file for specfem2d by using gmsh. I produced the mesh inside gmsh, I saved the mesh and then paste it to the MESH directory and then run the command python LibGmsh2Specfem_convert_Gmsh_to_Specfem2D_official.py model_shallow_v42 -t A -b A -r A -l A That produced the five files in the same directory. Please find attached the file I have been trying to run (model_shallow_v42). However, when running specfem2d gives error telling me about Jacobians, I have been trying for a couple of days and cannot get with the error, I checked the line loops and seems fine , I have resized the elements and with no success. I was hoping you could help me with this. Could you please check my file and indicate me what is wrong with it, and any solution. I would really appreciate your help, Thank you for your time and help, Sincerely, ------------------------------------------- Edwin Castillo -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: model_shallow_v42.geo Type: application/octet-stream Size: 5569 bytes Desc: not available URL: From thomas.spelce at intel.com Mon Jul 16 18:47:47 2018 From: thomas.spelce at intel.com (Spelce, Thomas) Date: Tue, 17 Jul 2018 01:47:47 +0000 Subject: [CIG-SEISMO] RAJA version of SW4lite Message-ID: <0D8F7D39-43FE-4E1A-8486-984EBA84E423@intel.com> I’d like to exercise the RAJA version of SW4lite. I see provisions for compiling in the Makefile and have built a RAJA library, but it doesn’t appear that the source enables and RAJA calls. Please advise how I can get a version that supports RAJA, ideally the lite version for performance testing using a simulator. Thanks, Tom Spelce -------------- next part -------------- An HTML attachment was scrubbed... URL: From komatitsch at lma.cnrs-mrs.fr Tue Jul 17 04:26:06 2018 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Tue, 17 Jul 2018 13:26:06 +0200 Subject: [CIG-SEISMO] Model error in mesh In-Reply-To: References: Message-ID: <604ee115-85f8-8e47-f30e-76f2b48cd489@lma.cnrs-mrs.fr> Dear Edwin, Thank you for your message. Let me forward it to my colleague Paul Cristini, who is a Gmsh + SPECFEM expert. Paul is currently traveling, thus he may need a couple of weeks to send you an answer. (he is the author of the LibGmsh2Specfem_convert_Gmsh_to_Specfem2D_official.py script) Best regards, Dimitri. On 07/17/2018 12:20 AM, ed cas wrote: > Dear developers and users, > > > I have been trying to produce a readable mesh file > for specfem2d by using gmsh. > > I produced the mesh inside gmsh, I saved the mesh and then paste it to > the MESH directory > > and then run the command > > python LibGmsh2Specfem_convert_Gmsh_to_Specfem2D_official.py > model_shallow_v42 -t A -b A -r A -l A > That produced the five files in the same directory. > > Please find attached the file I have been trying to run (model_shallow_v42). > > However, when running specfem2d gives error telling me about Jacobians, > I have been trying for a couple of days and cannot get with the error, I > checked the > line loops and seems  fine , I have resized the elements and with no > success. > > I was hoping you could help me with this. > Could you please check my file and indicate me what is wrong with it, > and any solution. > > I would really appreciate your help, > > Thank you for your time and help, > > Sincerely, > > ------------------------------------------- > Edwin Castillo > > > > _______________________________________________ > 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 petersson1 at llnl.gov Tue Jul 17 14:50:47 2018 From: petersson1 at llnl.gov (Petersson, Anders) Date: Tue, 17 Jul 2018 21:50:47 +0000 Subject: [CIG-SEISMO] RAJA version of SW4lite Message-ID: <27878CDA-38E1-4F3A-A33A-8E9E0908962D@llnl.gov> Hi Tom, Did you try the source code in the Raja branch of the sw4lite repository? Thanks, Anders On 7/17/18, 12:08 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org" wrote: RAJA version of SW4lite From maximilian.trapp at rub.de Wed Jul 25 06:40:05 2018 From: maximilian.trapp at rub.de (Max) Date: Wed, 25 Jul 2018 15:40:05 +0200 Subject: [CIG-SEISMO] Anisotropy in SpecFEM3D Message-ID: Dear all, I am using SpecFEM3D for a longer while now and need to implement anisotropy. Unfortunately, I am not very deep into coding and in the user manual there is not very much written about anisotropy. It is said that I have to call the subroutine by the command "call model_aniso(iflag_aniso, rho, vp, vs, & c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26, & c33,c34,c35,c36,c44,c45,c46,c55,c56,c66)". But where do I implement this? Where do I have to put it in the code? And is there an example on anisotropy somewhere? Best regards, Max -------------- next part -------------- An HTML attachment was scrubbed... URL: From komatitsch at lma.cnrs-mrs.fr Wed Jul 25 07:17:16 2018 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Wed, 25 Jul 2018 16:17:16 +0200 Subject: [CIG-SEISMO] Anisotropy in SpecFEM3D In-Reply-To: References: Message-ID: <0949b636-abc5-dc7c-63c0-af74ea450cc9@lma.cnrs-mrs.fr> Dear Max, The routine is in src/generate_databases/model_aniso.f90, which contains an example. Best regards, Dimitri. On 07/25/2018 03:40 PM, Max wrote: > Dear all, > > I am using SpecFEM3D for a longer while now and need to implement > anisotropy. Unfortunately, I am not very deep into coding and in the > user manual there is not very much written about anisotropy. > It is said that I have to call the subroutine by the command "call > model_aniso(iflag_aniso, rho, vp, vs, & > c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26, & > c33,c34,c35,c36,c44,c45,c46,c55,c56,c66)". > > But where do I implement this? Where do I have to put it in the code? > And is there an example on anisotropy somewhere? > > Best regards, > > Max > > > > _______________________________________________ > 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 maximilian.trapp at rub.de Wed Jul 25 09:48:52 2018 From: maximilian.trapp at rub.de (Max) Date: Wed, 25 Jul 2018 18:48:52 +0200 Subject: [CIG-SEISMO] Anisotropy in SpecFEM3D In-Reply-To: <0949b636-abc5-dc7c-63c0-af74ea450cc9@lma.cnrs-mrs.fr> References: <0949b636-abc5-dc7c-63c0-af74ea450cc9@lma.cnrs-mrs.fr> Message-ID: <1ce440df-1853-736d-f11f-3e749e14705d@rub.de> Dear Dimitri, thank you for your reply. I saw that file and tried to use it, but I think you have to include the 'call model_aniso'-command somewhere, don't you? If I mark the flag for anisotropy in the nummaterial_velocity_file, where do I enter the material parameters? Best regards, Max Am 25.07.2018 um 16:17 schrieb Dimitri Komatitsch: > > Dear Max, > > The routine is in src/generate_databases/model_aniso.f90, which > contains an example. > > Best regards, > Dimitri. > > On 07/25/2018 03:40 PM, Max wrote: >> Dear all, >> >> I am using SpecFEM3D for a longer while now and need to implement >> anisotropy. Unfortunately, I am not very deep into coding and in the >> user manual there is not very much written about anisotropy. >> It is said that I have to call the subroutine by the command "call >> model_aniso(iflag_aniso, rho, vp, vs, & >> c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26, & >> c33,c34,c35,c36,c44,c45,c46,c55,c56,c66)". >> >> But where do I implement this? Where do I have to put it in the code? >> And is there an example on anisotropy somewhere? >> >> Best regards, >> >> Max >> >> >> >> _______________________________________________ >> CIG-SEISMO mailing list >> CIG-SEISMO at geodynamics.org >> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo >> > From johannes.aichele at inserm.fr Fri Jul 27 10:10:22 2018 From: johannes.aichele at inserm.fr (johannes aichele) Date: Fri, 27 Jul 2018 19:10:22 +0200 Subject: [CIG-SEISMO] what is fi in get_poroelastic_velocities.f90 Message-ID: Dear Specfem experts, I am totally new to specfem and currently try to understand the get_poroelastic_velocities.f90 script. To my understanding it calculates the wave speeds in a poroelastic medium from the 7 Biot parameters. (mu_fr,phi, tort,rho_s,rho_f,eta_f,perm_xx) I understand furthermore  that f0 is the dominant frequency which I guess means the frequency for which the wave speed is analyzed. However I do not get the physical meaning of fi. I would be very grateful if someone could help me out. Kind regards, Johannes Aichele Johannes AICHELE Phd Student - Doctorant Laboratoire Applications Thérapeutiques des Ultrasons LabTAU, INSERM unité 1032 151 cours Albert Thomas 69424 LYON Cedex 03 France From komatitsch at lma.cnrs-mrs.fr Fri Jul 27 10:33:16 2018 From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch) Date: Fri, 27 Jul 2018 19:33:16 +0200 Subject: [CIG-SEISMO] what is fi in get_poroelastic_velocities.f90 In-Reply-To: References: Message-ID: Dear Johannes, You can refer to https://authors.library.caltech.edu/13230/1/MORgji08.pdf and https://onlinelibrary.wiley.com/doi/pdf/10.1111/j.1365-246X.2011.04985.x (you can also contact Christina, cc'ed) Best regards, Dimitri. On 07/27/2018 07:10 PM, johannes aichele wrote: > Dear Specfem experts, > > I am totally new to specfem and currently try to understand the > get_poroelastic_velocities.f90 script. > To my understanding it calculates the wave speeds in a poroelastic > medium from the 7 Biot parameters. > (mu_fr,phi, tort,rho_s,rho_f,eta_f,perm_xx) > > I understand furthermore  that f0 is the dominant frequency which I > guess means the frequency for which the wave speed is analyzed. > However I do not get the physical meaning of fi. > > I would be very grateful if someone could help me out. > > > Kind regards, Johannes Aichele > > Johannes AICHELE > Phd Student - Doctorant > Laboratoire Applications Thérapeutiques des Ultrasons > LabTAU, INSERM unité 1032 > 151 cours Albert Thomas > 69424 LYON Cedex 03 > France > > _______________________________________________ > 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 morency1 at llnl.gov Fri Jul 27 11:34:51 2018 From: morency1 at llnl.gov (Morency, Christina E.) Date: Fri, 27 Jul 2018 18:34:51 +0000 Subject: [CIG-SEISMO] what is fi in get_poroelastic_velocities.f90 In-Reply-To: References: Message-ID: <99AFE2B5-8E48-4B73-8335-F8D4F5AAAAB7@llnl.gov> Hi Dimitri, Thanks for the connection. J'espere que tout va bien! Johannes - I see another email from you. I will take this conversion offline shortly. Cheers, Chris On 7/27/18, 10:33 AM, "Dimitri Komatitsch" wrote: Dear Johannes, You can refer to https://authors.library.caltech.edu/13230/1/MORgji08.pdf and https://onlinelibrary.wiley.com/doi/pdf/10.1111/j.1365-246X.2011.04985.x (you can also contact Christina, cc'ed) Best regards, Dimitri. On 07/27/2018 07:10 PM, johannes aichele wrote: > Dear Specfem experts, > > I am totally new to specfem and currently try to understand the > get_poroelastic_velocities.f90 script. > To my understanding it calculates the wave speeds in a poroelastic > medium from the 7 Biot parameters. > (mu_fr,phi, tort,rho_s,rho_f,eta_f,perm_xx) > > I understand furthermore that f0 is the dominant frequency which I > guess means the frequency for which the wave speed is analyzed. > However I do not get the physical meaning of fi. > > I would be very grateful if someone could help me out. > > > Kind regards, Johannes Aichele > > Johannes AICHELE > Phd Student - Doctorant > Laboratoire Applications Thérapeutiques des Ultrasons > LabTAU, INSERM unité 1032 > 151 cours Albert Thomas > 69424 LYON Cedex 03 > France > > _______________________________________________ > 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