Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/604/A135 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-07
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/604/A135 Dynamics of ionized and neutral gas in M8 (Damiani+, 2017)
*================================================================================
*The Gaia-ESO Survey: dynamics of ionized and neutral gas in the
*Lagoon nebula (M8).
* Damiani F., Bonito R., Prisinzano L., Zwitter T., Bayo A., Kalari V.,
* Jimenez-Esteban F. M., Costado M. T., Jofre P., Randich S., Flaccomio E.,
* Lanzafame A. C., Lardo C., Morbidelli L., Zaggia S.
* <Astron. Astrophys. 604, A135 (2017)>
* =2017A&A...604A.135D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Best-fit parameters for Halpha, [NII], [SII]
integer*4 nr__
parameter (nr__=1177) ! Number of records
character*129 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*20 Id (nr__) ! Sky-position identifier (1)
real*8 RAdeg (nr__) ! (deg) Right Ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
real*4 RVHalpha (nr__) ! (km/s) ? Radial velocity of H-alpha line
real*4 sigHalpha (nr__) ! (km/s) ? Gaussian width of H-alpha line
real*8 NormHalpha (nr__) ! (km/s/min) ? Intensity of H-alpha line
* (in ADU*km/s/min)
real*4 RV_NII_6584(nr__) ! (km/s) ? Radial velocity of [NII] 6584
real*4 sig_NII_6584(nr__) ! (km/s) ? Gaussian width of [NII] 6584
real*8 Norm_NII_6584(nr__) ! (km/s/min) ? Intensity of [NII] 6584
* (in ADU*km/s/min)
real*4 RV_SII_6717(nr__) ! (km/s) ? Radial velocity of [SII] 6717
real*4 sig_SII_6717(nr__) ! (km/s) ? Gaussian width of [SII] 6717
real*8 Norm_SII_6717(nr__) ! (km/s/min) ? Intensity of [SII] 6717
* (in ADU*km/s/min)
real*4 RV_SII_6731(nr__) ! (km/s) ? Radial velocity of [SII] 6731
real*4 sig_SII_6731(nr__) ! (km/s) ? Gaussian width of [SII] 6731
real*8 Norm_SII_6731(nr__) ! (km/s/min) ? Intensity of [SII] 6731
* (in ADU*km/s/min)
*Note (1): Identifier from J2000 coordinates of sky/star position,
* SKY_HHMMSSss+DDMMSSs or HHMMSSss+DDMMSSs.
* Id's starting with 'SKY' are measurements from pure-sky fibers,
* all others are from faint-star fibers.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Multi-Gaussian best fits to NaI lines
integer*4 nr__1
parameter (nr__1=53) ! Number of records
character*118 ar__1 ! Full-size record
character*6 StarId (nr__1) ! Identifier from UVES setup/star Id (2)
character*16 Name (nr__1) ! Star identifier from J2000 coordinates
* (HHMMSSss+DDMMSSs)
real*4 RVComp1 (nr__1) ! (km/s) Radial velocity of component #1
real*4 sigComp1 (nr__1) ! (km/s) Gaussian width of component #1
real*8 NormComp1 (nr__1) ! (km/s) Intensity of component #1
real*4 RVComp2 (nr__1) ! (km/s) ? Radial velocity of component #2
real*4 sigComp2 (nr__1) ! (km/s) ? Gaussian width of component #2
real*8 NormComp2 (nr__1) ! (km/s) ? Intensity of component #2
real*4 RVComp3 (nr__1) ! (km/s) ? Radial velocity of component #3
real*4 sigComp3 (nr__1) ! (km/s) ? Gaussian width of component #3
real*8 NormComp3 (nr__1) ! (km/s) ? Intensity of component #3
real*4 RVComp4 (nr__1) ! (km/s) ? Radial velocity of component #4
real*4 sigComp4 (nr__1) ! (km/s) ? Gaussian width of component #4
real*4 NormComp4 (nr__1) ! (km/s) ? Intensity of component #4
real*4 RVComp5 (nr__1) ! (km/s) ? Radial velocity of component #5
real*4 sigComp5 (nr__1) ! (km/s) ? Gaussian width of component #5
real*4 NormComp5 (nr__1) ! (km/s) ? Intensity of component #5
*Note (2): Form of identifier is UUU-NN. UVES setup UUU is either 520 or 580nm;
* star number NN is sequential within the same UVES setup.
C=============================================================================
C Loading file 'table2.dat' ! Best-fit parameters for Halpha, [NII], [SII]
C Format for file interpretation
1 format(
+ A20,1X,F8.4,1X,F9.5,1X,F6.2,1X,F5.2,1X,F9.2,1X,F6.2,1X,F5.2,
+ 1X,F9.2,1X,F6.2,1X,F5.2,1X,F8.2,1X,F6.2,1X,F5.2,1X,F8.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,1177
read(1,'(A129)')ar__
read(ar__,1)
+ Id(i__),RAdeg(i__),DEdeg(i__),RVHalpha(i__),sigHalpha(i__),
+ NormHalpha(i__),RV_NII_6584(i__),sig_NII_6584(i__),
+ Norm_NII_6584(i__),RV_SII_6717(i__),sig_SII_6717(i__),
+ Norm_SII_6717(i__),RV_SII_6731(i__),sig_SII_6731(i__),
+ Norm_SII_6731(i__)
if(ar__(41:46) .EQ. '') RVHalpha(i__) = rNULL__
if(ar__(48:52) .EQ. '') sigHalpha(i__) = rNULL__
if(ar__(54:62) .EQ. '') NormHalpha(i__) = rNULL__
if(ar__(64:69) .EQ. '') RV_NII_6584(i__) = rNULL__
if(ar__(71:75) .EQ. '') sig_NII_6584(i__) = rNULL__
if(ar__(77:85) .EQ. '') Norm_NII_6584(i__) = rNULL__
if(ar__(87:92) .EQ. '') RV_SII_6717(i__) = rNULL__
if(ar__(94:98) .EQ. '') sig_SII_6717(i__) = rNULL__
if(ar__(100:107) .EQ. '') Norm_SII_6717(i__) = rNULL__
if(ar__(109:114) .EQ. '') RV_SII_6731(i__) = rNULL__
if(ar__(116:120) .EQ. '') sig_SII_6731(i__) = rNULL__
if(ar__(122:129) .EQ. '') Norm_SII_6731(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Id(i__),RAdeg(i__),DEdeg(i__),RVHalpha(i__),sigHalpha(i__),
+ NormHalpha(i__),RV_NII_6584(i__),sig_NII_6584(i__),
+ Norm_NII_6584(i__),RV_SII_6717(i__),sig_SII_6717(i__),
+ Norm_SII_6717(i__),RV_SII_6731(i__),sig_SII_6731(i__),
+ Norm_SII_6731(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Multi-Gaussian best fits to NaI lines
C Format for file interpretation
2 format(
+ A6,1X,A16,1X,F5.2,1X,F4.2,1X,F7.2,1X,F6.2,1X,F5.2,1X,F7.2,1X,
+ F6.2,1X,F4.2,1X,F7.2,1X,F6.2,1X,F4.2,1X,F5.2,1X,F6.2,1X,F4.2,
+ 1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,53
read(1,'(A118)')ar__1
read(ar__1,2)
+ StarId(i__),Name(i__),RVComp1(i__),sigComp1(i__),
+ NormComp1(i__),RVComp2(i__),sigComp2(i__),NormComp2(i__),
+ RVComp3(i__),sigComp3(i__),NormComp3(i__),RVComp4(i__),
+ sigComp4(i__),NormComp4(i__),RVComp5(i__),sigComp5(i__),
+ NormComp5(i__)
if(ar__1(44:49) .EQ. '') RVComp2(i__) = rNULL__
if(ar__1(51:55) .EQ. '') sigComp2(i__) = rNULL__
if(ar__1(57:63) .EQ. '') NormComp2(i__) = rNULL__
if(ar__1(65:70) .EQ. '') RVComp3(i__) = rNULL__
if(ar__1(72:75) .EQ. '') sigComp3(i__) = rNULL__
if(ar__1(77:83) .EQ. '') NormComp3(i__) = rNULL__
if(ar__1(85:90) .EQ. '') RVComp4(i__) = rNULL__
if(ar__1(92:95) .EQ. '') sigComp4(i__) = rNULL__
if(ar__1(97:101) .EQ. '') NormComp4(i__) = rNULL__
if(ar__1(103:108) .EQ. '') RVComp5(i__) = rNULL__
if(ar__1(110:113) .EQ. '') sigComp5(i__) = rNULL__
if(ar__1(115:118) .EQ. '') NormComp5(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ StarId(i__),Name(i__),RVComp1(i__),sigComp1(i__),
+ NormComp1(i__),RVComp2(i__),sigComp2(i__),NormComp2(i__),
+ RVComp3(i__),sigComp3(i__),NormComp3(i__),RVComp4(i__),
+ sigComp4(i__),NormComp4(i__),RVComp5(i__),sigComp5(i__),
+ NormComp5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end