Changeset 20797

Timestamp:

06/23/16 08:42:52 (9 years ago)

Author:

Mathieu Morlighem

Message:

CHG: minor

File:

• issm/trunk-jpl/examples/Pig/PigRegion.m (modified) (2 diffs)

issm/trunk-jpl/examples/Pig/PigRegion.m

@@ -8 +8 @@
 xmin    = ncreadatt(nsidc_vel,'/','xmin');
 ymax    = ncreadatt(nsidc_vel,'/','ymax');
-spacing = ncreadatt(nsidc_vel,'/','spacing');ls
+spacing = ncreadatt(nsidc_vel,'/','spacing');
 
-nx      = double(ncreadatt(nsidc_vel,'/','nx'));
-ny      = double(ncreadatt(nsidc_vel,'/','ny'));
-vx      = double(ncread(nsidc_vel,'vx'));
-vy      = double(ncread(nsidc_vel,'vy'));
+nx = double(ncreadatt(nsidc_vel,'/','nx'));
+ny = double(ncreadatt(nsidc_vel,'/','ny'));
+vx = double(ncread(nsidc_vel,'vx'));
+vy = double(ncread(nsidc_vel,'vy'));
 
 xmin = strtrim(xmin);  % this is a string, and we need to recover the double value
@@ -29 +29 @@
 
 %Limit the region to Pine Island
-posx=find(x<=-12.0e5 & x>=-18.0e5);
-x_pig=x(posx);
-posy=find(y<=1.0e5 & y>-4.0e5);
-y_pig=flipud(y(posy));
+posx  = find(x<=-12.0e5 & x>=-18.0e5);
+x_pig = x(posx);
+posy  = find(y<=1.0e5 & y>-4.0e5);
+y_pig = flipud(y(posy));
 
-vx_pig=flipud(vx(posx,posy)');
-vy_pig=flipud(vy(posx,posy)');
-vel_pig=sqrt(vx_pig.^2+vy_pig.^2);
+vx_pig  = flipud(vx(posx,posy)');
+vy_pig  = flipud(vy(posx,posy)');
+vel_pig = sqrt(vx_pig.^2+vy_pig.^2);
 
-imagesc(x_pig,y_pig,log(vel_pig+1))
-axis xy equal
+imagesc(x_pig,y_pig,log(vel_pig+1));
+axis xy equal tight;