info

Spectral Element Method for wave propagation and rupture dynamics.

 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 Program  S P E C F E M : start
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
                                                   
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 D a t e : 06 - 06 - 2007                                 T i m e  : 17:28:43
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 
***********************************************
*            I n p u t   p h a s e            *
***********************************************
 


 E c h o   S e t t i n g s
 =========================

     Echo during input phase. . . . . . . . . . . (Verbose(1)) = T
     Echo during init phase . . . . . . . . . . . (Verbose(2)) = T
     Echo during check phase. . . . . . . . . . . (Verbose(3)) = T
     Echo during run phase. . . . . . . . . . . . (Verbose(4)) = T
     Display frequency  . . . . . . . . . . . . .(ItSnapshots) = 5000
     First display. . . . . . . . . . . . . . . .(ItSnapshot1) = 0
     Basic info output frequency. . . . . . . . . . . (ItInfo) = 200



 G e n e r a l   P a r a m e t e r s
 ===================================

     Execution mode . . . . . . . . . . . . . . . (iexec) = solve     
     Number of nodes per edge . . . . . . . . . . .(ngll) = 9
     Number of d.o.f per node . . . . . . . . . . .(ndof) = 2
     Highest frequency to be resolved . . . . . . .(fmax) =   25.000E+00


 M e s h   G e n e r a t i o n
 =============================

     Method. . . . . . . . . . . . . . . . .(method) = CARTESIAN 
     Minimum X . . . . . . . . . . . . . . (xlim(1)) =  0.000E+00
     Maximum X . . . . . . . . . . . . . . (xlim(2)) =  4.000E+03
     Minimum Z . . . . . . . . . . . . . . (zlim(1)) = -2.000E+03
     Maximum Z . . . . . . . . . . . . . . (zlim(2)) =  0.000E+00
     Number of elements along X. . . . . .(nelem(1)) =    40
     Number of elements along Z. . . . . .(nelem(2)) =    20
     Cut by horizontal fault . . . . . . . .(faultx) = F


 T i me   i n t e g r a t i o n
 ==============================

     Scheme. . . . . . . . . . . . . .(kind) = leapfrog  
     Number of time steps. . . . . (NbSteps) = 3000
     Time step increment . . . . . . . .(Dt) =  500.000E-06
     Total simulation duration . (TotalTime) =    1.500E+00


 M a t e r i a l   s e t s :   2 D  e l a s t i c i t y
 ======================================================

     Number of material sets . . . . . . . . . =     1

     Material number . . . . . . . . . . (tag) =     1
     Type  . . . . . . . . . . . . . . .(mode) = ISOTR     
     P-wave velocity . . . . . . . . . . .(cp) =   3.200E+03
     S-wave velocity . . . . . . . . . . .(cs) =   1.848E+03
     Mass density. . . . . . . . . . . (denst) =   2.000E+03
     Poisson's ratio . . . . . . . . . (poiss) = 250.008E-03
     First Lame parameter Lambda . . . .(alam) =   6.827E+09
     Second Lame parameter Mu. . . . . . (amu) =   6.827E+09
     Bulk modulus K. . . . . . . . . . .(Kvol) =  11.378E+09
     Young's modulus E. . . . . . . . (young) =  17.066E+09


 B o u n d a r y   C o n d i t i o n s
 =====================================


     Boundary tag. . . . . . . . . . . . (tag) = 1
     Boundary condition. . . . . . . . .(kind) = ABSORB
     Type of absorbing boundary. . . .(stacey) = Stacey              
     Allow incident wave . . . . . . (let_wave) = T

     Boundary tag. . . . . . . . . . . . (tag) = 2
     Boundary condition. . . . . . . . .(kind) = ABSORB
     Type of absorbing boundary. . . .(stacey) = Stacey              
     Allow incident wave . . . . . . (let_wave) = T

     Boundary tag. . . . . . . . . . . . (tag) = 4
     Boundary condition. . . . . . . . .(kind) = ABSORB
     Type of absorbing boundary. . . .(stacey) = Stacey              
     Allow incident wave . . . . . . (let_wave) = T


 S o u r c e s
 =============

     X-position (meters). . . . .(coord(1)) =    1.500E+03
     Y-position (meters). . . . .(coord(2)) =  -50.000E+00
     Source time function . .(TimeFunction) = Ricker
     Fundamental frequency (Hz) . . . .(f0) =  10.000E+00
     Time delay (s) . . . . . . . . (onset) = 150.000E-03
     Multiplying factor . . . . . . (ampli) =   1.000E+00
     Source Type. . . . . . . . . . . . . . = Collocated Force
     If P-SV: counterclockwise angle / up . = 0.00


 R e c e i v e r s
 =================

     Number of receivers . . . . . . . . . . . . (number) = 2
     Subsampling for seismograms recording . . . .(isamp) = 1
     Field recorded. . . . . . . . . . . . . . . .(field) = D
     Axis of the seismogram plot . . . . . . . . .(irepr) = D


 S n a p s h o t   O u t p u t s
 ===============================

     Save results in PS file or not . . . . .(postscript) = T
     Save results in AVS file or not. . . . . . . . (avs) = F
     Save results in Visual3 file or not. . . . (visual3) = F
     Save results in binary file or not . . . . . . (bin) = T
     Selected fields :
       Displacement . . . . . . . . . . . . . . . . . . . = F
       Velocity . . . . . . . . . . . . . . . . . . . . . = T
       Acceleration . . . . . . . . . . . . . . . . . . . = F
       Strain . . . . . . . . . . . . . . . . . . . . . . = F
       Stress . . . . . . . . . . . . . . . . . . . . . . = F
     Selected components for PostScript snapshots :
       X  . . . . . . . . . . . . . . . . . . . . . . . . = F
       Y  . . . . . . . . . . . . . . . . . . . . . . . . = F
       Z  . . . . . . . . . . . . . . . . . . . . . . . . = F
       Amplitude  . . . . . . . . . . . . . . . . . . . . = T


 P o s t S c r i p t   O u t p u t s
 ===================================

     Plot mesh . . . . . . . . . . . . . . . . . . (mesh) = F
     Numbered mesh . . . . . . . . . . . . . . .(numbers) = F
     Background fill . . . . . . . . . . . . (background) = none      
     Subsampling for velocity model display  . (isubsamp) = 3
     Color display . . . . . . . . . . . . . . . .(color) = T
     Plot boundaries . . . . . . . . . . . . (boundaries) = T
     Plot symbols  . . . . . . . . . . . . . . .(symbols) = T
     Write legends . . . . . . . . . . . . . . . (legend) = T
     Plot vector fields  . . . . . . . . . . . .(vectors) = F
     Amplitude-Scaling . . . . . . . . . . . (ScaleField) = 0.00
     Interpolate vector field  . . . . . . . . (interpol) = T
     Points per edge for interpolation . . . (DisplayPts) = 3
 
