CALL WRITE_0D_R8( viscC2leith, INDEX_NONE,'viscC2leith =',
& viscC2leith, viscC4leith,
& viscC2leith, viscC2leithD,
_RL viscC2leith
& .OR. viscC2leith.NE.0. .OR. viscC4leith.NE.0.
IF ( (viscC2leith**2*grdVrt+viscC2leithD**2*grdDiv)
Alth2=sqrt(viscC2leith**2*grdVrt+viscC2leithD**2*grdDiv)*L3
Alth2=(viscC2leith*grdVrt+(viscC2leithD*grdDiv))*L3
IF ( viscC2leith.NE.0. .OR. viscAhGrid.NE.0.
& .OR.(viscC2leith.NE.0.)
viscC2leith = 0.D0