Help with ABAQUS SIGINI subroutine

Dear ABAQUS Users,

I am trying to assign the initial stresses for a region in my model using the SIGINI subroutine. I learned some fortran programming skills the last week and I managed to open the file where the stresses that I retrieve from an Abaqus analysis and assign these stresses in Array using fortran commands.

Now I am stuck with the writing of the SIGINI subroutine where I have the corresponding stresses that I store in my Array to the corresponding NOEL and NPT.

I wish to know if I have to include my fortran code that I use to retrieve the stresses and store them in the array. then I include the subroutine commands found in ABAQUS document where I include a Do loop statement to repeat the assignment for each integration point of the elements located in the region where I am interested to assign the stored stresses.

I will appreciate if anyone has similar code written in FORTRAN where I can learn the linking between the FORTRAN code to open and read the data in the external file (i.e. the external file includes the stresses generated as output by an ABAQUS analysis) where later the subroutine SIGINI statements should be added to assign the stored stresses for the second analysis.

Sorry for the long post but I wish to thank in advance everyone who has any light for such situation.

Best Regards

Sherif