Documentation of prec_GR_CT

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

Help text

print -dps2 PREC_CT_GR_t31_DOM.ps

Cross-Reference Information

Listing of script prec_GR_CT

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[pr1, pr2] = getnc('wgr_5093.nc', 'PRECC', 'PRECL');
pw = (pr1 + pr2) * 3600 * 24 * 1000;
[pr1, pr2] = getnc('cgr_5093.nc', 'PRECC', 'PRECL');
pc = (pr1 + pr2) * 3600 * 24 * 1000;
[lat, lon] = getll('wgr_5093.nc');
tind = [1 2 12];
pgr = squeeze(mean(pw(tind,:,:) - pc(tind,:,:)));
default_global; FRAME = [90 300 -40 40];
figure(1);
subplot(2,1,1)
     gcont(pgr, [-15:-1 1:15]);
     dc; box;
     title('GR\_DOM t31x15 Precipitation (WARM - COLD)')
     xlabel('Contour Interval 1 mm day^-^1')
[pr1, pr2] = getnc('w_cycle.nc', 'PRECC', 'PRECL');
pw = (pr1 + pr2) * 3600 * 24 * 1000;
[pr1, pr2] = getnc('c_cycle.nc', 'PRECC', 'PRECL');
pc = (pr1 + pr2) * 3600 * 24 * 1000;
tind = [2 3 4];
pgr = squeeze(mean(pw(tind,:,:) - pc(tind,:,:)));
subplot(2,1,2)
     gcont(pgr, [-20:2:-2 2:2:20]);
     dc; box;
     title('CT\_DOM t31x15 Precipitation (WARM - COLD)')
     xlabel('Contour Interval 2 mm day^-^1')
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/CT_GR

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data

[pr1, pr2] = getnc('wgr_SOM_hard.nc', 'PRECC', 'PRECL');
pw = (pr1 + pr2) * 3600 * 24 * 1000;
[pr1, pr2] = getnc('cgr_SOM_hard.nc', 'PRECC', 'PRECL');
pc = (pr1 + pr2) * 3600 * 24 * 1000;
pgr = squeeze(mean(pw(tind,:,:) - pc(tind,:,:)));
[lat, lon] = getll('wgr_5093.nc');

subplot(2,1,1)
     gcont(pgr, [-5:-1 1:5]);
     dc; box;
     title('GR\_SOM t31x15 Precipitation (WARM - COLD)')
     xlabel('Contour Interval 1 mm day^-^1')

[pr1, pr2] = getnc('w_SOM_cycle.nc', 'PRECC', 'PRECL');
pw = (pr1 + pr2) * 3600 * 24 * 1000;
[pr1, pr2] = getnc('c_SOM_cycle.nc', 'PRECC', 'PRECL');
pc = (pr1 + pr2) * 3600 * 24 * 1000;
pgr = squeeze(mean(pw - pc));

subplot(2,1,2)
     gcont(pgr, [-20:2:-2 2:2:20]);
     dc; box;
     title('CT\_SOM t31x15 Precipitation (WARM - COLD)')
     xlabel('Contour Interval 2 mm day^-^1')

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/CT_GR
%print -dps2 PREC_CT_GR_t31_SOM.ps

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data

[pr1, pr2] = getnc('wgr.nc', 'PRECC', 'PRECL');
pw = (pr1 + pr2) * 3600 * 24 * 1000;
[pr1, pr2] = getnc('cgr.nc', 'PRECC', 'PRECL');
pc = (pr1 + pr2) * 3600 * 24 * 1000;
pgr = squeeze(mean(pw - pc));
[lat, lon] = getll('wgr.nc');
ctlim = FRAME;
default_global

subplot(2,1,1)
     gcont(pgr, [-5:-1 1:5]);
     dc; box;
     title('GR\_DOM T42 Precipitation (WARM - COLD)')
     xlabel('Contour Interval 1 mm day^-^1')

[pr1, pr2] = getnc('wct.nc', 'PRECC', 'PRECL');
pw = (pr1 + pr2) * 3600 * 24 * 1000;
[pr1, pr2] = getnc('cct.nc', 'PRECC', 'PRECL');
pc = (pr1 + pr2) * 3600 * 24 * 1000;
pgr = squeeze(mean(pw - pc));

subplot(2,1,2)
     gcont(pgr, [-20:2:-2 2:2:20]);
     dc; box;
     title('CT\_DOM T42 Precipitation (WARM - COLD)')
     xlabel('Contour Interval 2 mm day^-^1')

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/CT_GR