@themos, depending on the Fortran processor, this may be a rather minor “editorial” comment or a significant issue with a data race condition that adversely impacts your coarray version: that is to avoid the use of SAVE statements and also DATA with its “implied save” behavior in Fortran parallel programs.
At first glance, it appears you can employ named constant arrays with both your tp table and hx character list. This comment may be of interest with this.