- Timestamp:
- 11/16/12 08:10:16 (12 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 13397-13398,13401,13407-13582,13584-13974
- Property svn:mergeinfo changed
-
issm/trunk/src
- Property svn:mergeinfo changed
-
issm/trunk/src/c/modules/PositiveDegreeDayx/PositiveDegreeDayx.cpp
r12706 r13975 35 35 IssmDouble PDup, PDCUT = 2.0; // PDcut: rain/snow cutoff temperature (C) 36 36 IssmDouble tstar; // monthly mean surface temp 37 37 38 38 IssmDouble* pdds=NULL; 39 39 IssmDouble* pds=NULL; 40 40 Element* element = NULL; 41 41 42 42 pdds=xNew<IssmDouble>(NPDMAX+1); 43 43 pds=xNew<IssmDouble>(NPDCMAX+1); 44 44 45 45 // initialize PDD (creation of a lookup table) 46 46 tstep = 0.1; … … 55 55 56 56 itm = reCast<int,IssmDouble>((2*siglim/DT + 1.5)); 57 57 58 58 if (itm >= NPDMAX){ 59 59 _printLine_("increase NPDMAX in massBalance.cpp"); … … 73 73 } 74 74 pdds[itm+1] = siglim + DT; 75 75 76 76 //*********compute PD(T) : snow/precip fraction. precipitation falling as snow 77 77 tstepc = 0.1; … … 100 100 pds[itm+1] = 0.; 101 101 // *******END initialize PDD 102 102 103 103 for(i=0;i<elements->Size();i++){ 104 element= (Element*)elements->GetObjectByOffset(i);104 element=dynamic_cast<Element*>(elements->GetObjectByOffset(i)); 105 105 element->PositiveDegreeDay(pdds,pds,signorm); 106 106 } … … 108 108 xDelete<IssmDouble>(pdds); 109 109 xDelete<IssmDouble>(pds); 110 110 111 111 }
Note:
See TracChangeset
for help on using the changeset viewer.