To do what Prof. Taylor suggests you would do something like the following, then you will have access to the nodal coordinates. This would be preferred to my other suggestion.
include 'eldata.h'
include 'pointer.h'
include 'sdata.h'
include 'comblk.h'
real (kind=8) :: xl(:,:)
xl(1:ndm,1:nel) => hr(np(44):np(44)+ndm*nel-1)
This will allow you to access xl( , ) in the 'usual way'. Note for this to work your comblk.h needs to be set up with the target attribute (the archive will need to be recompiled):
!-----[--+---------+---------+---------+---------+---------+---------+-]
! Modification log Date(dd/mm/year)
! 1. Increase dimension to 1024 to force loops on 03/10/2011
! long arrays
real (kind=8), target:: hr
integer , target:: mr
common /comblk/ hr(1024),mr(1024)