Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/701/A111 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Aug-16
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/A+A/701/A111 PDRs4All XVI. Orion Bar template spectra (Van De Putte+, 2025)
*================================================================================
*PDRs4All. XVI. Tracing aromatic infrared band characteristics in
*photodissociation region spectra with PAHFIT in the JWST era.
* Van De Putte D., Peeters E., Gordon K.D., Smith J.D.T., Lai T.S.-Y.,
* Maragkoudakis A., Schefter B., Sidhu A., Doshi D., Berne O., Cami J.,
* Boersma C., Dartois E., Habart E., Onaka T., Tielens A.G.G.M.
* <Astron. Astrophys. 701, A111 (2025)>
* =2025A&A...701A.111V (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! PAHFIT dust feature tuning results for
Atomic PDR DF3, and their mean
integer*4 nr__
parameter (nr__=63) ! Number of records
character*48 ar__ ! Full-size record
real*4 lambdaAtomic ! (um) ? Central wavelength fit in Atomic PDR
real*4 FWHMAtomic ! (um) ? FWHM fit in Atomic PDR
real*4 lambdaDF3 ! (um) ? Central wavelength fit in DF3
real*4 FWHMDF3 ! (um) ? FWHM fit in DF3
real*4 lambdamean ! (um) Mean central wavelength
real*4 FWHMmean ! (um) Mean FWHM
character*9 Selection ! Selected central wavelength and FWHM
* for use in PDR pack
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ngc7469.dat' ! NIRSpec + MRS aperture extractions of NGC7469
star forming ring
integer*4 nr__1
parameter (nr__1=13980) ! Number of records
character*263 ar__1 ! Full-size record
real*8 lambda ! (um) [0.97/28.7] Wavelength
real*8 Inu_SF1 ! (MJy/sr) SF1 aperture flux
real*8 e_Inu_SF1 ! (MJy/sr) SF1 aperture uncertainty
real*8 Inu_SF2 ! (MJy/sr) ? SF2 aperture flux
real*8 e_Inu_SF2 ! (MJy/sr) ? SF2 aperture uncertainty
real*8 Inu_SF3 ! (MJy/sr) ? SF3 aperture flux
real*8 e_Inu_SF3 ! (MJy/sr) ? SF3 aperture uncertainty
real*8 Inu_SF4 ! (MJy/sr) ? SF4 aperture flux
real*8 e_Inu_SF4 ! (MJy/sr) ? SF4 aperture uncertainty
real*8 Inu_SF5 ! (MJy/sr) ? SF5 aperture flux
real*8 e_Inu_SF5 ! (MJy/sr) ? SF5 aperture uncertainty
real*8 Inu_SF6 ! (MJy/sr) ? SF6 aperture flux
real*8 e_Inu_SF6 ! (MJy/sr) ? SF6 aperture uncertainty
real*8 Inu_SFR ! (MJy/sr) ? SFRing aperture flux
real*8 e_Inu_SFR ! (MJy/sr) ? SFRing aperture uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'orionbar.dat' ! NIRSpec + MRS template spectra of Orion Bar
(corrected version, 05-Nov-2025)
integer*4 nr__2
parameter (nr__2=20675) ! Number of records
character*191 ar__2 ! Full-size record
real*8 lambda_1 ! (um) [0.97/28.7] Wavelength
real*8 Inu_HII ! (MJy/sr) ? HII template spectrum flux
real*8 e_Inu_HII ! (MJy/sr) ? HII template spectrum uncertainty
real*8 Inu_Atm ! (MJy/sr) Atomic PDR template spectrum flux
real*8 e_Inu_Atm ! (MJy/sr) ? Atomic PDR template spectrum uncertainty
real*8 Inu_DF3 ! (MJy/sr) DF3 template spectrum flux
real*8 e_Inu_DF3 ! (MJy/sr) ? F3 template spectrum uncertainty
real*8 Inu_DF2 ! (MJy/sr) DF2 template spectrum flux
real*8 e_Inu_DF2 ! (MJy/sr) ? F2 template spectrum uncertainty
real*8 Inu_DF1 ! (MJy/sr) DF1 template spectrum flux
real*8 e_Inu_DF1 ! (MJy/sr) DF1 template spectrum uncertainty
C=============================================================================
C Loading file 'tableb1.dat' ! PAHFIT dust feature tuning results for
* Atomic PDR DF3, and their mean
C Format for file interpretation
1 format(F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,63
read(1,'(A48)')ar__
read(ar__,1)
+ lambdaAtomic,FWHMAtomic,lambdaDF3,FWHMDF3,lambdamean,FWHMmean,
+ Selection
if(ar__(1:6) .EQ. '') lambdaAtomic = rNULL__
if(ar__(8:12) .EQ. '') FWHMAtomic = rNULL__
if(ar__(14:19) .EQ. '') lambdaDF3 = rNULL__
if(ar__(21:25) .EQ. '') FWHMDF3 = rNULL__
c ..............Just test output...........
write(6,1)
+ lambdaAtomic,FWHMAtomic,lambdaDF3,FWHMDF3,lambdamean,FWHMmean,
+ Selection
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ngc7469.dat' ! NIRSpec + MRS aperture extractions of NGC7469
* star forming ring
C Format for file interpretation
2 format(
+ F11.8,4X,F14.8,6X,F12.8,4X,F14.8,6X,F12.8,4X,F14.8,6X,F12.8,
+ 4X,F14.8,6X,F12.8,4X,F14.8,6X,F12.8,4X,F14.8,6X,F12.8,4X,
+ F14.8,5X,F13.8)
C Effective file loading
open(unit=1,status='old',file=
+'ngc7469.dat')
write(6,*) '....Loading file: ngc7469.dat'
do i__=1,13980
read(1,'(A263)')ar__1
read(ar__1,2)
+ lambda,Inu_SF1,e_Inu_SF1,Inu_SF2,e_Inu_SF2,Inu_SF3,e_Inu_SF3,
+ Inu_SF4,e_Inu_SF4,Inu_SF5,e_Inu_SF5,Inu_SF6,e_Inu_SF6,Inu_SFR,
+ e_Inu_SFR
if(ar__1(52:65) .EQ. '') Inu_SF2 = rNULL__
if(ar__1(72:83) .EQ. '') e_Inu_SF2 = rNULL__
if(ar__1(88:101) .EQ. '') Inu_SF3 = rNULL__
if(ar__1(108:119) .EQ. '') e_Inu_SF3 = rNULL__
if(ar__1(124:137) .EQ. '') Inu_SF4 = rNULL__
if(ar__1(144:155) .EQ. '') e_Inu_SF4 = rNULL__
if(ar__1(160:173) .EQ. '') Inu_SF5 = rNULL__
if(ar__1(180:191) .EQ. '') e_Inu_SF5 = rNULL__
if(ar__1(196:209) .EQ. '') Inu_SF6 = rNULL__
if(ar__1(216:227) .EQ. '') e_Inu_SF6 = rNULL__
if(ar__1(232:245) .EQ. '') Inu_SFR = rNULL__
if(ar__1(251:263) .EQ. '') e_Inu_SFR = rNULL__
c ..............Just test output...........
write(6,2)
+ lambda,Inu_SF1,e_Inu_SF1,Inu_SF2,e_Inu_SF2,Inu_SF3,e_Inu_SF3,
+ Inu_SF4,e_Inu_SF4,Inu_SF5,e_Inu_SF5,Inu_SF6,e_Inu_SF6,Inu_SFR,
+ e_Inu_SFR
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'orionbar.dat' ! NIRSpec + MRS template spectra of Orion Bar
* (corrected version, 05-Nov-2025)
C Format for file interpretation
3 format(
+ F11.8,3X,F15.8,5X,F13.8,3X,F15.8,5X,F13.8,3X,F15.8,6X,F12.8,
+ 3X,F15.8,6X,F12.8,3X,F15.8,5X,F13.8)
C Effective file loading
open(unit=1,status='old',file=
+'orionbar.dat')
write(6,*) '....Loading file: orionbar.dat'
do i__=1,20675
read(1,'(A191)')ar__2
read(ar__2,3)
+ lambda_1,Inu_HII,e_Inu_HII,Inu_Atm,e_Inu_Atm,Inu_DF3,
+ e_Inu_DF3,Inu_DF2,e_Inu_DF2,Inu_DF1,e_Inu_DF1
if(ar__2(15:29) .EQ. '') Inu_HII = rNULL__
if(ar__2(35:47) .EQ. '') e_Inu_HII = rNULL__
if(ar__2(71:83) .EQ. '') e_Inu_Atm = rNULL__
if(ar__2(108:119) .EQ. '') e_Inu_DF3 = rNULL__
if(ar__2(144:155) .EQ. '') e_Inu_DF2 = rNULL__
c ..............Just test output...........
write(6,3)
+ lambda_1,Inu_HII,e_Inu_HII,Inu_Atm,e_Inu_Atm,Inu_DF3,
+ e_Inu_DF3,Inu_DF2,e_Inu_DF2,Inu_DF1,e_Inu_DF1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end