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

Listing of script plot_30s30n_terms

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

   [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 = [110 299 -30 30];
lag = 0;

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

for i = 1:6;
  tem = -1 * squeeze(uptb(i,:,:));
  subplot(6,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Phase = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(6,2,2)
  title([tit ' U'' dTbardx']);

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

%  Plot the vdtdy data

default_global;  XAX = lon_out; YAX = lat_out; FRAME = [110 299 -30 30];
lag = 0;

figure(2); figure_orient;
for i = 1:6;
  tem = -1 * squeeze(vbtp(i,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Phase = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(6,2,1)
  title([tit ' Vbar dT''dy']);

for i = 1:6;
  tem = -1 * squeeze(vptb(i,:,:));
  subplot(6,2,2*i);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Phase = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(6,2,2)
  title([tit ' V'' dTbardy']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 vdtdy_terms_80-270m_CPC1.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;

figure(3); figure_orient;
for i = 1:6;
  tem = -1 * squeeze(wbtp(i,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Phase = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(6,2,1)
  title([tit ' Wbar dT''dz']);

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

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 wdtdz_terms_80-270m_CPC1.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_ceof(hflx, pcs, nfrm);

%  Plot the HC and HFLX data

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

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

for i = 1:6;
  tem = squeeze(dhcdt(i,:,:));
  subplot(6,2,2*i-1);
    gcont(tem, clev);
    dc2(tem);
    ylabel(['Phase = ' 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_0-270m_CPC1.ps




cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
%save HP10_hctend_L1-7_yr101-1000.mat uptb ubtp vptb vbtp wptb wbtp lags tim lev nfrm
%save HP10_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 HP10_hctend_L1-7_yr101-1000.mat
uptb7 = uptb;
ubtp7 = ubtp; 
vptb7 = vptb; 
vbtp7 = vbtp; 
wptb7 = wptb; 
wbtp7 = wbtp;

load HP10_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;