Changeset 27439 for issm/trunk-jpl/src/c/modules/StochasticForcingx/StochasticForcingx.cpp

Timestamp:

11/28/22 09:05:24 (2 years ago)

Author:

vverjans

Message:

CHG moving computation of subglacial water pressure to Element

File:

• issm/trunk-jpl/src/c/modules/StochasticForcingx/StochasticForcingx.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/modules/StochasticForcingx/StochasticForcingx.cpp

@@ -207 +207 @@
                   IssmDouble p_water[numvertices];
                                                 element->GetInputValue(&dimensionid,StochasticForcingDefaultIdEnum);
-                                                Gauss* gauss=element->NewGauss();
-                                                Friction* friction = new Friction(element);
-                                                for(int i=0;i<numvertices;i++){
-                                                        gauss->GaussVertex(i);
-                                                        p_water_deterministic[i] = friction->SubglacialWaterPressure(gauss);
-                                                        p_water[i]               = p_water_deterministic[i] + noisefield[dimensionid]; 
-                                                }
-                                                element->AddInput(FrictionWaterPressureEnum,p_water,P1DGEnum);
-                                                delete gauss;
-                                                delete friction;
+                                                element->SubglacialWaterPressure(FrictionWaterPressureEnum);
+                  element->GetInputListOnVertices(&p_water_deterministic[0],FrictionWaterPressureEnum);
+                  for(int i=0;i<numvertices;i++) p_water[i] = p_water_deterministic[i] + noisefield[dimensionid];
+                  element->AddInput(FrictionWaterPressureEnum,p_water,P1DGEnum);
                                         }
                                         break;
@@ -228 +222 @@
                   IssmDouble p_water[numvertices];
                                                 element->GetInputValue(&dimensionid,StochasticForcingDefaultIdEnum);
-                                                Gauss* gauss=element->NewGauss();
-                                                Friction* friction = new Friction(element);
-                                                for(int i=0;i<numvertices;i++){
-                                                        gauss->GaussVertex(i);
-                                                        p_water_deterministic[i] = friction->SubglacialWaterPressure(gauss);
-                                                        p_water[i]               = p_water_deterministic[i] + noisefield[dimensionid]; 
-                                                }
-                                                element->AddInput(FrictionSchoofWaterPressureEnum,p_water,P1DGEnum);
-                                                delete gauss;
-                                                delete friction;
+                                                element->SubglacialWaterPressure(FrictionSchoofWaterPressureEnum);
+                  element->GetInputListOnVertices(&p_water_deterministic[0],FrictionSchoofWaterPressureEnum);
+                  for(int i=0;i<numvertices;i++) p_water[i] = p_water_deterministic[i] + noisefield[dimensionid];
+                  element->AddInput(FrictionSchoofWaterPressureEnum,p_water,P1DGEnum);
                                         }
                                         break;
@@ -249 +237 @@
                   IssmDouble p_water[numvertices];
                                                 element->GetInputValue(&dimensionid,StochasticForcingDefaultIdEnum);
-                                                Gauss* gauss=element->NewGauss();
-                                                Friction* friction = new Friction(element);
-                                                for(int i=0;i<numvertices;i++){
-                                                        gauss->GaussVertex(i);
-                                                        p_water_deterministic[i] = friction->SubglacialWaterPressure(gauss);
-                                                        p_water[i]               = p_water_deterministic[i] + noisefield[dimensionid]; 
-                                                }
+                                                element->SubglacialWaterPressure(FrictionCoulombWaterPressureEnum);
+                                                element->GetInputListOnVertices(&p_water_deterministic[0],FrictionCoulombWaterPressureEnum);
+                                                for(int i=0;i<numvertices;i++) p_water[i] = p_water_deterministic[i] + noisefield[dimensionid];
                                                 element->AddInput(FrictionCoulombWaterPressureEnum,p_water,P1DGEnum);
-                                                delete gauss;
-                                                delete friction;
                                         }
                                         break;