unadjustedPE.m

function[PEhot, Thot] = unadjustedPE
%% Returns a table of unadjusted PE values (mm/day) for temperatures over
% 26.5 C from Thornthwaite and Mather Table 5, page 226.
%
% [PE, T] = unadjustedPE
%
% ----- Outputs -----
%
% PE: A vector set of unadjusted PE values for hot temperatures
%
% T: The tempertaures corresponding to the PE values.

Thot = 26:.1:38;
PEhot = [...
      4.5 4.5 4.5 4.5 4.5 4.5 4.5 4.6 4.6 4.6 ...
      4.6 4.7 4.7 4.7 4.8 4.8 4.8 4.8 4.9 4.9 ...
      4.9 5.0 5.0 5.0 5.0 5.1 5.1 5.1 5.1 5.2 ...
      5.2 5.2 5.2 5.2 5.3 5.3 5.3 5.3 5.4 5.4 ...
      5.4 5.4 5.4 5.5 5.5 5.5 5.5 5.5 5.6 5.6 ...
      5.6 5.6 5.6 5.6 5.7 5.7 5.7 5.7 5.7 5.8 ...
      5.8 5.8 5.8 5.8 5.8 5.8 5.9 5.9 5.9 5.9 ...
      5.9 5.9 5.9 5.9 6.0 6.0 6.0 6.0 6.0 6.0 ...
      6.0 6.0 6.0 6.0 6.1 6.1 6.1 6.1 6.1 6.1 ...
      6.1 6.1 6.1 6.1 6.1 6.1 6.1 6.1 6.1 6.1 ...
      6.1 6.1 6.2 6.2 6.2 6.2 6.2 6.2 6.2 6.2 ...
      6.2 6.2 6.2 6.2 6.2 6.2 6.2 6.2 6.2 6.2 ...
      6.2 ];
  
end