NCL Home > Documentation > Functions > Meteorology, Crop

refevt_hargreaves_fao56

Use the Hargreaves ETo equation to derive reference evapotranspiration as described in FAO 56.

Available in version 6.4.0 and later.

Prototype

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl"

	function refevt_hargreaves_fao56 (
		tmin       : numeric,  
		tmax       : numeric,  
		rex        : numeric,  
		iounit [3] : integer   
	)

Arguments

tmin
tmax

Scalars or arrays of the same size, shape and size containing minimum and maximum temperatures. The units may be derees Celcius, Kelvin or Farenheit. See the iounit argument.

rex

A variable containing extraterrestrial radiation (mm/day. It must have the same dimension sizes as the tmin and tmax variables. See radext_fao56

iounit

An integer array of length 3 indicating the units of input arguments tmin, tmax, rex and the units of the returned argument.

  • iounit(0)=0 ; input tmin and tmax are in degrees C (degC).
  • iounit(0)=1 ; input tmin and tmax are in degrees K (degK).
  • iounit(0)=2 ; input tmin and tmax are in degrees F.
  • iounit(1)=0 ; input radext are in mm/day
  • iounit(1)=1 ; input radext are in MJ/(m2-day)
  • iounit(1)=2 ; input radext are in W/m2
  • iounit(2)=0 ; output units => mm/day
  • iounit(2)=1 ; output units => MJ/(m2-day)
  • iounit(2)=2 ; output units => W/m2

Return value

An array the same size as tmin containg the estimated evapotranspiration.

Description

Estimates evapotranspiration vis the simple "Hargreaves ETo" equation as described in the Food and Agriculture Organization (FAO) Irrigation and Drainage Paper 56 entitled: Crop evapotranspiration - Guidelines for computing crop water requirement . Specifically, see equation 52 of Chapter 3.

See Also

Crop & Evapotranspiration functions

Examples

Example 1: See Example 20 in Chapter 4 for Lyon, France in mid July. The FAO-56 example 20 returns an estimate of 5.0 mm/day using the Hargreaves method. For illustration, the refevt_turc estimates are shown.. 


  jday   = 196               ; 15 July
  lat    = 45.72             ; Lyon, France
  ra_0   = radext_fao56(jday, lat, 0)   ;  16.5463 mm/day 
  ra_1   = radext_fao56(jday, lat, 1)   ;  40.5546 MJ/(m2-day) 
  ra_2   = radext_fao56(jday, lat, 2)   ; 469.378  W/m2   

  print(ra_0 )          ; print with meta data
  print(ra_1 )
  print(ra_2 )

  tmin   = 14.8              ; degC
  tmax   = 26.6
                             ; HARGREAVES
  evtH_0 = refevt_hargreaves_fao56( tmin, tmax, ra_0, (/0,0,0/) ) ;   5.03 mm/day
  evtH_1 = refevt_hargreaves_fao56( tmin, tmax, ra_1, (/0,1,1/) ) ;  12.33 MJ/(m2-day) 
  evtH_2 = refevt_hargreaves_fao56( tmin, tmax, ra_2, (/0,2,2/) ) ; 141.92 W/m2   

  print("evtH_0="+evtH_0)   ;   5.03 mm/day         This agrees with Example 20, Chapter 4
  print("evtH_1="+evtH_1)   ;  12.33 MJ/(m2-day) 
  print("evtH_0="+evtH_2)   ; 141.92 W/m2   

  ---------------------------------------------------------------------
                            ; TURC (requires 'radsol')
  ---------------------------------------------------------------------
                                        ; http://www.lyon.climatemps.com/
  nsun   = 9.2                          ; 'observed' hrs/sun day for Lyon, Juy
  sunhrx = daylight_fao56(jday, lat)    ; max daylight/sun; hr per day
  print(sunhrx)                         ; 15.2  (matches web site)

  rs_0   = radsol_fao56(ra_0, sunhrx, nsun, (/0,0/), False)   ;    9.15 mm/day
  rs_1   = radsol_fao56(ra_1, sunhrx, nsun, (/1,1/), False)   ;   22.44 MJ/(m2-day)
  rs_2   = radsol_fao56(ra_2, sunhrx, nsun, (/2,2/), False)   ; 2596.72 W/m2

  print(rs_0 )                          ; print with meta data
  print("rs_1="+rs_1 )
  print("rs_2="+rs_2 )

  tmean     = (tmin+tmax)*0.5           ; Turc requires 'tmean'

  evturc_0  = refevt_turc( tmean, rs_0, (/0,0,0/) ) ;    7.87 mm/day
  evturc_1  = refevt_turc( tmean, rs_1, (/0,1,1/) ) ;   19.31 MJ/(m2-day) 
  evturc_2  = refevt_turc( tmean, rs_2, (/0,2,2/) ) ; 223.692 W/m2  

  print(evturc_0)               ;    7.87 mm/day    will print with meta data
  print("evturc_1="+evturc_1)   ;   19.31 MJ/(m2-day) 
  print("evturc_2="+evturc_2)   ;  223.69 W/m2