[cig-commits] r12705 - seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sat Aug 23 15:32:22 PDT 2008
Author: dkomati1
Date: 2008-08-23 15:32:22 -0700 (Sat, 23 Aug 2008)
New Revision: 12705
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/crustal_model.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
Log:
fixed a small bug: the code did not compile because of a "then" statement that I had forgotten
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/crustal_model.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/crustal_model.f90 2008-08-23 21:11:06 UTC (rev 12704)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/crustal_model.f90 2008-08-23 22:32:22 UTC (rev 12705)
@@ -334,7 +334,6 @@
implicit none
include "constants.h"
-
! argument variables
integer ierr
double precision rhtyp(NLAYERS_CRUST),thtp(NLAYERS_CRUST)
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90 2008-08-23 21:11:06 UTC (rev 12704)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90 2008-08-23 22:32:22 UTC (rev 12705)
@@ -1767,7 +1767,7 @@
Ufluidnorm = maxval(abs(displ_outer_core))
if(Usolidnorm > STABILITY_THRESHOLD .or. Usolidnorm < 0 .or. &
- Ufluidnorm > STABILITY_THRESHOLD .or. Ufluidnorm < 0) &
+ Ufluidnorm > STABILITY_THRESHOLD .or. Ufluidnorm < 0) then
print *,'stability problem in slice ',myrank,' in chunk ',ichunk
print *,'around latitude ',(PI/2.0d0 - ystore_crust_mantle(ibool_crust_mantle(1,1,1,1))) * 180.0d0/PI
print *,'and longitude ',zstore_crust_mantle(ibool_crust_mantle(1,1,1,1)) * 180.0d0/PI
More information about the cig-commits
mailing list