Changeset 23828

Timestamp:

04/07/19 20:24:49 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added check on bed

File:

• issm/trunk-jpl/src/m/classes/geometry.m (modified) (1 diff)

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

@@ -66 +66 @@
                                 if strcmp(solution,'TransientSolution') & md.transient.isgroundingline,
                                         md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+                                        if any(self.bed>self.base),
+                                                md = checkmessage(md,['base<bed on one or more vertex']);
+                                        end 
+                                        pos = find(md.mask.groundedice_levelset>0);
+                                        if any(abs(self.bed(pos)-self.base(pos))>10^-9),
+                                                md = checkmessage(md,['equality base=bed on grounded ice violated']);
+                                        end 
+                                        md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
                                 end
                         end