Actually, when I added a LC with ID 1 the solver starts, again without going through the subroutine.
When I change -1 to -2 on the boundary condition card and I add a LC with ID 2 the solver still gives the same error "Load Curve ID 1 not found".
So I am confused as to how that works.
Does the magnitude of the negative number mean anything to the solver?
Does the solver always look for LCID 1 when a negative number is applied?