I'm using your code as part of a larger code base. It works great for linear elastic experiments where there is only 1 time increment. However, when I do plastic experiments, the computation time to get the node COORD, U, and Element COORD, LE is very long due to the number of elements. What makes it way worse is the number of increments exacerbates this problem.
Is there a way to retrieve only the outputs from the LAST increment? This also will be useful so I can distinguish the output of the last increment from the last increment of the step before.