[cig-commits] r21177 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries

[elliott.sales.de.andrade at geodynamics.org]

Author: elliott.sales.de.andrade
Date: 2012-12-18 21:25:07 -0800 (Tue, 18 Dec 2012)
New Revision: 21177

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/combine_AVS_DX.f90
Log:
Print ocean percentage only when required.

The above/below information about oceans is only calculated when crustal
or topography is enabled. Trying to divide them when that's not true
results in a 0/0.



Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/combine_AVS_DX.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/combine_AVS_DX.f90	2012-12-17 00:33:59 UTC (rev 21176)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/auxiliaries/combine_AVS_DX.f90	2012-12-19 05:25:07 UTC (rev 21177)
@@ -1007,13 +1007,15 @@
   if(ivalue == 3) then
     print *
     print *,'elevation min, max = ',minval(elevation_sphere),maxval(elevation_sphere)
-    if(TOPOGRAPHY .or. CRUSTAL) print *,'elevation has been normalized for topography'
-    print *
+    if(TOPOGRAPHY .or. CRUSTAL) then
+      print *,'elevation has been normalized for topography'
+      print *
 
 ! print percentage of oceans at surface of the globe
+      print *
+      print *,'the oceans represent ',100. * below_zero / (above_zero + below_zero),' % of the surface of the mesh'
+    endif
     print *
-    print *,'the oceans represent ',100. * below_zero / (above_zero + below_zero),' % of the surface of the mesh'
-    print *
 
   endif