Documentation of plot_precip

Global Index (short | long) | Local contents | Local Index (short | long)

Cross-Reference Information

Listing of script plot_precip

clear
lims = [-0.1 360 -90 90];
tim = 101:1000;
lags = [0];
varn1= 'rnd';  tit1 = 'Mean Precipitation';  units1 = 'mm day^-^1';
varn2= 'ico';  tit2 = 'Mean Ice Concentration';  units2 = 'Percent';
varn1= 'evp';  tit1 = 'Mean Latent Heat Flux';  units1 = 'W m^-^2';
varn2= 'hfl';  tit2 = 'Mean Sensible Heat Flux';  units2 = 'W m^-^2';
varn1= 'rgn';  tit1 = 'Net Longwave';  units1 = 'W m^-^2';
varn2= 'sgn';  tit2 = 'Net Shortwave';  units2 = 'W m^-^2';
[lat, lon, depth, lm] = getll('temp', lims);
var1 = getflx(varn1, lims, tim);
var2 = getflx(varn2, lims, tim);
var1m = mean(var1);
var2m = mean(var2);
var1m = var1m * 
default_global; 
FRAME = [0 360 -90 90];
cint = 1; clev = [0:cint:15];
subplot(2,1,1)
  gcont(var1m, 0:15);
  dc
  title([tit1 ':  Years 101-1000']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units1]);
FRAME = [0 360 40 90];
cint = 0.1; clev = [-1:cint:1];
subplot(2,2,3);
  mcont(var2m, clev, 'stereo', [90 270]);
  title([tit2 ':  Years 101-1000']);
  xlabel(['Contour Interval:  ' num2str(100*cint) ' ' units2]);
  y = axis
FRAME = [0 360 -90 -40];
cint = 0.1; clev = [-1:cint:1];
subplot(2,2,4);
  mcont(var2m, clev, 'stereo', [-90 270]);
  setm(gca, 'FLatLimit', [40 90], 'MapLatLimit', [-40 -90]);
  title([tit2 ':  Years 101-1000']);
  xlabel(['Contour Interval:  ' num2str(100*cint) ' ' units2]);
cd /home/disk/tao/dvimont/matlab/CSIRO/Plots_Clim
default_global; 
FRAME = [0 360 -90 90];
cint = 25; clev = [-500:cint:500];
subplot(2,1,1)
  gcont(var1m, clev);
  dc
  title([tit1 ':  Years 101-1000']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units1]);
FRAME = [0 360 -90 90];
cint = 25; clev = [-500:cint:500];
subplot(2,1,2)
  gcont(-1*var2m, clev);
  dc
  title([tit2 ':  Years 101-1000']);
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units2]);
cd /home/disk/tao/dvimont/matlab/CSIRO/Plots_Clim
lh = getflx('evp', lims, tim);
lh = squeeze(mean(lh));
sh = getflx('hfl', lims, tim);
sh = squeeze(mean(sh));
lw = getflx('rgn', lims, tim);
lw = squeeze(mean(lw));
sw = getflx('sgn', lims, tim);
sw = squeeze(mean(sw));
net = getnc('heat', lims, 1, tim);
net = -1*squeeze(mean(net));
sumflx = lh + sh + lw - sw;
subplot(2,1,1);
  gcont(sumflx, clev);
  dc 
subplot(2,1,1);
  gcont(sumflx, clev);
  dc
  title('Sum of Mean Heat Flux Component:  Years 101-1000');
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units2]);
fadj = net - sumflx;
subplot(2,1,2);
  gcont(fadj, clev);
  dc
  title('Flux Adjustment');
  xlabel(['Contour Interval:  ' num2str(cint) ' ' units2]);
cd /home/disk/tao/dvimont/matlab/CSIRO/Plots_Clim