_RL sumcos
sumcos = 0. _d 0
sumcos = sumcos + cos(yc(i,j,bi,bj)*deg2rad)
if(sumcos.eq.0) sumcos=1
fctilemm = (fctilemm / sumcos)
if(sumcos.eq.0) sumcos=1.0