%Estimate the Hargreaves Samani ETo for Puerto Rico

%Load insolation and elevation files


%PLOTTING SCRIPT
res=0.0179856/2; %Pixel resolution
latlim=[17.8 18.7]; %Latitudes Limits
lonlim=[-67.3 -65.42]; %Longitudes Limits
[Lonm,Latm]=meshgrid(lonlim(1):res:lonlim(2),latlim(2):-res:latlim(1));
Latm = flipud(Latm); % For WorldMap Lats and Observations

%plot elevation
figure (1)
title('Elevation (m)')
worldmap(latlim,lonlim);
surfacem(Latm,Lonm,ele);
colorbar('SouthOutside')
saveas(gcf,'ELE.jpg')

%plot the solar insolation



%ESTIMATE, PLOT and SAVE Matlab file of the surface temperature for January.  Using regression
%equation for average temperature for January from Goyal et al., 1988. 




%ESTIMATE, PLOT and SAVE Matlab file for Reference ET using equatino 1 in Hargreaves and Samani
%(1985).  Be sure to convert the solar radiation from MJ/m2/day to mm of
%water.

