Changeset 25587 for issm/branches/trunk-larour-SLPS2020/src/m/plot/processmesh.m

Timestamp:

09/20/20 16:36:58 (4 years ago)

Author:

Eric.Larour

Message:

CHG: diverse

File:

• issm/branches/trunk-larour-SLPS2020/src/m/plot/processmesh.m (modified) (3 diffs)

issm/branches/trunk-larour-SLPS2020/src/m/plot/processmesh.m

@@ -24 +24 @@
 if strcmp(domaintype(md.mesh),'3Dsurface'),
         [x y z elements is2d isplanet] = processmesh(md.mesh,options);
-        if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
-                x=md.mesh.long;
-                y=md.mesh.lat;
+        if strcmpi(getfieldvalue(options,'coord','xy'),'latlon') | strcmpi(getfieldvalue(options,'coord','xy'),'latlong'),
+                x0=md.mesh.long;
+                y0=md.mesh.lat;
+                %add row at lat=90 and lat=-90
+                add=[(-180:.1:-1e-5)';(1e-5:.1:180)'];
+                nadd=length(add);
+                xextra=[add;add];
+                yextra=[90*ones(nadd,1); -90*ones(nadd,1)];
+                x=[x0;xextra];
+                y=[y0;yextra];
+
+                if strcmpi(getfieldvalue(options,'coordcent','atlantic'),'pacific'),
+                        pos=find(x>0);  x(pos)=-360+x(pos);
+                end
                 elements=delaunay(x,y);
-                z=md.mesh.lat; z(:)=0;
+                z=x; z(:)=0;
         end
         return;
@@ -42 +53 @@
 else
         elements=md.mesh.elements;
-        if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon') ,
+        if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon') &  ~strcmpi(getfieldvalue(options,'coord','xy'),'latlong') ,
                 x=md.mesh.x;
                 if isprop(md.mesh,'x2d'), x2d=md.mesh.x2d; end
@@ -50 +61 @@
                 x=md.mesh.long;
                 y=md.mesh.lat;
+                if strcmpi(getfieldvalue(options,'coordcent','atlantic'),'pacific'),
+                        pos=find(x>0);  x(pos)-360+x(pos);
+                end
         end
 end