real function evapol(tu,nu,tv,nv,c,rad,x,y)