itr = gencost_itracer(kgen)
& gencost_itracer,
gencost_itracer(k) = 1
& gencost_itracer(k)
itr=gencost_itracer(k)
& gencost_msk_pointer3d, gencost_itracer
integer gencost_itracer(NGENCOST)