[cig-commits] r8430 - in seismo/2D/SPECFEM2D/trunk: MAILLE90 SPECFEM90

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 15:45:39 -0800 (Fri, 07 Dec 2007)
New Revision: 8430

Modified:
   seismo/2D/SPECFEM2D/trunk/MAILLE90/Makefile
   seismo/2D/SPECFEM2D/trunk/SPECFEM90/Makefile
   seismo/2D/SPECFEM2D/trunk/SPECFEM90/compute_gradient_attenuation.f90
   seismo/2D/SPECFEM2D/trunk/SPECFEM90/specfem2D.f90
Log:
added support for g95 compiler in Makefile


Modified: seismo/2D/SPECFEM2D/trunk/MAILLE90/Makefile
===================================================================
--- seismo/2D/SPECFEM2D/trunk/MAILLE90/Makefile	2004-09-11 12:01:21 UTC (rev 8429)
+++ seismo/2D/SPECFEM2D/trunk/MAILLE90/Makefile	2007-12-07 23:45:39 UTC (rev 8430)
@@ -16,6 +16,12 @@
 #FLAGS=-O0 -implicitnone -warn stderrors -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -check bounds -C
 FLAGS=-fast -implicitnone -warn stderrors -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -check nobounds
 
+#
+# g95 (free f95 compiler from http://www.g95.org, still under development, but works)
+#
+#F90 = g95
+#FLAGS = -O
+
 # Dec Alpha
 #F90 = f90
 #FLAGS=-O0 -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -check nounderflow -check bounds -C

Modified: seismo/2D/SPECFEM2D/trunk/SPECFEM90/Makefile
===================================================================
--- seismo/2D/SPECFEM2D/trunk/SPECFEM90/Makefile	2004-09-11 12:01:21 UTC (rev 8429)
+++ seismo/2D/SPECFEM2D/trunk/SPECFEM90/Makefile	2007-12-07 23:45:39 UTC (rev 8430)
@@ -16,6 +16,12 @@
 #FLAGS=-O0 -implicitnone -warn stderrors -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -check bounds -C
 FLAGS=-fast -implicitnone -warn stderrors -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -check nobounds
 
+#
+# g95 (free f95 compiler from http://www.g95.org, still under development, but works)
+#
+#F90 = g95
+#FLAGS = -O
+
 # Dec Alpha
 #F90 = f90
 #FLAGS=-O0 -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -check nounderflow -check bounds -C

Modified: seismo/2D/SPECFEM2D/trunk/SPECFEM90/compute_gradient_attenuation.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/SPECFEM90/compute_gradient_attenuation.f90	2004-09-11 12:01:21 UTC (rev 8429)
+++ seismo/2D/SPECFEM2D/trunk/SPECFEM90/compute_gradient_attenuation.f90	2007-12-07 23:45:39 UTC (rev 8430)
@@ -11,7 +11,8 @@
 !
 !========================================================================
 
-  subroutine compute_gradient_attenuation(displ,duxdxl,duzdxl,duxdzl,duzdzl,xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,NSPEC,npoin)
+  subroutine compute_gradient_attenuation(displ,duxdxl,duzdxl,duxdzl,duzdzl, &
+         xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,NSPEC,npoin)
 
 ! compute Grad(displ) for attenuation
 

Modified: seismo/2D/SPECFEM2D/trunk/SPECFEM90/specfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/SPECFEM90/specfem2D.f90	2004-09-11 12:01:21 UTC (rev 8429)
+++ seismo/2D/SPECFEM2D/trunk/SPECFEM90/specfem2D.f90	2007-12-07 23:45:39 UTC (rev 8430)
@@ -736,9 +736,11 @@
 
 ! add the memory variables using the relaxed parameters (Carcione page 111)
 ! beware: there is a bug in Carcione's equation for sigma_zz
-    sigma_xx = sigma_xx + (lambdal_relaxed + mul_relaxed)* (e1_mech1(i,j,k) + e1_mech2(i,j,k)) + TWO * mul_relaxed * (e11_mech1(i,j,k) + e11_mech2(i,j,k))
+    sigma_xx = sigma_xx + (lambdal_relaxed + mul_relaxed)* &
+      (e1_mech1(i,j,k) + e1_mech2(i,j,k)) + TWO * mul_relaxed * (e11_mech1(i,j,k) + e11_mech2(i,j,k))
     sigma_xz = sigma_xz + mul_relaxed * (e13_mech1(i,j,k) + e13_mech2(i,j,k))
-    sigma_zz = sigma_zz + (lambdal_relaxed + mul_relaxed)* (e1_mech1(i,j,k) + e1_mech2(i,j,k)) - TWO * mul_relaxed * (e11_mech1(i,j,k) + e11_mech2(i,j,k))
+    sigma_zz = sigma_zz + (lambdal_relaxed + mul_relaxed)* &
+      (e1_mech1(i,j,k) + e1_mech2(i,j,k)) - TWO * mul_relaxed * (e11_mech1(i,j,k) + e11_mech2(i,j,k))
 
   else