CDS.headStuff2
FORTRAN Generation
Conversion of standardized ReadMe file for
catalog 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-19 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/MNRAS/389/1150 Spitzer photometry of globulars in 2 galaxies (Spitler+, 2008) *================================================================================ *Extending the baseline: Spitzer mid-infrared photometry of globular cluster *systems in the Centaurus A and Sombrero Galaxies. * Spitler L.R., Forbes D.A., Beasley M.A. * <Mon. Not. R. Astron. Soc., 389, 1150-1162 (2008)> * =2008MNRAS.389.1150S C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'n5128gcs.dat' ! NGC 5128 globular cluster photometry integer*4 nr__ parameter (nr__=415) ! Number of records character*142 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) integer*4 v_WHB2007_ (nr__) ! Woodley et al. (2007, Cat. J/AJ/134/494) * identification name (GCNNNN) 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 Umag (nr__) ! (mag) ?=99.99 U magnitude real*4 e_Umag (nr__) ! (mag) ?=9.99 U magnitude error real*4 Bmag (nr__) ! (mag) ?=99.99 B magnitude real*4 e_Bmag (nr__) ! (mag) ?=9.99 B magnitude error real*4 Vmag (nr__) ! (mag) ?=99.99 V magnitude real*4 e_Vmag (nr__) ! (mag) ?=9.99 V magnitude error real*4 Rmag (nr__) ! (mag) ?=99.99 R magnitude real*4 e_Rmag (nr__) ! (mag) ?=9.99 R magnitude error real*4 Imag (nr__) ! (mag) ?=99.99 I magnitude real*4 e_Imag (nr__) ! (mag) ?=9.99 I magnitude error real*4 Cmag (nr__) ! (mag) ?=99.99 C magnitude (Washington 421+/-32nm) real*4 e_Cmag (nr__) ! (mag) ?=9.99 C magnitude error real*4 Mmag (nr__) ! (mag) ?=99.99 M magnitude (Washington 802+/-26nm) real*4 e_Mmag (nr__) ! (mag) ?=9.99 M magnitude error real*4 T1mag (nr__) ! (mag) ?=99.99 T1 magnitude (Washington 633+/-80nm) real*4 e_T1mag (nr__) ! (mag) ?=9.99 T1 magnitude error real*4 v3_6mag (nr__) ! (mag) ?=99.99 [3.6] magnitude from Spitzer/IRAC real*4 e_3_6mag (nr__) ! (mag) ?=9.999 [3.6] magnitude error real*4 v4_5mag (nr__) ! (mag) ?=99.99 [4.5] magnitude from Spitzer/IRAC real*4 e_4_5mag (nr__) ! (mag) ?=9.999 [4.5] magnitude error c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'n4594gcs.dat' ! NGC 4594 globular cluster photometry integer*4 nr__1 parameter (nr__1=660) ! Number of records character*95 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*3 v_SLS2006_ (nr__1) ! Spitler et al. (2006, Cat. J/AJ/132/1593) * identification number (1) integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000) integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000) real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000) character*1 DE__1 (nr__1) ! Declination (J2000) integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000) integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000) real*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000) real*4 Bmag_1 (nr__1) ! (mag) ?=99.99 B magnitude real*4 e_Bmag_1 (nr__1) ! (mag) ?=9.99 B magnitude error real*4 Vmag_1 (nr__1) ! (mag) ?=99.99 V magnitude real*4 e_Vmag_1 (nr__1) ! (mag) ?=9.99 V magnitude error real*4 Rmag_1 (nr__1) ! (mag) ?=99.99 R magnitude real*4 e_Rmag_1 (nr__1) ! (mag) ?=9.99 R magnitude error real*4 v3_6mag_1 (nr__1) ! (mag) ?=99.99 [3.6] magnitude from Spitzer/IRAC real*4 e_3_6mag_1 (nr__1) ! (mag) ?=9.999 [3.6] magnitude error real*4 v4_5mag_1 (nr__1) ! (mag) ?=99.99 [4.5] magnitude from Spitzer/IRAC real*4 e_4_5mag_1 (nr__1) ! (mag) ?=9.999 [4.5] magnitude error real*4 Rhl (nr__1) ! (pc) ?=9.99 Half-light radius *Note (1): A ultra-compact dwarf galaxy is included in this catalogue with * name "ucd", identified as "NAME UCD 1" in Simbad. C============================================================================= C Loading file 'n5128gcs.dat' ! NGC 5128 globular cluster photometry C Format for file interpretation 1 format( + 3X,I3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X, + F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2, + 1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X, + F5.2,1X,F5.3,1X,F5.2,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'n5128gcs.dat') write(6,*) '....Loading file: n5128gcs.dat' do i__=1,415 read(1,'(A142)')ar__ read(ar__,1) + v_WHB2007_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),Umag(i__),e_Umag(i__),Bmag(i__),e_Bmag(i__), + Vmag(i__),e_Vmag(i__),Rmag(i__),e_Rmag(i__),Imag(i__), + e_Imag(i__),Cmag(i__),e_Cmag(i__),Mmag(i__),e_Mmag(i__), + T1mag(i__),e_T1mag(i__),v3_6mag(i__),e_3_6mag(i__), + v4_5mag(i__),e_4_5mag(i__) 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) + v_WHB2007_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),Umag(i__),e_Umag(i__),Bmag(i__),e_Bmag(i__), + Vmag(i__),e_Vmag(i__),Rmag(i__),e_Rmag(i__),Imag(i__), + e_Imag(i__),Cmag(i__),e_Cmag(i__),Mmag(i__),e_Mmag(i__), + T1mag(i__),e_T1mag(i__),v3_6mag(i__),e_3_6mag(i__), + v4_5mag(i__),e_4_5mag(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 'n4594gcs.dat' ! NGC 4594 globular cluster photometry C Format for file interpretation 2 format( + A3,4X,I2,1X,I2,1X,F5.2,3X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,1X,F4.2, + 1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.3,1X,F5.2,1X, + F5.3,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'n4594gcs.dat') write(6,*) '....Loading file: n4594gcs.dat' do i__=1,660 read(1,'(A95)')ar__1 read(ar__1,2) + v_SLS2006_(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__), + DEd_1(i__),DEm_1(i__),DEs_1(i__),Bmag_1(i__),e_Bmag_1(i__), + Vmag_1(i__),e_Vmag_1(i__),Rmag_1(i__),e_Rmag_1(i__), + v3_6mag_1(i__),e_3_6mag_1(i__),v4_5mag_1(i__),e_4_5mag_1(i__), + Rhl(i__) RAdeg_1(i__) = rNULL__ DEdeg_1(i__) = rNULL__ c Derive coordinates RAdeg_1 and DEdeg_1 from input data c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown) if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15. if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4. if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240. if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__) c ..............Just test output........... write(6,2) + v_SLS2006_(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__), + DEd_1(i__),DEm_1(i__),DEs_1(i__),Bmag_1(i__),e_Bmag_1(i__), + Vmag_1(i__),e_Vmag_1(i__),Rmag_1(i__),e_Rmag_1(i__), + v3_6mag_1(i__),e_3_6mag_1(i__),v4_5mag_1(i__),e_4_5mag_1(i__), + Rhl(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end