Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/625/A19 into FORTRAN code for loading all data files into arrays.
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-Oct-11
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/625/A19 33 GHz continuum map of M51 (Querejeta+, 2019)
*================================================================================
*Dense gas is not enough: environmental variations in the star formation
*efficiency of dense molecular gas at 100pc scales in M51.
* Querejeta M., Schinnerer E., Schruba A., Murphy E., Meidt S., Usero A.,
* Leroy A. K., Pety J., Bigiel F., Chevance M., Faesi C. M., Gallagher M.,
* Garcia-Burillo S., Glover S. C. O., Hygate A. P. S., Jimenez-Donaire M. J.,
* Kruijssen J. M. D., Momjian E., Rosolowsky E., Utomo D.
* <Astron. Astrophys. 625, A19 (2019)>
* =2019A&A...625A..19Q (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Measurements for apertures on HCN peaks
integer*4 nr__
parameter (nr__=27) ! Number of records
character*122 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*8 Name (nr__) ! Region in M51 (north, ring, south)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
integer*4 SigStar (nr__) ! (Msun/pc2) Stellar mass surface density
integer*4 e_SigStar (nr__) ! (Msun/pc2) Uncertainty stellar mass surface density
real*4 TIR (nr__) ! (10+6Lsun) Total Infrared Luminosity (TIR)
real*4 e_TIR (nr__) ! (10+6Lsun) Uncertainty TIR
real*4 F33GHz (nr__) ! (uJy) 33GHz Spectral Flux Density
real*4 e_F33GHz (nr__) ! (uJy) Uncertainty 33GHz Flux Density
real*4 tfrac (nr__) ! (%) Thermal fraction at 33GHz
real*4 e_tfrac (nr__) ! (%) Uncertainty thermal fraction 33GHz
real*4 SigSFR (nr__) ! (Msun/Myr/pc2) Star formation rate surface density
real*4 e_SigSFR (nr__) ! (Msun/Myr/pc2) Uncertainty SFR surface density
real*4 IHCN (nr__) ! (K.km/s) HCN(1-0) intensity
real*4 e_IHCN (nr__) ! (K.km/s) Uncertainty HCN(1-0) intensity
real*4 ICO (nr__) ! (K.km/s) ?=- CO(1-0) intensity
real*4 e_ICO (nr__) ! (K.km/s) ?=- Uncertainty CO(1-0) intensity
real*4 sigHCN (nr__) ! (km/s) ?=- HCN(1-0) velocity dispersion
real*4 e_sigHCN (nr__) ! (km/s) ?=- Uncertainty HCN(1-0) velocity
* dispersion
real*4 sigCO (nr__) ! (km/s) ?=- CO(1-0) velocity dispersion
real*4 e_sigCO (nr__) ! (km/s) ?=- Uncertainty CO(1-0) velocity
* dispersion
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Measurements for Nyquist-sampled apertures
integer*4 nr__1
parameter (nr__1=4280) ! Number of records
character*122 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*8 Name_1 (nr__1) ! Region in M51 (north, ring, south)
real*8 RAdeg_1 (nr__1) ! (deg) Right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination (J2000)
integer*4 SigStar_1 (nr__1) ! (Msun/pc2) Stellar mass surface density
integer*4 e_SigStar_1(nr__1) ! (Msun/pc2) Uncertainty stellar mass surface density
real*4 TIR_1 (nr__1) ! (10+6Lsun) Total Infrared Luminosity (TIR)
real*4 e_TIR_1 (nr__1) ! (10+6Lsun) Uncertainty TIR
real*4 F33GHz_1 (nr__1) ! (uJy) 33GHz Spectral Flux Density
real*4 e_F33GHz_1 (nr__1) ! (uJy) Uncertainty 33GHz Flux Density
real*4 tfrac_1 (nr__1) ! (%) Thermal fraction at 33GHz
real*4 e_tfrac_1 (nr__1) ! (%) Uncertainty thermal fraction 33GHz
real*4 SigSFR_1 (nr__1) ! (Msun/Myr/pc2) Star formation rate surface density
real*4 e_SigSFR_1 (nr__1) ! (Msun/Myr/pc2) Uncertainty SFR surface density
real*4 IHCN_1 (nr__1) ! (K.km/s) HCN(1-0) intensity
real*4 e_IHCN_1 (nr__1) ! (K.km/s) Uncertainty HCN(1-0) intensity
real*4 ICO_1 (nr__1) ! (K.km/s) ?=- CO(1-0) intensity
real*4 e_ICO_1 (nr__1) ! (K.km/s) ?=- Uncertainty CO(1-0) intensity
real*4 sigHCN_1 (nr__1) ! (km/s) ?=- HCN(1-0) velocity dispersion
real*4 e_sigHCN_1 (nr__1) ! (km/s) ?=- Uncertainty HCN(1-0) velocity
* dispersion
real*4 sigCO_1 (nr__1) ! (km/s) ?=- CO(1-0) velocity dispersion
real*4 e_sigCO_1 (nr__1) ! (km/s) ?=- Uncertainty CO(1-0) velocity
* dispersion
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! Information of fits datacube
integer*4 nr__2
parameter (nr__2=1) ! Number of records
character*88 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension of center (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination of center (J2000)
integer*4 Nx (nr__2) ! Number of pixels along X-axis
integer*4 Ny (nr__2) ! Number of pixels along Y-axis
character*19 Obs_date (nr__2) ! ("datime") Observation date
integer*4 size (nr__2) ! (Kibyte) Size of FITS file
character*8 FileName (nr__2) ! Name of FITS file, in subdirectory fits
character*27 Title (nr__2) ! Title of the FITS file
C=============================================================================
C Loading file 'tablea1.dat' ! Measurements for apertures on HCN peaks
C Format for file interpretation
1 format(
+ A8,1X,F8.4,1X,F8.5,1X,I4,1X,I4,1X,F5.1,1X,F4.1,1X,F5.1,1X,
+ F3.1,1X,F6.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,
+ 1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,27
read(1,'(A122)')ar__
read(ar__,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),SigStar(i__),e_SigStar(i__),
+ TIR(i__),e_TIR(i__),F33GHz(i__),e_F33GHz(i__),tfrac(i__),
+ e_tfrac(i__),SigSFR(i__),e_SigSFR(i__),IHCN(i__),e_IHCN(i__),
+ ICO(i__),e_ICO(i__),sigHCN(i__),e_sigHCN(i__),sigCO(i__),
+ e_sigCO(i__)
if (idig(ar__(93:97)).EQ.0) ICO(i__) = rNULL__
if (idig(ar__(99:102)).EQ.0) e_ICO(i__) = rNULL__
if (idig(ar__(104:107)).EQ.0) sigHCN(i__) = rNULL__
if (idig(ar__(109:111)).EQ.0) e_sigHCN(i__) = rNULL__
if (idig(ar__(113:117)).EQ.0) sigCO(i__) = rNULL__
if (idig(ar__(119:122)).EQ.0) e_sigCO(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),SigStar(i__),e_SigStar(i__),
+ TIR(i__),e_TIR(i__),F33GHz(i__),e_F33GHz(i__),tfrac(i__),
+ e_tfrac(i__),SigSFR(i__),e_SigSFR(i__),IHCN(i__),e_IHCN(i__),
+ ICO(i__),e_ICO(i__),sigHCN(i__),e_sigHCN(i__),sigCO(i__),
+ e_sigCO(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Measurements for Nyquist-sampled apertures
C Format for file interpretation
2 format(
+ A8,1X,F8.4,1X,F8.5,1X,I4,1X,I4,1X,F5.1,1X,F4.1,1X,F5.1,1X,
+ F3.1,1X,F6.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,
+ 1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,4280
read(1,'(A122)')ar__1
read(ar__1,2)
+ Name_1(i__),RAdeg_1(i__),DEdeg_1(i__),SigStar_1(i__),
+ e_SigStar_1(i__),TIR_1(i__),e_TIR_1(i__),F33GHz_1(i__),
+ e_F33GHz_1(i__),tfrac_1(i__),e_tfrac_1(i__),SigSFR_1(i__),
+ e_SigSFR_1(i__),IHCN_1(i__),e_IHCN_1(i__),ICO_1(i__),
+ e_ICO_1(i__),sigHCN_1(i__),e_sigHCN_1(i__),sigCO_1(i__),
+ e_sigCO_1(i__)
if (idig(ar__1(93:97)).EQ.0) ICO_1(i__) = rNULL__
if (idig(ar__1(99:102)).EQ.0) e_ICO_1(i__) = rNULL__
if (idig(ar__1(104:107)).EQ.0) sigHCN_1(i__) = rNULL__
if (idig(ar__1(109:111)).EQ.0) e_sigHCN_1(i__) = rNULL__
if (idig(ar__1(113:117)).EQ.0) sigCO_1(i__) = rNULL__
if (idig(ar__1(119:122)).EQ.0) e_sigCO_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),RAdeg_1(i__),DEdeg_1(i__),SigStar_1(i__),
+ e_SigStar_1(i__),TIR_1(i__),e_TIR_1(i__),F33GHz_1(i__),
+ e_F33GHz_1(i__),tfrac_1(i__),e_tfrac_1(i__),SigSFR_1(i__),
+ e_SigSFR_1(i__),IHCN_1(i__),e_IHCN_1(i__),ICO_1(i__),
+ e_ICO_1(i__),sigHCN_1(i__),e_sigHCN_1(i__),sigCO_1(i__),
+ e_sigCO_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'list.dat' ! Information of fits datacube
C Format for file interpretation
3 format(F9.5,F9.5,1X,I3,1X,I3,1X,A19,1X,I4,1X,A8,1X,A27)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,1
read(1,'(A88)')ar__2
read(ar__2,3)
+ RAdeg_2(i__),DEdeg_2(i__),Nx(i__),Ny(i__),Obs_date(i__),
+ size(i__),FileName(i__),Title(i__)
c ..............Just test output...........
write(6,3)
+ RAdeg_2(i__),DEdeg_2(i__),Nx(i__),Ny(i__),Obs_date(i__),
+ size(i__),FileName(i__),Title(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end