fctiletot = 0.
fctiletot = fctiletot
_GLOBAL_SUM_RL( fctiletot , myThid )
& + ((fctiletot/sumtot+1.D0)/wmean_hflux)**2
& fctiletot/sumtot,
_RL fctiletot
& + ( (fctiletot/sumtot)/wmean_sflux )**2