I think the key challenge in my code is calling these in line functions:
rr(τ) = 1/((1+prstp)^(tstep*(τ-1)))
I am calling it in the following way:
@NLconstraint(dice_model,
UTILITY == sum((C[i]^(1-elasmu)*l[i]/(1-elasmu)) * rr(1.0) for i in 1:100))