%** this script produces maps of surface
% temerature in both marine and terrestrial locations
%
% globalclimatology.mat  relevant variables
% Ta        monthly mean temperature in degrees C, 
%           0.5 degree gridded data  360x720x12
%           layer 1: January; layer 2: February; etc.  
% Ta_ann    mean annual surface air temperature, degrees C, from Ta              
% Tlat      grid latitudes                
% Tlon      grid longitudes
%
% references
% Legates, D. R. and C. J. Willmott, 1990, Mean Seasonal 
%   and Spatial Variability Global Surface Air Temperature, 
%   Theoretical and Applied Climatology , 41, 11-21.


load globalclimatology
load coast
 
%* SET map month and figure numbers
mapT=1;     % select index number of month to plot
FN=15;      % starting figure number
 
%* ranges for color scales
Ctmin=min(min(min(Ta)));                 
Ctmax=max(max(max(Ta)));

 
%* map plot window
wind=[min(Tlon) max(Tlon) min(Tlat) max(Tlat)]; 
 
%********
% surface air temperature flat map for specified month
%********
 
figure(FN)                      % open figure window
clf
caxis([Ctmin Ctmax])            % set color range for this plot
axis equal
axis(wind)                      % define axes for this plot
hold on
pcolor(Tlon, Tlat, Ta(:,:,mapT))   % simple Plate Carree projection
shading flat
xlabel('longitude')
ylabel('latitude')
title('mean monthly surface air temperature (degrees C)')
text(190, 90, int2str(mapT))
colorbar('horiz')
hold on
plot(long, lat, 'k-')