%* map global precipitation data set from Tyndall Center
%
%  precip_monthly.mat   30-year mean monthly precip. over land surface
%                       1961 to 1990; millimeters/day; 0.5 degree graticule
%                       360x720x12
%                       layer 1: January; layer 2: February; etc.
%  plat                 grid latitudes
%  plon                 grid longitudes
%

load globalclimatology
load coast

%* SET map month and figure numbers
mapM=1;     % select index number of month to plot
FNZ=10;     % starting figure number, program uses FNZ and FNZ+1
 
%* ranges for color scales
Cmin=min(min(min(precip)));                 
Cmax=max(max(max(precip)));

 
%* map plot window
wind=[min(plon) max(plon) min(plat) max(plat)]; 
 
 
%********
% world precipitation flat map for specified month
%********
 
figure(FNZ)                            %open figure window
clf
caxis([Cmin Cmax])
axis equal
axis(wind)                             % define axes for this plot
hold on
pcolor(plon, plat, precip(:,:,mapM))   % simple Plate Carree projection
shading flat
xlabel('longitude')
ylabel('latitude')
title('1961 to 1990 mean monthly precipitation (mm/day)')
text(190, 90, int2str(mapM))
colorbar('horiz')
hold on
plot(long, lat, 'r-')