Usage
NetRad(lat, Jday, Tx, Tn, albedo = 0.18, forest = 0, slope = 0,
aspect = 0, airtemp = (Tn+Tx)/2, cloudiness = "Estimate",
surfemissivity = 0.97, surftemp = (Tn+Tx)/2, units = "kJm2d",
AEparams=list(vp=NULL, opt="linear"))
Arguments
lat
latitude. Default is radians, but will automatically convert from degrees if value is larger than pi/2 or less than -pi/2
forest
Forest canopy or shade cover (0-1) [-]
aspect
Ground aspect [radians from north]
airtemp
Average air temp [C]
cloudiness
Cloudiness percentage (0-1) [-]
Default value will estimate cloudiness based on daily temperature difference.
surfemissivity
Surface emissivity [-]
surftemp
Surface temperature [C]
units
Units of the result. Defaults to "kJm2d" which returns result in kJ/m2/d, the other option is "Wm2", which will make output in W/m2
AEparams
Atmospheric Emissivity options. Defaults to linear approximation. To use Brutsaert option, include vapor pressure values (kPa)