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