Changeset 28101

Timestamp:

02/15/24 11:47:28 (14 months ago)

Author:

Mathieu Morlighem

Message:

BUG: check that observations are in chronological order

Location:

issm/trunk-jpl/src

Files:

• c/classes/Cfsurfacesquaretransient.cpp (modified) (1 diff)
• m/classes/cfsurfacesquaretransient.m (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Cfsurfacesquaretransient.cpp

@@ -52 +52 @@
         xMemCpy<bool>(this->passedflags,in_passedflags,this->num_datatimes);
 
+        #ifdef _ISSM_DEBUG_ 
+        for(int i=0;i<this->num_datatimes-1;i++){
+                if(this->datatimes[i+1]<=this->datatimes[i]){
+                        _error_("time series is not in chronological order");
+                }
+        }
+        #endif
+
         this->J = in_J;
 }

issm/trunk-jpl/src/m/classes/cfsurfacesquaretransient.m

@@ -57 +57 @@
                         end
                         md = checkfield(md,'fieldname','self.definitionstring','field',self.definitionstring,'values',OutputdefinitionStringArray);
-                        md = checkfield(md,'fieldname','self.observations','field',self.observations,'size',[md.mesh.numberofvertices+1 NaN],'NaN',1,'Inf',1);
+                        md = checkfield(md,'fieldname','self.observations','field',self.observations,'size',[md.mesh.numberofvertices+1 NaN],'NaN',1,'Inf',1,'timeseries',1);
                         md = checkfield(md,'fieldname','self.weights','field',self.weights,'size',[md.mesh.numberofvertices+1 NaN],'NaN',1,'Inf',1);