Changeset 21429

Timestamp:

12/02/16 16:02:08 (8 years ago)

Author:

langchar

Message:

CHG: flipud dz after computing levels, not before

File:

• issm/trunk-jpl/src/m/plot/plotgemb.m (modified) (5 diffs)

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

@@ -20 +20 @@
                         deltat=results(i+1).time-time;
                 else
-                        deltat=time-results(i-1).time;
+                        deltat=time-results(i-1).time;                  
                 end
 
                 %figure out number of levels: 
-                dz=flipud(results(i).SmbDz(element,:)') ;
+                dz=results(i).SmbDz(element,:);
                 if numlevels==-1,
                         nlevels=length(dz); 
@@ -30 +30 @@
                         nlevels=numlevels;
                 end
-                dz=dz(1:nlevels);
+        dz=flipud(dz(1:nlevels)') ;
 
                 %retrieve values: 
@@ -38 +38 @@
                 %build vertical values: 
                 nz=length(dz); 
-                %for j=1:nz, 
-                %       patch([time-deltat/2 ;time+deltat/2; time+deltat/2; time-deltat/2],[z0; z0; z0+dz(j); z0+dz(j)],[T(j);T(j);T(j);T(j)]);
-                %       z0=z0+dz(j);
-                %end
 
                 xA=(time-deltat/2)*ones(nz,1);
@@ -53 +49 @@
                 zD=zeros(nz,1);
                 
+        
                 for j=1:nz, 
                         zA(j)=z0;
@@ -59 +56 @@
                         zD(j)=z0+dz(j);
                         z0=z0+dz(j);
-                end
+        end       
+    
 
-                patch([xA,xB,xC,xD]',[zA,zB,zC,zD]',[T,T,T,T]');
+                patch([xA,xB,xC,xD]',[zA,zB,zC,zD]',[T,T,T,T]','EdgeColor','none');
 
                 if i>=maxstep,
                         break;
                 end
-        end
+    end
+