Ignore:
Timestamp:
10/08/14 02:06:27 (10 years ago)
Author:
jbondzio
Message:

CHG: use solutionsequence_thermal_nonlinear for enthalpy

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/solutionsequences/convergence.cpp

    r15771 r18589  
    2323        IssmDouble nF;
    2424        IssmDouble solver_residue,res;
     25        int analysis_type;
    2526
    2627        /*convergence options*/
     
    4142        /*get convergence options*/
    4243        parameters->FindParam(&eps_res,StressbalanceRestolEnum);
    43         parameters->FindParam(&eps_rel,StressbalanceReltolEnum);
     44        // get analysis type
     45        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
     46        if(analysis_type==StressbalanceAnalysisEnum)
     47                parameters->FindParam(&eps_rel,StressbalanceReltolEnum);
     48        else if(analysis_type==EnthalpyAnalysisEnum)
     49                parameters->FindParam(&eps_rel,ThermalReltolEnum);
    4450        parameters->FindParam(&eps_abs,StressbalanceAbstolEnum);
    4551        parameters->FindParam(&yts,ConstantsYtsEnum);
Note: See TracChangeset for help on using the changeset viewer.