Global Index (short | long) | Local contents | Local Index (short | long)
Look at CT runs:
This script calls | |
---|---|
clear if 0 % cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct_5093/som_cycle/data cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct_5093/cycle/data varn1 = 'Z3'; varn2 = 'AIRT'; varn3 = 'V'; varn4 = 'U';% var5 = 'VV'; var6 = 'U'; varn7 = 'PS' ind = 5:18; for i = 1:19; i if i < 10; filin = ['000' num2str(i) '.nc']; else filin = ['00' num2str(i) '.nc']; end nc = netcdf(filin, 'nowrite'); v1 = squeeze(nc{varn1}(:,ind,:,:)); v2 = squeeze(nc{varn2}(:,ind,:,:)); v3 = squeeze(nc{varn3}(:,ind,:,:)); v4 = squeeze(nc{varn4}(:,ind,:,:)); % v5 = squeeze(nc{var5}(:,ind,:,:)); % v6 = squeeze(nc{var6}(:,ind,:,:)); v7 = squeeze(nc{varn7}(:,:,:)); nc = close(nc); if i == 1; tind = [1:2]; else tind = 12*(i-2) + 2 + [1:12]; end var1(tind,:,:,:) = v1; var2(tind,:,:,:) = v2; var3(tind,:,:,:) = v3; var4(tind,:,:,:) = v4; ps(tind,:,:) = v7; end clear v1 v2 v3 v4 v5 v6 v7 for i = 1:9; indc = 24*(i-1)+[1:6]; indw = 24*(i-1)+[13:18]; var1w(i,:,:,:) = mean(var1(indw,:,:,:)); var1c(i,:,:,:) = mean(var1(indc,:,:,:)); var2w(i,:,:,:) = mean(var2(indw,:,:,:)); var2c(i,:,:,:) = mean(var2(indc,:,:,:)); var3w(i,:,:,:) = mean(var3(indw,:,:,:)); var3c(i,:,:,:) = mean(var3(indc,:,:,:)); var4w(i,:,:,:) = mean(var4(indw,:,:,:)); var4c(i,:,:,:) = mean(var4(indc,:,:,:)); pw(i,:,:) = mean(ps(indw,:,:)); pc(i,:,:) = mean(ps(indc,:,:)); end clear var1 var2 var3 var4 ps cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data [lat, lon, lev, hyam, hybm, P0] = getll('wgr_SOM_hard.nc'); hyam = hyam(ind); hybm = hybm(ind); get_global; default_global; newlev = [150 250 300 500 850 970 1000]; z3w = atlev(var1w, newlev, pw, hyam, hybm, P0); z3c = atlev(var1c, newlev, pc, hyam, hybm, P0); airtw = atlev(var2w, newlev, pw, hyam, hybm, P0); airtc = atlev(var2c, newlev, pc, hyam, hybm, P0); vw = atlev(var3w, newlev, pw, hyam, hybm, P0); vc = atlev(var3c, newlev, pc, hyam, hybm, P0); uw = atlev(var4w, newlev, pw, hyam, hybm, P0); uc = atlev(var4c, newlev, pc, hyam, hybm, P0); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data %save uvzt_ctdom.mat z3w z3c airtw airtc vw vc uw uc newlev lat lon end if 1 cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data load uvzt_ctsom.mat z3w1 = z3w; z3c1 = z3c; airtw1 = airtw; airtc1 = airtc; vw1 = vw; vc1 = vc; uw1 = uw; uc1 = uc; load uvzt_ctdom.mat end get_global; default_global; lev = 4; varn = 'z3'; tit1 = ['NDJFMA CT\_SOM: ' num2str(newlev(lev)) 'mb HEIGHT (WCT - CCT)/2']; tit2 = ['NDJFMA CT\_DOM: ' num2str(newlev(lev)) 'mb HEIGHT (WCT - CCT)/2']; tit3 = ['NDJFMA CT: ' num2str(newlev(lev)) 'mb HEIGHT']; figure(1); FRAME = [0 360 -90 90]; cint = 10; clev = [-90:cint:90]; xl = ' m'; lev = 4; varn = 'u'; tit1 = ['NDJFMA CT\_SOM: ' num2str(newlev(lev)) 'mb U (WCT - CCT)/2']; tit2 = ['NDJFMA CT\_DOM: ' num2str(newlev(lev)) 'mb U (WCT - CCT)/2']; tit3 = ['NDJFMA CT: ' num2str(newlev(lev)) 'mb U']; figure(1); FRAME = [0 360 -90 90]; cint = 1; clev = [-9:cint:9]; xl = ' m s^-^1'; eval(['var1 = squeeze(' varn 'w1(:,lev,:,:));']); eval(['var2 = squeeze(' varn 'c1(:,lev,:,:));']); [tstat, dof] = ttest(var1, var2); score = tscore(dof, 2.5); top = squeeze(mean(var1 - var2))/2; sp(1) if ismap(gca); clma; end; mcont(top, clev); greyshd(abs(tstat), score); title(tit1); xlabel(['Contour Interval: ' num2str(cint) xl]); eval(['var1 = squeeze(' varn 'w(:,lev,:,:));']); eval(['var2 = squeeze(' varn 'c(:,lev,:,:));']); [tstat, dof] = ttest(var1, var2); score = tscore(dof, 2.5); top = squeeze(mean(var1 - var2))/2; sp(2) if ismap(gca); clma; end; mcont(top, clev); greyshd(abs(tstat), score); title(tit2); xlabel(['Contour Interval: ' num2str(cint) xl]); figure(2); figure_orient; i = 1; i = 2; eval(['var1 = squeeze(' varn 'w1(:,lev,:,:)) - squeeze(' varn 'c1(:,lev,:,:));']); eval(['var2 = squeeze(' varn 'w(:,lev,:,:)) - squeeze(' varn 'c(:,lev,:,:));']); [tstat, dof] = ttest(var1, var2); score = tscore(dof, 2.5); top = squeeze(mean(var1 - var2))/2; sp(i) if ismap(gca); clma; end; mcont(top, clev); greyshd(abs(tstat), score); title([tit3 ': 0.5*((WSOM-CSOM)-(WDOM-CDOM))']) xlabel(['Contour Interval: ' num2str(cint) xl]);