Changeset 13975 for issm/trunk/src/c/modules/PositiveDegreeDayx/PositiveDegreeDayx.cpp

Timestamp:

11/16/12 08:10:16 (12 years ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 13974

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/c/modules/PositiveDegreeDayx/PositiveDegreeDayx.cpp (modified) (5 diffs)

issm/trunk/src/c/modules/PositiveDegreeDayx/PositiveDegreeDayx.cpp

@@ -35 +35 @@
   IssmDouble PDup, PDCUT = 2.0;    // PDcut: rain/snow cutoff temperature (C)
   IssmDouble tstar; // monthly mean surface temp
-  
+
   IssmDouble* pdds=NULL; 
   IssmDouble* pds=NULL; 
   Element* element = NULL;
-  
+
   pdds=xNew<IssmDouble>(NPDMAX+1); 
   pds=xNew<IssmDouble>(NPDCMAX+1); 
-  
+
   // initialize PDD (creation of a lookup table)
   tstep = 0.1;
@@ -55 +55 @@
 
   itm = reCast<int,IssmDouble>((2*siglim/DT + 1.5));
-  
+
   if (itm >= NPDMAX){
     _printLine_("increase NPDMAX in massBalance.cpp");
@@ -73 +73 @@
   }
   pdds[itm+1] = siglim + DT;
-  
+
   //*********compute PD(T) : snow/precip fraction. precipitation falling as snow
   tstepc = 0.1;
@@ -100 +100 @@
   pds[itm+1] = 0.;
   //     *******END initialize PDD
-  
+
   for(i=0;i<elements->Size();i++){
-    element=(Element*)elements->GetObjectByOffset(i);
+    element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
     element->PositiveDegreeDay(pdds,pds,signorm);
   }
@@ -108 +108 @@
   xDelete<IssmDouble>(pdds);
   xDelete<IssmDouble>(pds);
-  
+
 }