theEnergy = 0.
theEnergy = theEnergy + theEng(n)
CALL MON_OUT_RL(mon_var,theEnergy, mon_sufx(0), myThid)
_RL theEng(nLatBnd), theEnergy
theEnergy = -rhos*Lfresh*theMeanG
theEnergy = theEnergy -rhoi*(theMean1+theMean2)*theVolG/2
CALL MON_OUT_RL(mon_var, theEnergy, mon_sufx(0), myThid)
_RL theMean1, theMean2, theEnergy