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

[gurnis at geodynamics.org]

Author: gurnis
Date: 2008-05-05 12:17:31 -0700 (Mon, 05 May 2008)
New Revision: 11917

Modified:
   mc/3D/CitcomS/trunk/lib/Regional_read_input_from_files.c
Log:
Bug Fix



Modified: mc/3D/CitcomS/trunk/lib/Regional_read_input_from_files.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Regional_read_input_from_files.c	2008-05-05 18:23:53 UTC (rev 11916)
+++ mc/3D/CitcomS/trunk/lib/Regional_read_input_from_files.c	2008-05-05 19:17:31 UTC (rev 11917)
@@ -366,6 +366,8 @@
 
       for(i=1;i<=nnn;i++)   {
          fscanf(fp1,"%f",&(TB1[i]));
+         /* if( E->parallel.me == 0)  
+        fprintf(stderr, "\nINSIDE regional_read_input_files_for_timesteps TB1=%f %d\n",TB1[i],i); */
          if (pos_age) {
              fscanf(fp2,"%f",&(TB2[i]));
          }
@@ -381,7 +383,7 @@
 		if (pos_age) { /* positive ages - we must interpolate */
                     E->sphere.cap[1].TB[1][nodel] = (TB1[nodeg] + (TB2[nodeg]-TB1[nodeg])/(newage2-newage1)*(age-newage1));
                     E->sphere.cap[1].TB[2][nodel] = (TB1[nodeg] + (TB2[nodeg]-TB1[nodeg])/(newage2-newage1)*(age-newage1));
-                    E->sphere.cap[1].TB[2][nodel] = (TB1[nodeg] + (TB2[nodeg]-TB1[nodeg])/(newage2-newage1)*(age-newage1));
+                    E->sphere.cap[1].TB[3][nodel] = (TB1[nodeg] + (TB2[nodeg]-TB1[nodeg])/(newage2-newage1)*(age-newage1));
 		}
 		else { /* negative ages - don't do the interpolation */
                     E->sphere.cap[1].TB[1][nodel] = TB1[nodeg];