slaoffset_sum = 0. _d 0
slaoffset_sum=slaoffset_sum+RA(i,j,bi,bj)
_GLOBAL_SUM_RL( slaoffset_sum , mythid )
if (slaoffset_sum .eq. 0.0) then
& slaoffset,slaoffset_sum,slaoffset_weight
slaoffset = slaoffset / slaoffset_sum
& ( slaoffset_sum.GT.0.25 _d 0 * globalArea ) ) THEN
_RL slaoffset_sum