Changeset 27294

Timestamp:

10/05/22 08:54:44 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added CFL for hydro models

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• Penta.cpp (modified) (3 diffs)
• Tria.cpp (modified) (3 diffs)

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

@@ -3901 +3901 @@
 
         /*intermediary: */
+        bool       ishydro;
         IssmDouble C;
         IssmDouble xyz_list[NUMVERTICES][3];
@@ -3906 +3907 @@
         /*get CFL coefficient:*/
         this->parameters->FindParam(&C,TimesteppingCflCoefficientEnum);
+        this->parameters->FindParam(&ishydro,TransientIshydrologyEnum);
 
         /*Get for Vx and Vy, the max of abs value: */
@@ -3939 +3941 @@
         /*CFL criterion: */
         IssmDouble dt = C/(maxabsvx/dx+maxabsvy/dy+maxabsvz/dz + 1.e-18);
+
+        /*Check hydro timestep also and take the minimum*/
+        if(ishydro){
+                if(IsOnBase()){
+                        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();
+
+                        /*CFL criterion: */
+                        IssmDouble dt2 = C/(maxabsvx/dx+maxabsvy/dy + 1.e-18);
+                        if(dt2<dt) dt = dt2;
+                }
+        }
 
         return dt;

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

@@ -5253 +5253 @@
 
         /*intermediary: */
+        bool       ishydro;
         IssmDouble C;
         IssmDouble xyz_list[NUMVERTICES][3];
@@ -5258 +5259 @@
         /*get CFL coefficient:*/
         this->parameters->FindParam(&C,TimesteppingCflCoefficientEnum);
+        this->parameters->FindParam(&ishydro,TransientIshydrologyEnum);
 
         /*Get for Vx and Vy, the max of abs value: */
@@ -5284 +5286 @@
         /*CFL criterion: */
         IssmDouble dt = C/(maxabsvx/dx+maxabsvy/dy + 1.e-18);
+
+        /*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();
+
+                /*CFL criterion: */
+                IssmDouble dt2 = C/(maxabsvx/dx+maxabsvy/dy + 1.e-18);
+                if(dt2<dt) dt = dt2;
+        }
 
         return dt;