%  Routine for defining the Aleutian Low.  
%
%  Note that you will have to download some .m files from:
%  http://www.aos.wisc.edu/~dvimont/matlab
%  
%  Also, you will need to have the NetCDF toolboxes installed for MATLAB,
%  which is currently being done on the 14th floor computers.  I've
%  deposited the SST and SLP data in ~dvimont/matlab/Data on the 14th
%  floor computers if you would like to do this yourself.
%

clear

%  Start by loading SST data.  Define region:
lims = [120 270 20 75];
lev = 1;
tim = get_time(1948, 2003, 1948);

%  Load the data - note that this uses getnc2, which is avialable via the
%  website above.  You need the netcdf toolboxes to use this, so you can
%  also load the data from ~dvimont/matlab/Data

cd /Volumes/Data/NCEP/monthly   %  You will have to replace this line
[slp, lat, lon] = getnc2('slp.mon.mean.nc', 'slp', lims, 1, tim);

%  Now, I'll save the data to ~dvimont/matlab/Data (don't do this)
%
%  cd ~dvimont/matlab/Data
%  save slp.aleut.mat slp lat lon lims lev tim
%
%  If you can't run the getnc2 command above, then load the sst data from
%  my home account.

%  Remove the annual cycle
[slp, slpc] = remove_anncyc(slp);

%  Weight data by sqrt(cos(lat)):
slpw = cosweight(slp, lat);

%  Now, perform EOF analysis
[lam, lds, pcs, per] = eof_dan(slpw);

%  We'll define the first PC as the Aleutian Low.

%  Now, generate regression map onto that PC
[regmap, corrmap] = regress_eof(slp, pcs(:,1));

%  And finally, plot it up nicely.  Again, you'll need a number of
%  routines from the website above.

%  Set up figure window
figure_tall(1); clf;

%  Define globals - YAX, XAX, and FRAME.  These are used by my mapping
%  routines.  
global_latlon(lat, lon, lims);

%  Now, set up subplot in the current window.
global_axes(6, 2.5, .5, .5, 1.5);
subplot2(1,1); cla;

%  Set up a map axis for the mapping toolbox - this makes pretty plots.
map_axis('giso');
map_contour_pn(regmap, 0.5, 'zero');

%  OK!  finish it up:
fill_landmap('under');
setm(gca, 'MLineLocation', 10, 'PLineLocation', 15);
gridm on; tightmap;
xlabel('Contour:  0.5 mb', 'fontsize', 12);
title('\bf Aleutian Low regression map', 'fontsize', 14);

%  Print this to a file
cd ~/Class/aos575/Homework/HW4
print -depsc2 -painters Aleut_regmap.eps