% CalCOFI objective maps. - Manu

load CalCOFI.mat salt temp TempSeasM SaltSeasM spice SpiceSeasM ...
                 lon lat z mask3d h year month timemask
		     
% CalCOFI data from  2/1949 to 9/2000 OBJECTIVE MAPS

% Time dependent monthly values
% temp(x,y,z,t)  TEMPERATURE      
% salt(x,y,z,t)  SALINITY
% spice(x,y,z,t) SPICE   
% mask3d(x,y,z)  LANDMASK
% h(x,y)         TOPOTGRAPHY

% Seasonal Means
% TempSeasM(x,y,z,t=1:12)  TEMPERATURE      
% SaltSeasM(x,y,z,t=1:12)  SALINITY
% SpiceSeasM(x,y,z,t=1:12) SPICE


% lon LONGITUDE (2D array)
% lat LATITUDE  (2D array)
% z   DEPTHS    (1D array)

% year(t)     YEAR
% month(t)    MONTH
% timemask(t) IS NaN for those time when OA map does not have 
%             sufficient data coverage

% OA scales: dx, dy = 100 km ; dt = 4 months.

% Reference to OA map and errors is:
%  Di Lorenzo et al., 2005: The warming of the California Current: 
%  Dynamics and ecosystem implications. 
%  Journal of Physical Oceanography, 35 (3), 336-362.
%  e-version: http://o3d.org/manu/papers/DiLorenzo_jpo2003.pdf


% Horizontal plot for SST on JULY 1997
figure;
in=find(year == 1997 & month == 7);
subplot(2,1,1)
rnt_contourfill(lon,lat,temp(:,:,1,in).*mask3d(:,:,1),20);
title 'JULY 1997 SST'
colorbar
subplot(2,1,2)
in=find(year == 1997 & month == 7);
rnt_contourfill(lon,lat,(temp(:,:,1,in)-TempSeasM(:,:,1,7)).*mask3d(:,:,1),20);
title 'JULY 1997 SST ANOMALY'
colorbar
print -djpeg100 SSTHorizontal.jpg

% Section along LINE 90 (note: line 90 is j=8 on this grid)
figure;
j=8;
in=find(year == 1997 & month == 7);
% create 3D mask
c1=repmat(lon(:,j),[1 length(z)]);
c2=repmat(z, [1 size(lon,1)])';

subplot(2,1,1)
tsect=(temp(:,:,:,in)).*mask3d;
tsect=squeeze(tsect(:,j,:));
rnt_contourfill(c1,c2,tsect,20); colorbar;shading flat
title 'JULY 1997 TEMP'
subplot(2,1,2)
tano=(temp(:,:,:,in)-TempSeasM(:,:,:,7)).*mask3d;
tano=squeeze(tano(:,j,:));
rnt_contourfill(c1,c2,tano,20); colorbar;shading flat
title 'JULY 1997 TEMP ANOMALY'
print -djpeg100 TEMPVerticalSection.jpg