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