Etp: Computes Potential Evapo-Transpiration (PET)

In r2dRue: 2d Rain Use Efficience model

Description

PET (mm) is computed from mean maximum, mean and mean minimum air temperatures in Celsius degree and extra-terrestrial solar radiation (cal/cm2/day) using the Hargreaves-Samani equation.

batchPet () provides a batch process to compute PET over sets of images, typically time-series, specified as raster groups.

Usage

petHgsm(Tmin,Tmax,Tmed,Rad,month)
batchPetHgsm(outFl, monthIni, Tmin, Tmed, Tmax, Rad, ...)

Arguments

monthIni	The month order-number of the first image in the time-series. An integer between 1 and 12.
month	The target month order-number to compute PET An integer between 1 and 12.
Tmin	In petHgsm, a SpatialGridDataFrame with mean minimum air temperature in Celsius degrees. In batchpetHgsm, a raster group with the file names containing mean minimum air temperature.
Tmax	In petHgsm a SpatialGridDataFrame with mean maximum air temperature in Celsius degrees. In batchpetHgsm, a raster group with the file names containing mean maximum air temperature.
Tmed	In petHgsm a SpatialGridDataFrame with mean air temperature (Celsius degree). In batchpetHgsm, a raster group with the file names containing mean air temperature.
Rad	In petHgsm a SpatialGridDataFrame with extra-terrestial solar radiation (cal/cm2/day). In batchpetHgsm a raster group with the file names containing extra-terrestial solar radiation.
outFl	A raster group with file names to save the results
...	Any unmatched parameter will be passed to the writeGDAL routine.

Value

petHgsm returns a spatialGridDataFrame representing the potential evapotranspiration. batchEtp writes to disc the computed rasters as SpatialGridDataFrame objects.

Author(s)

Gabriel del Barrio, Marieta E. SanJuan & Alberto Ruiz

References

Estimating potential evapotranspiration. Hargreaves, G.H., Samani, Z.A. Journal of the Irrigation & Drainage Division - ASCE Volume 108, Issue IR3, 1982, Pages 225-230

Examples

#not run
## Read  the temperature and radiation data 
#tmin=readGDAL('Tmin_Ene1990')
#tmax=readGDAL('Tmax_Ene1990')
#tmea=readGDAL('Tmea_Ene1990')
#rad=readGDAL('SolarRad_Ene1990')
## Calc the PET for Ene 1990
#PetEne=petHgsm(tmin,tmax,tmea,rad, 1)
#writeGDAL(PetEne,'PetEne1990')
## Or to calculate and write to disk, in a batch process, a set of PETs
## create the Input file list
#Flmin=c('Tmin_Ene1990','Tmin_Feb1990', 'Tmin_Mar1990')
#Flmax=c('Tmax_Ene1990','Tmax_Feb1990', 'Tmax_Mar1990')
#Flmea=c('Tmea_Ene1990','Tmea_Feb1990', 'Tmea_Mar1990')
#Flrad=c('Rad_Ene','Rad_Feb', 'Rad_Mar')
#PETFileNames=c('PET_Ene1990','PET_Feb1990', 'PET_Mar1990')
##define the date of the first image
#iniDate=c(1,1990)
#batchPET(iniDate, Flmin, Flmax, Flmea, Flrad, PETFileNames)

r2dRue documentation built on May 1, 2019, 10:10 p.m.