Changeset 12867

Timestamp:

08/02/12 08:46:28 (13 years ago)

Author:

lemorzad

Message:

delta18o update

Location:

issm/trunk-jpl/src/c/shared/Elements

Files:

• ComputeDelta18oTemperaturePrecipitation.cpp (modified) (1 diff)
• PddSurfaceMassBalance.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/shared/Elements/ComputeDelta18oTemperaturePrecipitation.cpp

@@ -15 +15 @@
   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);
+
   for (int imonth = 0; imonth<12; imonth++){  
-    monthlytemperaturestmp[imonth] = glacialindex*TemperaturesLgm[imonth] + (1-glacialindex)*TemperaturesPresentday[imonth];
-    monthlyprectmp[imonth] = 1.5*pow(2,((monthlytemperaturestmp[imonth]-0)/10)); //equation from rob's paper
+    monthlytemperaturestmp[imonth] = TemperaturesPresentday[imonth];//glacialindex*TemperaturesLgm[imonth] + (1-glacialindex)*TemperaturesPresentday[imonth];
+    //monthlyprectmp[imonth] = 1.5*pow(2,((monthlytemperaturestmp[imonth]-273.15-0)/10)); //equation from rob's paper
+    monthlyprectmp[imonth] = PrecipitationsPresentday[imonth];
+
     /*Assign output pointer*/
     *(monthlytemperaturesout+imonth) = monthlytemperaturestmp[imonth];
     *(monthlyprecout+imonth) = monthlyprectmp[imonth];
   }
+  // printf(" tempera %f\n",monthlytemperaturestmp[1]);
 }

issm/trunk-jpl/src/c/shared/Elements/PddSurfaceMassBalance.cpp

@@ -61 +61 @@
   
   sconv=(rho_water/rho_ice)/12.; //rhow_rain/rhoi / 12 months
-  
+
   /*PDD constant*/
   siglim = 2.5*signorm;