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