Documentation of plot_30s30n_terms

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

  Get ubar_dHCprimedx and vbar_dHCprimedy

Cross-Reference Information

This script calls
dc2 default_global figure_orient frame gcont get_dhcdt get_ubar_gradtprime get_uprime_gradtbar get_wbar_dtprimedz get_wprime_dtbardz getll getnc regress_eof subplot

Listing of script plot_30s30n_terms


clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat
cd ~/matlab/CSIRO/Heat/Old_routines
nfrm = -8:8; lags = nfrm;
tim = 101:1000;
lev = 1:3;
lims = [106 312 -31 31];
tit = '80-270m LP9';
tit = '0-270m LP9';

   [ubtp, vbtp, lat_out, lon_out, depth_out] = ...
           get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);

   [uptb, vptb, lat_out, lon_out, depth_out] = ...
           get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);

%  Plot the data

[lat_out, lon_out, depth_out] = getll('u', lims);
default_global;  XAX = lon_out; YAX = lat_out; FRAME = [105 299 -30 30];
lag = 0;

tem1 = -1*ubtp;
tem1 = -1*uptb;
figure(1); figure_orient(1);
cint = 0.5; clev = [-10:cint:-cint cint:cint:10];
for i = 1:6;
  tem = -1 * squeeze(tem1(i+2,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+2))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
  tem = -1 * squeeze(tem1(i+8,:,:));
  subplot(6,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+8))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
for i = 1:2;
  subplot(6,2,i)
    title([tit ' Ubar dT''dx']);
    title([tit ' U'' dTbardx']);
end

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
%print -dps2 ubtp_terms_0-270m_LP9.ps
%print -dps2 uptb_terms_0-270m_LP9.ps

%  Plot the vdtdy data

default_global;  XAX = lon_out; YAX = lat_out; FRAME = [110 210 -20 20];
lag = 0;
figure(2); figure_orient(1);

tem1 = -1*vbtp;
tem1 = -1*vptb;
cint = 0.5; clev = [-10:cint:-cint cint:cint:10];
for i = 1:6;
  tem = -1 * squeeze(tem1(i+2,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+2))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
  tem = -1 * squeeze(tem1(i+8,:,:));
  subplot(6,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+8))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
for i = 1:2;
  subplot(6,2,i)
    title([tit ' Vbar dT''dy']);
    title([tit ' V'' dTbardy']);
end

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
%print -dps2 vbtp_terms_0-270m_LP9.ps
%print -dps2 vptb_terms_0-270m_LP9.ps

%  Get the wdtdz terms

cd ~/matlab/CSIRO/Heat/Old_routines

   [wbtp, lat_out, lon_out] = ...
        get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
   [wptb, lat_out, lon_out] = ...
        get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);

%  Plot the wdtdz data

[lat_out, lon_out, depth_out] = getll('wl', lims);
default_global;  XAX = lon_out; YAX = lat_out; FRAME = [110 299 -30 30];
lag = 0;

default_global;  XAX = lon_out; YAX = lat_out; FRAME = [110 299 -30 30];
lag = 0;
figure(3); figure_orient(1);

tem1 = -1*wbtp;
tem1 = -1*wptb;
cint = 0.5; clev = [-10:cint:-cint cint:cint:10];
for i = 1:6;
  tem = -1 * squeeze(tem1(i+2,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+2))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
  tem = -1 * squeeze(tem1(i+8,:,:));
  subplot(6,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+8))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
for i = 1:2;
  subplot(6,2,i)
    title([tit ' Wbar dT''dz']);
    title([tit ' W'' dTbardz']);
end

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
%print -dps2 wbtp_terms_0-270m_LP9.ps
%print -dps2 wptb_terms_0-270m_LP9.ps



%  Plot the heat flux and actual heat content terms

hflx = getnc('heat', lims, lev, tim);
[lat, lon, depth] = getll('heat', lims);

cd ~/matlab/CSIRO/Heat
[dhcdt, lat_out, lon_out, tim2] = get_dhcdt(pcs, lims, nfrm, tim, lev);

hflx_reg = regress_eof(hflx, pcs, nfrm);

%  Plot the HC and HFLX data

default_global;  XAX = lon; YAX = lat; FRAME = [110 299 -30 30];
lag = 0;

