Changeset 27449

Timestamp:

12/01/22 21:18:18 (2 years ago)

Author:

youngmc3

Message:

CHG: fixed RignotMeltParameterization for 3D

File:

• issm/trunk-jpl/src/c/classes/Elements/Element.cpp (modified) (1 diff)

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

@@ -3922 +3922 @@
 
    /* Start looping on the number of vertices: */
-   GaussTria gauss;
+        Gauss* gauss=this->NewGauss();
    for(int iv=0;iv<numvertices;iv++){
-      gauss.GaussVertex(iv);
+                gauss->GaussVertex(iv);
 
       /* Get variables */
-      bed_input->GetInputValue(&bed,&gauss);
-      qsg_input->GetInputValue(&qsg,&gauss);
-                TF_input->GetInputValue(&TF,&gauss);
+      bed_input->GetInputValue(&bed,gauss);
+      qsg_input->GetInputValue(&qsg,gauss);
+                TF_input->GetInputValue(&TF,gauss);
 
       if(basin_icefront_area[basinid]==0.) meltrates[iv]=0.;