_RL leithD2fac, leithD4fac
leithD2fac=(viscC2leithD/pi)**6
leithD2fac=(viscC2leithD/pi)**3
leithD2fac=0. _d 0
& SQRT(leith2fac*grdVrt+leithD2fac*grdDiv)*L3
& SQRT(leithD2fac*grdDiv)*L3
viscAh_Dlth(i,j)=(leith2fac*grdVrt+(leithD2fac*grdDiv))*L3
viscAh_DlthD(i,j)=((leithD2fac*grdDiv))*L3
viscAh_ZLth(i,j)=(leith2fac*grdVrt+(leithD2fac*grdDiv))*L3
viscAh_ZLthD(i,j)=(leithD2fac*grdDiv)*L3