Changeset 110

Timestamp:

04/28/09 15:42:49 (16 years ago)

Author:

Eric.Larour

Message:

Had wiped off viscosity overshoot from helene

File:

• issm/trunk/src/c/ModelProcessorx/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp (modified) (3 diffs)

issm/trunk/src/c/ModelProcessorx/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp

@@ -66 +66 @@
         double tria_meanvel;/*!scaling ratio for velocities*/
         double tria_epsvel; /*!minimum velocity to avoid infinite velocity ratios*/
+        double tria_viscosity_overshoot; 
 
         /*matice constructor input: */
@@ -274 +275 @@
                         tria_epsvel=model->epsvel;
 
+                        /*viscosity_overshoot*/
+                        tria_viscosity_overshoot=model->viscosity_overshoot;
+
                         /*Create tria element using its constructor:*/
-                        tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel);
+                        tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot);
 
                         /*Add tria element to elements dataset: */
@@ -571 +575 @@
                 node_x[1]=model->y[i];
                 node_x[2]=model->z[i];
+
                 
                 node_onbed=(int)model->gridonbed[i];
-                node_onsurface=(int)model->gridonsurface[i];
-                
+                node_onsurface=(int)model->gridonsurface[i];    
                 if (strcmp(model->meshtype,"3d")==0){
                         if (isnan(model->uppernodes[i])){