eta_qrN = 0. _d 0
eta_qrN = eta_qrN + eta_qrNtile
_GLOBAL_SUM_R8(eta_qrN, myThid)
cgBeta = eta_qrN/eta_qrNM1
eta_qrNM1 = eta_qrN
alpha = eta_qrN/alpha
_RL eta_qrN,eta_qrNtile
_RL eta_qrN, eta_qrNtile