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

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-08-21 15:41:31 -0700 (Tue, 21 Aug 2007)
New Revision: 7861

Modified:
   mc/3D/CitcomS/trunk/lib/Pan_problem_misc_functions.c
   mc/3D/CitcomS/trunk/lib/Phase_change.c
Log:
Disable phase change in compressible flow

Modified: mc/3D/CitcomS/trunk/lib/Pan_problem_misc_functions.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Pan_problem_misc_functions.c	2007-08-21 22:03:46 UTC (rev 7860)
+++ mc/3D/CitcomS/trunk/lib/Pan_problem_misc_functions.c	2007-08-21 22:41:31 UTC (rev 7861)
@@ -141,6 +141,7 @@
     double temp,temp2;
     void remove_horiz_ave2(struct All_variables*, double**);
 
+    /* Rayleigh number */
     temp = E->control.Atemp;
 
     for(m=1;m<=E->sphere.caps_per_proc;m++)
@@ -157,6 +158,7 @@
     if(E->control.tracer &&
        (E->composition.ichemical_buoyancy)) {
         for(j=0;j<E->composition.ncomp;j++) {
+            /* TODO: how to scale chemical buoyancy wrt reference density? */
             temp2 = E->composition.buoyancy_ratio[j] * temp;
             for(m=1;m<=E->sphere.caps_per_proc;m++)
                 for(i=1;i<=E->lmesh.nno;i++)
@@ -164,6 +166,7 @@
         }
     }
 
+    /* TODO: phase change function doesn't know about reference density */
     phase_change_apply_410(E, buoy);
     phase_change_apply_670(E, buoy);
     phase_change_apply_cmb(E, buoy);

Modified: mc/3D/CitcomS/trunk/lib/Phase_change.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Phase_change.c	2007-08-21 22:03:46 UTC (rev 7860)
+++ mc/3D/CitcomS/trunk/lib/Phase_change.c	2007-08-21 22:41:31 UTC (rev 7861)
@@ -164,6 +164,12 @@
   pt5 = 0.5;
   one = 1.0;
 
+  /* TODO: compute pressure with reference density */
+  /* disable phase change in compressible flow for now */
+  if(E->control.inv_gruneisen != 0) {
+      myerror(E, "Error: phase change is not implemented in compressible flow.");
+  }
+
   for(m=1;m<=E->sphere.caps_per_proc;m++)     {
     for(i=1;i<=E->lmesh.nno;i++)  {
       e_pressure = (E->sphere.ro-E->sx[m][3][i]) - depth