_RL zice
zice = 0.0005 _d 0
rdn=xkar/(log(zref/zice))
q100=ssq+qstar*(dlog(100. _d 0/zice)-psx100)
_RL rdn,ren,rhn,zice,zref
zice=.0005 _d 0