***********************************************
*   I n i t i a l i z a t i o n   p h a s e   *
***********************************************
 
 
     Defining the FEM mesh ...... [OK]
     Saving node coordinates in file MeshNodesCoord_sem2d.tab ...... [OK]
     Saving element connectivity in file ElmtNodes_sem2d.tab ...... [OK]
 
 S p e c t r a l   e l e m e n t s   g r i d
 ===========================================
 
     Numbering GLL points ...... [OK]
     Total number of GLL points. . . . . . . = 51681
 
     Saving element/node table in binary file ibool_sem2d.dat ...... [OK]
     Defining nodes coordinates ...... [OK]
 
     Saving the grid coordinates (coord) in a text file ...... [OK]
     Saving the grid coordinates (coord) in a binary file ...... [OK]
 
 M a t e r i a l   p r o p e r t i e s
 =====================================
 
     Translating input velocity model ...... [OK]
     Defining elasticity work arrays ...... [OK]
     Exporting model ...... [OK]
 
 M e s h   p r o p e r t i e s
 =============================
 
     Checking mesh ...... [OK]
    Max mesh size =   18.156E+00
    Min mesh size =    5.012E+00
    Ratio max/min =    3.622E+00
 
    RESOLUTION: nodes per min wavelength =    5.912E+00
 
     Dump PostScript Resolution_sem2d.ps ...... [OK]
     Dump PostScript Stability_sem2d.ps ...... [OK]
 
 T i m e   s o l v e r
 =====================
 
    Time step (secs)      =  500.000E-06
    Number of time steps  = 3000
    Total duration (secs) =    1.500E+00
    Courant number        =  319.227E-03
 
    STABILITY:  CFL number               =  319.227E-03
     Initializing kinematic fields ...... [OK]
       Max displ =    0.000E+00
       Max veloc =    0.000E+00
     Building the mass matrix ...... [OK]
     Defining boundary conditions ...... [OK]
     Initializing receivers ...

 R e c e i v e r s
 =================

 Receivers have been relocated to the nearest GLL node

 Receiver  x-requested  z-requested   x-obtained   z-obtained   distance

        1    2.200E+03    0.000E+00    2.200E+03    0.000E+00    0.000E+00
        2    2.700E+03    0.000E+00    2.700E+03    0.000E+00    0.000E+00

  Maximum distance between asked and real =   0.000E+00
 
  Sampling rate (Hz)        =    2.000E+03
  Sampling timestep (secs)  =  500.000E-06
  Total number of samples   = 3001
  Number of receivers       = 2
 
... [OK]
     Initializing sources ...

 S o u r c e s
 =============

 Sources have been relocated to the nearest GLL node

   Source  x-requested  z-requested   x-obtained   z-obtained     distance

        1    1.500E+03  -50.000E+00    1.500E+03  -50.000E+00    0.000E+00

  Maximum distance between requested and real =   0.000E+00
... [OK]
Timestep #       0  t =   0.000E+00  vmax =   0.000E+00  dmax =   0.000E+00

Snapshot at timestep = 0
     Dump PostScript va_000_sem2d.ps ...... [OK]
 
 
***********************************************
*           S o l v e r   p h a s e           *
***********************************************
 

---  CPU TIME ESTIMATES (in seconds) :
  CPU time for initialization . . 868.053E-03
  CPU time per timestep . . . . .  40.003E-03
  Total solver CPU time . . . . . 120.009E+00
                   (mins) . . . .   2.000E+00
                   (hours). . . .  33.336E-03

Timestep #     200  t = 100.000E-03  vmax = 886.553E-12  dmax =  13.701E-12
Timestep #     400  t = 200.000E-03  vmax =   1.643E-09  dmax =  24.377E-12
Timestep #     600  t = 300.000E-03  vmax =   1.338E-09  dmax =  19.253E-12
Timestep #     800  t = 400.000E-03  vmax =   1.318E-09  dmax =  17.920E-12
Timestep #    1000  t = 500.000E-03  vmax =   1.405E-09  dmax =  17.545E-12
Timestep #    1200  t = 600.000E-03  vmax =   1.406E-09  dmax =  18.063E-12
Timestep #    1400  t = 700.000E-03  vmax =   1.428E-09  dmax =  18.358E-12
Timestep #    1600  t = 800.000E-03  vmax =   1.448E-09  dmax =  18.333E-12
Timestep #    1800  t = 900.000E-03  vmax =   1.384E-09  dmax =  18.334E-12
Timestep #    2000  t =   1.000E+00  vmax =   1.398E-09  dmax =  18.348E-12
Timestep #    2200  t =   1.100E+00  vmax =   1.444E-09  dmax =  18.633E-12
Timestep #    2400  t =   1.200E+00  vmax =   1.449E-09  dmax =  18.477E-12
Timestep #    2600  t =   1.300E+00  vmax =   1.423E-09  dmax =  18.351E-12
Timestep #    2800  t =   1.400E+00  vmax =   1.349E-09  dmax =  18.187E-12
Timestep #    3000  t =   1.500E+00  vmax =   1.446E-09  dmax =  18.419E-12


---  CPU TIME INFORMATION (in seconds) :
  CPU time for initialization . . 868.053E-03
  CPU time per timestep . . . . .  41.097E-03
  Total solver CPU time . . . . . 123.292E+00
                   (mins) . . . .   2.055E+00
                   (hours). . . .  34.248E-03

 Storing sismos data (SEP format) ...


 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 Program  S P E C F E M :  end
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 Lamb's problem                                    
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 D a t e : 06 - 06 - 2007                                 T i m e  : 17:30:47
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------