Conversion of standardized ReadMe file for
file /./ftp/cats/J/AZh/74/63 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-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/AZh/74/63 Cluster analysis of pulsar parameters (Malov+, 1997) *================================================================================ *Cluster analysis of pulsar parameters * Malov I.F., Malov O.I. * <Astron. Zh. 74, 63 (1997)> * =1997AZh....74...63M (SIMBAD/NED BibCode) * =1997ARep...41...56M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Pulsar catalogue integer*4 nr__ parameter (nr__=218) ! Number of records character*100 ar__ ! Full-size record integer*4 Seq (nr__) ! [001-218] Running number character*9 PSR (nr__) ! Name of pulsar real*4 Per (nr__) ! (s) Period real*8 Pdot (nr__) ! (fs/s) Period derivative real*8 w50 (nr__) ! (ms) ? Profile width at 50% of peak real*4 logL (nr__) ! ([10-7W]) Log of integrated radio luminosity real*4 logB (nr__) ! ([gauss]) Log of magnetic field in gauss at the surface real*4 logAge (nr__) ! ([yr]) Log of characteristic age in years real*4 v_z_ (nr__) ! (kpc) Distance from the galactic plane real*4 alf2 (nr__) ! ? Spectral index of the linear part of * pulsar spectrum real*4 numax (nr__) ! (GHz) ? Frequency of the maximum of the spectrum real*4 nucut (nr__) ! (GHz) ? Frequency of the cutoff in the spectrum real*4 logBlc (nr__) ! ([gauss]) Magnetic field in gauss at the light cylinder real*4 beta (nr__) ! (deg) ? Angle between rotation and magnetic axes real*4 C (nr__) ! ? Maximum derivative of linear polarization * position angle real*4 logeta (nr__) ! Log of efficiency of the transformation of * rotation energy into radio emission * (always negative) C============================================================================= C Loading file 'table1.dat' ! Pulsar catalogue C Format for file interpretation 1 format( + I3,1X,A9,1X,F5.3,1X,F9.5,1X,F8.4,1X,F5.2,3X,F5.2,1X,F4.2,1X, + F4.2,2X,F3.1,1X,F4.2,1X,F4.2,2X,F4.2,1X,F4.1,1X,F5.2,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,218 read(1,'(A100)')ar__ read(ar__,1) + Seq(i__),PSR(i__),Per(i__),Pdot(i__),w50(i__),logL(i__), + logB(i__),logAge(i__),v_z_(i__),alf2(i__),numax(i__), + nucut(i__),logBlc(i__),beta(i__),C(i__),logeta(i__) if(ar__(31:38) .EQ. '') w50(i__) = rNULL__ if(ar__(65:67) .EQ. '') alf2(i__) = rNULL__ if(ar__(69:72) .EQ. '') numax(i__) = rNULL__ if(ar__(74:77) .EQ. '') nucut(i__) = rNULL__ if(ar__(85:88) .EQ. '') beta(i__) = rNULL__ if(ar__(90:94) .EQ. '') C(i__) = rNULL__ c ..............Just test output........... write(6,1) + Seq(i__),PSR(i__),Per(i__),Pdot(i__),w50(i__),logL(i__), + logB(i__),logAge(i__),v_z_(i__),alf2(i__),numax(i__), + nucut(i__),logBlc(i__),beta(i__),C(i__),logeta(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end