O viscArDefault, diffKrTDefault, diffKrSDefault,
diffKrS = diffKrSDefault
_RL diffKrSDefault
diffKrSDefault = 0. _d -3