source: issm/oecreview/Archive/20545-21336/ISSM-20981-20982.diff@ 21337

Index: ../trunk-jpl/src/c/cores/movingfront_core.cpp
===================================================================
--- ../trunk-jpl/src/c/cores/movingfront_core.cpp       (revision 20981)
+++ ../trunk-jpl/src/c/cores/movingfront_core.cpp       (revision 20982)
@@ -13,7 +13,7 @@
 
        /* intermediaries */
        bool save_results,isstressbalance,ismasstransport,isthermal,isenthalpy,islevelset,ismovingfront;
-       int domaintype, num_extrapol_vars, index;
+       int  domaintype, num_extrapol_vars, index,reinit_frequency,step;
        int* extrapol_vars=NULL;
        Analysis  *analysis=NULL;
 
@@ -24,13 +24,16 @@
        femmodel->parameters->FindParam(&ismasstransport,TransientIsmasstransportEnum);
        femmodel->parameters->FindParam(&isthermal,TransientIsthermalEnum);
        femmodel->parameters->FindParam(&ismovingfront,TransientIsmovingfrontEnum);
+       femmodel->parameters->FindParam(&reinit_frequency,LevelsetReinitFrequencyEnum);
+       femmodel->parameters->FindParam(&step,StepEnum);
        if(isthermal && domaintype==Domain3DEnum) femmodel->parameters->FindParam(&isenthalpy,ThermalIsenthalpyEnum);
 
        if(!ismovingfront) return;
 
        /* start the work from here */
        Calvingx(femmodel);
-       if(VerboseSolution()) _printf0_("   computing levelset transport\n");
+       if(VerboseSolution()) _printf0_("   computing level set transport\n");
+
        /* smoothen slope of lsf for computation of normal on ice domain*/
        levelsetfunctionslope_core(femmodel);
 
@@ -66,6 +69,12 @@
        analysis->Core(femmodel);
        delete analysis;
 
+       /*Reset levelset if needed*/
+       if(reinit_frequency && (step%reinit_frequency==0)){
+               if(VerboseSolution()) _printf0_("   reinitializing level set\n");
+               femmodel->ResetLevelset();
+       }
+
        /* update vertices included for next calculation */
        GetMaskOfIceVerticesLSMx(femmodel);