Thread safety fix.
General WARNING about lines like the following:

integer, dimension(:), pointer :: nodes => null()

The initialisation expression '=> null()' has the, usually unintended, consequence of adding the save status to the variable, i.e. it becomes a static variable whose value (in this case association status) is retained between calls. The initialisation expression is in fact only "evaluated" in the first call. Do not use initialisation expressions for local variables, unless this is what you want - in which case it would be nice to explicitly add the save attribute.