Stability Functions and Pardiso

Dear Guru:

I have been using PARDISO for a while with a structural analysis program from Harrison 1973 Fortran.  Updated of course.

I have no problems with it in the standard solve form. 

I am now looking at high compression members which require a reformulation of the stiffness matrix at each stage using STABILITY FUNCTIONS. These now work and the problems solve in reasonable time. 

I now started a problem - dnl.inp. which is an analysis of a 1 km arch, although that is not relevant. If I just use self weight of the beams it solves nicely, but I started to add the mass of the trucks on the bridge, in this case using a simple method just to develop the model - the 1.1 in the equation increases the deadload. The actual code is in the ELEMENTS.F90 file in the attached solution folder. 

do i = 1,n
        write(*,104)i,nodeloads(i,2),-(nodemass(i)*gr),(-(nodemass(i)*gr) + nodeloads(i,2))
104     Format(' Node : ',I4,' Applied Load :: ',F15.3, ' Self Weight :: ',F15.3, ' Total Load :: ',f15.3)
        nodeloads(i,2) = ZERO - ((nodemass(i)*gr)*1.1)
    end do

If I set 1.1 to 1.0 it runs nicely - about 1.01 it spits out a PARDISO -2 error code at about the 5th iteration.  I cannot see in the help files how to solve for memory size problem. 

Help Please - to run just type dnl at the input. 

John