Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/121/1192 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-Mar-29 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/AJ/121/1192 Proper-motion measurements with the VLA (McGary+, 2001) *================================================================================ *Proper-motion measurements with the VLA. *I. Wide-field imaging and pulse-gating techniques. * McGary R.S., Brisken W.F., Fruchter A.S., Goss W.M., Thorsett S.E. * <Astron. J. 121, 1192 (2001)> * =2001AJ....121.1192M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Reference sources for proper motions integer*4 nr__ parameter (nr__=44) ! Number of records character*74 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*8 Name (nr__) ! Pulsar name, Bhhmm+dd (1) real*4 F20 (nr__) ! (mJy) The 20cm flux density of reference real*4 e_F20 (nr__) ! (mJy) Uncertainty in F20 integer*4 RAh (nr__) ! (h) Right Ascension of reference (J2000) integer*4 RAm (nr__) ! (min) Right Ascension of reference (J2000) real*4 RAs (nr__) ! (s) Right Ascension of reference (J2000) real*4 e_RAs (nr__) ! (s) Uncertainty in RAs character*1 DE_ (nr__) ! Sign of the reference Declination (J2000) integer*4 DEd (nr__) ! (deg) Declination of reference (J2000) integer*4 DEm (nr__) ! (arcmin) Declination of reference (J2000) real*4 DEs (nr__) ! (arcsec) Declination of reference (J2000) real*4 e_DEs (nr__) ! (arcsec) Uncertainty in DEs integer*4 dRA (nr__) ! (arcsec) RA offset from pulsar position integer*4 dDE (nr__) ! (arcsec) Dec offset from pulsar position integer*4 pmRA (nr__) ! (mas/yr) ? Proper Motion along RA integer*4 e_pmRA (nr__) ! (mas/yr) ? Uncertainty in pmRA integer*4 pmDE (nr__) ! (mas/yr) ? Proper Motion along Dec integer*4 e_pmDE (nr__) ! (mas/yr) ? Uncertainty in pmDE *Note (1): The objects are identified in Simbad as [MBF2001] JHHMMSS.ss+DDMMSS.s c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Proper motions of B0919+06, B1237+25 and B1937-26 integer*4 nr__1 parameter (nr__1=3) ! Number of records character*69 ar__1 ! Full-size record character*8 Name_1 (nr__1) ! Pulsar Name real*4 F20_1 (nr__1) ! (mJy) The 20cm flux density real*4 e_F20_1 (nr__1) ! (mJy) Uncertainty in F20 real*4 GF20 (nr__1) ! (mJy) Gated 20cm flux density real*4 e_GF20 (nr__1) ! (mJy) Uncertainty in GF20 integer*4 Epoch1 (nr__1) ! (yr) Epoch of observation (1900+) character*1 n_Epoch1 (nr__1) ! [a] a: gated observation integer*4 Epoch2 (nr__1) ! (yr) Epoch of observation (1900+) character*1 n_Epoch2 (nr__1) ! [a] a: gated observation integer*4 Epoch3 (nr__1) ! (yr) Epoch of observation (1900+) character*1 n_Epoch3 (nr__1) ! [a] a: gated observation integer*4 Epoch4 (nr__1) ! (yr) ? Epoch of observation (1900+) character*1 n_Epoch4 (nr__1) ! [a] a: gated observation integer*4 Num (nr__1) ! Number of reference objects used real*4 pmRA_1 (nr__1) ! (mas/yr) Proper motion along RA real*4 e_pmRA_1 (nr__1) ! (mas/yr) Uncertainty in pmRA real*4 pmDE_1 (nr__1) ! (mas/yr) Proper motion along Dec real*4 e_pmDE_1 (nr__1) ! (mas/yr) Uncertainty in pmDE real*4 Cov (nr__1) ! Covariance of RA & Dec proper motion * components C============================================================================= C Loading file 'table1.dat' ! Reference sources for proper motions C Format for file interpretation 1 format( + A8,1X,F4.1,1X,F3.1,1X,I2,1X,I2,1X,F5.2,1X,F3.2,1X,A1,I2,1X,I2, + 1X,F4.1,1X,F2.1,1X,I5,1X,I5,1X,I3,1X,I2,1X,I2,1X,I2) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,44 read(1,'(A74)')ar__ read(ar__,1) + Name(i__),F20(i__),e_F20(i__),RAh(i__),RAm(i__),RAs(i__), + e_RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),e_DEs(i__), + dRA(i__),dDE(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__) if(ar__(62:64) .EQ. '') pmRA(i__) = iNULL__ if(ar__(66:67) .EQ. '') e_pmRA(i__) = iNULL__ if(ar__(69:70) .EQ. '') pmDE(i__) = iNULL__ if(ar__(72:73) .EQ. '') e_pmDE(i__) = iNULL__ 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__),F20(i__),e_F20(i__),RAh(i__),RAm(i__),RAs(i__), + e_RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),e_DEs(i__), + dRA(i__),dDE(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_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 'table2.dat' ! Proper motions of B0919+06, B1237+25 and B1937-26 C Format for file interpretation 2 format( + A8,1X,F3.1,1X,F3.1,1X,F5.1,1X,F3.1,1X,I2,A1,1X,I2,A1,1X,I2,A1, + 1X,I2,A1,1X,I2,1X,F6.1,1X,F3.1,1X,F4.1,1X,F3.1,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,3 read(1,'(A69)')ar__1 read(ar__1,2) + Name_1(i__),F20_1(i__),e_F20_1(i__),GF20(i__),e_GF20(i__), + Epoch1(i__),n_Epoch1(i__),Epoch2(i__),n_Epoch2(i__), + Epoch3(i__),n_Epoch3(i__),Epoch4(i__),n_Epoch4(i__),Num(i__), + pmRA_1(i__),e_pmRA_1(i__),pmDE_1(i__),e_pmDE_1(i__),Cov(i__) if(ar__1(40:41) .EQ. '') Epoch4(i__) = iNULL__ c ..............Just test output........... write(6,2) + Name_1(i__),F20_1(i__),e_F20_1(i__),GF20(i__),e_GF20(i__), + Epoch1(i__),n_Epoch1(i__),Epoch2(i__),n_Epoch2(i__), + Epoch3(i__),n_Epoch3(i__),Epoch4(i__),n_Epoch4(i__),Num(i__), + pmRA_1(i__),e_pmRA_1(i__),pmDE_1(i__),e_pmDE_1(i__),Cov(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end