Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/shared/Elements/ComputeDelta18oTemperaturePrecipitation.cpp
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/shared/Elements/ComputeDelta18oTemperaturePrecipitation.cpp	(revision 12922)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/shared/Elements/ComputeDelta18oTemperaturePrecipitation.cpp	(revision 12923)
@@ -15,15 +15,9 @@
   IssmDouble delta18oLapseRate=-6.2*pow(10.,-3);
   IssmDouble glacialindex; // used to vary present day temperature
   int imonth=0;
-  // printf("Delta18oTime %f\n",Delta18oTime);
-  // printf("Delta18oPresent %f\n",Delta18oPresent);
-  // printf("Delta18oLgm %f\n",Delta18oLgm);
 
-  glacialindex = (Delta18oTime-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceTime-Delta18oSurfacePresent))
-    /(Delta18oLgm-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceLgm-Delta18oSurfacePresent));
-  // printf(" monthlytemperatureseeeeee %f\n",TemperaturesPresentday[1]);
-  // printf(" monthlypreceeeeeeee %f\n",TemperaturesLgm[1]);
-  // printf(" glacialindex %f\n",glacialindex);
+  glacialindex = 0;//(Delta18oTime-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceTime-Delta18oSurfacePresent))
+  //  /(Delta18oLgm-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceLgm-Delta18oSurfacePresent));
 
   for (int imonth = 0; imonth<12; imonth++){  
     monthlytemperaturestmp[imonth] = glacialindex*TemperaturesLgm[imonth] + (1-glacialindex)*TemperaturesPresentday[imonth];