Changeset 17357


Ignore:
Timestamp:
02/26/14 12:12:12 (11 years ago)
Author:
bdef
Message:

CHG: minor cleanup

File:
1 edited

Legend:

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

    r17356 r17357  
    3434
    3535        bool       sedconverged,eplconverged,hydroconverged;
    36         bool       transfered;
    3736        bool       isefficientlayer;
    3837        int        constraints_converged;
     
    4039        int        sedcount,eplcount,hydrocount;
    4140        int        hydro_maxiter;
    42         IssmDouble sediment_kmax,time;
     41        IssmDouble sediment_kmax;
    4342        IssmDouble eps_hyd;
    4443        IssmDouble ndu_sed,nu_sed;
    4544        IssmDouble ndu_epl,nu_epl;
    46         IssmDouble EplConv;
    4745
    4846        /*Recover parameters: */
     
    5149        femmodel->parameters->FindParam(&hydro_maxiter,HydrologydcMaxIterEnum);
    5250        femmodel->parameters->FindParam(&eps_hyd,HydrologydcRelTolEnum);
    53         femmodel->parameters->FindParam(&time,TimeEnum);
    5451
    5552        hydrocount=1;
     
    7471        for(;;){
    7572
    76                 EplConv=1.0;
    7773                sedcount=1;
    7874                eplcount=1;
     
    169165                        InputUpdateFromConstantx(femmodel,false,ConvergedEnum);
    170166                        femmodel->parameters->SetParam(HydrologyEfficientEnum,HydrologyLayerEnum);
    171                         EplConv=1.0;
    172 
    173167
    174168                        for(;;){
     
    220214                                if (ndu_epl==0.0 && nu_epl==0.0) nu_epl=1.0e-6; /*Hacking the case where the EPL is used but empty*/
    221215                                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;
    223216                                if((ndu_epl/nu_epl)<eps_hyd*10.) eplconverged=true;
    224217                                if (eplcount>=hydro_maxiter){
    225218                                        _error_("   maximum number of EPL iterations (" << hydro_maxiter << ") exceeded");
    226219                                }
    227                                
    228220                                eplcount++;
    229221                               
Note: See TracChangeset for help on using the changeset viewer.