& mult_smooth_ic,
_RL mult_smooth_ic
& mult_smooth_ic, mult_smooth_bc,
mult_smooth_ic = 0. _d 0
& +mult_smooth_ic*objf_temp0smoo(bi,bj) )
& +mult_smooth_ic*objf_salt0smoo(bi,bj) )