[cig-commits] r19562 - seismo/2D/SPECFEM2D/trunk/src/specfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-02-01 14:50:01 -0800 (Wed, 01 Feb 2012)
New Revision: 19562

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
Log:
added a warning about materials with a negative Poisson's ratio


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90	2012-02-01 21:57:23 UTC (rev 19561)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90	2012-02-01 22:50:01 UTC (rev 19562)
@@ -281,6 +281,14 @@
            ! material can be acoustic (fluid) or elastic (solid)
            if(poroelastcoef(2,1,n) > TINYVAL) then    ! elastic
               write(IOUT,200) n,cp,cs,density(1),poisson,lambda,mu,kappa,young,QKappa,Qmu
+              if(poisson < 0.d0) then
+                write(IOUT,*)
+                write(IOUT,*) 'Materials with a negative Poisson''s ratio can exist,'
+                write(IOUT,*) 'see e.g. R. Lakes, "Science" vol. 235, p. 1038-1040 (1987),'
+                write(IOUT,*) 'but are extremely rare.'
+                write(IOUT,*) 'Hope you know what you are doing...'
+                write(IOUT,*)
+              endif
            else                                       ! acoustic
               write(IOUT,300) n,cp,density(1),kappa,QKappa,Qmu
            endif
@@ -290,6 +298,14 @@
            ! material is poroelastic (solid/fluid)
            write(iout,500) n,sqrt(cpIsquare),sqrt(cpIIsquare),sqrt(cssquare)
            write(iout,600) density(1),poisson_s,lambda_s,mu_s,kappa_s,young_s
+           if(poisson_s < 0.d0) then
+             write(IOUT,*)
+             write(IOUT,*) 'Materials with a negative Poisson''s ratio can exist,'
+             write(IOUT,*) 'see e.g. R. Lakes, "Science" vol. 235, p. 1038-1040 (1987),'
+             write(IOUT,*) 'but are extremely rare.'
+             write(IOUT,*) 'Hope you know what you are doing...'
+             write(IOUT,*) 
+           endif
            write(iout,700) density(2),kappa_f,eta_f
            write(iout,800) lambda_fr,mu_fr,kappa_fr,porosity_array(n),tortuosity_array(n),&
                 permeability(1,n),permeability(2,n),permeability(3,n),Qmu
@@ -313,11 +329,11 @@
        'P-wave velocity. . . . . . . . . . . (cp) =',1pe15.8,/5x, &
        'S-wave velocity. . . . . . . . . . . (cs) =',1pe15.8,/5x, &
        'Mass density. . . . . . . . . . (density) =',1pe15.8,/5x, &
-       'Poisson''s ratio . . . . . . . .(poisson) =',1pe15.8,/5x, &
+       'Poisson''s ratio. . . . . . . . .(poisson) =',1pe15.8,/5x, &
        'First Lame parameter Lambda. . . (lambda) =',1pe15.8,/5x, &
        'Second Lame parameter Mu. . . . . . .(mu) =',1pe15.8,/5x, &
        'Bulk modulus Kappa . . . . . . . .(kappa) =',1pe15.8,/5x, &
-       'Young''s modulus E . . . . . . . .(young) =',1pe15.8,/5x, &
+       'Young''s modulus E. . . . . . . . .(young) =',1pe15.8,/5x, &
        'QKappa_attenuation .  . . . . . .(QKappa) =',1pe15.8,/5x, &
        'Qmu_attenuation . . . . . . . . . . (Qmu) =',1pe15.8)
 
@@ -354,11 +370,11 @@
 600 format(//5x,'-------------------------------',/5x, &
        '-- Solid phase properties --',/5x, &
        'Mass density. . . . . . . . . . (density_s) =',1pe15.8,/5x, &
-       'Poisson''s ratio. . . . . . . . (poisson_s) =',1pe15.8,/5x, &
+       'Poisson''s ratio . . . . . . . . (poisson_s) =',1pe15.8,/5x, &
        'First Lame parameter Lambda. . . (lambda_s) =',1pe15.8,/5x, &
        'Second Lame parameter Mu. . . . . . .(mu_s) =',1pe15.8,/5x, &
        'Solid bulk modulus Kappa . . . . .(kappa_s) =',1pe15.8,/5x, &
-       'Young''s modulus E. . . . . . .  .(young_s) =',1pe15.8)
+       'Young''s modulus E . . . . . . .  .(young_s) =',1pe15.8)
 
 700 format(//5x,'-------------------------------',/5x, &
        '-- Fluid phase properties --',/5x, &
@@ -385,7 +401,7 @@
          'H. . . . . . . . . .=',1pe15.8,/5x, &
          'C. . . . . . . . . .=',1pe15.8,/5x, &
          'M. . . . . . . . . .=',1pe15.8,/5x, &
-         'characteristic freq =',1pe15.8)
+         'Characteristic freq =',1pe15.8)
 
   end subroutine gmat01
 

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90	2012-02-01 21:57:23 UTC (rev 19561)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90	2012-02-01 22:50:01 UTC (rev 19562)
@@ -300,7 +300,7 @@
   'Add Bielak conditions . . . .(add_Bielak_conditions) = ',l6/5x, &
   'Assign external model . . . .(assign_external_model) = ',l6/5x, &
   'Read external SEP file . . .(READ_EXTERNAL_SEP_FILE) = ',l6/5x, &
-  'Turn attenuation on or off. . .(ATTENUATION_VISCOELASTIC_SOLID) = ',l6/5x, &
+  'Attenuation on/off .(ATTENUATION_VISCOELASTIC_SOLID) = ',l6/5x, &
   'Save grid in external file or not. . . (output_grid) = ',l6/5x, &
   'Save a file with total energy or not.(output_energy) = ',l6)