Changeset 2723

Timestamp:

12/11/09 09:07:19 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added check for balancedvelocities solution

File:

• issm/trunk/src/m/classes/public/ismodelselfconsistent.m (modified) (1 diff)

issm/trunk/src/m/classes/public/ismodelselfconsistent.m

@@ -327 +327 @@
         if any(md.spcthickness(find(md.gridonboundary))~=1),
                 error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcthickness']);
+        end
+end
+
+%BALANCEDVELOCITIES
+if md.analysis_type==BalancedvelocitiesAnalysisEnum
+
+        %VELOCITIES MELTING AND ACCUMULATION
+        fields={'vx','vy','accumulation','melting'};
+        checksize(md,fields,[md.numberofgrids 1]);
+        checknan(md,fields);
+
+        %SPC
+        if any(md.spcvelocity(find(md.gridonboundary),[1:2])~=1),
+                error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcvelocity']);
         end
 end