%** this script plots zonal mean climatology data
%
% zonalclimatology.mat 
% annual arrays: 180x1    monthly arrays: 180x12
% Taal_zm       zonal mean annual land surface temperature, degrees C        
% Tal_zm        zonal mean monthly land surface temperature, degrees C
% SST_zm        zonal mean annual sea surface temperature, degrees C
% maplat        latitudes corresponding to Taal, Tal, SST zonal data
% precip_zm     zonal mean monthly precipitation, mm/day
% pzmlat        latitudes corresponding to precip_zm



load zonalclimatology

% strings to plot 12 months with different colors and line types
% and create a legend
pcs=['r- '; 'g- '; 'c- '; 'm- '; 'b- '; 'k- '; 'r-.'; 'g-.'; 'c-.'; 
    'm-.'; 'b-.'; 'k-.'];
lgs=['Jan';'Feb';'Mar';'Apr';'May';'Jun';'Jul';'Aug';'Sep';'Nov';'Dec']

%* plotting windows for Ta and precip
windt=[-90 90 floor(min(min(Tal_zm))) ceil(max(max(Tal_zm)))];
windp=[-90 90 0 ceil(max(max(precip_zm)))];


%***** 
%  plot 12 months of zonal mean temperature over land
%*****
figure(201)
clf
axis(windt)
hold on
for n=1:size(Tal_zm)*[0 1]'
    plot(maplat, Tal_zm(:,n), pcs(n,:))
end
grid on
legend(lgs);
title('zonal mean surface air temperature over land')
xlabel('latitude')
ylabel('temperature (degrees C)')


%*****
%  plot 12 months of zonal mean precipitation over land
%*****
figure(202)
clf
axis(windp)
hold on
for n=1:size(precipl_zm)*[0 1]'
    plot(maplat, precipl_zm(:,n), pcs(n,:))
end
grid on
legend(lgs);
title('zonal mean precipitation rate over land')
xlabel('latitude')
ylabel('precipitation rate (mm/day)')


%***** 
%  plot 12 months of zonal mean temperature over land and ocean
%*****
figure(203)
clf
axis(windt)
hold on
for n=1:size(Ta_zm)*[0 1]'
    plot(maplat, Ta_zm(:,n), pcs(n,:))
end
grid on
legend(lgs);
title('zonal mean surface air temperature')
xlabel('latitude')
ylabel('temperature (degrees C)')


%*****
%  plot 12 months of zonal mean precipitation over land and ocean
%*****
figure(204)
clf
axis(windp)
hold on
for n=1:size(precip_zm)*[0 1]'
    plot(maplat, precip_zm(:,n), pcs(n,:))
end
grid on
legend(lgs);
title('zonal mean precipitation rate')
xlabel('latitude')
ylabel('precipitation rate (mm/day)')