Changeset 17176

Timestamp:

01/27/14 16:29:23 (11 years ago)

Author:

bdef

Message:

CHG: Removing hydrological relaxation that mess up with the accuracy of the results

Location:

issm/trunk-jpl/src/c

Files:

• analyses/HydrologyDCEfficientAnalysis.cpp (modified) (1 diff)
• classes/Elements/Tria.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp

@@ -334 +334 @@
 
         /*Get previous water head*/
-        basalelement->GetInputListOnNodes(&eplOldHeads[0],EplHeadOldEnum);
+        /* basalelement->GetInputListOnNodes(&eplOldHeads[0],EplHeadEnum); */
         
-        for(i=0;i<numnodes;i++) {
-                eplHeads[i] = eplOldHeads[i]+0.8*(eplHeads[i]-eplOldHeads[i]);
-        }
+        /* for(i=0;i<numnodes;i++) { */
+        /*      eplHeads[i] = eplOldHeads[i]+0.8*(eplHeads[i]-eplOldHeads[i]); */
+        /* } */
         /*Add input to the element: */
         element->AddBasalInput(EplHeadEnum,eplHeads,P1Enum);

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -4812 +4812 @@
                                         
                                         /*Relaxation stuff*/
-                                        wh_trans=preceding_transfer[i]+0.8*(wh_trans-preceding_transfer[i]);
+                                        //                                      wh_trans=preceding_transfer[i]+0.8*(wh_trans-preceding_transfer[i]);
                                         
                                         /*Assign output pointer*/
@@ -4977 +4977 @@
                                         
                                 /*Relaxation stuff*/
-                                thickness[i] = preceding_thickness[i]+0.8*(thickness[i]-preceding_thickness[i]);
+                                //                              thickness[i] = preceding_thickness[i]+0.8*(thickness[i]-preceding_thickness[i]);
                         }
                 }