_RL L2, L3, L5, L2rdt, L4rdt, recip_dt
L5 = (L2*L3)
& SQRT(leith4fac*grdVrt+leithD4fac*grdDiv)*L5
& SQRT(leithD4fac*grdDiv)*L5
viscA4_Dlth(i,j)=(leith4fac*grdVrt+(leithD4fac*grdDiv))*L5
viscA4_DlthD(i,j)=((leithD4fac*grdDiv))*L5
viscA4_ZLth(i,j)=(leith4fac*grdVrt+(leithD4fac*grdDiv))*L5
viscA4_ZLthD(i,j)=(leithD4fac*grdDiv)*L5