[cig-commits] r8420 - seismo/2D/SPECFEM2D/trunk/MAILLE90

walter at geodynamics.org walter at geodynamics.org
Fri Dec 7 15:44:50 PST 2007


Author: walter
Date: 2007-12-07 15:44:49 -0800 (Fri, 07 Dec 2007)
New Revision: 8420

Added:
   seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_aniso_ok
   seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_iso_ok
Log:
added anisotropic Par file for SEM 2D Dimitri for tests in Rouen


Added: seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_aniso_ok
===================================================================
--- seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_aniso_ok	2004-05-21 00:40:31 UTC (rev 8419)
+++ seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_aniso_ok	2007-12-07 23:44:49 UTC (rev 8420)
@@ -0,0 +1,81 @@
+# ----------------------------------------------------------------
+#
+#    This is the parameter file
+#    Put variable names first and actual value after 15th column
+#
+# ----------------------------------------------------------------
+<-           ->
+#
+# File names and path for different outputs
+#
+title         =Test pour INSA Rouen full anisotropy
+topofile      =topo_rouen_aniso.dat
+interffile    =none
+#
+# geometry of the model (origin lower-left corner = 0,0) and mesh description
+#
+xmin          =0.0d0         ! abscissa of left side of the model
+xmax          =0.06d0        ! abscissa of right side of the model
+nx            =80            ! number of elements along X
+nz            =80            ! number of elements along Z
+ngnod         =4             ! nb noeuds de controle pour blocs generes (4 ou 9)
+ratio         =0.967741      ! ratio pour separation en deux zones
+initialfield  =.false.       ! use a plane wave as source or not
+ireadmodel    =.false.       ! use a plane wave as source or not
+#
+# absorbing boundaries parameters
+#
+absorbhaut    =.false.       ! Absorbing boundary active or not
+absorbbas     =.false.
+absorbgauche  =.false.
+absorbdroite  =.false.
+#
+# time step parameters
+#
+nt            =1400          ! nb total de pas de temps
+dt            =10.d-9        ! valeur du pas de temps
+#
+# source parameters
+#
+isource_surf  =.false.       ! sources dans le volume ou a la surface
+xs            =0.03d0        ! source location x in meters
+zs            =0.03d0        ! source location z in meters
+f0            =1.3d6         ! central source frequency (Hz)
+t0            =0.0000008d0   ! time delay of the source in seconds
+source_type   =1             ! source type : force=1 or explosion=2
+angle         =0.            ! angle of the source (for a force only)
+factor        =1.d6          ! amplification factor
+#
+# receiver line parameters
+#
+ienreg_surf   =.false.       ! enregistrement dans le volume ou a la surface
+isismostype   =1             ! record 1=displacement 2=velocity 3=acceleration
+nrec          =11            ! number of receivers
+xdeb          =0.            ! first receiver x in meters
+zdeb          =0.048d0       ! first receiver z in meters
+xfin          =0.06d0        ! last receiver x in meters
+zfin          =0.048d0       ! last receiver z in meters
+anglerec      =0.d0          ! angle to rotate the components at the receivers
+#
+# display parameters
+#
+itaff         =100           ! display frequency in time steps
+ivecttype     =1             ! display 1=displacement 2=velocity 3=acceleration
+cutvect       =1.            ! amplitude min affichee en % pour vector plots
+imeshvect     =.true.        ! display mesh on vector plots or not
+imodelvect    =.false.       ! display velocity model on vector plots or not
+iboundvect    =.true.        ! display boundary conditions on vector plots
+interpol      =.true.        ! interpolation of the display or not
+iptsdisp      =6             ! nb of points for interpolation of the display
+isubsamp      =2             ! subsampling of color snapshots
+ignuplot      =.false.       ! generate a GNUPLOT file for the grid
+ioutputgrid   =.false.       ! save the grid in a text file or not
+#
+# velocity and density model (nx,nz)
+#
+nbmodels      =2             ! nb de modeles differents (rho,vp,vs)
+1 0 8920.d0 2500.d0 1443.375d0 0 0
+2 0 8920.d0 2500.d0 1443.375d0 0 0
+nbzone        =2             ! nb of zones and model number for each zone
+1 80 1 80 1
+3 51 3 38 2

