Changeset 26458

Timestamp:

09/23/21 10:55:59 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixed basal velocities if 2d model with SSA

File:

• issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp (modified) (1 diff)

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

@@ -144 +144 @@
         iomodel->FetchDataToInput(inputs,elements,"md.initialization.watercolumn",HydrologySheetThicknessEnum);
         iomodel->FetchDataToInput(inputs,elements,"md.initialization.hydraulic_potential",HydraulicPotentialEnum);
-        iomodel->FetchDataToInput(inputs,elements,"md.initialization.vx",VxEnum);
-        iomodel->FetchDataToInput(inputs,elements,"md.initialization.vy",VyEnum);
-        iomodel->FindConstant(&frictionlaw,"md.friction.law");
+        if iomodel->domaintype==Domain2DhorizontalEnum){
+                iomodel->FetchDataToInput(inputs,elements,"md.initialization.vx",VxBaseEnum);
+                iomodel->FetchDataToInput(inputs,elements,"md.initialization.vy",VyBaseEnum);
+                iomodel->FindConstant(&frictionlaw,"md.friction.law");
+        }
+        else{
+                iomodel->FetchDataToInput(inputs,elements,"md.initialization.vx",VxEnum);
+                iomodel->FetchDataToInput(inputs,elements,"md.initialization.vy",VyEnum);
+                iomodel->FindConstant(&frictionlaw,"md.friction.law");
+        }
 
         /*Friction law variables*/