Changeset 12748 for issm/trunk-jpl/src/c/solutions/transient_core.cpp

Timestamp:

07/26/12 10:08:15 (13 years ago)

Author:

lemorzad

Message:

addind delta18o temperature and precipitation methode

File:

• issm/trunk-jpl/src/c/solutions/transient_core.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/solutions/transient_core.cpp

@@ -24 +24 @@
         /*parameters: */
         IssmDouble starttime,finaltime,dt,yts;
-        bool   isdiagnostic,isprognostic,isthermal,isgroundingline,isenthalpy;
+        bool   isdiagnostic,isprognostic,isthermal,isgroundingline,isenthalpy,ispdd,isdelta18o;
         bool   save_results,dakota_analysis;
         bool   time_adapt=false;
@@ -53 +53 @@
         femmodel->parameters->FindParam(&numoutputs,TransientNumRequestedOutputsEnum);
         if(numoutputs) femmodel->parameters->FindParam(&requested_outputs,&numoutputs,TransientRequestedOutputsEnum);
+        femmodel->parameters->FindParam(&ispdd,SurfaceforcingsIspddEnum);
+        femmodel->parameters->FindParam(&isdelta18o,SurfaceforcingsIsdelta18oEnum);
 
         /*initialize: */
@@ -106 +108 @@
                         #endif
                 }
+                
+                if(ispdd){
+                        if(isdelta18o){
+                                if(VerboseSolution()) _pprintLine_("   call Delta18oParametrization module");
+                                Delta18oParameterizationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
+                        } 
+                        if(VerboseSolution()) _pprintLine_("   call positive degree day module");
+                        PositiveDegreeDayx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
+                }
 
                 if(isdiagnostic){