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