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