% globaltopo.mat    Earth surface elevation, meters
%   topoE       present-day 180 x 360 matrix, 1 degree grid
%   maplat      grid latitude 
%   maplon      grid longitude
%
% 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


load globaltopo
load globalclimatology

goodies=find(topoE>geoidE);
Tplot=Ta_ann(1:2:360,1:2:720);

figure(100)
clf
plot(Tplot(goodies), topoE(goodies), 'b.')
hold on
title('mean annual atmospheric temperature lapse rate')
xlabel('temperature (C)')
ylabel('altitude (m)')