Changeset 17357
- Timestamp:
- 02/26/14 12:12:12 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
r17356 r17357 34 34 35 35 bool sedconverged,eplconverged,hydroconverged; 36 bool transfered;37 36 bool isefficientlayer; 38 37 int constraints_converged; … … 40 39 int sedcount,eplcount,hydrocount; 41 40 int hydro_maxiter; 42 IssmDouble sediment_kmax ,time;41 IssmDouble sediment_kmax; 43 42 IssmDouble eps_hyd; 44 43 IssmDouble ndu_sed,nu_sed; 45 44 IssmDouble ndu_epl,nu_epl; 46 IssmDouble EplConv;47 45 48 46 /*Recover parameters: */ … … 51 49 femmodel->parameters->FindParam(&hydro_maxiter,HydrologydcMaxIterEnum); 52 50 femmodel->parameters->FindParam(&eps_hyd,HydrologydcRelTolEnum); 53 femmodel->parameters->FindParam(&time,TimeEnum);54 51 55 52 hydrocount=1; … … 74 71 for(;;){ 75 72 76 EplConv=1.0;77 73 sedcount=1; 78 74 eplcount=1; … … 169 165 InputUpdateFromConstantx(femmodel,false,ConvergedEnum); 170 166 femmodel->parameters->SetParam(HydrologyEfficientEnum,HydrologyLayerEnum); 171 EplConv=1.0;172 173 167 174 168 for(;;){ … … 220 214 if (ndu_epl==0.0 && nu_epl==0.0) nu_epl=1.0e-6; /*Hacking the case where the EPL is used but empty*/ 221 215 if(VerboseConvergence()) _printf0_(setw(50) << left << " Inner EPL Convergence criterion:" << ndu_epl/nu_epl*100 << "%, aiming lower than " << eps_hyd*10*100 << " %\n"); 222 EplConv=ndu_epl/nu_epl;223 216 if((ndu_epl/nu_epl)<eps_hyd*10.) eplconverged=true; 224 217 if (eplcount>=hydro_maxiter){ 225 218 _error_(" maximum number of EPL iterations (" << hydro_maxiter << ") exceeded"); 226 219 } 227 228 220 eplcount++; 229 221
Note:
See TracChangeset
for help on using the changeset viewer.