Changeset 12866

Timestamp:

08/02/12 08:45:54 (13 years ago)

Author:

lemorzad

Message:

delta18o update

File:

• issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp (modified) (6 diffs)

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

@@ -931 +931 @@
                         input2->GetInputValue(&TemperaturesLgm[iv][month],gauss,month/12.*yts);
                         input3->GetInputValue(&PrecipitationsPresentday[iv][month],gauss,month/12.*yts);
-                        monthlyprec[iv][month]=monthlyprec[iv][month]*yts; // convertion to m/yr
+                        PrecipitationsPresentday[iv][month]=PrecipitationsPresentday[iv][month]/yts; // converion in m/sec
                 }
         }
         /*Recover delta18o and Delta18oSurface at present day, lgm and at time t*/
-        this->parameters->FindParam(&Delta18oPresent,SurfaceforcingsDelta18oEnum,finaltime*yts);
-        this->parameters->FindParam(&Delta18oLgm,SurfaceforcingsDelta18oEnum,(finaltime-21000)*yts);
-        this->parameters->FindParam(&Delta18oTime,SurfaceforcingsDelta18oEnum,time*yts);
-        this->parameters->FindParam(&Delta18oSurfacePresent,SurfaceforcingsDelta18oSurfaceEnum,finaltime*yts);
-        this->parameters->FindParam(&Delta18oSurfaceLgm,SurfaceforcingsDelta18oSurfaceEnum,(finaltime-21000)*yts);
-        this->parameters->FindParam(&Delta18oSurfaceTime,SurfaceforcingsDelta18oSurfaceEnum,time*yts);
-
+        this->parameters->FindParam(&Delta18oPresent,SurfaceforcingsDelta18oEnum,finaltime);
+        this->parameters->FindParam(&Delta18oLgm,SurfaceforcingsDelta18oEnum,(finaltime-(21000*yts)));
+        this->parameters->FindParam(&Delta18oTime,SurfaceforcingsDelta18oEnum,time);
+        this->parameters->FindParam(&Delta18oSurfacePresent,SurfaceforcingsDelta18oSurfaceEnum,finaltime);
+        this->parameters->FindParam(&Delta18oSurfaceLgm,SurfaceforcingsDelta18oSurfaceEnum,(finaltime-(21000*yts)));
+        this->parameters->FindParam(&Delta18oSurfaceTime,SurfaceforcingsDelta18oSurfaceEnum,time);
+        
         /*Compute the temperature and precipitation*/
         for(int iv=0;iv<NUMVERTICES;iv++){
@@ -950 +950 @@
                                         &monthlytemperatures[iv][0], &monthlyprec[iv][0]);
         }
+        //if(id==1)     printf(" monthlytemperatures1 %f\n",monthlytemperatures[iv][0]);
+        //if(id==1)     printf(" monthlyprec1 %f\n",monthlyprec[iv][0]*yts);
         /*Update inputs*/ 
         TransientInput* NewTemperatureInput = new TransientInput(SurfaceforcingsMonthlytemperaturesEnum);
         TransientInput* NewPrecipitationInput = new TransientInput(SurfaceforcingsPrecipitationEnum);
         for (int imonth=0;imonth<12;imonth++) {
-                for(int iv=0;iv<NUMVERTICES;iv++) {
-                        TriaP1Input* newmonthinput1 = new TriaP1Input(SurfaceforcingsMonthlytemperaturesEnum,&monthlytemperatures[iv][imonth]);
-                        NewTemperatureInput->AddTimeInput(newmonthinput1,imonth/12.*yts);
-                        TriaP1Input* newmonthinput2 = new TriaP1Input(SurfaceforcingsPrecipitationEnum,&monthlyprec[iv][imonth]);
-                        NewPrecipitationInput->AddTimeInput(newmonthinput2,imonth/12.*yts);
-                }
+                        TriaP1Input* newmonthinput1 = new TriaP1Input(SurfaceforcingsMonthlytemperaturesEnum,&TemperaturesPresentday[0][imonth]);
+                        //TriaP1Input* newmonthinput1 = new TriaP1Input(SurfaceforcingsMonthlytemperaturesEnum,&monthlytemperatures[0][imonth]);
+                        NewTemperatureInput->AddTimeInput(newmonthinput1,time+imonth/12.*yts);
+                        TriaP1Input* newmonthinput2 = new TriaP1Input(SurfaceforcingsPrecipitationEnum,&monthlyprec[0][imonth]);
+                        NewPrecipitationInput->AddTimeInput(newmonthinput2,time+imonth/12.*yts);
         }
         this->inputs->AddInput(NewTemperatureInput);
@@ -1385 +1386 @@
 
         int    i;
-        Input *input = NULL;
+        Input *input = NULL;    
 
         /*Go through all the input objects, and find the one corresponding to enum_type, if it exists: */
@@ -2154 +2155 @@
    IssmDouble agd[NUMVERTICES];             // surface mass balance
    IssmDouble monthlytemperatures[NUMVERTICES][12],monthlyprec[NUMVERTICES][12];
-   IssmDouble h[NUMVERTICES],s[NUMVERTICES]; // ,b
+   IssmDouble h[NUMVERTICES],s[NUMVERTICES];
    IssmDouble rho_water,rho_ice;
 
@@ -2170 +2171 @@
        monthlytemperatures[iv][month]=monthlytemperatures[iv][month]-273.15; // conversion from Kelvin to celcius
        input2->GetInputValue(&monthlyprec[iv][month],gauss,time+month/12.*yts);
-       monthlyprec[iv][month]=monthlyprec[iv][month]*yts; // convertion to m/yr
+       monthlyprec[iv][month]=monthlyprec[iv][month]*yts; // convertion in m/y
      }
    }
@@ -2181 +2182 @@
   rho_ice=matpar->GetRhoIce();
   rho_water=matpar->GetRhoFreshwater();
+
+  //if(id==1) printf(" monthlytemperatures %f\n",monthlytemperatures[1][0]);
+  //if(id==1) printf("monthlyprec %f\n",monthlyprec[1][0]);
 
    /*measure the surface mass balance*/