Changeset 19008


Ignore:
Timestamp:
01/13/15 15:24:20 (10 years ago)
Author:
lemorzad
Message:

update comments on deltaO18 pdd method

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/shared/Elements/ComputeDelta18oTemperaturePrecipitation.cpp

    r18950 r19008  
    1717  IssmDouble glacialindex; // used to vary present day temperature
    1818
    19   // CURENTLY THE TEMPERATURE COMPUTED HERE ARE NOT USED.
    20   // Tdiff is used isnstead in PddSurfaceMassBalance.cpp to get the interpolated temperature
    21 
    2219  glacialindex = (Delta18oTime-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceTime-Delta18oSurfacePresent))
    2320    /(Delta18oLgm-Delta18oPresent-delta18oLapseRate*(Delta18oSurfaceLgm-Delta18oSurfacePresent)); // Tarasov 2004 paper
     
    2522  for (int imonth = 0; imonth<12; imonth++){
    2623    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 paper
    2824    monthlyprectmp[imonth] = PrecipitationsPresentday[imonth];
    2925
     
    3228    *(monthlyprecout+imonth) = monthlyprectmp[imonth];
    3329  }
    34   // printf(" tempera %f\n",monthlytemperaturestmp[1]);
    3530}
Note: See TracChangeset for help on using the changeset viewer.