Global Index (short | long) | Local contents | Local Index (short | long)
cnum = 10; [clevt, cintt] = get_contour_int(top, cnum); [clevb, cintb] = get_contour_int(bot, cnum);
This script calls | |
---|---|
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data [lat1, lon1, lev1, hyam, hybm, P0] = getll('wct_globsst_t31.nc'); [lat2, lon2, lev2] = getll('wct_globsst_t31.nc'); [whgt1, wprc1, wprl1, wu1, wv1, wps1, wpsl1] = ... getnc('wct_globsst_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ... 'PSL'); [chgt1, cprc1, cprl1, cu1, cv1, cps1, cpsl1] = ... getnc('cct_globsst_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ... 'PSL'); [whgt2, wprc2, wprl2, wu2, wv2, wps2, wpsl2] = ... getnc('wct_TropPac_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ... 'PSL'); [chgt2, cprc2, cprl2, cu2, cv2, cps2, cpsl2] = ... getnc('cct_TropPac_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ... 'PSL'); whgt1 = squeeze(mean(whgt1)); chgt1 = squeeze(mean(chgt1)); whgt2 = squeeze(mean(whgt2)); chgt2 = squeeze(mean(chgt2)); wprc1 = squeeze(mean(wprc1)); cprc1 = squeeze(mean(cprc1)); wprc2 = squeeze(mean(wprc2)); cprc2 = squeeze(mean(cprc2)); wprl1 = squeeze(mean(wprl1)); cprl1 = squeeze(mean(cprl1)); wprl2 = squeeze(mean(wprl2)); cprl2 = squeeze(mean(cprl2)); wu1 = squeeze(mean(wu1)); cu1 = squeeze(mean(cu1)); wu2 = squeeze(mean(wu2)); cu2 = squeeze(mean(cu2)); wv1 = squeeze(mean(wv1)); cv1 = squeeze(mean(cv1)); wv2 = squeeze(mean(wv2)); cv2 = squeeze(mean(cv2)); wps1 = squeeze(mean(wps1)); cps1 = squeeze(mean(cps1)); wps2 = squeeze(mean(wps2)); cps2 = squeeze(mean(cps2)); wpsl1 = squeeze(mean(wpsl1)); cpsl1 = squeeze(mean(cpsl1)); wpsl2 = squeeze(mean(wpsl2)); cpsl2 = squeeze(mean(cpsl2)); get_global FRAME = [0 360 -90 90]; top = squeeze((wprc1 + wprl1) - (cprc1 + cprl1)) * (3600 * 24 * 1000) / 2; bot = squeeze((wprc2 + wprl2) - (cprc2 + cprl2)) * (3600 * 24 * 1000) / 2; units = ['mm (day)^-^1']; head = ['Precipitation: WARM - COLD']; FRAME = [0 360 -40 40]; cnum = 10; top = squeeze(wps1 - cps1) / 200; bot = squeeze(wps2 - cps2) / 200; units = ['mb']; head = ['Surface Pressure: WARM - COLD']; FRAME = [0 360 -90 90]; top = squeeze(wpsl1 - cpsl1) / 200; bot = squeeze(wpsl2 - cpsl2) / 200; units = ['mb']; head = ['SLP: WARM - COLD']; FRAME = [0 360 -90 90]; num = 18; top = squeeze(wu1(num,:,:) - cu1(num,:,:)) / 2; bot = squeeze(wu2(num,:,:) - cu2(num,:,:)) / 2; units = ['M/S']; head = [num2str(10 * round(lev1(num) / 10)) 'mb Zonal Wind: WARM - COLD']; FRAME = [0 360 -90 90]; cnum = 12; num = 17; divw1 = sph_div(squeeze(wu1(num,:,:)), squeeze(wv1(num,:,:)), lon1, lat1); divc1 = sph_div(squeeze(cu1(num,:,:)), squeeze(cv1(num,:,:)), lon1, lat1); divw2 = sph_div(squeeze(wu2(num,:,:)), squeeze(wv2(num,:,:)), lon2, lat2); divc2 = sph_div(squeeze(cu2(num,:,:)), squeeze(cv2(num,:,:)), lon2, lat2); top = (divw1 - divc1) / 2; bot = (divw2 - divc2) / 2; head = [num2str(10 * round(lev1(num) / 10)) 'mb Zonal Wind: WARM - COLD']; units = ['S^-^1']; cnum = 10; FRAME = [0 360 -40 40]; num = 8; curlw1 = sph_curl(squeeze(wu1(num,:,:)), squeeze(wv1(num,:,:)), lon1, lat1); curlc1 = sph_curl(squeeze(cu1(num,:,:)), squeeze(cv1(num,:,:)), lon1, lat1); curlw2 = sph_curl(squeeze(wu2(num,:,:)), squeeze(wv2(num,:,:)), lon2, lat2); curlc2 = sph_curl(squeeze(cu2(num,:,:)), squeeze(cv2(num,:,:)), lon2, lat2); top = (curlw1 - curlc1) / 2; bot = (curlw2 - curlc2) / 2; head = [num2str(10 * round(lev1(num) / 10)) 'mb Relative Vort.: WARM - COLD']; units = ['S^-^1']; cnum = 12; FRAME = [0 360 -90 90]; cnum = 12; [clevt, cintt] = get_contour_int(top, cnum); [clevb, cintb] = get_contour_int(top, cnum); figure(1) sp(1) XAX = lon1; YAX = lat1; mcont(top, clevt, 'giso', [0 180]); title([head ', GLOBAL SST']); xlabel(['Units: ' num2str(cintt) ' ' units]); sp(2) XAX = lon2; YAX = lat2; mcont(bot, clevb, 'giso', [0 180]); title([head ', TPAC SST']); xlabel(['Units: ' num2str(cintb) ' ' units]); cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots num = 11;cintt=10;cintb=10; num = 8;cintt=10;cintb=20; top = squeeze(whgt1(num,:,:) - chgt1(num,:,:)) / 2; bot = squeeze(whgt2(num,:,:) - chgt2(num,:,:)) / 2; units = ['M']; head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: WARM - COLD']; FRAME = [0 360 -90 90]; num = 11; top = squeeze(whgt1(num,:,:) - whgt2(num,:,:)) / 2; bot = squeeze(chgt1(num,:,:) - chgt2(num,:,:)) / 2; units = ['M']; head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: GLOB - TPAC']; FRAME = [0 360 -90 90]; num = 8; top = squeeze(wu1(num,:,:) - cu1(num,:,:)) / 2; bot = squeeze(wu2(num,:,:) - cu2(num,:,:)) / 2; units = ['M/S']; head = [num2str(10 * round(lev1(num) / 10)) 'mb Zonal Wind: WARM - COLD']; FRAME = [0 360 -90 90]; cintt = 2; cintb = 2; num = 10; top = squeeze(wu1(num,:,:) - wu2(num,:,:)) / 2; bot = squeeze(cu1(num,:,:) - cu2(num,:,:)) / 2; units = ['M/S']; head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: GLOB - TPAC']; FRAME = [0 360 -90 90]; clevt = [-100:cintt:100]; clevb = [-100:cintb:100]; figure(1) sp(1) XAX = lon1; YAX = lat1; mcont(top, clevt, 'stereo', [90 270]); title([head ', GLOBAL SST']); xlabel(['Units: ' num2str(cintt) ' ' units]); sp(2) XAX = lon2; YAX = lat2; mcont(bot, clevb, 'stereo', [90 270]); title([head ', TPAC SST']); xlabel(['Units: ' num2str(cintb) ' ' units]); cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data [wuu1, wvv1, wuv1] = getnc('wct_globsst_t31.nc', ... 'UU', 'VV', 'UV'); [cuu1, cvv1, cuv1] = getnc('cct_globsst_t31.nc', ... 'UU', 'VV', 'UV'); [wuu2, wvv2, wuv2] = getnc('wct_TropPac_t31.nc', ... 'UU', 'VV', 'UV'); [cuu2, cvv2, cuv2] = getnc('cct_TropPac_t31.nc', ... 'UU', 'VV', 'UV'); wuu1 = squeeze(mean(wuu1)); wvv1 = squeeze(mean(wvv1)); wuv1 = squeeze(mean(wuv1)); cuu1 = squeeze(mean(cuu1)); cvv1 = squeeze(mean(cvv1)); cuv1 = squeeze(mean(cuv1)); wuu2 = squeeze(mean(wuu2)); wvv2 = squeeze(mean(wvv2)); wuv2 = squeeze(mean(wuv2)); cuu2 = squeeze(mean(cuu2)); cvv2 = squeeze(mean(cvv2)); cuvfigure(1) sp(1) XAX = lon1; YAX = lat1; gquiv(top1, top2, 2); % mquiv(top1, top2, 1, 'stereo', [90 270]); title([head1 ', TPAC SST']); xlabel(['Units: ' units]); dc sp(2) XAX = lon2; YAX = lat2; gquiv(bot1, bot2, 2); % mquiv(bot1, bot2, 1, 'stereo', [90 270]); title([head2 ', TPAC SST']); xlabel(['Units: ' units]); dc 2 = squeeze(mean(cuv2)); weke1 = (wuu1 + wvv1 - wu1.^2 - wv1.^2); ceke1 = (cuu1 + cvv1 - cu1.^2 - cv1.^2); weke2 = (wuu2 + wvv2 - wu2.^2 - wv2.^2); ceke2 = (cuu2 + cvv2 - cu2.^2 - cv2.^2); num = 11; top = squeeze(weke1(num,:,:) - ceke1(num,:,:)) / 2; bot = squeeze(weke2(num,:,:) - ceke2(num,:,:)) / 2; units = ['M^2/S^2']; head = [num2str(10 * round(lev1(num) / 10)) 'mb EKE: WARM - COLD']; FRAME = [0 360 -90 90]; num = 11; top = squeeze(weke1(num,:,:) - weke2(num,:,:)) / 2; bot = squeeze(ceke1(num,:,:) - ceke2(num,:,:)) / 2; units = ['M^2/S^2']; head = [num2str(10 * round(lev1(num) / 10)) 'mb EKE: WARM - COLD']; FRAME = [0 360 -90 90]; cnum = 12; [clevt, cintt] = get_contour_int(top, cnum); [clevb, cintb] = get_contour_int(bot, cnum); figure(1) sp(1) XAX = lon1; YAX = lat1; mcont(top, clevt, 'stereo', [90 270]); title([head ', GLOBAL SST']); xlabel(['Units: ' num2str(cintt) ' ' units]); sp(2) XAX = lon2; YAX = lat2; mcont(bot, clevb, 'stereo', [90 270]); title([head ', TPAC SST']); xlabel(['Units: ' num2str(cintb) ' ' units]); wep1_1 = wvv1 - wv1.^2 - wuu1 + wu1.^2; wep1_2 = wu1.*wv1-wuv1; cep1_1 = cvv1 - cv1.^2 - cuu1 + cu1.^2; cep1_2 = cu1.*cv1-cuv1; wep2_1 = wvv2 - wv2.^2 - wuu2 + wu2.^2; wep2_2 = wu2.*wv2-wuv2; cep2_1 = cvv2 - cv2.^2 - cuu2 + cu2.^2; cep2_2 = cu2.*cv2-cuv2; clf num = 11; top1 = squeeze(wep1_1(num,:,:)-cep1_1(num,:,:))/2; top2 = squeeze(wep1_2(num,:,:)-cep1_2(num,:,:))/2; bot1 = squeeze(wep2_1(num,:,:)-cep2_1(num,:,:))/2; bot2 = squeeze(wep2_2(num,:,:)-cep2_2(num,:,:))/2; units = ['M^2/S^2']; head = [num2str(10 * round(lev1(num) / 10)) 'mb E-VECT: WARM - COLD']; FRAME = [120 270 0 90]; num = 8; top1 = squeeze(wep1_1(num,:,:)); top2 = squeeze(wep1_2(num,:,:)); bot1 = squeeze(cep1_1(num,:,:)); bot2 = squeeze(cep1_2(num,:,:)); units = ['M^2/S^2']; head1 = [num2str(10 * round(lev1(num) / 10)) 'mb E-VECT: WARM']; head2 = [num2str(10 * round(lev1(num) / 10)) 'mb E-VECT: COLD']; FRAME = [0 150 -20 80]; figure(1) sp(1) XAX = lon1; YAX = lat1; gquiv(top1, top2, 2); % mquiv(top1, top2, 1, 'stereo', [90 270]); title([head1 ', TPAC SST']); xlabel(['Units: ' units]); dc sp(2) XAX = lon2; YAX = lat2; gquiv(bot1, bot2, 2); % mquiv(bot1, bot2, 1, 'stereo', [90 270]); title([head2 ', TPAC SST']); xlabel(['Units: ' units]); dc