[cig-commits] r14838 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2009-05-01 15:18:00 -0700 (Fri, 01 May 2009)
New Revision: 14838

Modified:
   mc/3D/CitcomS/trunk/lib/Mineral_physics_models.c
Log:
Depth is measured from r=1.0, not outer radius, because in coupled model, the outer radius can be < 1.0.

Modified: mc/3D/CitcomS/trunk/lib/Mineral_physics_models.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Mineral_physics_models.c	2009-05-01 22:17:27 UTC (rev 14837)
+++ mc/3D/CitcomS/trunk/lib/Mineral_physics_models.c	2009-05-01 22:18:00 UTC (rev 14838)
@@ -170,7 +170,7 @@
     double *rhor, *vpr, *vsr, *depthkm;
     double d, d2, dT, dC, drho, dvp, dvs;
 
-    /* compute horizontal average, if not done yet */
+    /* compute horizontal average */
     if(!E->output.horiz_avg)
         compute_horiz_avg(E);
 
@@ -183,7 +183,7 @@
 
     for(nz=1; nz<=E->lmesh.noz; nz++) {
         get_prem(E->sx[m][3][nz], &vpr[nz], &vsr[nz], &rhor[nz]);
-        depthkm[nz] = (E->sphere.ro - E->sx[m][3][nz]) * E->data.radius_km;
+        depthkm[nz] = (1.0 - E->sx[m][3][nz]) * E->data.radius_km;
     }
 
     /* deviation from the reference */
@@ -278,7 +278,7 @@
         break;
 
     case 100:
-        /* user-defined initial temperature goes here */
+        /* user-defined mineral physics model goes here */
         fprintf(stderr,"Need user definition for mineral physics model: 'mineral_physics_model=%d'\n",
                 E->control.mineral_physics_model);
         parallel_process_termination();