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