objf_bottomdrag(bi,bj) = objf_bottomdrag(bi,bj) + fctile
& objf_bottomdrag,
_RL objf_bottomdrag(nsx,nsy)
& + mult_bottomdrag * objf_bottomdrag(bi,bj)
f_bottomdrag = f_bottomdrag + objf_bottomdrag(bi,bj)
objf_bottomdrag(bi,bj) = 0. _d 0