Changeset 20795
- Timestamp:
- 06/22/16 21:20:29 (9 years ago)
- Location:
- issm/trunk-jpl/examples/Pig
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/Pig/PigRegion.m
r20730 r20795 8 8 xmin = ncreadatt(nsidc_vel,'/','xmin'); 9 9 ymax = ncreadatt(nsidc_vel,'/','ymax'); 10 spacing = ncreadatt(nsidc_vel,'/','spacing'); 10 spacing = ncreadatt(nsidc_vel,'/','spacing');ls 11 11 12 nx = double(ncreadatt(nsidc_vel,'/','nx')); 12 13 ny = double(ncreadatt(nsidc_vel,'/','ny')); … … 31 32 x_pig=x(posx); 32 33 posy=find(y<=1.0e5 & y>-4.0e5); 33 y_pig= y(posy);34 y_pig=flipud(y(posy)); 34 35 35 36 vx_pig=flipud(vx(posx,posy)'); … … 37 38 vel_pig=sqrt(vx_pig.^2+vy_pig.^2); 38 39 39 imagesc(x_pig,y_pig,vel_pig) 40 imagesc(x_pig,y_pig,log(vel_pig+1)) 41 axis xy equal -
issm/trunk-jpl/examples/Pig/runme.m
r20764 r20795 15 15 %ploting 16 16 plotmodel(md,'data','mesh') 17 return; 17 18 18 19 % Load Velocities … … 144 145 md = loadmodel('./Models/PIG_Control_drag'); 145 146 146 plotmodel(md,' nlines',2,'ncols',2,'unit#all','km','axis#all','equal',...147 plotmodel(md,'axis#all','equal',... 147 148 'xlim#all',[min(md.mesh.x) max(md.mesh.x)]/10^3,... 148 149 'ylim#all',[min(md.mesh.y) max(md.mesh.y)]/10^3,... … … 181 182 surf=find(mdHO.mesh.vertexonsurface); 182 183 183 plotmodel(mdHO,'nlines',3,'ncols',2,' unit#all','km','axis#all','equal',...184 plotmodel(mdHO,'nlines',3,'ncols',2,'axis#all','equal',... 184 185 'xlim#all',[min(mdHO.mesh.x) max(mdHO.mesh.x)]/10^3,... 185 186 'ylim#all',[min(mdHO.mesh.y) max(mdHO.mesh.y)]/10^3,...
Note:
See TracChangeset
for help on using the changeset viewer.