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