Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 12384)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 12385)
@@ -2306,4 +2306,6 @@
    double pddtj[NUMVERTICES], hmx2;
 
+   //printf("ok1p\n");
+
    /*Recover monthly temperatures and precipitation*/
    Input*     input=inputs->GetInput(SurfaceforcingsMonthlytemperaturesEnum); _assert_(input);
@@ -2322,7 +2324,9 @@
      }
    } 
-   /*Recover info at the vertices: */
+   /*Recover info a the vertices: */
    GetInputListOnVertices(&h[0],ThicknessEnum);
    GetInputListOnVertices(&s[0],SurfaceEnum);
+
+   //printf("ok2p\n");
 
    /*Get material parameters :*/
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 12384)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 12385)
@@ -2136,6 +2136,8 @@
    double fsupT=0.5,  fsupndd=0.5;  // Tsurf mode factors for supice
    double pddtj[NUMVERTICES], hmx2;
-
-   if(!IsOnBed()) return;
+   
+   //printf("ok1t\n");
+   //if(!IsOnBed()) return;
+   //printf("ok2t\n");
 
    /*Recover monthly temperatures and precipitation*/
@@ -2168,6 +2170,6 @@
    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;