& , objf_smrsst
_RL objf_smrsst (nsx,nsy)
& objf_smrsst,objf_smrsss
& objf_smrsst, num_smrsst,
& + mult_smrsst * objf_smrsst(bi,bj)
f_smrsst = f_smrsst + objf_smrsst(bi,bj)
objf_smrsst(bi,bj) = 0. _d 0