Documentation of look_dutdx

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

Help text

  Load CPCs

Cross-Reference Information

Listing of script look_dutdx

clear
lims = [105 310 -63 63];
lev = 1:7;
tim = 101:550;
nfrm = 6;
[lat, lon, depth] = getll('u', lims);
biff2 = 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 /home/disk/tao/dvimont/matlab/CSIRO/Heat

[dubartdx, dtbarudx, lon1, lat1] = get_dutdx(nfrm, pcs, lims, lev, tim);
[dvbartdy, dtbarvdy, lon2, lat2] = get_dvtdy(nfrm, pcs, lims, lev, tim);

%  Plot the data

lag = 0;
for biff = 4:4;
  if biff == 1;
    get_global; FRAME = [110 299 -60 60]; XAX = lon1; YAX = lat1;
    tem1 = dubartdx;
    tit2 = ['Convergence of Ubar*T'';  '];
    ptit = 'dubarHCdx';
  elseif biff == 2;
    get_global; FRAME = [110 299 -60 60]; XAX = lon1; YAX = lat1;
    tem1 = dtbarudx;
    tit2 = ['Convergence of U''*Tbar;  '];
    ptit = 'dHCbarudx';
  elseif biff == 3;
    get_global; FRAME = [110 299 -60 60]; XAX = lon2; YAX = lat2;
    tem1 = dvbartdy;
    tit2 = ['Convergence of Vbar*T'';  '];
    ptit = 'dvbarHCdy';
  elseif biff == 4;
    get_global; FRAME = [110 299 -60 60]; XAX = lon2; YAX = lat2;
    tem1 = dtbarvdy;
    tit2 = ['Convergence of V''*Tbar;  '];
    ptit = 'dHCbarvdy';
  end

  figure(biff); figure_orient;
  cint = 10; clev = [-500:cint:-cint cint:cint:500];
  for i = 1:nfrm
    tem = squeeze(tem1(i,:,:));
    subplot(3,2,i);
      gcont(tem, clev);
      dc2(tem);
      title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
      box on
    if i > 4;
      xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
    end
  end
  subplot(3,2,3);
    ylabel([tit ':  ' tit2 ' Depth = 0:' num2str(depth(max(lev))) 'm']);

  cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Plot_HCadv
%  eval(['print -dps2 ' ptit '_yr1.ps']);

end