Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/399/425 into FORTRAN code for reading data files line by line.
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-Apr-17 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/MNRAS/399/425 (De)excitation rates of ortho-NH3 and para-NH3 (Faure+, 2009) *================================================================================ *On the robustness of the ammonia thermometer. * Maret S., Faure A., Scifoni E., Wiesenfeld L. * <Mon. Not. R. Astron. Soc. 399, 425 (2009)> * =2009MNRAS.399..425M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'onh3lev.dat' ! Energy levels of ortho-NH3, in cm^-1^ integer*4 nr__ parameter (nr__=6) ! Number of records character*37 ar__ ! Full-size record integer*4 Level ! [1,10] level number real*8 E ! (cm-1) level energy real*8 F ! statistical weight character*9 N ! quantum numbers c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'pnh3lev.dat' ! Energy levels of para-NH3, in cm^-1^ integer*4 nr__1 parameter (nr__1=10) ! Number of records character*37 ar__1 ! Full-size record integer*4 Level_1 ! [1,10] level number real*8 E_1 ! (cm-1) level energy real*8 F_1 ! statistical weight character*9 N_1 ! quantum numbers c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'onh3ph2.dat' ! De-excitation rates (in cm^3^/s) of ortho-NH3 in collision with para-H2 integer*4 nr__2 parameter (nr__2=15) ! Number of records character*208 ar__2 ! Full-size record integer*4 Trans ! [1,45] Transition number integer*4 i ! [2,10] Initial level integer*4 f_2 ! [1,10] Final level real*4 k5 ! (cm3/s) Rate coefficient at 5K real*4 k10 ! (cm3/s) Rate coefficient at 10K real*4 k15 ! (cm3/s) Rate coefficient at 15K real*4 k20 ! (cm3/s) Rate coefficient at 20K real*4 k25 ! (cm3/s) Rate coefficient at 25K real*4 k30 ! (cm3/s) Rate coefficient at 30K real*4 k35 ! (cm3/s) Rate coefficient at 35K real*4 k40 ! (cm3/s) Rate coefficient at 40K real*4 k45 ! (cm3/s) Rate coefficient at 45K real*4 k50 ! (cm3/s) Rate coefficient at 50K real*4 k55 ! (cm3/s) Rate coefficient at 55K real*4 k60 ! (cm3/s) Rate coefficient at 60K real*4 k65 ! (cm3/s) Rate coefficient at 65K real*4 k70 ! (cm3/s) Rate coefficient at 70K real*4 k75 ! (cm3/s) Rate coefficient at 75K real*4 k80 ! (cm3/s) Rate coefficient at 80K real*4 k85 ! (cm3/s) Rate coefficient at 85K real*4 k90 ! (cm3/s) Rate coefficient at 90K real*4 k95 ! (cm3/s) Rate coefficient at 95K real*4 k100 ! (cm3/s) Rate coefficient at 100K c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'pnh3ph2.dat' ! De-excitation rates (in cm^3^/s) of para-NH3 in collision with para-H2 integer*4 nr__3 parameter (nr__3=45) ! Number of records character*208 ar__3 ! Full-size record integer*4 Trans_1 ! [1,45] Transition number integer*4 i_1 ! [2,10] Initial level integer*4 f_3 ! [1,10] Final level real*4 k5_1 ! (cm3/s) Rate coefficient at 5K real*4 k10_1 ! (cm3/s) Rate coefficient at 10K real*4 k15_1 ! (cm3/s) Rate coefficient at 15K real*4 k20_1 ! (cm3/s) Rate coefficient at 20K real*4 k25_1 ! (cm3/s) Rate coefficient at 25K real*4 k30_1 ! (cm3/s) Rate coefficient at 30K real*4 k35_1 ! (cm3/s) Rate coefficient at 35K real*4 k40_1 ! (cm3/s) Rate coefficient at 40K real*4 k45_1 ! (cm3/s) Rate coefficient at 45K real*4 k50_1 ! (cm3/s) Rate coefficient at 50K real*4 k55_1 ! (cm3/s) Rate coefficient at 55K real*4 k60_1 ! (cm3/s) Rate coefficient at 60K real*4 k65_1 ! (cm3/s) Rate coefficient at 65K real*4 k70_1 ! (cm3/s) Rate coefficient at 70K real*4 k75_1 ! (cm3/s) Rate coefficient at 75K real*4 k80_1 ! (cm3/s) Rate coefficient at 80K real*4 k85_1 ! (cm3/s) Rate coefficient at 85K real*4 k90_1 ! (cm3/s) Rate coefficient at 90K real*4 k95_1 ! (cm3/s) Rate coefficient at 95K real*4 k100_1 ! (cm3/s) Rate coefficient at 100K C============================================================================= C Loading file 'onh3lev.dat' ! Energy levels of ortho-NH3, in cm^-1^ C Format for file interpretation 1 format(I5,F16.9,F7.1,A9) C Effective file loading open(unit=1,status='old',file= +'onh3lev.dat') write(6,*) '....Loading file: onh3lev.dat' do i__=1,6 read(1,'(A37)')ar__ read(ar__,1)Level,E,F,N c ..............Just test output........... write(6,1)Level,E,F,N c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'pnh3lev.dat' ! Energy levels of para-NH3, in cm^-1^ C Format for file interpretation 2 format(I5,F16.9,F7.1,A9) C Effective file loading open(unit=1,status='old',file= +'pnh3lev.dat') write(6,*) '....Loading file: pnh3lev.dat' do i__=1,10 read(1,'(A37)')ar__1 read(ar__1,2)Level_1,E_1,F_1,N_1 c ..............Just test output........... write(6,2)Level_1,E_1,F_1,N_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'onh3ph2.dat' ! De-excitation rates (in cm^3^/s) of ortho-NH3 * in collision with para-H2 C Format for file interpretation 3 format( + I2,1X,I2,1X,I2,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X, + E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3, + 1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3) C Effective file loading open(unit=1,status='old',file= +'onh3ph2.dat') write(6,*) '....Loading file: onh3ph2.dat' do i__=1,15 read(1,'(A208)')ar__2 read(ar__2,3) + Trans,i,f_2,k5,k10,k15,k20,k25,k30,k35,k40,k45,k50,k55,k60, + k65,k70,k75,k80,k85,k90,k95,k100 c ..............Just test output........... write(6,3) + Trans,i,f_2,k5,k10,k15,k20,k25,k30,k35,k40,k45,k50,k55,k60, + k65,k70,k75,k80,k85,k90,k95,k100 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'pnh3ph2.dat' ! De-excitation rates (in cm^3^/s) of para-NH3 in * collision with para-H2 C Format for file interpretation 4 format( + I2,1X,I2,1X,I2,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X, + E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3, + 1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3,1X,E9.3) C Effective file loading open(unit=1,status='old',file= +'pnh3ph2.dat') write(6,*) '....Loading file: pnh3ph2.dat' do i__=1,45 read(1,'(A208)')ar__3 read(ar__3,4) + Trans_1,i_1,f_3,k5_1,k10_1,k15_1,k20_1,k25_1,k30_1,k35_1, + k40_1,k45_1,k50_1,k55_1,k60_1,k65_1,k70_1,k75_1,k80_1,k85_1, + k90_1,k95_1,k100_1 c ..............Just test output........... write(6,4) + Trans_1,i_1,f_3,k5_1,k10_1,k15_1,k20_1,k25_1,k30_1,k35_1, + k40_1,k45_1,k50_1,k55_1,k60_1,k65_1,k70_1,k75_1,k80_1,k85_1, + k90_1,k95_1,k100_1 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end