& lambdaThetaZonRelax,
_RL lambdaThetaZonRelax
& -lambdaThetaZonRelax*maskC(i,j,kSurface,bi,bj)*
& tauThetaZonRelax, lambdaThetaZonRelax,
lambdaThetaZonRelax = 0.
lambdaThetaZonRelax = 1./tauThetaZonRelax