source: issm/oecreview/Archive/24307-24683/ISSM-24483-24484.diff@ 24684

../trunk-jpl/src/c/classes/FemModel.cpp

@@ -5338 +5338 @@
                 hmaxvertices_serial=xNew<IssmDouble>(numberofvertices);
                 for(int i=0;i<numberofvertices;i++) hmaxvertices_serial[i]=NAN;
                 /*Fill hmaxvertices*/
+                if(this->amrbamg->thicknesserror_threshold>0)   this->GethmaxVerticesFromEstimators(hmaxvertices_serial,ThicknessErrorEstimatorEnum);
+                if(this->amrbamg->deviatoricerror_threshold>0)  this->GethmaxVerticesFromEstimators(hmaxvertices_serial,DeviatoricStressErrorEstimatorEnum);
                 if(this->amrbamg->groundingline_distance>0)             this->GethmaxVerticesFromZeroLevelSetDistance(hmaxvertices_serial,MaskGroundediceLevelsetEnum);
                 if(this->amrbamg->icefront_distance>0)                          this->GethmaxVerticesFromZeroLevelSetDistance(hmaxvertices_serial,MaskIceLevelsetEnum);
-                if(this->amrbamg->thicknesserror_threshold>0)   this->GethmaxVerticesFromEstimators(hmaxvertices_serial,ThicknessErrorEstimatorEnum);
-                if(this->amrbamg->deviatoricerror_threshold>0)  this->GethmaxVerticesFromEstimators(hmaxvertices_serial,DeviatoricStressErrorEstimatorEnum);
         }
 
         if(my_rank==0){
@@ -5521 +5521 @@
         IssmDouble L1,L2,L3;
         int vid,v1,v2,v3;
         bool refine;
+        
 
         /*Fill variables*/
         switch(errorestimator_type){
@@ -5542 +5543 @@
         }
         if(!error_elements) _error_("error_elements is NULL!\n");
         if(groupthreshold<DBL_EPSILON) _error_("group threshold is too small!");
+        
+        /*Get mesh*/
+        this->GetMesh(&index,&x,&y,&numberofvertices,&numberofelements);
+        if(numberofelements<0) _error_("number of elements is negative!\n");
+        if(numberofvertices<0) _error_("number of vertices is negative!\n");
+        maxlength               = xNew<IssmDouble>(numberofelements);
+        error_vertices  = xNewZeroInit<IssmDouble>(numberofvertices);
 
         /*Find the max of the estimators if it was not provided*/
         if(maxerror<DBL_EPSILON){
@@ -5552 +5560 @@
         }
         }
 
-        /*Get mesh*/
-        this->GetMesh(&index,&x,&y,&numberofvertices,&numberofelements);
-        maxlength               = xNew<IssmDouble>(numberofelements);
-        error_vertices  = xNewZeroInit<IssmDouble>(numberofvertices);
-
         /*Fill error_vertices (this is the sum of all elements connected to the vertex)*/
         for(int i=0;i<numberofelements;i++){
                 v1=index[i*elementswidth+0]-1;//Matlab to C indexing