Changeset 20564
- Timestamp:
- 04/26/16 15:48:16 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/plot/plotgemb.m
r19586 r20564 38 38 %build vertical values: 39 39 nz=length(dz); 40 %for j=1:nz, 41 % 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)]); 42 % z0=z0+dz(j); 43 %end 44 45 xA=(time-deltat/2)*ones(nz,1); 46 xB=(time+deltat/2)*ones(nz,1); 47 xC=(time+deltat/2)*ones(nz,1); 48 xD=(time-deltat/2)*ones(nz,1); 49 50 zA=zeros(nz,1); 51 zB=zeros(nz,1); 52 zC=zeros(nz,1); 53 zD=zeros(nz,1); 54 40 55 for j=1:nz, 41 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)]); 56 zA(j)=z0; 57 zB(j)=z0; 58 zC(j)=z0+dz(j); 59 zD(j)=z0+dz(j); 42 60 z0=z0+dz(j); 43 61 end 62 63 patch([xA,xB,xC,xD]',[zA,zB,zC,zD]',[T,T,T,T]'); 64 44 65 if i>=maxstep, 45 66 break;
Note:
See TracChangeset
for help on using the changeset viewer.