Documentation of wind_CT_GR

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

Cross-Reference Information

Listing of script wind_CT_GR

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[uw, vw] = getnc('wgr_5093.nc', 18, 'U', 'V');
[uc, vc] = getnc('cgr_5093.nc', 18, 'U', 'V');
[lat, lon] = getll('wgr_5093.nc');
tind = [1 2 3 4 11 12];
u = squeeze(mean(uw(tind,:,:,:) - uc(tind,:,:,:)));
v = squeeze(mean(vw(tind,:,:,:) - vc(tind,:,:,:)));
ctlim = [90 300 -40 40];
default_global;
figure(1)
subplot(2,1,1)
     [h, a, c] = gquiv(u, v, 1, 1);
     dc; box;
     title('GR\_DOM t31x15 Surface Winds (WARM - COLD)')
     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
[uw, vw] = getnc('w_cycle.nc', 18, 'U', 'V');
[uc, vc] = getnc('c_cycle.nc', 18, 'U', 'V');
[lat, lon] = getll('wgr_5093.nc');
tind = [1 2 3 4 11 12];
u = squeeze(mean(uw - uc));
v = squeeze(mean(vw - vc));
ctlim = [90 300 -40 40];
default_global;
figure(1)
subplot(2,1,2)
     [h, a, c] = gquiv(u, v, 1, 1);
     dc; box;
     title('CT\_DOM t31x15 Surface Winds (WARM - COLD)')
     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/CT_GR
print -dps2 SURF_WIND_CT_GR_t31_DOM.ps
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[uw, vw] = getnc('wgr_SOM_hard.nc', 18, 'U', 'V');
[uc, vc] = getnc('cgr_SOM_hard.nc', 18, 'U', 'V');
tind = [1 2 3 4 11 12];
u = squeeze(mean(uw(tind,:,:,:) - uc(tind,:,:,:)));
v = squeeze(mean(vw(tind,:,:,:) - vc(tind,:,:,:)));
ctlim = [90 300 -40 40];
default_global;
figure(1)
subplot(2,1,1)
     [h, a, c] = gquiv(u, v, 1, 1);
     dc; box;
     title('GR\_SOM t31x15 Surface Winds (WARM - COLD)')
     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
[uw, vw] = getnc('w_SOM_cycle.nc', 18, 'U', 'V');
[uc, vc] = getnc('c_SOM_cycle.nc', 18, 'U', 'V');
[lat, lon] = getll('wgr_5093.nc');
tind = [1 2 3 4 11 12];
u = squeeze(mean(uw - uc));
v = squeeze(mean(vw - vc));
ctlim = [90 300 -40 40];
default_global;
figure(1)
subplot(2,1,2)
     [h, a, c] = gquiv(u, v, 1, 1);
     dc; box;
     title('CT\_SOM t31x15 Surface Winds (WARM - COLD)')
     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/CT_GR
print -dps2 SURF_WIND_CT_GR_t31_SOM.ps
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[uw, vw] = getnc('wgr.nc', 18, 'U', 'V');
[uc, vc] = getnc('cgr.nc', 18, 'U', 'V');
[lat, lon] = getll('wgr.nc');
u = squeeze(mean(uw - uc));
v = squeeze(mean(vw - vc));
ctlim = [90 300 -40 40];
default_global;
figure(1)
subplot(2,1,1)
     [h, a, c] = gquiv(u, v, 1, 2);
     dc; box;
     title('GR\_DOM T42 Surface Winds (WARM - COLD)')
     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
[uw, vw] = getnc('wct.nc', 18, 'U', 'V');
[uc, vc] = getnc('cct.nc', 18, 'U', 'V');
u = squeeze(mean(uw - uc));
v = squeeze(mean(vw - vc));
ctlim = [90 300 -40 40];
default_global;
figure(1)
subplot(2,1,2)
     [h, a, c] = gquiv(u, v, 1, 2);
     dc; box;
     title('CT\_DOM T42 Surface Winds (WARM - COLD)')
     xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/CT_GR
print -dps2 SURF_WIND_CT_GR_T42_DOM.ps