Global Index (short | long) | Local contents | Local Index (short | long)
save eddy_interped_grsom.mat vtw vtc kew kec newlev lat lon
This script calls | |
---|---|
clear if 0; cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/som_5093/data var1 = 'VT'; var2 = 'AIRT'; var3 = 'V'; var4 = 'UU'; var5 = 'VV'; var6 = 'U'; var7 = 'PS' ind = 5:18; nc = netcdf('EDDY_time_wgr.nc', 'nowrite'); v1 = squeeze(nc{var1}(:,ind,:,:)); v2 = squeeze(nc{var2}(:,ind,:,:)); v3 = squeeze(nc{var3}(:,ind,:,:)); v4 = squeeze(nc{var4}(:,ind,:,:)); v5 = squeeze(nc{var5}(:,ind,:,:)); v6 = squeeze(nc{var6}(:,ind,:,:)); v7 = squeeze(nc{var7}(:,:,:)); nc = close(nc); psw(:,:,:) = v7(:,:,:); e1w = v1-v2.*v3; e2w = v4+v5-v6.^2-v3.^2; clear v1 v2 v3 v4 v5 v6 v7 nc = netcdf('EDDY_time_cgr.nc', 'nowrite'); v1 = squeeze(nc{var1}(:,ind,:,:)); v2 = squeeze(nc{var2}(:,ind,:,:)); v3 = squeeze(nc{var3}(:,ind,:,:)); v4 = squeeze(nc{var4}(:,ind,:,:)); v5 = squeeze(nc{var5}(:,ind,:,:)); v6 = squeeze(nc{var6}(:,ind,:,:)); v7 = squeeze(nc{var7}(:,:,:)); nc = close(nc); e1c = v1-v2.*v3; e2c = v4+v5-v6.^2-v3.^2; psc(:,:,:) = v7(:,:,:); clear v1 v2 v3 v4 v5 v6 v7 for i = 1:10; ind = 12*(i-1)+[1:6]; vtw(i,:,:,:) = mean(e1w(ind,:,:,:)); vtc(i,:,:,:) = mean(e1c(ind,:,:,:)); kew(i,:,:,:) = mean(e2w(ind,:,:,:)); kec(i,:,:,:) = mean(e2c(ind,:,:,:)); pw(i,:,:) = mean(psw(ind,:,:)); pc(i,:,:) = mean(psc(ind,:,:)); end ind = 5:18; 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]; vtw = atlev(vtw, newlev, pw, hyam, hybm, P0); vtc = atlev(vtc, newlev, pc, hyam, hybm, P0); kew = atlev(kew, newlev, pw, hyam, hybm, P0); kec = atlev(kec, newlev, pc, hyam, hybm, P0); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data else cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data load eddy_interped_grsom.mat % load eddy_interped_grdom.mat end get_global; default_global; levt = 5; levb = 3; [tstatt, dof] = ttest(vtw(:,levt,:,:), vtc(:,levt,:,:)); scoret = tscore(dof, 2.5); top = squeeze(mean(vtw(:,levt,:,:) - vtc(:,levt,:,:)))/2; cintt = 1; clevt = [-5:cintt:5]; [tstatb, dof] = ttest(kew(:,levb,:,:), kec(:,levb,:,:)); scoreb = tscore(dof, 2.5); bot = squeeze(mean(kew(:,levb,:,:) - kec(:,levb,:,:)))/2; cintb = 15; clevb = [-90:cintb:90]; tit1 = ['NDJFMAM GR Run: ' num2str(newlev(levt)) 'mb v''T'', (WGR - CGR)/2']; tit2 = ['NDJFMAM GR Run: ' num2str(newlev(levb)) 'mb EKE, (WGR - CGR)/2']; clabt = [' K m s^-^1']; clabb = [' m^2 s^-^2']; lev = 5; top = squeeze(mean(vtw(:,lev,:,:))); bot = squeeze(mean(vtc(:,lev,:,:))); cintt = 7.5; clevt = [-60:cintt:60]; cintb = cintt; clevb = clevt; clabt = [' K m s^-^1']; clabb = clabt; tit1 = ['WGR ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM']; tit2 = ['CGR ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM']; lev = 3; top = squeeze(mean(kew(:,lev,:,:))); bot = squeeze(mean(kec(:,lev,:,:))); cintt = 100; clevt = [-600:cintt:600]; cintb = cintt; clevb = clevt; clabt = [' K m s^-^1']; clabb = clabt; tit1 = ['WGR ' num2str(newlev(lev)) 'mb EKE, NDJFMAM']; tit2 = ['CGR ' num2str(newlev(lev)) 'mb EKE, NDJFMAM']; figure(1); FRAME = [0 360 -90 90]; sp(1) if ismap(gca); clma; end; mcont(top, clevt); title(tit1); xlabel(['Contour Interval: ' num2str(cintt) clabt]) % greyshd(abs(tstatt), scoret); sp(2) if ismap(gca); clma; end; mcont(bot, clevb); title(tit2); xlabel(['Contour Interval: ' num2str(cintb) clabb]) % greyshd(abs(tstatb), scoreb); % Look at CT runs: 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 var1 = 'VT'; var2 = 'AIRT'; var3 = 'V'; var4 = 'UU'; var5 = 'VV'; var6 = 'U'; var7 = '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{var1}(:,ind,:,:)); v2 = squeeze(nc{var2}(:,ind,:,:)); v3 = squeeze(nc{var3}(:,ind,:,:)); v4 = squeeze(nc{var4}(:,ind,:,:)); v5 = squeeze(nc{var5}(:,ind,:,:)); v6 = squeeze(nc{var6}(:,ind,:,:)); v7 = squeeze(nc{var7}(:,:,:)); nc = close(nc); if i == 1; tind = [1:2]; else tind = 12*(i-2) + 2 + [1:12]; end e1(tind,:,:,:) = v1-v2.*v3; e2(tind,:,:,:) = v4+v5-v6.^2-v3.^2; 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]; vtw(i,:,:,:) = mean(e1(indw,:,:,:)); vtc(i,:,:,:) = mean(e1(indc,:,:,:)); kew(i,:,:,:) = mean(e2(indw,:,:,:)); kec(i,:,:,:) = mean(e2(indc,:,:,:)); pw(i,:,:) = mean(ps(indw,:,:)); pc(i,:,:) = mean(ps(indc,:,:)); end 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]; vtw = atlev(vtw, newlev, pw, hyam, hybm, P0); vtc = atlev(vtc, newlev, pc, hyam, hybm, P0); kew = atlev(kew, newlev, pw, hyam, hybm, P0); kec = atlev(kec, newlev, pc, hyam, hybm, P0); cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data %save eddy_interped_ctdom.mat vtw vtc kew kec newlev lat lon else cd /home/disk/hayes2/dvimont/ccm/ccm3.6/matlab_data load eddy_interped_ctsom.mat % load eddy_interped_ctdom.mat end get_global; default_global; levt = 5; levb = 3; [tstatt, dof] = ttest(vtw(:,levt,:,:), vtc(:,levt,:,:)); scoret = tscore(dof, 2.5); top = squeeze(mean(vtw(:,levt,:,:) - vtc(:,levt,:,:)))/2; cintt = 1; clevt = [-5:cintt:5]; [tstatb, dof] = ttest(kew(:,levb,:,:), kec(:,levb,:,:)); scoreb = tscore(dof, 2.5); bot = squeeze(mean(kew(:,levb,:,:) - kec(:,levb,:,:)))/2; cintb = 15; clevb = [-90:cintb:90]; tit1 = ['NDJFMAM CT Run: ' num2str(newlev(levt)) 'mb v''T'', (WCT - CCT)/2']; tit2 = ['NDJFMAM CT Run: ' num2str(newlev(levb)) 'mb EKE, (WCT - CCT)/2']; clabt = [' K m s^-^1']; clabb = [' m^2 s^-^2']; lev = 5; top = squeeze(mean(vtw(:,lev,:,:))); bot = squeeze(mean(vtc(:,lev,:,:))); cintt = 7.5; clevt = [-60:cintt:60]; cintb = cintt; clevb = clevt; clabt = [' K m s^-^1']; clabb = clabt; tit1 = ['WCT ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM']; tit2 = ['CCT ' num2str(newlev(lev)) 'mb v''T'', NDJFMAM']; lev = 3; top = squeeze(mean(kew(:,lev,:,:))); bot = squeeze(mean(kec(:,lev,:,:))); cintt = 100; clevt = [-600:cintt:600]; cintb = cintt; clevb = clevt; clabt = [' K m s^-^1']; clabb = clabt; tit1 = ['WCT ' num2str(newlev(lev)) 'mb EKE, NDJFMAM']; tit2 = ['CCT ' num2str(newlev(lev)) 'mb EKE, NDJFMAM']; figure(1); FRAME = [0 360 -90 90]; sp(1) if ismap(gca); clma; end; mcont(top, clevt); title(tit1); xlabel(['Contour Interval: ' num2str(cintt) clabt]) % greyshd(abs(tstatt), scoret); sp(2) if ismap(gca); clma; end; mcont(bot, clevb); title(tit2); xlabel(['Contour Interval: ' num2str(cintb) clabb]) % greyshd(abs(tstatb), scoreb);