%
%       map_erbe.m
%
%   ....This is a simple script to plot ERBE data contained in erbe.mat.
%
%       Written by Von P. Walden
%                  6 April 1998

%
%   ....Loads ERBE data.
%
load erbe.mat

%
%   ....Creates a legend for ERBE maps.
%
erbe_legend = [1/2.5 87.5 0];

%
%   ....Create map figure with Mollweide projection.
%
figure
axesm mollweid

%
%   ....Plot shortwave radiation
%
meshm(sw,erbe_legend)

%
%   ....Add continental outline.
%
load coast
plotm(lat,long,'k')

%
%   ....Adds a title and colorbar.
%
title('ERBE global data set from July 1986, NOAA-9, Shortwave Radiation')
colorbar('horiz')
text(3.1, -2.6,'W m-2')

%
%   ....Repeat above commands for different variable.
%
figure
axesm mollweid
meshm(lw,erbe_legend)
plotm(lat,long,'k')
title('ERBE global data set from July 1986, NOAA-9, Longwave Radiation')
colorbar('horiz')
text(3.1, -2.6,'W m-2')

%
%   ....Repeat above commands for different variable.
%
figure
axesm mollweid
meshm(net,erbe_legend)
plotm(lat,long,'k')
title('ERBE global data set from July 1986, NOAA-9, Net Radiation')
colorbar('horiz')
text(3.1, -2.6,'W m-2')

