FORTRAN Generation
(/./ftp/cats/J/ApJ/416/276)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/416/276 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/ApJ/416/276  Eq. of State in Magnetic Neutron Stars    (Rognvaldsson+ 1993)
*================================================================================
*Thomas-Fermi calculations of atoms and matter in magnetic neutron
*stars: effects of higher Landau bands
*   Rognvaldsson O.E., Fushiki I., Gudmundsson E.H., Pethick C.J.,
*   Yngvason J.
*   <Astrophys. J., 416, 276 (1993)>
*   =1993ApJ...416..276R
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b10_0.dat'	! 1 million (10^6^) Tesla

      integer*4 nr__
      parameter (nr__=40)	! Number of records
      character*68 ar__   	! Full-size record

      real*4        density    (nr__) ! (kg/m3) density of matter
      real*4        radius     (nr__) ! (m) radius of cell
      real*4        pressure   (nr__) ! (Pa) pressure
      real*4        eatom      (nr__) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot    (nr__) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12       (nr__) ! (m) depth.(surface gravity)/(10^12 m/s^2)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b10_5.dat'	! 3 million (10^6.5^) Tesla

      integer*4 nr__1
      parameter (nr__1=46)	! Number of records
      character*68 ar__1  	! Full-size record

      real*4        density_1  (nr__1) ! (kg/m3) density of matter
      real*4        radius_1   (nr__1) ! (m) radius of cell
      real*4        pressure_1 (nr__1) ! (Pa) pressure
      real*4        eatom_1    (nr__1) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot_1  (nr__1) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12_1     (nr__1) ! (m) depth.(surface gravity)/(10^12 m/s^2)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b11_0.dat'	! 10 million (10^7^) Tesla

      integer*4 nr__2
      parameter (nr__2=41)	! Number of records
      character*68 ar__2  	! Full-size record

      real*4        density_2  (nr__2) ! (kg/m3) density of matter
      real*4        radius_2   (nr__2) ! (m) radius of cell
      real*4        pressure_2 (nr__2) ! (Pa) pressure
      real*4        eatom_2    (nr__2) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot_2  (nr__2) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12_2     (nr__2) ! (m) depth.(surface gravity)/(10^12 m/s^2)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b11_5.dat'	! 30 million (10^7.5^) Tesla

      integer*4 nr__3
      parameter (nr__3=42)	! Number of records
      character*68 ar__3  	! Full-size record

      real*4        density_3  (nr__3) ! (kg/m3) density of matter
      real*4        radius_3   (nr__3) ! (m) radius of cell
      real*4        pressure_3 (nr__3) ! (Pa) pressure
      real*4        eatom_3    (nr__3) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot_3  (nr__3) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12_3     (nr__3) ! (m) depth.(surface gravity)/(10^12 m/s^2)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b12_0.dat'	! 100 million (10^8^) Tesla

      integer*4 nr__4
      parameter (nr__4=30)	! Number of records
      character*68 ar__4  	! Full-size record

      real*4        density_4  (nr__4) ! (kg/m3) density of matter
      real*4        radius_4   (nr__4) ! (m) radius of cell
      real*4        pressure_4 (nr__4) ! (Pa) pressure
      real*4        eatom_4    (nr__4) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot_4  (nr__4) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12_4     (nr__4) ! (m) depth.(surface gravity)/(10^12 m/s^2)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b12_5.dat'	! 300 million (10^8.5^) Tesla

      integer*4 nr__5
      parameter (nr__5=27)	! Number of records
      character*68 ar__5  	! Full-size record

      real*4        density_5  (nr__5) ! (kg/m3) density of matter
      real*4        radius_5   (nr__5) ! (m) radius of cell
      real*4        pressure_5 (nr__5) ! (Pa) pressure
      real*4        eatom_5    (nr__5) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot_5  (nr__5) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12_5     (nr__5) ! (m) depth.(surface gravity)/(10^12 m/s^2)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'b13_0.dat'	! 1000 million (10^9^) Tesla

      integer*4 nr__6
      parameter (nr__6=20)	! Number of records
      character*68 ar__6  	! Full-size record

      real*4        density_6  (nr__6) ! (kg/m3) density of matter
      real*4        radius_6   (nr__6) ! (m) radius of cell
      real*4        pressure_6 (nr__6) ! (Pa) pressure
      real*4        eatom_6    (nr__6) ! (keV) total energy of atom (in 1.602E-16J)
      real*4        chempot_6  (nr__6) ! (keV) chemical potential (in 1.602E-16J)
      real*8        zg12_6     (nr__6) ! (m) depth.(surface gravity)/(10^12 m/s^2)

