Global Index (short | long) | Local contents | Local Index (short | long)
Load GR pattern and time series
This script calls | |
---|---|
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data load ct_gr.mat; % Get CCM data: cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data filin = 'wgr_TropPac.nc'; nc = netcdf(filin, 'nowrite'); hgtw = nc{'Z3'}(1:3,:,:,:); psw = nc{'PS'}(1:3,:,:); tsw = nc{'TS'}(1:3,:,:); nc = close(nc); filin = 'cgr_TropPac.nc'; nc = netcdf(filin, 'nowrite'); hgtc = nc{'Z3'}(1:3,:,:,:); psc = nc{'PS'}(1:3,:,:); tsc = nc{'TS'}(1:3,:,:); lat = nc{'lat'}(:); lon = nc{'lon'}(:); hybm = nc{'hybm'}(:); hyam = nc{'hyam'}(:); levs = nc{'lev'}(:); P0 = nc{'P0'}(:); nc = close(nc); % Get NMC data: cd /home/disk/tao/data/nmc.reanalysis/monthly filin = 'hgt.mon.mean.nc'; nc = netcdf(filin, 'nowrite'); nmclev = nc{'level'}(:); nmclat = nc{'lat'}(:); nmclon = nc{'lon'}(:); nc = close(nc); hgtw_nmc = atlev(hgtw, nmclev, psw, hyam, hybm, P0); hgtc_nmc = atlev(hgtc, nmclev, psc, hyam, hybm, P0); cd /home/disk/tao/dvimont/matlab/CCM/GR filin = 'hgt.mon.mean.nc'; var_nam = 'hgt'; lev = 1000; [hgt1000gr, hgt1000cl] = reg_gr(filin, var_nam, lev); lev = 500; [hgt500gr, hgt500cl] = reg_gr(filin, var_nam, lev); lev = 250; [hgt250gr, hgt250cl] = reg_gr(filin, var_nam, lev); %lev = 200; %[hgt200gr, hgt200cl] = reg_gr(filin, var_nam, lev); %lev = 850; %[hgt850gr, hgt850cl] = reg_gr(filin, var_nam, lev); cd /home/disk/tao/dvimont/matlab/CCM/GR/DATA %save hgtlp.mat hgt1000gr hgt1000cl hgt500gr hgt500cl hgt250gr hgt250cl hgtw_nmc hgtc_nmc %load hgt.mat load hgtlp.mat cd /home/disk/tao/dvimont/matlab/CCM/GR [lat, lon] = getll('wgr.nc'); [nmclat, nmclon] = getll('hgt.mon.mean.nc'); [nmclev] = getnc('hgt.mon.mean.nc', [1:17], 'level'); define_global XAX = nmclon; YAX = nmclat; FRAME = [0 360 -90 90]; top = hgt250gr; lev = 250; cint = 5; figure(1) sp(1) mcont(-1 * top, [-50:cint:50], [90 270]); title(['NMC: NH ' num2str(lev) 'mb HGT Regressed on GR']); xlabel(['Contour Interval: ' num2str(cint) 'm']) sp(2) mcont(-1 * top, [-50:cint:50], [-90 270]); title(['NMC: SH ' num2str(lev) 'mb HGT Regressed on GR']); xlabel(['Contour Interval: ' num2str(cint) 'm']) cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots XAX = lon; YAX = lat; lev = 250; cint = 5; lind = find(nmclev == lev); top = squeeze(mean(hgtw_nmc(:,lind,:,:) - hgtc_nmc(:,lind,:,:)))/2; figure(2) sp(1) mcont(top, [-500:cint:500], [90 270]); title(['GR Run: NH ' num2str(lev) 'mb HGT (WARM - COLD)/2']); xlabel(['Contour Interval: ' num2str(cint) 'm']) sp(2) mcont(1 * top, [-500:cint:500], [-90 270]); title(['GR Run: SH ' num2str(lev) 'mb HGT (WARM - COLD)/2']); xlabel(['Contour Interval: ' num2str(cint) 'm']) cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots % Check TS out: figure(2) orient landscape gcont(squeeze(mean(tsw - tsc)/2), [-2:.25:2]); dc % Look at individual maps of hgtw_nmc and hgtc_nmc lev = 500; cint = 50; lind = find(nmclev == lev); top = squeeze(mean(hgtw_nmc(:,lind,:,:))); top = top - (mean2(top'))' * ones(1, size(top, 2)); bot = squeeze(mean(hgtc_nmc(:,lind,:,:))); bot = bot - (mean2(bot'))' * ones(1, size(bot, 2)); figure(1) sp(1) mcont(top, [-500:cint:500], [90 270]); title(['GR Run: NH ' num2str(lev) 'mb HGT (WARM - ZONAL\_MEAN)']); xlabel(['Contour Interval: ' num2str(2*cint) 'm']) sp(2) mcont(bot, [-500:cint:500], [90 270]); title(['GR Run: NH ' num2str(lev) 'mb HGT (COLD - ZONAL\_MEAN)']); xlabel(['Contour Interval: ' num2str(2*cint) 'm']) cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots