& objf_atempm, objf_atempsmoo,
& objf_atemp, objf_atempm, objf_atempsmoo,
_RL objf_atempsmoo(nsx,nsy)
& +mult_smooth_bc*objf_atempsmoo(bi,bj) )
f_atempsmoo = f_atempsmoo + objf_atempsmoo(bi,bj)