Yes of course, we’ll wait.
I also have a modern Fortran reimplementation of FFTPACK here: hfsolver/fourier.f90 at master · certik/hfsolver · GitHub, but that should go into a separate package. For fortran-lang/fftpack
I would like to preserve the original algorithm, for benchmarking and reference.