Changeset 16570

Timestamp:

10/28/13 17:04:21 (11 years ago)

Author:

bdef

Message:

Hacking the convergence criterion when the sediment layer is empty

File:

• issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp

@@ -166 +166 @@
                         delete ug_sed_main_iter;
                         if (xIsNan<IssmDouble>(ndu_sed) || xIsNan<IssmDouble>(nu_sed)) _error_("Sed convergence criterion is NaN!");
+                        if (ndu_sed==0.0 && nu_sed==0.0) nu_sed=1.0e-6; /*Hacking the case where the Sediment is used but empty*/
                         if (!xIsNan<IssmDouble>(eps_hyd)){
                                 if (!isefficientlayer){