[cig-commits] r12697 - seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Thu Aug 21 16:25:16 PDT 2008
Author: dkomati1
Date: 2008-08-21 16:25:16 -0700 (Thu, 21 Aug 2008)
New Revision: 12697
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.F90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
Log:
added an error message in case of a stability problem: print the chunk number, slice number, latitude and longitude at which the problem occurred
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.F90 2008-08-21 21:20:21 UTC (rev 12696)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.F90 2008-08-21 23:25:16 UTC (rev 12697)
@@ -258,8 +258,8 @@
endif
! get the orientation of the seismometer
- thetan=(90.0d0+stdip)*PI/180.0d0
- phin=stazi*PI/180.0d0
+ thetan = (90.0d0+stdip)*PI/180.0d0
+ phin = stazi*PI/180.0d0
! we use the same convention as in Harvard normal modes for the orientation
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90 2008-08-21 21:20:21 UTC (rev 12696)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90 2008-08-21 23:25:16 UTC (rev 12697)
@@ -1766,6 +1766,13 @@
Ufluidnorm = maxval(abs(displ_outer_core))
+ if(Usolidnorm > STABILITY_THRESHOLD .or. Usolidnorm < 0 .or. &
+ Ufluidnorm > STABILITY_THRESHOLD .or. Ufluidnorm < 0) &
+ 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
+ endif
+
! compute the maximum of the maxima for all the slices using an MPI reduction
#ifdef USE_MPI
call MPI_REDUCE(Usolidnorm,Usolidnorm_all,1,CUSTOM_MPI_TYPE,MPI_MAX,0, &
More information about the cig-commits
mailing list