Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/504/853 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-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/A+A/504/853 Molecular lines in 5 massive dense cores (Herpin, 2009) *================================================================================ *S-bearing molecules in massive dense cores * Herpin F., Marseille M., Wakelam V., Bontemps S., Lis D.C. * <Astron. Astrophys. 504, 853 (2009)> * =2009A&A...504..853H C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Observed line emissions parameters for W43MM1 integer*4 nr__ parameter (nr__=25) ! Number of records character*115 ar__ ! Full-size record character*9 Species ! Species character*17 Trans ! Observed line transition real*4 dV ! (km/s) ? Linewidth real*4 e_dV ! (km/s) ? rms uncertainty on dV character*1 u_dV ! [?] Uncertainty flag on dV real*4 Tmb ! (K) ? Main beam temperature character*1 u_Tmb ! [?] Uncertainty flag on Tmb integer*4 e_Tmb ! (mK) ? rms uncertainty on Tmb character*1 l_Flux ! Limit flag on Flux real*4 Flux ! (K.km/s) ? Integrated flux real*4 e_Flux ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux ! [?] Uncertainty flag on Flux real*4 dV3s ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s ! [?] Uncertainty flag on dV3s real*4 dV1_2 ! (km/s) ? half-power linewidth character*1 u_dV1_2 ! [?] Uncertainty flag on dV1/2 character*1 Flag ! [abcdfg] Flag on profile (1) character*14 Com ! Comments character*13 FileName ! Name of the fits file in subdirectory fits *Note (1): Flags as follows: * a = self-absorption at dV_source * b = line blend with H_2_CO line from image band * c = unidentified line at 34.87km/s (216.716GHz) * d = + other line at 29.1km/s (135.78223GHz), * Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17 * f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s * g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Observed line emissions parameters for IRAS18264-1152 integer*4 nr__1 parameter (nr__1=25) ! Number of records character*115 ar__1 ! Full-size record character*9 Species_1 ! Species character*17 Trans_1 ! Observed line transition real*4 dV_1 ! (km/s) ? Linewidth real*4 e_dV_1 ! (km/s) ? rms uncertainty on dV character*1 u_dV_1 ! [?] Uncertainty flag on dV real*4 Tmb_1 ! (K) ? Main beam temperature character*1 u_Tmb_1 ! [?] Uncertainty flag on Tmb integer*4 e_Tmb_1 ! (mK) ? rms uncertainty on Tmb character*1 l_Flux_1 ! Limit flag on Flux real*4 Flux_1 ! (K.km/s) ? Integrated flux real*4 e_Flux_1 ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux_1 ! [?] Uncertainty flag on Flux real*4 dV3s_1 ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s_1 ! [?] Uncertainty flag on dV3s real*4 dV1_2_1 ! (km/s) ? half-power linewidth character*1 u_dV1_2_1 ! [?] Uncertainty flag on dV1/2 character*1 Flag_1 ! [abcdfg] Flag on profile (1) character*14 Com_1 ! Comments character*13 FileName_1 ! Name of the fits file in subdirectory fits *Note (1): Flags as follows: * a = self-absorption at dV_source * b = line blend with H_2_CO line from image band * c = unidentified line at 34.87km/s (216.716GHz) * d = + other line at 29.1km/s (135.78223GHz), * Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17 * f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s * g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Observed line emissions parameters for IRAS05358+3543 integer*4 nr__2 parameter (nr__2=25) ! Number of records character*115 ar__2 ! Full-size record character*9 Species_2 ! Species character*17 Trans_2 ! Observed line transition real*4 dV_2 ! (km/s) ? Linewidth real*4 e_dV_2 ! (km/s) ? rms uncertainty on dV character*1 u_dV_2 ! [?] Uncertainty flag on dV real*4 Tmb_2 ! (K) ? Main beam temperature character*1 u_Tmb_2 ! [?] Uncertainty flag on Tmb integer*4 e_Tmb_2 ! (mK) ? rms uncertainty on Tmb character*1 l_Flux_2 ! Limit flag on Flux real*4 Flux_2 ! (K.km/s) ? Integrated flux real*4 e_Flux_2 ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux_2 ! [?] Uncertainty flag on Flux real*4 dV3s_2 ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s_2 ! [?] Uncertainty flag on dV3s real*4 dV1_2_2 ! (km/s) ? half-power linewidth character*1 u_dV1_2_2 ! [?] Uncertainty flag on dV1/2 character*1 Flag_2 ! [abcdfg] Flag on profile (1) character*14 Com_2 ! Comments character*13 FileName_2 ! Name of the fits file in subdirectory fits *Note (1): Flags as follows: * a = self-absorption at dV_source * b = line blend with H_2_CO line from image band * c = unidentified line at 34.87km/s (216.716GHz) * d = + other line at 29.1km/s (135.78223GHz), * Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17 * f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s * g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Observed line emissions parameters for IRAS18162-2048 integer*4 nr__3 parameter (nr__3=25) ! Number of records character*115 ar__3 ! Full-size record character*9 Species_3 ! Species character*17 Trans_3 ! Observed line transition real*4 dV_3 ! (km/s) ? Linewidth real*4 e_dV_3 ! (km/s) ? rms uncertainty on dV character*1 u_dV_3 ! [?] Uncertainty flag on dV real*4 Tmb_3 ! (K) ? Main beam temperature character*1 u_Tmb_3 ! [?] Uncertainty flag on Tmb integer*4 e_Tmb_3 ! (mK) ? rms uncertainty on Tmb character*1 l_Flux_3 ! Limit flag on Flux real*4 Flux_3 ! (K.km/s) ? Integrated flux real*4 e_Flux_3 ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux_3 ! [?] Uncertainty flag on Flux real*4 dV3s_3 ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s_3 ! [?] Uncertainty flag on dV3s real*4 dV1_2_3 ! (km/s) ? half-power linewidth character*1 u_dV1_2_3 ! [?] Uncertainty flag on dV1/2 character*1 Flag_3 ! [abcdfg] Flag on profile (1) character*14 Com_3 ! Comments character*13 FileName_3 ! Name of the fits file in subdirectory fits *Note (1): Flags as follows: * a = self-absorption at dV_source * b = line blend with H_2_CO line from image band * c = unidentified line at 34.87km/s (216.716GHz) * d = + other line at 29.1km/s (135.78223GHz), * Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17 * f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s * g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! Opacities derived from modeling the molecular emission of the sources integer*4 nr__4 parameter (nr__4=24) ! Number of records character*55 ar__4 ! Full-size record character*9 Species_4 ! Species character*17 Trans_4 ! Observed line transition real*4 Opac1 ! ? Line opacity for source W43MM1 character*2 n_Opac1 ! [ND] ND for no data real*4 Opac2 ! Line opacity for source IRAS18264-1152 real*4 Opac3 ! ? Line opacity for source IRAS05358+3543 character*2 n_Opac3 ! [ND] ND for no data real*4 Opac4 ! ? Line opacity for source IRAS18162-2048 character*2 n_Opac4 ! [ND] ND for no data C============================================================================= C Loading file 'table6.dat' ! Observed line emissions parameters for W43MM1 C Format for file interpretation 1 format( + A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3, + 1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,25 read(1,'(A115)')ar__ read(ar__,1) + Species,Trans,dV,e_dV,u_dV,Tmb,u_Tmb,e_Tmb,l_Flux,Flux,e_Flux, + u_e_Flux,dV3s,u_dV3s,dV1_2,u_dV1_2,Flag,Com,FileName if(ar__(30:33) .EQ. '') dV = rNULL__ if(ar__(35:38) .EQ. '') e_dV = rNULL__ if(ar__(42:45) .EQ. '') Tmb = rNULL__ if(ar__(48:50) .EQ. '') e_Tmb = iNULL__ if(ar__(54:59) .EQ. '') Flux = rNULL__ if(ar__(61:64) .EQ. '') e_Flux = rNULL__ if(ar__(69:73) .EQ. '') dV3s = rNULL__ if(ar__(77:80) .EQ. '') dV1_2 = rNULL__ c ..............Just test output........... write(6,1) + Species,Trans,dV,e_dV,u_dV,Tmb,u_Tmb,e_Tmb,l_Flux,Flux,e_Flux, + u_e_Flux,dV3s,u_dV3s,dV1_2,u_dV1_2,Flag,Com,FileName c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Observed line emissions parameters for IRAS18264-1152 C Format for file interpretation 2 format( + A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3, + 1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,25 read(1,'(A115)')ar__1 read(ar__1,2) + Species_1,Trans_1,dV_1,e_dV_1,u_dV_1,Tmb_1,u_Tmb_1,e_Tmb_1, + l_Flux_1,Flux_1,e_Flux_1,u_e_Flux_1,dV3s_1,u_dV3s_1,dV1_2_1, + u_dV1_2_1,Flag_1,Com_1,FileName_1 if(ar__1(30:33) .EQ. '') dV_1 = rNULL__ if(ar__1(35:38) .EQ. '') e_dV_1 = rNULL__ if(ar__1(42:45) .EQ. '') Tmb_1 = rNULL__ if(ar__1(48:50) .EQ. '') e_Tmb_1 = iNULL__ if(ar__1(54:59) .EQ. '') Flux_1 = rNULL__ if(ar__1(61:64) .EQ. '') e_Flux_1 = rNULL__ if(ar__1(69:73) .EQ. '') dV3s_1 = rNULL__ if(ar__1(77:80) .EQ. '') dV1_2_1 = rNULL__ c ..............Just test output........... write(6,2) + Species_1,Trans_1,dV_1,e_dV_1,u_dV_1,Tmb_1,u_Tmb_1,e_Tmb_1, + l_Flux_1,Flux_1,e_Flux_1,u_e_Flux_1,dV3s_1,u_dV3s_1,dV1_2_1, + u_dV1_2_1,Flag_1,Com_1,FileName_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Observed line emissions parameters for IRAS05358+3543 C Format for file interpretation 3 format( + A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3, + 1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,25 read(1,'(A115)')ar__2 read(ar__2,3) + Species_2,Trans_2,dV_2,e_dV_2,u_dV_2,Tmb_2,u_Tmb_2,e_Tmb_2, + l_Flux_2,Flux_2,e_Flux_2,u_e_Flux_2,dV3s_2,u_dV3s_2,dV1_2_2, + u_dV1_2_2,Flag_2,Com_2,FileName_2 if(ar__2(30:33) .EQ. '') dV_2 = rNULL__ if(ar__2(35:38) .EQ. '') e_dV_2 = rNULL__ if(ar__2(42:45) .EQ. '') Tmb_2 = rNULL__ if(ar__2(48:50) .EQ. '') e_Tmb_2 = iNULL__ if(ar__2(54:59) .EQ. '') Flux_2 = rNULL__ if(ar__2(61:64) .EQ. '') e_Flux_2 = rNULL__ if(ar__2(69:73) .EQ. '') dV3s_2 = rNULL__ if(ar__2(77:80) .EQ. '') dV1_2_2 = rNULL__ c ..............Just test output........... write(6,3) + Species_2,Trans_2,dV_2,e_dV_2,u_dV_2,Tmb_2,u_Tmb_2,e_Tmb_2, + l_Flux_2,Flux_2,e_Flux_2,u_e_Flux_2,dV3s_2,u_dV3s_2,dV1_2_2, + u_dV1_2_2,Flag_2,Com_2,FileName_2 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! Observed line emissions parameters for IRAS18162-2048 C Format for file interpretation 4 format( + A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3, + 1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,25 read(1,'(A115)')ar__3 read(ar__3,4) + Species_3,Trans_3,dV_3,e_dV_3,u_dV_3,Tmb_3,u_Tmb_3,e_Tmb_3, + l_Flux_3,Flux_3,e_Flux_3,u_e_Flux_3,dV3s_3,u_dV3s_3,dV1_2_3, + u_dV1_2_3,Flag_3,Com_3,FileName_3 if(ar__3(30:33) .EQ. '') dV_3 = rNULL__ if(ar__3(35:38) .EQ. '') e_dV_3 = rNULL__ if(ar__3(42:45) .EQ. '') Tmb_3 = rNULL__ if(ar__3(48:50) .EQ. '') e_Tmb_3 = iNULL__ if(ar__3(54:59) .EQ. '') Flux_3 = rNULL__ if(ar__3(61:64) .EQ. '') e_Flux_3 = rNULL__ if(ar__3(69:73) .EQ. '') dV3s_3 = rNULL__ if(ar__3(77:80) .EQ. '') dV1_2_3 = rNULL__ c ..............Just test output........... write(6,4) + Species_3,Trans_3,dV_3,e_dV_3,u_dV_3,Tmb_3,u_Tmb_3,e_Tmb_3, + l_Flux_3,Flux_3,e_Flux_3,u_e_Flux_3,dV3s_3,u_dV3s_3,dV1_2_3, + u_dV1_2_3,Flag_3,Com_3,FileName_3 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! Opacities derived from modeling the molecular * emission of the sources C Format for file interpretation 5 format(A9,1X,A17,1X,F6.2,A2,1X,F4.2,1X,F4.2,A2,1X,F4.2,A2) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,24 read(1,'(A55)')ar__4 read(ar__4,5) + Species_4,Trans_4,Opac1,n_Opac1,Opac2,Opac3,n_Opac3,Opac4, + n_Opac4 if(ar__4(29:34) .EQ. '') Opac1 = rNULL__ if(ar__4(43:46) .EQ. '') Opac3 = rNULL__ if(ar__4(50:53) .EQ. '') Opac4 = rNULL__ c ..............Just test output........... write(6,5) + Species_4,Trans_4,Opac1,n_Opac1,Opac2,Opac3,n_Opac3,Opac4, + n_Opac4 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end