& objf_thsice
_RL objf_thsice (nsx,nsy)
& + mult_thsice * objf_thsice(bi,bj)
f_thsice = f_thsice + objf_thsice(bi,bj)
objf_thsice(bi,bj) = 0. _d 0
objf_thsice(bi,bj) = objf_thsice(bi,bj) +