Changeset 19137

Timestamp:

02/20/15 09:27:26 (10 years ago)

Author:

lemorzad

Message:

PDD temperature surface forcing modif in Tria and Penta

Location:

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

Files:

• Penta.cpp (modified) (4 diffs)
• Tria.cpp (modified) (5 diffs)

issm/trunk-jpl/src/c/classes/Elements/Penta.cpp ¶

@@ -2166 +2166 @@
    IssmDouble agd[NUMVERTICES];             // surface mass balance
    IssmDouble monthlytemperatures[NUMVERTICES][12],monthlyprec[NUMVERTICES][12];
+   IssmDouble yearlytemperatures[NUMVERTICES]={0.};
    IssmDouble tmp[NUMVERTICES];
    IssmDouble h[NUMVERTICES],s[NUMVERTICES];
    IssmDouble rho_water,rho_ice,desfac,s0p,s0t,rlaps,rlapslgm;
    IssmDouble PfacTime,TdiffTime,sealevTime;
+   IssmDouble mavg=1./12.; //factor for monthly average
 
    /*Get material parameters :*/
@@ -2194 +2196 @@
        gauss->GaussVertex(iv);
        input->GetInputValue(&monthlytemperatures[iv][month],gauss,time+month/12.*yts);
-       monthlytemperatures[iv][month]=monthlytemperatures[iv][month]-273.15; // conversion from Kelvin to celcius
+       yearlytemperatures[iv]=yearlytemperatures[iv]+monthlytemperatures[iv][month]*mavg; // Has to be in Kelvin 
+       monthlytemperatures[iv][month]=monthlytemperatures[iv][month]-273.15; // conversion from Kelvin to celcius for PDD module
        input2->GetInputValue(&monthlyprec[iv][month],gauss,time+month/12.*yts);
      }
@@ -2241 +2244 @@
 
 
+   this->inputs->AddInput(new PentaInput(TemperatureEnum,&yearlytemperatures[0],P1Enum));
    this->inputs->AddInput(new PentaInput(SurfaceforcingsMassBalanceEnum,&agd[0],P1Enum));
    // this->inputs->AddInput(NewTemperatureInput);
@@ -2246 +2250 @@
    // //this->inputs->AddInput(new PentaVertexInput(ThermalSpcTemperatureEnum,&Tsurf[0]));
     this->InputExtrude(SurfaceforcingsMassBalanceEnum,-1);
+    this->InputExtrude(TemperatureEnum,-1);
    // this->InputExtrude(SurfaceforcingsMonthlytemperaturesEnum,-1);
    // this->InputExtrude(SurfaceforcingsPrecipitationEnum,-1);

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp ¶

@@ -2574 +2574 @@
    IssmDouble agd[NUMVERTICES];             // surface mass balance
    IssmDouble monthlytemperatures[NUMVERTICES][12],monthlyprec[NUMVERTICES][12];
+   IssmDouble yearlytemperatures[NUMVERTICES]={0.};
    IssmDouble tmp[NUMVERTICES];
    IssmDouble h[NUMVERTICES],s[NUMVERTICES];
    IssmDouble rho_water,rho_ice,desfac,s0p,s0t,rlaps,rlapslgm;
    IssmDouble PfacTime,TdiffTime,sealevTime;
+   IssmDouble mavg=1./12.; //factor for monthly average
    
    /*Get material parameters :*/
@@ -2590 +2592 @@
   rlapslgm=matpar->GetMaterialParameter(SurfaceforcingsRlapslgmEnum);
 
-   /*Recover monthly temperatures and precipitation and Present day and LGm ones*/
+   /*Recover monthly temperatures and precipitation and compute the yearly mean temperatures*/
    Input*     input=inputs->GetInput(SurfaceforcingsMonthlytemperaturesEnum); _assert_(input);
    Input*     input2=inputs->GetInput(SurfaceforcingsPrecipitationEnum); _assert_(input2);
@@ -2598 +2600 @@
    this->parameters->FindParam(&yts,ConstantsYtsEnum);
 
-
    for(int month=0;month<12;month++) {
      for(int iv=0;iv<NUMVERTICES;iv++) {
        gauss->GaussVertex(iv);
        input->GetInputValue(&monthlytemperatures[iv][month],gauss,time+month/12.*yts);
-       monthlytemperatures[iv][month]=monthlytemperatures[iv][month]-273.15; // conversion from Kelvin to celcius
+       yearlytemperatures[iv]=yearlytemperatures[iv]+monthlytemperatures[iv][month]*mavg; // Has to be in Kelvin 
+       monthlytemperatures[iv][month]=monthlytemperatures[iv][month]-273.15; // conversion from Kelvin to celcius for PDD module
        input2->GetInputValue(&monthlyprec[iv][month],gauss,time+month/12.*yts);
      }
-   }
+   }     
 
   /*Recover Pfac, Tdiff and sealev at time t:
@@ -2626 +2628 @@
    /*measure the surface mass balance*/
   for (int iv = 0; iv<NUMVERTICES; iv++){
-     agd[iv]=PddSurfaceMassBalance(&monthlytemperatures[iv][0], &monthlyprec[iv][0],
+    agd[iv]=PddSurfaceMassBalance(&monthlytemperatures[iv][0], &monthlyprec[iv][0],
                                   pdds, pds, signorm, yts, h[iv], s[iv],
                                   desfac, s0t, s0p,rlaps,rlapslgm,TdiffTime,sealevTime,
                                   rho_water,rho_ice);
-   }
+  }
 
    /*Update inputs*/    
@@ -2648 +2650 @@
 
 
-   this->inputs->AddInput(new TriaInput(SurfaceforcingsMassBalanceEnum,&agd[0],P1Enum));
+  this->inputs->AddInput(new TriaInput(TemperatureEnum,&yearlytemperatures[0],P1Enum));
+  this->inputs->AddInput(new TriaInput(SurfaceforcingsMassBalanceEnum,&agd[0],P1Enum));
    // this->inputs->AddInput(NewTemperatureInput);
    // this->inputs->AddInput(NewPrecipitationInput);