Changeset 19008
- Timestamp:
- 01/13/15 15:24:20 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/shared/Elements/ComputeDelta18oTemperaturePrecipitation.cpp
r18950 r19008 17 17 IssmDouble glacialindex; // used to vary present day temperature 18 18 19 // CURENTLY THE TEMPERATURE COMPUTED HERE ARE NOT USED.20 // Tdiff is used isnstead in PddSurfaceMassBalance.cpp to get the interpolated temperature21 22 19 glacialindex = (Delta18oTime-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceTime-Delta18oSurfacePresent)) 23 20 /(Delta18oLgm-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceLgm-Delta18oSurfacePresent)); // Tarasov 2004 paper … … 25 22 for (int imonth = 0; imonth<12; imonth++){ 26 23 monthlytemperaturestmp[imonth] = glacialindex*TemperaturesLgm[imonth] + (1.-glacialindex)*TemperaturesPresentday[imonth]; 27 //monthlyprectmp[imonth] = 1.5*pow(2,((monthlytemperaturestmp[imonth]-273.15-0)/10)); //equation from rob's paper28 24 monthlyprectmp[imonth] = PrecipitationsPresentday[imonth]; 29 25 … … 32 28 *(monthlyprecout+imonth) = monthlyprectmp[imonth]; 33 29 } 34 // printf(" tempera %f\n",monthlytemperaturestmp[1]);35 30 }
Note:
See TracChangeset
for help on using the changeset viewer.