Documentation of make_Z500_plots

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

Cross-Reference Information

Listing of script make_Z500_plots

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct_t31_ense/data
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[zw, pw] = getnc('wgr_ice.nc', 'Z3', 'PS');
[zc, pc] = getnc('cgr_ice.nc', 'Z3', 'PS');
[zw, pw] = getnc('wgr_SOM_hard.nc', 'Z3', 'PS');
[zc, pc] = getnc('cgr_SOM_hard.nc', 'Z3', 'PS');
[zw, pw] = getnc('wct_ense.nc', 'Z3', 'PS');
[zc, pc] = getnc('cct_ense.nc', 'Z3', 'PS');
[lat, lon, lev, hyam, hybm, P0] = getll('cct_ense.nc');
get_global
default_global
if size(zw,1) == 12;
  zw = zw([1 2 3 4 11 12], :, :, :);
  pw = pw([1 2 3 4 11 12], :, :);
  zc = zc([1 2 3 4 11 12], :, :, :);
  pc = pc([1 2 3 4 11 12], :, :);
end
lev2 = 500;
zw2c = atlev(zw(:,5:18,:,:), lev2, pw, hyam(5:18), hybm(5:18), P0);
zc2c = atlev(zc(:,5:18,:,:), lev2, pc, hyam(5:18), hybm(5:18), P0);
top = squeeze(mean(zw2c-zc2c))/2;
figure(1);
sp(1)
     if ismap(gca); clma; end;
     mcont(top, [-100:20:100]);
     title('CT\_ENSEMBLE Run:  DJF 500mb Height (WARM - COLD)/2');
     xlabel('Contour Interval 20m');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
zw = getnc('wgr_SOM_hard.nc', 11, 'Z3');
zc = getnc('cgr_SOM_hard.nc', 11, 'Z3');
zw = getnc('wgr_ice.nc', 11, 'Z3');
zc = getnc('cgr_ice.nc', 11, 'Z3');
zw = squeeze(zw);
zc = squeeze(zc);
top = squeeze(mean(zw2-zc2))/2;
sp(2)
     if ismap(gca); clma; end;
     mcont(top, [-100:5:100]);
     title('GR\_SOM\_Hard Run:  500mb Height (WARM - COLD)/2');
     xlabel('Contour Interval 5m');
cd /home/disk/tao/dvimont/pub/david
FRAME = [-.1 360 -90 90];
[xk, yk] = keep_var(FRAME, lon, lat);
lat2 = lat(yk); lon2 = lon(xk);
XAX = lon2; YAX = lat2; var2 = top(yk, xk);