C=============================================================================

C  Loading file 'b10_0.dat'	! 1 million (10^6^) Tesla

C  Format for file interpretation

    1 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b10_0.dat')
      write(6,*) '....Loading file: b10_0.dat'
      do i__=1,40
        read(1,'(A68)')ar__
        read(ar__,1)
     +  density(i__),radius(i__),pressure(i__),eatom(i__),
     +  chempot(i__),zg12(i__)
c    ..............Just test output...........
        write(6,1)
     +  density(i__),radius(i__),pressure(i__),eatom(i__),
     +  chempot(i__),zg12(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'b10_5.dat'	! 3 million (10^6.5^) Tesla

C  Format for file interpretation

    2 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b10_5.dat')
      write(6,*) '....Loading file: b10_5.dat'
      do i__=1,46
        read(1,'(A68)')ar__1
        read(ar__1,2)
     +  density_1(i__),radius_1(i__),pressure_1(i__),eatom_1(i__),
     +  chempot_1(i__),zg12_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  density_1(i__),radius_1(i__),pressure_1(i__),eatom_1(i__),
     +  chempot_1(i__),zg12_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'b11_0.dat'	! 10 million (10^7^) Tesla

C  Format for file interpretation

    3 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b11_0.dat')
      write(6,*) '....Loading file: b11_0.dat'
      do i__=1,41
        read(1,'(A68)')ar__2
        read(ar__2,3)
     +  density_2(i__),radius_2(i__),pressure_2(i__),eatom_2(i__),
     +  chempot_2(i__),zg12_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  density_2(i__),radius_2(i__),pressure_2(i__),eatom_2(i__),
     +  chempot_2(i__),zg12_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'b11_5.dat'	! 30 million (10^7.5^) Tesla

C  Format for file interpretation

    4 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b11_5.dat')
      write(6,*) '....Loading file: b11_5.dat'
      do i__=1,42
        read(1,'(A68)')ar__3
        read(ar__3,4)
     +  density_3(i__),radius_3(i__),pressure_3(i__),eatom_3(i__),
     +  chempot_3(i__),zg12_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  density_3(i__),radius_3(i__),pressure_3(i__),eatom_3(i__),
     +  chempot_3(i__),zg12_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'b12_0.dat'	! 100 million (10^8^) Tesla

C  Format for file interpretation

    5 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b12_0.dat')
      write(6,*) '....Loading file: b12_0.dat'
      do i__=1,30
        read(1,'(A68)')ar__4
        read(ar__4,5)
     +  density_4(i__),radius_4(i__),pressure_4(i__),eatom_4(i__),
     +  chempot_4(i__),zg12_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  density_4(i__),radius_4(i__),pressure_4(i__),eatom_4(i__),
     +  chempot_4(i__),zg12_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'b12_5.dat'	! 300 million (10^8.5^) Tesla

C  Format for file interpretation

    6 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b12_5.dat')
      write(6,*) '....Loading file: b12_5.dat'
      do i__=1,27
        read(1,'(A68)')ar__5
        read(ar__5,6)
     +  density_5(i__),radius_5(i__),pressure_5(i__),eatom_5(i__),
     +  chempot_5(i__),zg12_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  density_5(i__),radius_5(i__),pressure_5(i__),eatom_5(i__),
     +  chempot_5(i__),zg12_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'b13_0.dat'	! 1000 million (10^9^) Tesla

C  Format for file interpretation

    7 format(E8.2,1X,E11.5,1X,E11.5,1X,E12.5,1X,E12.5,F10.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b13_0.dat')
      write(6,*) '....Loading file: b13_0.dat'
      do i__=1,20
        read(1,'(A68)')ar__6
        read(ar__6,7)
     +  density_6(i__),radius_6(i__),pressure_6(i__),eatom_6(i__),
     +  chempot_6(i__),zg12_6(i__)
c    ..............Just test output...........
        write(6,7)
     +  density_6(i__),radius_6(i__),pressure_6(i__),eatom_6(i__),
     +  chempot_6(i__),zg12_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end