Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/504/853 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-28 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 (nr__) ! Species character*17 Trans (nr__) ! Observed line transition real*4 dV (nr__) ! (km/s) ? Linewidth real*4 e_dV (nr__) ! (km/s) ? rms uncertainty on dV character*1 u_dV (nr__) ! [?] Uncertainty flag on dV real*4 Tmb (nr__) ! (K) ? Main beam temperature character*1 u_Tmb (nr__) ! [?] Uncertainty flag on Tmb integer*4 e_Tmb (nr__) ! (mK) ? rms uncertainty on Tmb character*1 l_Flux (nr__) ! Limit flag on Flux real*4 Flux (nr__) ! (K.km/s) ? Integrated flux real*4 e_Flux (nr__) ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux (nr__) ! [?] Uncertainty flag on Flux real*4 dV3s (nr__) ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s (nr__) ! [?] Uncertainty flag on dV3s real*4 dV1_2 (nr__) ! (km/s) ? half-power linewidth character*1 u_dV1_2 (nr__) ! [?] Uncertainty flag on dV1/2 character*1 Flag (nr__) ! [abcdfg] Flag on profile (1) character*14 Com (nr__) ! Comments character*13 FileName (nr__) ! 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 (nr__1) ! Species character*17 Trans_1 (nr__1) ! Observed line transition real*4 dV_1 (nr__1) ! (km/s) ? Linewidth real*4 e_dV_1 (nr__1) ! (km/s) ? rms uncertainty on dV character*1 u_dV_1 (nr__1) ! [?] Uncertainty flag on dV real*4 Tmb_1 (nr__1) ! (K) ? Main beam temperature character*1 u_Tmb_1 (nr__1) ! [?] Uncertainty flag on Tmb integer*4 e_Tmb_1 (nr__1) ! (mK) ? rms uncertainty on Tmb character*1 l_Flux_1 (nr__1) ! Limit flag on Flux real*4 Flux_1 (nr__1) ! (K.km/s) ? Integrated flux real*4 e_Flux_1 (nr__1) ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux_1 (nr__1) ! [?] Uncertainty flag on Flux real*4 dV3s_1 (nr__1) ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s_1 (nr__1) ! [?] Uncertainty flag on dV3s real*4 dV1_2_1 (nr__1) ! (km/s) ? half-power linewidth character*1 u_dV1_2_1 (nr__1) ! [?] Uncertainty flag on dV1/2 character*1 Flag_1 (nr__1) ! [abcdfg] Flag on profile (1) character*14 Com_1 (nr__1) ! Comments character*13 FileName_1 (nr__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 (nr__2) ! Species character*17 Trans_2 (nr__2) ! Observed line transition real*4 dV_2 (nr__2) ! (km/s) ? Linewidth real*4 e_dV_2 (nr__2) ! (km/s) ? rms uncertainty on dV character*1 u_dV_2 (nr__2) ! [?] Uncertainty flag on dV real*4 Tmb_2 (nr__2) ! (K) ? Main beam temperature character*1 u_Tmb_2 (nr__2) ! [?] Uncertainty flag on Tmb integer*4 e_Tmb_2 (nr__2) ! (mK) ? rms uncertainty on Tmb character*1 l_Flux_2 (nr__2) ! Limit flag on Flux real*4 Flux_2 (nr__2) ! (K.km/s) ? Integrated flux real*4 e_Flux_2 (nr__2) ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux_2 (nr__2) ! [?] Uncertainty flag on Flux real*4 dV3s_2 (nr__2) ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s_2 (nr__2) ! [?] Uncertainty flag on dV3s real*4 dV1_2_2 (nr__2) ! (km/s) ? half-power linewidth character*1 u_dV1_2_2 (nr__2) ! [?] Uncertainty flag on dV1/2 character*1 Flag_2 (nr__2) ! [abcdfg] Flag on profile (1) character*14 Com_2 (nr__2) ! Comments character*13 FileName_2 (nr__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 (nr__3) ! Species character*17 Trans_3 (nr__3) ! Observed line transition real*4 dV_3 (nr__3) ! (km/s) ? Linewidth real*4 e_dV_3 (nr__3) ! (km/s) ? rms uncertainty on dV character*1 u_dV_3 (nr__3) ! [?] Uncertainty flag on dV real*4 Tmb_3 (nr__3) ! (K) ? Main beam temperature character*1 u_Tmb_3 (nr__3) ! [?] Uncertainty flag on Tmb integer*4 e_Tmb_3 (nr__3) ! (mK) ? rms uncertainty on Tmb character*1 l_Flux_3 (nr__3) ! Limit flag on Flux real*4 Flux_3 (nr__3) ! (K.km/s) ? Integrated flux real*4 e_Flux_3 (nr__3) ! (K.km/s) ? rms uncertainty on Flux character*1 u_e_Flux_3 (nr__3) ! [?] Uncertainty flag on Flux real*4 dV3s_3 (nr__3) ! (km/s) ? linewidth at 3 sigma character*1 u_dV3s_3 (nr__3) ! [?] Uncertainty flag on dV3s real*4 dV1_2_3 (nr__3) ! (km/s) ? half-power linewidth character*1 u_dV1_2_3 (nr__3) ! [?] Uncertainty flag on dV1/2 character*1 Flag_3 (nr__3) ! [abcdfg] Flag on profile (1) character*14 Com_3 (nr__3) ! Comments character*13 FileName_3 (nr__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 (nr__4) ! Species character*17 Trans_4 (nr__4) ! Observed line transition real*4 Opac1 (nr__4) ! ? Line opacity for source W43MM1 character*2 n_Opac1 (nr__4) ! [ND] ND for no data real*4 Opac2 (nr__4) ! Line opacity for source IRAS18264-1152 real*4 Opac3 (nr__4) ! ? Line opacity for source IRAS05358+3543 character*2 n_Opac3 (nr__4) ! [ND] ND for no data real*4 Opac4 (nr__4) ! ? Line opacity for source IRAS18162-2048 character*2 n_Opac4 (nr__4) ! [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(i__),Trans(i__),dV(i__),e_dV(i__),u_dV(i__),Tmb(i__), + u_Tmb(i__),e_Tmb(i__),l_Flux(i__),Flux(i__),e_Flux(i__), + u_e_Flux(i__),dV3s(i__),u_dV3s(i__),dV1_2(i__),u_dV1_2(i__), + Flag(i__),Com(i__),FileName(i__) if(ar__(30:33) .EQ. '') dV(i__) = rNULL__ if(ar__(35:38) .EQ. '') e_dV(i__) = rNULL__ if(ar__(42:45) .EQ. '') Tmb(i__) = rNULL__ if(ar__(48:50) .EQ. '') e_Tmb(i__) = iNULL__ if(ar__(54:59) .EQ. '') Flux(i__) = rNULL__ if(ar__(61:64) .EQ. '') e_Flux(i__) = rNULL__ if(ar__(69:73) .EQ. '') dV3s(i__) = rNULL__ if(ar__(77:80) .EQ. '') dV1_2(i__) = rNULL__ c ..............Just test output........... write(6,1) + Species(i__),Trans(i__),dV(i__),e_dV(i__),u_dV(i__),Tmb(i__), + u_Tmb(i__),e_Tmb(i__),l_Flux(i__),Flux(i__),e_Flux(i__), + u_e_Flux(i__),dV3s(i__),u_dV3s(i__),dV1_2(i__),u_dV1_2(i__), + Flag(i__),Com(i__),FileName(i__) 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(i__),Trans_1(i__),dV_1(i__),e_dV_1(i__),u_dV_1(i__), + Tmb_1(i__),u_Tmb_1(i__),e_Tmb_1(i__),l_Flux_1(i__), + Flux_1(i__),e_Flux_1(i__),u_e_Flux_1(i__),dV3s_1(i__), + u_dV3s_1(i__),dV1_2_1(i__),u_dV1_2_1(i__),Flag_1(i__), + Com_1(i__),FileName_1(i__) if(ar__1(30:33) .EQ. '') dV_1(i__) = rNULL__ if(ar__1(35:38) .EQ. '') e_dV_1(i__) = rNULL__ if(ar__1(42:45) .EQ. '') Tmb_1(i__) = rNULL__ if(ar__1(48:50) .EQ. '') e_Tmb_1(i__) = iNULL__ if(ar__1(54:59) .EQ. '') Flux_1(i__) = rNULL__ if(ar__1(61:64) .EQ. '') e_Flux_1(i__) = rNULL__ if(ar__1(69:73) .EQ. '') dV3s_1(i__) = rNULL__ if(ar__1(77:80) .EQ. '') dV1_2_1(i__) = rNULL__ c ..............Just test output........... write(6,2) + Species_1(i__),Trans_1(i__),dV_1(i__),e_dV_1(i__),u_dV_1(i__), + Tmb_1(i__),u_Tmb_1(i__),e_Tmb_1(i__),l_Flux_1(i__), + Flux_1(i__),e_Flux_1(i__),u_e_Flux_1(i__),dV3s_1(i__), + u_dV3s_1(i__),dV1_2_1(i__),u_dV1_2_1(i__),Flag_1(i__), + Com_1(i__),FileName_1(i__) 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(i__),Trans_2(i__),dV_2(i__),e_dV_2(i__),u_dV_2(i__), + Tmb_2(i__),u_Tmb_2(i__),e_Tmb_2(i__),l_Flux_2(i__), + Flux_2(i__),e_Flux_2(i__),u_e_Flux_2(i__),dV3s_2(i__), + u_dV3s_2(i__),dV1_2_2(i__),u_dV1_2_2(i__),Flag_2(i__), + Com_2(i__),FileName_2(i__) if(ar__2(30:33) .EQ. '') dV_2(i__) = rNULL__ if(ar__2(35:38) .EQ. '') e_dV_2(i__) = rNULL__ if(ar__2(42:45) .EQ. '') Tmb_2(i__) = rNULL__ if(ar__2(48:50) .EQ. '') e_Tmb_2(i__) = iNULL__ if(ar__2(54:59) .EQ. '') Flux_2(i__) = rNULL__ if(ar__2(61:64) .EQ. '') e_Flux_2(i__) = rNULL__ if(ar__2(69:73) .EQ. '') dV3s_2(i__) = rNULL__ if(ar__2(77:80) .EQ. '') dV1_2_2(i__) = rNULL__ c ..............Just test output........... write(6,3) + Species_2(i__),Trans_2(i__),dV_2(i__),e_dV_2(i__),u_dV_2(i__), + Tmb_2(i__),u_Tmb_2(i__),e_Tmb_2(i__),l_Flux_2(i__), + Flux_2(i__),e_Flux_2(i__),u_e_Flux_2(i__),dV3s_2(i__), + u_dV3s_2(i__),dV1_2_2(i__),u_dV1_2_2(i__),Flag_2(i__), + Com_2(i__),FileName_2(i__) 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(i__),Trans_3(i__),dV_3(i__),e_dV_3(i__),u_dV_3(i__), + Tmb_3(i__),u_Tmb_3(i__),e_Tmb_3(i__),l_Flux_3(i__), + Flux_3(i__),e_Flux_3(i__),u_e_Flux_3(i__),dV3s_3(i__), + u_dV3s_3(i__),dV1_2_3(i__),u_dV1_2_3(i__),Flag_3(i__), + Com_3(i__),FileName_3(i__) if(ar__3(30:33) .EQ. '') dV_3(i__) = rNULL__ if(ar__3(35:38) .EQ. '') e_dV_3(i__) = rNULL__ if(ar__3(42:45) .EQ. '') Tmb_3(i__) = rNULL__ if(ar__3(48:50) .EQ. '') e_Tmb_3(i__) = iNULL__ if(ar__3(54:59) .EQ. '') Flux_3(i__) = rNULL__ if(ar__3(61:64) .EQ. '') e_Flux_3(i__) = rNULL__ if(ar__3(69:73) .EQ. '') dV3s_3(i__) = rNULL__ if(ar__3(77:80) .EQ. '') dV1_2_3(i__) = rNULL__ c ..............Just test output........... write(6,4) + Species_3(i__),Trans_3(i__),dV_3(i__),e_dV_3(i__),u_dV_3(i__), + Tmb_3(i__),u_Tmb_3(i__),e_Tmb_3(i__),l_Flux_3(i__), + Flux_3(i__),e_Flux_3(i__),u_e_Flux_3(i__),dV3s_3(i__), + u_dV3s_3(i__),dV1_2_3(i__),u_dV1_2_3(i__),Flag_3(i__), + Com_3(i__),FileName_3(i__) 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(i__),Trans_4(i__),Opac1(i__),n_Opac1(i__), + Opac2(i__),Opac3(i__),n_Opac3(i__),Opac4(i__),n_Opac4(i__) if(ar__4(29:34) .EQ. '') Opac1(i__) = rNULL__ if(ar__4(43:46) .EQ. '') Opac3(i__) = rNULL__ if(ar__4(50:53) .EQ. '') Opac4(i__) = rNULL__ c ..............Just test output........... write(6,5) + Species_4(i__),Trans_4(i__),Opac1(i__),n_Opac1(i__), + Opac2(i__),Opac3(i__),n_Opac3(i__),Opac4(i__),n_Opac4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end