& mult_shifwflx,
_RL mult_shifwflx
& + mult_shifwflx * objf_shifwflx(bi,bj)
mult_shifwflx = 0. _d 0
& , mult_shifwflx, wshifwflx0, shifwflx_errfile