[cig-commits] r12714 - mc/2D/ConMan/trunk/src

[sdk at geodynamics.org]

Author: sdk
Date: 2008-08-26 11:30:08 -0700 (Tue, 26 Aug 2008)
New Revision: 12714

Modified:
   mc/2D/ConMan/trunk/src/geoid.F
Log:
essageadded explicit instructions for changing code for Blankenback problem

Modified: mc/2D/ConMan/trunk/src/geoid.F
===================================================================
--- mc/2D/ConMan/trunk/src/geoid.F	2008-08-26 13:33:48 UTC (rev 12713)
+++ mc/2D/ConMan/trunk/src/geoid.F	2008-08-26 18:30:08 UTC (rev 12714)
@@ -26,9 +26,14 @@
       real    geoid1(NWAVE),rtzzt(NWAVE),rtzzb(NWAVE), geoid2(NWAVE)
 c
 c constants for dimensionalization (geoid and topography in meters)
-c  the values used for the Blankenback et al. Benchmark 1a
-c   (Ra=1.0e4, constant viscosity)
 c
+      pi = 3.14159272e0
+      aratio = x(1,numnp)
+c
+c  The values used for the Blankenback et al. Benchmark 1a, 1b, 1c
+c  remove the comments from the following eight lines of code and
+c  comment out the lines below for the Blankenbach benchmarks
+c
 c     d = 1.0e6
 c     g = 1.0e1
 c     rho = 4.0e3
@@ -38,6 +43,14 @@
 c     alfa = 2.5e-5
 c     deltaT = 1.0e3
 c
+c  The values below are more representative of whole mantle convection
+c  values.  Notice that this combination of values assumes a specific
+c  values of the Rayleigh number.   If the input value of the Rayleigh
+c  number differs, then these values will not give meaningful geoid
+c  and topography values.
+c
+c  These lines should be commented out for the Blankenbach problems.
+c
       d = 3.000e6
       g = 1.0e1
       rho = 4.0e3
@@ -46,9 +59,9 @@
       bigG = 6.673e-11
       alfa = 2.0e-5
       deltaT = 2.000e3
-      pi = 3.14159272e0
-      aratio = x(1,numnp)
 c
+c  Below here nothing needs to be changed.
+c
 c if the problem uses reflecting boundary conditions, this is correct
 c  we also need to reflect the solution about the right boundary
 c  to get the correct symmetry.
@@ -197,4 +210,4 @@
 c
 1001  format(4e15.7)
       return
-      end
\ No newline at end of file
+      end