& cosPower, cAdjFreq,
_RL cosPower
& cosPower, viscAstrain, viscAtension,
IF (cosPower.NE.0.) THEN
& **cosPower
cosPower = 0.