Documentation of np_vs_eq_ceof_look

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

Help text

load LPregmap_3deof_totpac.mat

Cross-Reference Information

Listing of script np_vs_eq_ceof_look

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data
load LPregmap_3deof_12.5to62.5pac.mat
regmapnp = regmap;
load LPregmap_3deof_-30to30pac.mat
regmapeq = regmap;
tem = FRAME;
get_global; default_global; FRAME = tem;
ntim = size(pcs, 1);
[nlev, nlat, nlon] = size(clim);

cd /home/disk/hayes2/dvimont/csiro/data

filin = ['temp_A_L1-10.nc'];
nc = netcdf(filin, 'nowrite');
  lind = 1:10;
  depth = nc{'depth'}(lind);
nc = close(nc);

%  Look at various phases of these patterns

figure(1); clf; figure_landscape;
lag = 0; lg = lag*pi/180;
num = 1; lind = 1;
nfrm = 4;

cint = [.1 0 .1 0 .1 0 .05 0 .025]; 
lind = 1;
for fig = 1:2; figure(fig); figure_landscape;
for i = 1:nfrm;

  if i < 3;
    regmap = regmapnp;
  else
    regmap = regmapeq;
  end

  wgt = conj(exp(j * ((i-1) * pi/nfrm + lg) ));% * ones(npt, 1));  
  tem = squeeze(real(wgt .* regmap(num,lind,:,:)));

     subplot(2,2,i);
       if ismap(gca); clma; end
       mcont(tem, [-1:cint(lind):1], 'giso', [0 mean(FRAME(1:2))]);
       title(['LP:  <TEMP, CPC' num2str(num) '>, Phase = ' num2str(((i-1)*180)/nfrm+lag) ...
              ', Depth = ' num2str(depth(lind)/100) 'm'])
       xlabel(['Contour Interval:  ' num2str(cint(lind)) ' K std^-^1'])

end

%cd /home/disk/tao/dvimont/matlab/CSIRO/Plots2
%eval(['print -dps2 LP_CEOF3D_lev' num2str(lind) '_12.5to62.5pac.ps']);

end