Documentation of plot_wdtdz_terms

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

Help text

  Load PCS:

Cross-Reference Information

Listing of script plot_wdtdz_terms

clear

for biff2 = 2:3;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
  load LP10_L1-7_CEOF.mat; tit = 'Lowpass Filtered Data ( > 10 Years )';
  ptit = 'LP10';
elseif biff2 == 2;
  load HP8_L1-7_CEOF.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
  ptit = 'HP8';
elseif biff2 == 3;
  load RAW_L1-7_CEOF.mat; tit = 'Unfiltered Data';
  ptit = 'RAW';
end

cd ~/matlab/CSIRO/Heat/Old_routines

lims = [108 302 -62 62];
lev = 1:7;
tim = 147:550;
nfrm = 6;

%  Start with wbar dtprime dz

   [wdhdz, lat_out, lon_out, depth_out] = ...
           get_wbar_dtprimedz(pcs(tim-100,:), lims, nfrm, tim, lev);

%  Plot the data

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

figure(1); figure_orient;
cint = 1; clev = [-50:cint:-cint cint:cint:50];
for i = 1:nfrm
  tem = -1 * squeeze(wdhdz(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem)
    title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i > 4;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3);
  ylabel([tit ':  Vertical Advection of HC'' by Wbar;  '...
          'Depth = 0:' num2str(depth_out(max(lev))) 'm']);

if 1;
  cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv
  eval(['print -dps2 ' ptit '_wbar_dHCprimedz_yr1.ps'])
end



%  Next, do wprime dtbar dz

cd ~/matlab/CSIRO/Heat/Old_routines
   [wdhdz, lat_out, lon_out, depth_out] = ...
           get_wprime_dtbardz(pcs(tim-100,:), lims, nfrm, tim, lev);

%  Plot the data

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

figure(2); figure_orient;
cint = 1; clev = [-50:cint:-cint cint:cint:50];
for i = 1:nfrm
  tem = -1 * squeeze(wdhdz(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem)
    title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
  if i > 4;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3);
  ylabel([tit ':  Vertical Advection of HCbar by W'';  '...
          'Depth = 0:' num2str(depth_out(max(lev))) 'm']);

if 1;
  cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Old_routines/Plot_HCadv
  eval(['print -dps2 ' ptit '_wprime_dHCbardz_yr1.ps']);
end

end