_RL recip_rhoConstFresh
recip_rhoConstFresh = 1. _d 0/rhoConstFresh
evap(i,j,bi,bj) = -recip_rhoConstFresh*tau(i,j)*qstar(i,j)
& *recip_lhEvap*recip_rhoConstFresh
_RL recip_rhoConstFresh, recip_lhEvap
recip_rhoConstFresh = 1./rhoConstFresh