O viscArDefault, diffKrTDefault, diffKrSDefault,
IF ( viscAr .EQ. UNSET_RL ) viscAr = viscArDefault
_RL viscArDefault
viscArDefault = 0.D-3