[issm-svn] r12385 - issm/trunk/src/c/objects/Elements

[lemorzad at issm.ess.uci.edu]

Author: lemorzad
Date: 2012-06-07 10:50:31 -0700 (Thu, 07 Jun 2012)
New Revision: 12385

Modified:
   issm/trunk/src/c/objects/Elements/Penta.cpp
   issm/trunk/src/c/objects/Elements/Tria.cpp
Log:
debuging pdd

Modified: issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- issm/trunk/src/c/objects/Elements/Penta.cpp	2012-06-07 17:35:20 UTC (rev 12384)
+++ issm/trunk/src/c/objects/Elements/Penta.cpp	2012-06-07 17:50:31 UTC (rev 12385)
@@ -2305,6 +2305,8 @@
    double fsupT=0.5,  fsupndd=0.5;  // Tsurf mode factors for supice
    double pddtj[NUMVERTICES], hmx2;
 
+   //printf("ok1p\n");
+
    /*Recover monthly temperatures and precipitation*/
    Input*     input=inputs->GetInput(SurfaceforcingsMonthlytemperaturesEnum); _assert_(input);
    Input*     input2=inputs->GetInput(SurfaceforcingsPrecipitationEnum); _assert_(input2);
@@ -2321,10 +2323,12 @@
        monthlyprec[iv][month]=monthlyprec[iv][month]*yts; // convertion to m/yr
      }
    } 
-   /*Recover info at the vertices: */
+   /*Recover info a the vertices: */
    GetInputListOnVertices(&h[0],ThicknessEnum);
    GetInputListOnVertices(&s[0],SurfaceEnum);
 
+   //printf("ok2p\n");
+
    /*Get material parameters :*/
    rho_ice=matpar->GetRhoIce();
    rho_water=matpar->GetRhoFreshwater();

Modified: issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- issm/trunk/src/c/objects/Elements/Tria.cpp	2012-06-07 17:35:20 UTC (rev 12384)
+++ issm/trunk/src/c/objects/Elements/Tria.cpp	2012-06-07 17:50:31 UTC (rev 12385)
@@ -2135,9 +2135,11 @@
    double supice,supcap,diffndd;
    double fsupT=0.5,  fsupndd=0.5;  // Tsurf mode factors for supice
    double pddtj[NUMVERTICES], hmx2;
+   
+   //printf("ok1t\n");
+   //if(!IsOnBed()) return;
+   //printf("ok2t\n");
 
-   if(!IsOnBed()) return;
-
    /*Recover monthly temperatures and precipitation*/
    Input*     input=inputs->GetInput(SurfaceforcingsMonthlytemperaturesEnum); _assert_(input);
    Input*     input2=inputs->GetInput(SurfaceforcingsPrecipitationEnum); _assert_(input2);
@@ -2167,8 +2169,8 @@
      /*PDD constant*/
    siglim = 2.5*signorm; 
    siglimc = 2.5*signormc;
-   siglim0 =  siglim/DT + 0.5;
-   siglim0c =  siglimc/DT + 0.5;
+   siglim0 = siglim/DT + 0.5;
+   siglim0c = siglimc/DT + 0.5;
    PDup = siglimc+PDCUT;
    
    // seasonal loop