source: issm/oecreview/Archive/16554-17801/ISSM-17221-17222.diff@ 17802

../trunk-jpl/src/c/shared/Enum/EnumDefinitions.h

@@ -678 +678 @@
         /*Levelset related enums (will be moved to appropriate place when finished){{{*/
         TransientIslevelsetEnum,
         ExtrapolationVariableEnum,
+        IceMaskNodeActivationEnum,
         /*}}}*/
         MaximumNumberOfDefinitionsEnum
 };

../trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp

@@ -637 +637 @@
                 case LliboutryDuvalEnum : return "LliboutryDuval";
                 case TransientIslevelsetEnum : return "TransientIslevelset";
                 case ExtrapolationVariableEnum : return "ExtrapolationVariable";
+                case IceMaskNodeActivationEnum : return "IceMaskNodeActivation";
                 case MaximumNumberOfDefinitionsEnum : return "MaximumNumberOfDefinitions";
                 default : return "unknown";
 

../trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp

@@ -652 +652 @@
               else if (strcmp(name,"LliboutryDuval")==0) return LliboutryDuvalEnum;
               else if (strcmp(name,"TransientIslevelset")==0) return TransientIslevelsetEnum;
               else if (strcmp(name,"ExtrapolationVariable")==0) return ExtrapolationVariableEnum;
+              else if (strcmp(name,"IceMaskNodeActivation")==0) return IceMaskNodeActivationEnum;
               else if (strcmp(name,"MaximumNumberOfDefinitions")==0) return MaximumNumberOfDefinitionsEnum;
          else stage=7;
    }

../trunk-jpl/src/m/classes/mask.m

@@ -38 +38 @@
                 function marshall(obj,md,fid) % {{{
                         WriteData(fid,'object',obj,'fieldname','groundedice_levelset','format','DoubleMat','mattype',1);
                         WriteData(fid,'object',obj,'fieldname','ice_levelset','format','DoubleMat','mattype',1);
+
+                        % get mask of vertices of elements with ice
+                        isice=md.mask.ice_levelset<=0.;
+                        vlist = zeros(md.mesh.numberofvertices,1);
+                        pos=find(sum(isice(md.mesh.elements),2)>0);
+                        vlist(md.mesh.elements(pos,:))=1;
+                        WriteData(fid,'data',vlist,'enum',IceMaskNodeActivationEnum(),'format','DoubleMat','mattype',1);
+
                 end % }}}
         end
 end

../trunk-jpl/src/m/classes/mask.py

@@ +1 @@
+import numpy
 from fielddisplay import fielddisplay
 from EnumDefinitions import *
 from checkfield import *
@@ -42 +43 @@
         def marshall(self,md,fid):    # {{{
                 WriteData(fid,'object',self,'fieldname','groundedice_levelset','format','DoubleMat','mattype',1)
                 WriteData(fid,'object',self,'fieldname','ice_levelset','format','DoubleMat','mattype',1)
+
+                # get mask of vertices of elements with ice
+                isice=numpy.array(md.mask.ice_levelset<=0.,int)
+                vlist = numpy.zeros((md.mesh.numberofvertices,1), dtype=int)
+                pos=numpy.nonzero(numpy.sum(isice[md.mesh.elements-1],axis=1))[0]
+                vlist[md.mesh.elements[pos,:]-1]=1
+                WriteData(fid,'data',vlist,'enum',IceMaskNodeActivationEnum(),'format','DoubleMat','mattype',1);
         # }}}

../trunk-jpl/src/m/enum/EnumDefinitions.py

@@ -629 +629 @@
 def LliboutryDuvalEnum(): return StringToEnum("LliboutryDuval")[0]
 def TransientIslevelsetEnum(): return StringToEnum("TransientIslevelset")[0]
 def ExtrapolationVariableEnum(): return StringToEnum("ExtrapolationVariable")[0]
+def IceMaskNodeActivationEnum(): return StringToEnum("IceMaskNodeActivation")[0]
 def MaximumNumberOfDefinitionsEnum(): return StringToEnum("MaximumNumberOfDefinitions")[0]