Changeset 25621

Timestamp:

10/01/20 13:55:55 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added depthaverage option

Location:

issm/trunk-jpl/src/m/plot

Files:

• colormaps/demmap.m (modified) (1 diff)
• processdata.m (modified) (1 diff)
• processmesh.m (modified) (2 diffs)

issm/trunk-jpl/src/m/plot/colormaps/demmap.m

@@ -12 +12 @@
 %Input checks
 if nargin<3,
-        help landseacolor
+        help demmap
         error('3 or 4 arguments necessary');
 elseif nargin>4,
-        help landseacolor
+        help demmap
         error('3 or 4 arguments necessary');
 end

issm/trunk-jpl/src/m/plot/processdata.m

@@ -205 +205 @@
 
 %layer projection? 
-if getfieldvalue(options,'layer',0)>=1,
+if getfieldvalue(options,'layer',0)>=1
         data=project2d(md,data,getfieldvalue(options,'layer')); %project onto 2d mesh
+end
+if getfieldvalue(options,'depthaverage',0)
+        data=depthaverage(md,data)); %project onto 2d mesh
 end
 

issm/trunk-jpl/src/m/plot/processmesh.m

@@ -64 +64 @@
 
 %is it a 2d plot?
-if md.mesh.dimension()==2 | getfieldvalue(options,'layer',0)>=1,
+if md.mesh.dimension()==2 || getfieldvalue(options,'layer',0)>=1 || getfieldvalue(options,'depthaverage',0)
         is2d=1;
 else
@@ -71 +71 @@
 
 %layer projection? 
-if getfieldvalue(options,'layer',0)>=1,
+if getfieldvalue(options,'layer',0)>=1 || getfieldvalue(options,'depthaverage',0)
         if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
                 error('processmesh error message: cannot work with 3D meshes for now');