Changeset 12385

Timestamp:

06/07/12 10:50:31 (13 years ago)

Author:

lemorzad

Message:

debuging pdd

Location:

issm/trunk/src/c/objects/Elements

Files:

• Penta.cpp (modified) (2 diffs)
• Tria.cpp (modified) (2 diffs)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -2306 +2306 @@
    double pddtj[NUMVERTICES], hmx2;
 
+   //printf("ok1p\n");
+
    /*Recover monthly temperatures and precipitation*/
    Input*     input=inputs->GetInput(SurfaceforcingsMonthlytemperaturesEnum); _assert_(input);
@@ -2322 +2324 @@
      }
    } 
-   /*Recover info at the vertices: */
+   /*Recover info a the vertices: */
    GetInputListOnVertices(&h[0],ThicknessEnum);
    GetInputListOnVertices(&s[0],SurfaceEnum);
+
+   //printf("ok2p\n");
 
    /*Get material parameters :*/

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -2136 +2136 @@
    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 +2170 @@
    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;