Documentation of plot_60s60n_terms

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

Help text

  Load PCS:

Cross-Reference Information

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

Listing of script plot_60s60n_terms


clear

biff2 = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
%  load LP10_L1-7_CEOF.mat; tit = 'LP10';
  load BP10-60yr_L1-7_EOF.mat; tit = 'BP10-60';
elseif biff2 == 2;
  load HP8_L1-7_CEOF.mat; tit = 'HP8';
elseif biff2 == 3;
  load RAW_L1-7_CEOF.mat; tit = 'RAW';
end

cd ~/matlab/CSIRO/Heat/Old_routines

nfrm = [-6:5];%6;
tim = 101:550;
lev = 1:3; lev2 = 1:7
lims = [106 312 -63 63];

%  Get ubar_dHCprimedx and vbar_dHCprimedy

   [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);

   [ubtp2, vbtp2, lat_out, lon_out, depth_out] = ...
           get_ubar_gradtprime(pcs, lims, nfrm, tim, lev2);

   [uptb2, vptb2, lat_out, lon_out, depth_out] = ...
           get_uprime_gradtbar(pcs, lims, nfrm, tim, lev2);

ubtp = ubtp2 - ubtp;
uptb = uptb2 - uptb;
vbtp = vbtp2 - vbtp;
vptb = vptb2 - vptb;

%  Plot the data

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

figure(1); figure_orient;
cint = 0.5; clev = [-10:cint:-cint cint:cint:10];
for i = 1:6;
  tem = -1 * squeeze(ubtp(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated Ubar dT''dx']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_ubtp_terms_80-270m.ps

for i = 1:6;
  tem = -1 * squeeze(uptb(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated U'' dTbardx']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_uptb_terms_80-270m.ps


%  Plot the vdtdy data

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

figure(2); figure_orient;
for i = 1:6;
  tem = -1 * squeeze(vbtp(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated Vbar dT''dy']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_vbtp_terms_80-270m.ps

for i = 1:6;
  tem = -1 * squeeze(vptb(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated V'' dTbardy']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_vptb_terms_80-270m.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);

   [wbtp2, lat_out2, lon_out] = ...
        get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev2);
   [wptb2, lat_out2, lon_out] = ...
        get_wprime_dtbardz(pcs, lims, nfrm, tim, lev2);

wbtp = wbtp2 - wbtp;
wptb = wptb2 - wptb;

%  Plot the wdtdz data

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

figure(3); figure_orient;
for i = 1:6;
  tem = -1 * squeeze(wbtp(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated Wbar dT''dy']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_wbtp_terms_80-270m.ps

for i = 1:6;
  tem = -1 * squeeze(wptb(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated W'' dTbardy']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_wptb_terms_80-270m.ps


%  Plot the heat flux and actual heat content terms

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

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat
[heat, lath, lonh, depthh, middepth] = getheat(lev, tim, lims);
[heat2, lath, lonh, depthh, middepth] = getheat(lev2, tim, lims);
heat = heat2-heat;

hflx_reg = regress_eof(hflx, pcs, nfrm);
heat_reg = regress_eof(heat, pcs, nfrm);

%  Plot the HC and HFLX data

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

figure(4); figure_orient;
for i = 1:6;
  tem = squeeze(hflx_reg(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated Heat Flux (Pos Down)']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_hflx_80-270m.ps

figure(3); figure_orient(2);
cint = 0.25; clev = [-10:cint:10];
for i = 1:6;
  tem = 1e-8*squeeze(heat_reg(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' x 10^8 J m^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated Heat Content']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_HC_80-270m.ps


%  Plot sfc temp and zonal wind

[temp, taux] = getnc('temp', 'taux', lims, 8, tim);
[lat, lon, depth] = getll('temp', lims);
[latx, lonx, depthx] = getll('taux', lims);

[treg1, treg] = regress_eof(temp, pcs, nfrm);
xreg = regress_eof(taux, pcs, nfrm);

%  Plot the HC and HFLX data

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

for biff = 1:3;
cint = 0.025; clev = [-1:cint:1];
figure(biff); figure_orient;
for i = [1:6]+6*(biff-1);
  tem = squeeze(treg1(i,:,:));
  subplot(3,2,rem(i-1,6)+1);
    gcont(tem, clev);
    dc2(tem);
    title(['Lag ' num2str(nfrm(i))]);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);
  end
end
end

subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated L1 Temperature']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_temp_80-270m.ps

default_global;  XAX = lonx; YAX = latx; FRAME = [110 299 -60 60];
cint = 1; clev = [-10:cint:10];
for i = 1:6;
  tem = 100*squeeze(xreg(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' x 10^2 dynes cm^-^2']);
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated SFC Zonal Wind Stress']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_taux_80-270m.ps


%  Plot ocean velocities

[u, v] = getnc('u', 'v', lims, lev2, tim);
[lat, lon, depth] = getll('u', lims);
[latw, lonw, depthw] = getll('wl', lims);
dz = diff([0; depthw(lev2)]);

ureg = regress_eof(u, pcs, lags);
vreg = regress_eof(v, pcs, lags);
szu = size(u);

ureg = dz' * reshape(shiftdim(ureg, 1), [szu(2) nfrm*prod(szu(3:4))]) ./ sum(dz);
vreg = dz' * reshape(shiftdim(vreg, 1), [szu(2) nfrm*prod(szu(3:4))]) ./ sum(dz);

ureg = reshape(ureg, [szu(3:4) nfrm]);
ureg = shiftdim(ureg, 2);
vreg = reshape(vreg, [szu(3:4) nfrm]);
vreg = shiftdim(vreg, 2);

ureg3 = ureg - ureg2;
vreg3 = vreg - vreg2;

%  Plot u, v

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

figure(2); figure_orient;
for i = 1:6;
  utem = squeeze(ureg(i,:,:));
  vtem = squeeze(vreg(i,:,:));
  subplot(3,2,i);
    gquiv(utem, vtem, -0.125, 2, 'cm s^-^1');
%    gcont(utem, [-1:.025:1]);
    dc2(utem);
%    ylabel(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
    title(['Lag = ' num2str(lags(i))]);
  if i < 5;
    xlabel('');
  end
end
subplot(3,2,3)
  ylabel([tit ':  80-270m Integrated Ocean Velocities']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv/LP10_ttend
%print -dps2 LP10_oce_vel_80-270m.ps