Global Index (short | long) | Local contents | Local Index (short | long)
Plot out patterns regressed on LP filtered GR time series:
This script calls | |
---|---|
clear if 0; ind = 11; cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data var1 = getnc('Z3_wgr.nc', ind, 'Z3'); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data var2 = getnc('Z3_cgr.nc', ind, 'Z3'); [ntim, nlev, nlat, nlon] = size(var1); [lat, lon] = getll('Z3_cgr.nc'); end if 0; cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/wgr/data [a, b] = getnc('PREC_wgr.nc', 'PRECL', 'PRECC'); var1 = (a+b)*(1200*24*3600); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/gr_5093/cgr/data [a, b] = getnc('PREC_cgr.nc', 'PRECL', 'PRECC'); var2 = (a+b)*(1200*24*3600); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data [a, b] = getnc('PREC_clim.nc', 'PRECL', 'PRECC'); var3 = (a+b)*(1200*24*3600); [ntim, nlat, nlon] = size(var1); [lat, lon] = getll('Z3_cgr.nc'); end get_global default_global win = sort([2:12:ntim 3:12:ntim 4:12:ntim]); cint = 7.5; var1win = squeeze(var1(win,:,:,:)); var2win = squeeze(var2(win,:,:,:)); [tstat, dof] = ttest(var1win, var2win); score = tscore(dof, 2.5); top = squeeze(mean(var1win - var2win))/2; figure(2) sp(1); clma mcont(top, [-90:cint:110], [90 270]); greyshd(abs(tstat), score); title('500mb Z3: (WGR\_NOICE - CGR\_NOICE) / 2'); xlabel(['Contour Interval: ' num2str(cint) ' m']); sp(2); clma mcont(top, [-90:cint:110], [-90 270]); greyshd(abs(tstat), score); title('500mb Z3: (WGR\_NOICE - CGR\_NOICE) / 2'); xlabel(['Contour Interval: ' num2str(cint) ' m']); cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots figure(2) FRAME = [40 300 -45 45]; cint = 1; var1win = squeeze(var1(win,:,:,:)); var2win = squeeze(var3(win,:,:,:)); [tstat, dof] = ttest(var1win, var2win); score = tscore(dof, 2.5); top = squeeze(mean(var1win - var2win))/1; sp(1); clma mcont(top, [-10:cint:10], 'giso'); greyshd(abs(tstat), score); title('PRECIP: (WGR\_NOICE - CGR\_NOICE) / 2'); xlabel(['Contour Interval: ' num2str(cint) ' mm (day)^-^1']); var1win = squeeze(var2(win,:,:,:)); var2win = squeeze(var3(win,:,:,:)); [tstat, dof] = ttest(var1win, var2win); score = tscore(dof, 2.5); top = squeeze(mean(var1win - var2win))/1; sp(2); clma mcont(top, [-10:cint:10], 'giso'); greyshd(abs(tstat), score); title('PRECIP: (WGR\_NOICE - CGR\_NOICE) / 2'); xlabel(['Contour Interval: ' num2str(cint) ' mm (day)^-^1']); get_global cd /home/disk/tao/data/nmc.reanalysis/monthly filin = 'hgt.mon.mean.nc'; [lat, lon] = getll(filin); default_global cd /home/disk/tao/dvimont/matlab/CCM/GR filin = 'hgt.mon.mean.nc'; var = 'hgt'; [pat500, clim] = reg_gr(filin, var, 250); figure(1) sp(1) cla mcont(pat500, [-100:10:100]); title('NMC 250mb HEIGHT REGRESSED ON GR') xlabel('Contour Interval: 10 m (std)^-^1') sp(2) cla mcont(pat500, [-100:10:100], [-90 270]); title('NMC 250mb HEIGHT REGRESSED ON GR') xlabel('Contour Interval: 10 m (std)^-^1') cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots