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

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-01-16 10:43:28 -0800 (Tue, 16 Jan 2007)
New Revision: 5793

Modified:
   mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c
Log:
When "tic_method=1", the temperature is isothermal (== bottom b.c.) with a
uniformly cold plate (thickness specified by 'half_space_age')


Modified: mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c	2007-01-16 17:41:55 UTC (rev 5792)
+++ mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c	2007-01-16 18:43:28 UTC (rev 5793)
@@ -387,6 +387,16 @@
 
   }
   else if (E->convection.tic_method == 1) {
+      /* set up a top thermal boundary layer */
+      for(m=1;m<=E->sphere.caps_per_proc;m++)
+          for(i=1;i<=noy;i++)
+              for(j=1;j<=nox;j++)
+                  for(k=1;k<=noz;k++) {
+                      node=k+(j-1)*noz+(i-1)*nox*noz;
+                      r1=E->sx[m][3][node];
+                      temp = 0.2*(E->sphere.ro-r1) * 0.5/sqrt(E->convection.half_space_age/E->data.scalet);
+                      E->T[m][node] = E->control.TBCbotval*erf(temp);
+                  }
 
   }
   else if (E->convection.tic_method == 2) {