Changeset 27296

Timestamp:

10/05/22 10:33:09 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: time stepping may be looking for vxhydro/vyhydro

File:

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

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

@@ -5289 +5289 @@
         /*Check hydro timestep also and take the minimum*/
         if(ishydro){
-                Input* vx_input = this->GetInput(HydrologyWaterVxEnum); _assert_(vx_input);
-                Input* vy_input = this->GetInput(HydrologyWaterVyEnum); _assert_(vy_input);
-                IssmDouble maxabsvx = vx_input->GetInputMaxAbs();
-                IssmDouble maxabsvy = vy_input->GetInputMaxAbs();
+                Input* vx_input = this->GetInput(HydrologyWaterVxEnum);
+                Input* vy_input = this->GetInput(HydrologyWaterVyEnum);
+                IssmDouble maxabsvx = 0.03;
+                IssmDouble maxabsvy = 0.03;
+                if(vx_input) maxabsvx = vx_input->GetInputMaxAbs();
+                if(vy_input) maxabsvy = vy_input->GetInputMaxAbs();
 
                 /*CFL criterion: */