sem2dsolve.out

Spectral Element Method for wave propagation and rupture dynamics.

 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 Program  S P E C F E M : start
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 Anisotropic test                                  
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------
 D a t e : 16 - 11 - 2006                                 T i m e  : 12:07:10
 -------------------------------------------------------------------------------
 -------------------------------------------------------------------------------


 E c h o   C o n t r o l   c a r d
 ==================================

     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) = 200
     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   C o n t r o l   c a r d
 ==================================

     Execution mode . . . . . . . . . . . . . . . (iexec) =1
             ==  0     data check only                     
             ==  1     resolution                          
     Number of nodes per edge . . . . . . . . . . .(ngll) = 6
     Number of d.o.f per node . . . . . . . . . . .(ndof) = 2
     Highest frequency to be resolved . . . . . . .(fmax) =  400.000E+03


 I t e r a t i o n   i n f o s 
 =============================

     Time integration scheme . . . . .(kind) = leapfrog  
     Number of time steps. . . . . (NbSteps) = 2000
     Time step increment . . . . . . . .(Dt) =   50.000E-09
     Total simulation duration . (TotalTime) =  100.000E-06


 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 . . . . . . . . . =    2


     ------------------------
     Material set number . . . . . . . . . . . =    1
     -- Transverse anisotropic material --
     c11 coefficient (Pascal). . . . . . (c11) = 1.65000000E+11
     c13 coefficient (Pascal). . . . . . (c13) = 5.00000000E+10
     c33 coefficient (Pascal). . . . . . (c33) = 6.20000000E+10
     c44 coefficient (Pascal). . . . . . (c44) = 3.96000000E+10
     Mass density. . . . . . . . . . . (denst) = 7.10000000E+03
     Velocity of qP along vertical axis. . . . = 2.95506250E+03
     Velocity of qP along horizontal axis. . . = 4.82072988E+03
     Velocity of qSV along vertical axis . . . = 2.36166568E+03
     Velocity of qSV along horizontal axis . . = 2.36166568E+03


     ------------------------
     Material set number . . . . . . . . . . . =    2
     -- Isotropic material --
     P-wave velocity. . . . . . . . . . . (cp) = 4.82072988E+03
     S-wave velocity. . . . . . . . . . . (cs) = 2.36166568E+03
     Mass density. . . . . . . . . . . (denst) = 7.10000000E+03
     Poisson's ratio . . . . . . . . . (poiss) = 3.42105263E-01
     First Lame parameter Lambda. . . . (alam) = 8.58000000E+10
     Second Lame parameter Mu. . . . . . (amu) = 3.96000000E+10
     Bulk modulus K . . . . . . . . . . (Kvol) = 1.12200000E+11
     Young's modulus E . . . . . . . . (young) = 1.06294737E+11


 S o u r c e  F u n c t i o n s
 ============================

     X-position (meters). . . . . . . . . . = -20.000E-03
     Y-position (meters). . . . . . . . . . =   0.000E+00
     Function Type. . . . . . . . . . . . . = Ricker
     Fundamental frequency (Hz) . . . . . . = 170.000E+03
     Time delay (s) . . . . . . . . . . . . =   6.000E-06
     Multiplying factor . . . . . . . . . . =  10.000E+12
     Source Type. . . . . . . . . . . . . . = Collocated Force
     If P-SV: counterclockwise angle / up . =180.00


 R e c e i v e r s   c o n t r o l   c a r d
 ==================================

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


 P l o t   C o n t r o l   c a r d
 ==================================

     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 . . . . . . . . . . . . . . . . . . . = T
       Velocity . . . . . . . . . . . . . . . . . . . . . = F
       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   C o n t r o l   c a r d
 ==================================

     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) = F
     Points per edge for interpolation . . . (DisplayPts) = 3
     Defining the FEM mesh ...... [OK]


 G L L   g r i d
 ===============

     Listing macro-vertex connectivity ...... [OK]
     Numbering GLL points ...
       Processed       100006  nodes
       Processed       200011  nodes
       Processed       300011  nodes
       Processed       400016  nodes
     ... [OK]
 
     Nodes of the global mesh  . . . . = 423801
     Interior points . . . . . . . . . = 270400
     Edge points (without corners) . . = 136240
     Corner points . . . . . . . . . . = 17161
 
     Saving node index table (ibool) in a binary file ...... [OK]
     Defining shape functions ...... [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]
     Defining boundary topologies ...... [OK]
     Translating input velocity model ...... [OK]
     Defining elasticity work arrays ...... [OK]
     Exporting model ...... [OK]
     Checking mesh ...... [OK]
 
===================================================
=          Mesh properties checkboard             =
===================================================
 
    Max mesh size =    1.426E-03
    Min mesh size =  587.362E-06
    Ratio max/min =    2.428E+00