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