_RL smag2fac, smag4fac
smag2fac=(viscC2smag/pi)**2
smag2fac=0. _d 0
viscAh_DSmg(i,j)=smag2fac*viscAh_DSmg(i,j)
viscAh_ZSmg(i,j)=smag2fac*viscAh_ZSmg(i,j)