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