_RL gammaTmoleT, gammaTmoleS, gammaTurb, gammaTurbConst
gammaTurb = LOG( ustarSq * shiZetaN * etastar**2
gammaTurb = gammaTurbConst
& ustar/(gammaTurb + gammaTmoleT) )
& ustar/(gammaTurb + gammaTmoleS) )