NCL Home > Documentation > Functions > Crop

crop_water_need

Compute "crop water need", the quantity of water needed to attain maximum crop productivity.

Available in version 6.4.0 and later.

Prototype

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl"  ; This library is automatically loaded
                                                      ; from NCL V6.5.0 onward.
                                                      ; No need for user to explicitly load.

	function crop_water_need (
		pet         : numeric,  
		act_evaptr  : numeric,  
		opt     [1] : logical   
	)

Arguments

pet

Potential evapotranspiration. See: thornthwaite

act_evaptr

Actual evapotranspiration: the amount of water that is actually removed from a surface due to the processes of evaporation and transpiration. Same units, size, shape and type as pet.

opt

Currently not used. Set to False.

Return value

An array with the same dimensionality and type as pet containing the desired quantity.

Description

Here, this is simply:

    crop_water_need = pet - act_evaptr
All appropriate metadata is returned.

References: Actual, Reference and Potential Evapotranspiration

   Jeff Andresen, and Aaron Pollyea, Michigan State University Extension
   What is evapotranspiration and why it matters
   WWW: http://msue.anr.msu.edu/news/what_is_evapotranspiration_and_why_it_matters

   C. Brouwer and M. Heibloem (1986)
   Irrigation Water Management: Irrigation Water Needs
   WWW: http://www.fao.org/docrep/s2022e/s2022e00.htm
   
   Pidwirny, M. (2006).
   Actual and Potential Evapotranspiration
   Fundamentals of Physical Geography, 2nd Edition.
   WWW: http://www.physicalgeography.net/fundamentals/8j.html

   Irmak, S. and D.Z. Haman (2014; most recent)
   Evapotranspiration: Potential or Reference?
   WWW: https://edis.ifas.ufl.edu/ae256

See Also

Crop & Evapotranspiration functions

Examples

Example 1: This is a bogus example. The actual evapotranspiration values are randomly created. The actual evapotranspiration may exceed those estimated by the thornthwaite function.

No examples are currently available for this function.

   tmp           = (/23.3, 21.1, 19.6, 17.2, 12.6, 10.9 \
                    ,10.0, 11.0, 13.0, 15.8, 17.8, 20.1/)
   tmp@units     = "degC"
   tmp@long_name = "Sellers, Physical Climatology; 1969, p173"

   latAsp        = -38.0                                 ; Latitude of Aspendale, Australia
   pet           = thornthwaite(tmp, latAsp, False, 0)
   pet@long_name = "PET: Aspendale, Australia"
   pet@units     = "mm/month"                            ; ie, monthly total


   act_evap      = random_uniform(0.75,1.1,dimsizes(tmp))*pet
   cwn           = crop_water_need(pet, act_evap, False)
   printVarSummary(cwn)

   
   print("")
   print(pet+"  "+act_evap+"  "+cwn)

The edited output is:


Variable: cwn
Type: float
Total Size: 48 bytes
            12 values
Number of Dimensions: 1
Dimensions and sizes:	[12]
Coordinates: 
Number Of Attributes: 3
  _FillValue :	9.96921e+36
  long_name :	crop water need: PET-ActualEvap
  units :	mm/month



        PET      ACTUAL     DIFF
(0)    133.829  115.535    18.2947
(1)	96.9449  83.7155   13.2294
(2)	86.9764  80.2518   6.72452
(3)	61.0302  51.9401   9.09012
(4)	34.4376  30.0831   4.35456
(5)	25.0275  20.1192   4.90823
(6)	22.9359  21.7472   1.18868
(7)	29.1031  30.1873  -1.08415
(8)	41.1918  35.5866   5.60525
(9)	64.795   56.3086   8.48633
(10)	82.0458  87.2403  -5.19447
(11)   106.838   92.715   14.1229