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