Added: seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_iso_ok
===================================================================
--- seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_iso_ok	2004-05-21 00:40:31 UTC (rev 8419)
+++ seismo/2D/SPECFEM2D/trunk/MAILLE90/Par_rouen_iso_ok	2007-12-07 23:44:49 UTC (rev 8420)
@@ -0,0 +1,81 @@
+# ----------------------------------------------------------------
+#
+#    This is the parameter file
+#    Put variable names first and actual value after 15th column
+#
+# ----------------------------------------------------------------
+<-           ->
+#
+# File names and path for different outputs
+#
+title         =Test pour INSA Rouen
+topofile      =topo_rouen.dat
+interffile    =none
+#
+# geometry of the model (origin lower-left corner = 0,0) and mesh description
+#
+xmin          =0.0d0         ! abscissa of left side of the model
+xmax          =4000.d0       ! abscissa of right side of the model
+nx            =80            ! number of elements along X
+nz            =60            ! number of elements along Z
+ngnod         =4             ! nb noeuds de controle pour blocs generes (4 ou 9)
+ratio         =0.967741      ! ratio pour separation en deux zones
+initialfield  =.false.       ! use a plane wave as source or not
+ireadmodel    =.false.       ! use a plane wave as source or not
+#
+# absorbing boundaries parameters
+#
+absorbhaut    =.false.       ! Absorbing boundary active or not
+absorbbas     =.true.
+absorbgauche  =.true.
+absorbdroite  =.true.
+#
+# time step parameters
+#
+nt            =1600          ! nb total de pas de temps
+dt            =1.d-3         ! valeur du pas de temps
+#
+# source parameters
+#
+isource_surf  =.false.       ! sources dans le volume ou a la surface
+xs            =2000.         ! source location x in meters
+zs            =1500.         ! source location z in meters
+f0            =10.0          ! central source frequency (Hz)
+t0            =0.11          ! time delay of the source in seconds
+source_type   =1             ! source type : force=1 or explosion=2
+angle         =0.            ! angle of the source (for a force only)
+factor        =1.d10         ! amplification factor
+#
+# receiver line parameters
+#
+ienreg_surf   =.false.       ! enregistrement dans le volume ou a la surface
+isismostype   =2             ! record 1=displacement 2=velocity 3=acceleration
+nrec          =11            ! number of receivers
+xdeb          =0.            ! first receiver x in meters
+zdeb          =2200.         ! first receiver z in meters
+xfin          =4000.         ! last receiver x in meters
+zfin          =2200.         ! last receiver z in meters
+anglerec      =0.d0          ! angle to rotate the components at the receivers
+#
+# display parameters
+#
+itaff         =100           ! display frequency in time steps
+ivecttype     =2             ! display 1=displacement 2=velocity 3=acceleration
+cutvect       =1.            ! amplitude min affichee en % pour vector plots
+imeshvect     =.true.        ! display mesh on vector plots or not
+imodelvect    =.false.       ! display velocity model on vector plots or not
+iboundvect    =.true.        ! display boundary conditions on vector plots
+interpol      =.true.        ! interpolation of the display or not
+iptsdisp      =6             ! nb of points for interpolation of the display
+isubsamp      =2             ! subsampling of color snapshots
+ignuplot      =.false.       ! generate a GNUPLOT file for the grid
+ioutputgrid   =.false.       ! save the grid in a text file or not
+#
+# velocity and density model (nx,nz)
+#
+nbmodels      =2             ! nb de modeles differents (rho,vp,vs)
+1 0 2200.d0 2500.d0 1443.375d0 0 0
+2 0 2200.d0 2500.d0 1443.375d0 0 0
+nbzone        =2             ! nb of zones and model number for each zone
+1 80 1 60 1
+3 51 3 38 2



More information about the cig-commits mailing list