Changeset 22323 for issm/trunk-jpl/src/m/classes/mesh3dprisms.m

Timestamp:

01/03/18 21:02:33 (7 years ago)

Author:

seroussi

Message:

CHG: adding area correction to compute volumes and areas

File:

• issm/trunk-jpl/src/m/classes/mesh3dprisms.m (modified) (4 diffs)

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

@@ -17 +17 @@
                 long                        = NaN;
                 epsg                        = 0;
+                scale_factor                = NaN;
 
                 vertexonbase                = NaN;
@@ -120 +121 @@
                                 md = checkmessage(md,'md.mesh.z is not consistent with md.geometry.thickness, you changed the geometry after extrusion !!');
                         end
+                        if numel(md.mesh.scale_factor)>1,
+                                md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+                        end
                 end % }}}
                 function disp(self) % {{{
@@ -161 +165 @@
                         fielddisplay(self,'long','vertices longitude [degrees]');
                         fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
+                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
                 end % }}}
                 function marshall(self,prefix,md,fid) % {{{
@@ -181 +186 @@
                         WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofvertices2d','format','Integer');
                         WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofelements2d','format','Integer');
+                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
                 end % }}}
                 function type = domaintype(self) % {{{