Changeset 19987 for issm/trunk-jpl/src/m/classes/maskpsl.m

Timestamp:

01/22/16 11:01:02 (9 years ago)

Author:

Eric.Larour

Message:

CHG:

• upload legendre coefficients if they are precomputed.
• new land mask for the psl mask class.

File:

• issm/trunk-jpl/src/m/classes/maskpsl.m (modified) (5 diffs)

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

@@ -9 +9 @@
                 ice_levelset         = NaN;
                 ocean_levelset = NaN;
+                land_levelset = NaN;
         end
         methods (Static)
@@ -27 +28 @@
                         self.ice_levelset=project3d(md,'vector',self.ice_levelset,'type','node');
                         self.ocean_levelset=project3d(md,'vector',self.ocean_levelset,'type','node');
+                        self.land_levelset=project3d(md,'vector',self.land_levelset,'type','node');
                 end % }}}
                 function self = mask(varargin) % {{{
@@ -44 +46 @@
                         md = checkfield(md,'fieldname','mask.ice_levelset'        ,'size',[md.mesh.numberofvertices 1]);
                         md = checkfield(md,'fieldname','mask.ocean_levelset','size',[md.mesh.numberofvertices 1]);
+                        md = checkfield(md,'fieldname','mask.land_levelset','size',[md.mesh.numberofvertices 1]);
                         isice=(md.mask.ice_levelset<=0);
                         if sum(isice)==0,
@@ -62 +65 @@
                         fielddisplay(self,'ice_levelset','presence of ice if < 0, icefront position if = 0, no ice if > 0');
                         fielddisplay(self,'ocean_levelset','is the vertex on the ocean ? yes if = 1, no if = 0');
+                        fielddisplay(self,'land_levelset','is the vertex on the land ? yes if = 1, no if = 0');
                 end % }}}
                 function marshall(self,md,fid) % {{{
@@ -67 +71 @@
                         WriteData(fid,'object',self,'class','mask','fieldname','ice_levelset','format','DoubleMat','mattype',1);
                         WriteData(fid,'object',self,'class','mask','fieldname','ocean_levelset','format','DoubleMat','mattype',1);
+                        WriteData(fid,'object',self,'class','mask','fieldname','land_levelset','format','DoubleMat','mattype',1);
 
                         % get mask of vertices of elements with ice