Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/592/A18 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 2026-Oct-07
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/A+A/592/A18 Molecule and grain abundances (Marchand+, 2016)
*================================================================================
*Chemical solver to compute molecule and grain abundances and non-ideal MHD
*resistivities in prestellar core-collapse calculations.
* Marchand P., Masson J., Chabrier G., Hennebelle P., Commercon B., Vaytet N.
* <Astron. Astrophys. 592, A18 (2016)>
* =2016A&A...592A..18M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'abund.dat' ! Table of chemical abundances
integer*4 nr__
parameter (nr__=120000) ! Number of records
character*804 ar__ ! Full-size record
real*4 rho (nr__) ! (cm-3) Density
real*4 T (nr__) ! (K) Temperature
real*4 xi (nr__) ! (s-1) Ionisation rate
real*4 e (nr__) ! Abundance of electrons
real*4 Mg (nr__) ! Abundance of ion metal (Mg, etc.)
real*4 HCO_ (nr__) ! Abundance of molecular ion HCO^+^
real*4 H3_ (nr__) ! Abundance of molecular ion H_3_^+^
real*4 H_ (nr__) ! Abundance of molecular ion H^+^
real*4 C_ (nr__) ! Abundance of molecular ion C^+^
real*4 He_ (nr__) ! Abundance of molecular ion He^+^
real*4 K_ (nr__) ! Abundance of molecular ion K^+^
real*4 Na_ (nr__) ! Abundance of molecular ion Na^+^
real*4 grains_1 (nr__) ! Abundances of grains+ for bin 1
real*4 grains_1_1 (nr__) ! Abundances of grains- for bin 1
real*4 grains01 (nr__) ! Abundances of grains0 for bin 1
real*4 grains_2 (nr__) ! Abundances of grains+ for bin 2
real*4 grains_2_1 (nr__) ! Abundances of grains- for bin 2
real*4 grains02 (nr__) ! Abundances of grains0 for bin 2
real*4 grains_3 (nr__) ! Abundances of grains+ for bin 3
real*4 grains_3_1 (nr__) ! Abundances of grains- for bin 3
real*4 grains03 (nr__) ! Abundances of grains0 for bin 3
real*4 grains_4 (nr__) ! Abundances of grains+ for bin 4
real*4 grains_4_1 (nr__) ! Abundances of grains- for bin 4
real*4 grains04 (nr__) ! Abundances of grains0 for bin 4
real*4 grains_5 (nr__) ! Abundances of grains+ for bin 5
real*4 grains_5_1 (nr__) ! Abundances of grains- for bin 5
real*4 grains05 (nr__) ! Abundances of grains0 for bin 5
real*4 grains_t (nr__) ! Total abundance of grains+
real*4 grains_t_1 (nr__) ! Total abundance of grains-
real*4 grains0t (nr__) ! Total abundance of grains0
real*4 Tint (nr__) ! (s) Integration time
C=============================================================================
C Loading file 'abund.dat' ! Table of chemical abundances
C Format for file interpretation
1 format(
+ 1X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,
+ 2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,
+ 2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,
+ 2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,
+ 2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,2X,E24.17,
+ 2X,E23.17)
C Effective file loading
open(unit=1,status='old',file=
+'abund.dat')
write(6,*) '....Loading file: abund.dat'
do i__=1,120000
read(1,'(A804)')ar__
read(ar__,1)
+ rho(i__),T(i__),xi(i__),e(i__),Mg(i__),HCO_(i__),H3_(i__),
+ H_(i__),C_(i__),He_(i__),K_(i__),Na_(i__),grains_1(i__),
+ grains_1_1(i__),grains01(i__),grains_2(i__),grains_2_1(i__),
+ grains02(i__),grains_3(i__),grains_3_1(i__),grains03(i__),
+ grains_4(i__),grains_4_1(i__),grains04(i__),grains_5(i__),
+ grains_5_1(i__),grains05(i__),grains_t(i__),grains_t_1(i__),
+ grains0t(i__),Tint(i__)
c ..............Just test output...........
write(6,1)
+ rho(i__),T(i__),xi(i__),e(i__),Mg(i__),HCO_(i__),H3_(i__),
+ H_(i__),C_(i__),He_(i__),K_(i__),Na_(i__),grains_1(i__),
+ grains_1_1(i__),grains01(i__),grains_2(i__),grains_2_1(i__),
+ grains02(i__),grains_3(i__),grains_3_1(i__),grains03(i__),
+ grains_4(i__),grains_4_1(i__),grains04(i__),grains_5(i__),
+ grains_5_1(i__),grains05(i__),grains_t(i__),grains_t_1(i__),
+ grains0t(i__),Tint(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end