## DESCRIPTION

i.evapo.time_integration integrates ETa in time following a reference ET (typically) from a set of meteorological stations dataset. Inputs:

- ETa images
- ETa images DOY (Day of Year)
- ETo images
- ETo DOYmin as a single value

Method:

1. each ETa pixel is divided by the same day ETo and become ETrF
2. each ETrF pixel is multiplied by the ETo sum for the representative days
3. Sum all n temporal [ETrF*ETo_sum] pixels to make a summed(ET) in [DOYmin;DOYmax]

representative days calculation: let assume i belongs to range [DOYmin;DOYmax]

    DOYbeforeETa[i] = ( DOYofETa[i] - DOYofETa[i-1] ) / 2
    DOYafterETa[i] = ( DOYofETa[i+1] - DOYofETa[i] ) / 2

## NOTES

ETo images preparation: If you only have one meteorological station data, the easiest way is:

    n=0
    for ETo_val in Eto[1] Eto[2] ...
    do
    	r.mapcalc "eto$n = $ETo_val"
    	`expr n = n + 1'
    done

with Eto[1], Eto[2], etc being a simple copy and paste from your data file of all ETo values separated by an empty space from each other. If you have several meteorological stations data, then you need to grid them, Thiessen polygons or interpolation for each day. For multi-year calculations, just continue incrementing DOY values above 366, it will continue working, up to maximum input of 400 satellite images.

## SEE ALSO

i.eb.eta, i.evapo.potrad, i.evapo.senay, r.surf.idw, r.surf.idw2, r.resamp.interp

## AUTHORS

Yann Chemin, International Rice Research Institute, The Philippines

Last changed: $Date$