& mult_bottomdrag,
_RL mult_bottomdrag
& + mult_bottomdrag * objf_bottomdrag(bi,bj)
& mult_edtau, mult_bottomdrag,
mult_bottomdrag = 0. _d 0