theta2PE = theta2PE + th2peTile
_GLOBAL_SUM_R8(theta2PE,myThid)
theta2PE = theta2PE / theArea
theta2PE = theta2PE * rhoConst*recip_horiVertRatio
CALL MON_OUT_RL('_T2PE',theta2PE, mon_foot_mean ,myThid)
_RL vT_Mean, vS_Mean, vT_Heat, theta2PE
theta2PE = 0.