Documentation of ml_coup_slp_eof

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

Help text

load ML_ANN_sst.mat; slp = sst2;

Cross-Reference Information

Listing of script ml_coup_slp_eof

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load ML_ANN_slp.mat;

%lims = [110 285 20 90];  % SST - NHEM
%lims = [100 255 20 90];  % NPAC
%lims = [100 255 20 75];  % NPAC2
%lims = [-0.1 360 20 90];  % SLP - NHEM
%lims = [-0.1 360 -90 -20];  % SLP - SHEM
%lims = [120 250 20 75];  % NPAC3
lims = [110 260 20 90];  % NPAC4

%slp = getflx('psl', lims, 101:1000);
[lat, lon] = getll('temp', lims);

[xk, yk] = keep_var(lims, lon, lat);

slp = slp(:, yk, xk);
slp = detrend(slp);

for i = 0:2; i
  if i == 0;
    slp2 = slp;
  elseif i == 1;
    [b1, a1] = butter(9, 2/10);
    slp2 = slp - filtfilt(b1, a1, slp);
  elseif i == 2;
    [b2, a2] = butter(9, 2/9);
    slp2 = filtfilt(b2, a2, slp);
  end
  lat2 = lat(yk);
  lon2 = lon(xk);
  
  slp2 = cosweight(slp2, lat2);
  [ntim, nlat, nlon] = size(slp2);
  slp2 = reshape(slp2, ntim, nlat*nlon);
  
  if i == 0; 
    [rlam, rlds, rpcs, rper] = eof_dan(slp2);
  elseif i == 1;
    [hlam, hlds, hpcs, hper] = eof_dan(slp2);
  elseif i == 2;
    [llam, llds, lpcs, lper] = eof_dan(slp2);
  end
end

lat = lat2; lon = lon2;
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs/data
%save ml_slp_eof_npac.mat lims lat lon rlam rlds rpcs rper ...
%save coup_slp_eof_npac.mat lims lat lon rlam rlds rpcs rper ...
%save ml_sst_eof_npac.mat lims lat lon rlam rlds rpcs rper ...
%save ml_slp_eof_nhem.mat lims lat lon rlam rlds rpcs rper ...
%save coup_slp_eof_shem.mat lims lat lon rlam rlds rpcs rper ...
%save ml_slp_eof_npac2.mat lims lat lon rlam rlds rpcs rper ...
%save coup_slp_eof_npac2.mat lims lat lon rlam rlds rpcs rper ...
%save ml_slp_eof_npac3.mat lims lat lon rlam rlds rpcs rper ...
%save coup_slp_eof_npac4.mat lims lat lon rlam rlds rpcs rper ...
save ml_slp_eof_npac4.mat lims lat lon rlam rlds rpcs rper ...
   hlam hlds hpcs hper llam llds lpcs lper




%  Take a look

nlat = length(lat);
nlon = length(lon);

cint = 0.5;

figure(1); fo; clf;
default_global;
subplot(3,2,1);
  mcont3(reshape(lds1(:,1), nlat, nlon), cint, 'stereo', [90 270]);
  drawmapm; gridm on; axis_limits;
subplot(3,2,2);
  mcont3(reshape(lds1(:,2), nlat, nlon), cint, 'stereo', [90 270]);
  drawmapm; gridm on; axis_limits;
subplot(3,2,3);
  mcont3(reshape(lds1(:,3), nlat, nlon), cint, 'stereo', [90 270]);
  drawmapm; gridm on; axis_limits;
subplot(3,2,4);
  mcont3(reshape(lds1(:,4), nlat, nlon), cint, 'stereo', [90 270]);
  drawmapm; gridm on; axis_limits;
subplot(3,2,5);
  mcont3(reshape(llds(:,1), nlat, nlon), cint, 'stereo', [90 270]);
  drawmapm; gridm on; axis_limits;
subplot(3,2,6);
  mcont3(reshape(llds(:,2), nlat, nlon), cint, 'stereo', [90 270]);
  drawmapm; gridm on; axis_limits;