cint = 0.5; clev = [-10:cint:-cint cint:cint:10];
figure(4); figure_orient;
for i = 1:6;
  tem = squeeze(hflx_reg(i+2,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+2))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(6,2,2)
  title(['LP9 Heat Flux (Pos Down)']);

for i = 1:6;
  tem = squeeze(dhcdt(i,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(6,2,1)
  title([tit ' 80-270m HC Tendency']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 HFLX_dHCdt_80-270m_CPC1.ps




cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
save LP9_hctend_L1-3_yr101-1000.mat uptb ubtp vptb vbtp wptb wbtp...
     dhcdt3 hflx_reg lags tim lev nfrm
%save LP10_hctend_L1-3_yr101-1000.mat uptb ubtp vptb vbtp wptb wbtp lags tim lev nfrm


cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_hctend_L1-7_yr101-1000.mat
uptb7 = uptb;
ubtp7 = ubtp; 
vptb7 = vptb; 
vbtp7 = vbtp; 
wptb7 = wptb; 
wbtp7 = wbtp;

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_hctend_L1-3_yr101-1000.mat
uptb3 = uptb;
ubtp3 = ubtp; 
vptb3 = vptb; 
vbtp3 = vbtp; 
wptb3 = wptb; 
wbtp3 = wbtp;


uptb = uptb7 - uptb3;
ubtp = ubtp7 - ubtp3; 
vptb = vptb7 - vptb3; 
vbtp = vbtp7 - vbtp3; 
wptb = wptb7 - wptb3; 
wbtp = wbtp7 - wbtp3;


dhcdt7 = dhcdt;
dhcdt3 = dhcdt;
dhcdt = dhcdt7 - dhcdt3;

%  Get residual

[nfrm, nlat, nlon] = size(hflx_reg);
[lat, lon] = getll('u', lims);
[latw, lonw] = getll('wl', lims);

for i = 1:nfrm;
  tem1(i,:,:) = interp2(lonw, latw, squeeze(dhcdt7(i,:,:)), lon', lat);
  tem2(i,:,:) = interp2(lonw, latw, squeeze(hflx_reg(i,:,:)), lon', lat);
  tem3(i,:,:) = interp2(lonw, latw, squeeze(wptb7(i,:,:)), lon', lat);
  tem4(i,:,:) = interp2(lonw, latw, squeeze(wbtp7(i,:,:)), lon', lat);
end
dhcdt7 = tem1;
hflx_reg = tem2;
wptb7 = tem3;
wbtp7 = tem4;
res7 = dhcdt7 + (ubtp7 + vbtp7 + wbtp7 + uptb7 + vptb7 + wptb7 - hflx_reg);

for i = 1:nfrm;
  tem1(i,:,:) = interp2(lonw, latw, squeeze(dhcdt3(i,:,:)), lon', lat);
%  tem2(i,:,:) = interp2(lonw, latw, squeeze(hflx_reg(i,:,:)), lon', lat);
  tem3(i,:,:) = interp2(lonw, latw, squeeze(wptb3(i,:,:)), lon', lat);
  tem4(i,:,:) = interp2(lonw, latw, squeeze(wbtp3(i,:,:)), lon', lat);
end
dhcdt3 = tem1;
%hflx_reg = tem2;
wptb3 = tem3;
wbtp3 = tem4;
res3 = dhcdt3 + (ubtp3 + vbtp3 + wbtp3 + uptb3 + vptb3 + wptb3 - hflx_reg);

res47 = res7 - res3;

default_global;  XAX = lon; YAX = lat; FRAME = [110 299 -30 30];
lag = 0;

cint = 0.5; clev = [-5:cint:-cint cint:cint:5];
figure(3); figure_orient;
for i = 1:7;
  tem = squeeze(res47(i,:,:));
  subplot(7,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i))]);
  if i == 7;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
  tem = squeeze(res47(i+7,:,:));
  subplot(7,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i+7))]);
  if i == 7;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(7,2,1)
  title(['0-80m LP10 Residual']);

for i = 1:7;
  tem = squeeze(res47(i,:,:));
  subplot(7,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Lag = ' num2str(lags(i))]);
  if i == 7;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(7,2,2)
  title(['80-270m HP10 Residual']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 Residual_0-80m_80-270m.ps