Changeset 27363

Timestamp:

11/08/22 16:20:30 (2 years ago)

Author:

sehrenfe

Message:

CHG:Fix typo for moving GL

File:

• issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp (modified) (6 diffs)

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

@@ -418 +418 @@
                 b_input->GetInputValue(&b,gauss);
                 H_input->GetInputValue(&H,gauss);
+
+                /*Hard code B*/
+                B = Cuffey(273.15-2);
 
                 /*Get norm of gradient of hydraulic potential and make sure it is >0*/
@@ -523 +526 @@
                 b_input->GetInputValue(&b,gauss);
                 H_input->GetInputValue(&H,gauss);
+                
+                /*Hard code B*/
+                B = Cuffey(273.15-2);
 
                 /*Get basal velocity*/
@@ -674 +680 @@
                 element->GetInputListOnNodes(&mask[0],MaskOceanLevelsetEnum);
                 element->GetInputListOnNodes(&bed[0],BaseEnum);
-                element->GetInputListOnNodes(&thickness[0],BaseEnum);
+                element->GetInputListOnNodes(&thickness[0],ThicknessEnum);
                 element->GetInputListOnNodes(&ls_active[0],HydrologyMaskNodeActivationEnum);
 
@@ -779 +785 @@
                 iceLS_input->GetInputValue(&iceLS,gauss);
 
+                /*Hard code B*/
+                B = Cuffey(273.15-2);
+
                 /*Set sheet thickness to zero if floating or no ice*/
                 if(oceanLS<0. || iceLS>0.){
@@ -812 +821 @@
         }
         }
-
-        /*Force floating ice to have zero sheet thickness*/
-        /*if(!element->IsAllGrounded() || !element->IsIceInElement()){
-                for(int iv=0;iv<numvertices;iv++) h_new[iv] = 0.;
-        }*/
 
         element->AddInput(HydrologySheetThicknessEnum,h_new,P1Enum);
@@ -878 +882 @@
                 p_i = rho_ice*g*H;
 
-                /*Copmute overburden potential*/
+                /*Compute overburden potential*/
                 phi_0 = phi_m + p_i;