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