Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/525/A90 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 2024-Apr-20
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/525/A90 Thick disc vertical properties (Katz+, 2011)
*================================================================================
*Probing the Galactic thick disc vertical properties and interfaces.
* Katz D., Soubiran C., Cayrel R., Barbuy B., Friel E., Bienayme O.,
* Perrin M.-N.
* <Astron. Astrophys. 525, A90 (2011)>
* =2011A&A...525A..90K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Program stars coordinates and magnitudes
integer*4 nr__
parameter (nr__=346) ! Number of records
character*101 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*10 Name (nr__) ! Designation (FMN-NNNNNN;
* <[KSC2011] FMN-NNNNNN> in Simbad)
character*7 Spectro (nr__) ! Spectrograph: CARELEC, CFHT or BOTH
character*7 ExpID (nr__) ! Id of the exposure (could be several
* exposures)
integer*4 RAh (nr__) ! (h) Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Right ascension (J2000)
real*4 RAs (nr__) ! (s) Right ascension (J2000)
character*1 DE_ (nr__) ! Declination sign (J2000)
integer*4 DEd (nr__) ! (deg) Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Declination (J2000)
real*4 Vmag (nr__) ! (mag) Photographic V apparent magnitude
real*4 B_V (nr__) ! (mag) ? B-V photographic colour index
real*4 V_I (nr__) ! (mag) ? V-I colour index
real*4 Jmag (nr__) ! (mag) 2MASS J apparent magnitude
real*4 Hmag (nr__) ! (mag) 2MASS H apparent magnitude
real*4 Ksmag (nr__) ! (mag) ? 2MASS Ks apparent magnitude
real*4 pmRA (nr__) ! (mas/yr) Proper motion along right ascension
real*4 pmDE (nr__) ! (mas/yr) Proper motion along declination
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Program stars stellar parameters
integer*4 nr__1
parameter (nr__1=346) ! Number of records
character*45 ar__1 ! Full-size record
character*10 Name_1 (nr__1) ! Designation
integer*4 Teff (nr__1) ! (K) Effective temperature
real*4 logg (nr__1) ! ([cm/s2]) Logarithm of the surface gravity
real*4 v_Fe_H_ (nr__1) ! ([Sun]) Metallicity
real*4 VMAG_1 (nr__1) ! (mag) Absolute visual magnitude
real*8 RV (nr__1) ! (km/s) ? Radial velocity
integer*4 Teff1 (nr__1) ! (K) ? Effective temperature from SDSS photometry
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! Program stars metallicities, positions, velocities
integer*4 nr__2
parameter (nr__2=346) ! Number of records
character*80 ar__2 ! Full-size record
character*10 Name_2 (nr__2) ! Designation
real*4 v_Fe_H__1 (nr__2) ! ([Sun]) Metallicity
real*8 d (nr__2) ! (pc) Distance (from the Sun)
real*8 X (nr__2) ! (pc) Galactic position X
real*8 Y (nr__2) ! (pc) Galactic position Y
real*8 Z (nr__2) ! (pc) Galactic position Z
real*8 Uvel (nr__2) ! (km/s) ? Galactic velocity U (1)
real*8 Vvel (nr__2) ! (km/s) ? Galactic velocity V
real*8 Wvel (nr__2) ! (km/s) ? Galactic velocity W
real*8 U_W (nr__2) ! (km/s) ? U(cos47) - W(sin47)
*Note (1): U positive toward the Galactic center
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Elodie library reference stars
integer*4 nr__3
parameter (nr__3=1398) ! Number of records
character*35 ar__3 ! Full-size record
character*9 Name_3 (nr__3) ! Designation
integer*4 Teff_1 (nr__3) ! (K) Effective temperature (1)
real*4 logg_1 (nr__3) ! ([cm/s2]) Logarithm of the surface gravity (1)
real*4 v_Fe_H__2 (nr__3) ! ([Sun]) ? Metallicity (1)
real*4 VMAG_2 (nr__3) ! (mag) ? Absolute visual magnitude (1)
*Note (1): The stellar parameters have been derived from the PASTEL
* catalogue: Soubiran et al., 2010, Cat. B/pastel .
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb2.dat' ! CARELEC validation stars
integer*4 nr__4
parameter (nr__4=45) ! Number of records
character*35 ar__4 ! Full-size record
character*9 Name_4 (nr__4) ! Designation
integer*4 Teff_2 (nr__4) ! (K) Effective temperature (1)
real*4 logg_2 (nr__4) ! ([cm/s2]) Logarithm of the surface gravity (1)
real*4 v_Fe_H__3 (nr__4) ! ([Sun]) ? Metallicity (1)
real*4 VMAG_3 (nr__4) ! (mag) ? Absolute visual magnitude (1)
*Note (1): The stellar parameters have been derived from the PASTEL
* catalogue: Soubiran et al., 2010, Cat. B/pastel .
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb3.dat' ! CFHT validation stars
integer*4 nr__5
parameter (nr__5=46) ! Number of records
character*39 ar__5 ! Full-size record
character*9 Name_5 (nr__5) ! Designation
integer*4 Nspe (nr__5) ! Number of spectra recorded
integer*4 Teff_3 (nr__5) ! (K) Effective temperature (1)
real*4 logg_3 (nr__5) ! ([cm/s2]) Logarithm of the surface gravity (1)
real*4 v_Fe_H__4 (nr__5) ! ([Sun]) Metallicity (1)
real*4 Vmag_4 (nr__5) ! (mag) Apparent visual magnitude (1)
real*4 B_V_1 (nr__5) ! (mag) ? B-V colour index (1)
*Note (1): Stellar parameters for CFHT observations of validation stars
* of the globular clusters:
* * M5: identification names, V and B-V are from Buonanno et al.
* (1981MNRAS.196..435B); stellar parameters Teff, logg, are interpolated
* within values determined by Ramirez & Cohen (2003AJ....125..224R) for
* stars from the tip of the RGB down to the TO.
* * M71: identification names, V and B-V are from Arp & Hartwick
* (1971ApJ...167..499A). Stellar parameters are interpolated within
* values determined by Cohen et al. (2001AJ....122.1420C) for stars from
* the tip of the RGB down to the TO, for which a metallicity [Fe/H]=-0.7
* was found. Stars indicated by KC, S are on the numbering system of
* Cudworth (1985AJ.....90...65C).
* * M67: identifications are from Johnson & Sandage (1955ApJ...121..616J)
* and V and B-V are from Montgomery et al. (1993, Cat. J/AJ/106/181);
* stellar parameters are interpolated in calculations by VandenBerg &
* Clem (2003AJ....126..778V), where Teff values are from calibrations
* by Alonso et al. (1999A&AS..140..261A) and Sekiguchi & Fukugita
* (2000AJ....120.1072S) for a metallicity [Fe/H]=0.0 (An et al.,
* 2009ApJ...700..523A).
* * M92: Identification names, V and B-V are from Sandage & Walker
* (1966ApJ...143..313S) and Buonanno et al. (1983A&AS...53....1B), and
* F means Fagerholm. Stellar parameters are interpolated within values
* by Sneden et al. (2000AJ....120.1351S).
C=============================================================================
C Loading file 'tablea1.dat' ! Program stars coordinates and magnitudes
C Format for file interpretation
1 format(
+ A10,1X,A7,1X,A7,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F5.2,1X,F4.2,1X,F4.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.1,1X,F6.1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,346
read(1,'(A101)')ar__
read(ar__,1)
+ Name(i__),Spectro(i__),ExpID(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Vmag(i__),B_V(i__),
+ V_I(i__),Jmag(i__),Hmag(i__),Ksmag(i__),pmRA(i__),pmDE(i__)
if(ar__(58:61) .EQ. '') B_V(i__) = rNULL__
if(ar__(63:66) .EQ. '') V_I(i__) = rNULL__
if(ar__(82:87) .EQ. '') Ksmag(i__) = rNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),Spectro(i__),ExpID(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Vmag(i__),B_V(i__),
+ V_I(i__),Jmag(i__),Hmag(i__),Ksmag(i__),pmRA(i__),pmDE(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Program stars stellar parameters
C Format for file interpretation
2 format(A10,1X,I4,1X,F4.2,1X,F5.2,1X,F5.2,1X,F7.2,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,346
read(1,'(A45)')ar__1
read(ar__1,2)
+ Name_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),VMAG_1(i__),
+ RV(i__),Teff1(i__)
if(ar__1(34:40) .EQ. '') RV(i__) = rNULL__
if(ar__1(42:45) .EQ. '') Teff1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),VMAG_1(i__),
+ RV(i__),Teff1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea3.dat' ! Program stars metallicities, positions, velocities
C Format for file interpretation
3 format(
+ A10,1X,F5.2,1X,F7.1,1X,F7.1,1X,F7.1,1X,F7.1,1X,F7.1,1X,F7.1,
+ 1X,F7.1,1X,F7.1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,346
read(1,'(A80)')ar__2
read(ar__2,3)
+ Name_2(i__),v_Fe_H__1(i__),d(i__),X(i__),Y(i__),Z(i__),
+ Uvel(i__),Vvel(i__),Wvel(i__),U_W(i__)
if(ar__2(50:56) .EQ. '') Uvel(i__) = rNULL__
if(ar__2(58:64) .EQ. '') Vvel(i__) = rNULL__
if(ar__2(66:72) .EQ. '') Wvel(i__) = rNULL__
if(ar__2(74:80) .EQ. '') U_W(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),v_Fe_H__1(i__),d(i__),X(i__),Y(i__),Z(i__),
+ Uvel(i__),Vvel(i__),Wvel(i__),U_W(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Elodie library reference stars
C Format for file interpretation
4 format(A9,2X,I5,1X,F5.2,1X,F5.2,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,1398
read(1,'(A35)')ar__3
read(ar__3,4)
+ Name_3(i__),Teff_1(i__),logg_1(i__),v_Fe_H__2(i__),VMAG_2(i__)
if(ar__3(24:28) .EQ. '') v_Fe_H__2(i__) = rNULL__
if(ar__3(30:35) .EQ. '') VMAG_2(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),Teff_1(i__),logg_1(i__),v_Fe_H__2(i__),VMAG_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb2.dat' ! CARELEC validation stars
C Format for file interpretation
5 format(A9,2X,I5,1X,F5.2,1X,F5.2,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'tableb2.dat')
write(6,*) '....Loading file: tableb2.dat'
do i__=1,45
read(1,'(A35)')ar__4
read(ar__4,5)
+ Name_4(i__),Teff_2(i__),logg_2(i__),v_Fe_H__3(i__),VMAG_3(i__)
if(ar__4(24:28) .EQ. '') v_Fe_H__3(i__) = rNULL__
if(ar__4(30:35) .EQ. '') VMAG_3(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4(i__),Teff_2(i__),logg_2(i__),v_Fe_H__3(i__),VMAG_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb3.dat' ! CFHT validation stars
C Format for file interpretation
6 format(A9,2X,I1,1X,I4,1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb3.dat')
write(6,*) '....Loading file: tableb3.dat'
do i__=1,46
read(1,'(A39)')ar__5
read(ar__5,6)
+ Name_5(i__),Nspe(i__),Teff_3(i__),logg_3(i__),v_Fe_H__4(i__),
+ Vmag_4(i__),B_V_1(i__)
if(ar__5(36:39) .EQ. '') B_V_1(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Name_5(i__),Nspe(i__),Teff_3(i__),logg_3(i__),v_Fe_H__4(i__),
+ Vmag_4(i__),B_V_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end