Changeset 1267

Timestamp:

07/08/09 20:06:22 (15 years ago)

Author:

Eric.Larour

Message:

Added final checks to geography2.m

File:

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

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

@@ -121 +121 @@
 end
 
+
+%some final checks: 
+%check that no grid thinks it's on an ice shelf or ice sheet, and lies actually in the middle of the water.
+gridsgrounded=find(~gridonwater);
+lengthconnectivity=size(md.nodeconnectivity,2);
+groundedcounters=md.nodeconnectivity(gridsgrounded,lengthconnectivity);
+groundedconnectivity=md.nodeconnectivity(gridsgrounded,1:lengthconnectivity-1);
+pos=find(groundedconnectivity);
+groundedconnectivity(pos)=elementonwater(groundedconnectivity(pos));
+groundedsum=sum(groundedconnectivity,2);
+errorflags=find(groundedsum==groundedcounters);
+errorgrids=gridsgrounded(errorflags);
+
+gridonwater(errorgrids)=1;
+gridonicesheet(errorgrids)=0;
+gridoniceshelf(errorgrids)=0;
+
 %Return: 
 md.gridoniceshelf=gridoniceshelf;