Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/814/133 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-Apr-16
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/814/133 Herschel Galactic plane survey of [NII] (Goldsmith+, 2015)
*================================================================================
*Herschel Galactic plane survey of [NII] fine structure emission.
* Goldsmith P.F., Yildiz U.A., Langer W.D., Pineda J.L.
* <Astrophys. J., 814, 133 (2015)>
* =2015ApJ...814..133G (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Observed PACS Integrated Intensities for
[N II] 122 and 205 {mu}m Transitions are
Tabulated Together with the Inferred Electron
Volume Densities and N^+^ Column Densities
integer*4 nr__
parameter (nr__=70) ! Number of records
character*91 ar__ ! Full-size record
C Position composed of: GLON GLAT
character*10 ID (nr__) ! Line of sight (LOS) label (GLLL.l+B.b)
real*8 GLON (nr__) ! (deg) Galactic longitude
real*4 GLAT (nr__) ! (deg) Galactic latitude
integer*8 OBSID (nr__) ! Herschel OBSID
integer*4 R (nr__) ! [0/3]? Detection status (1)
character*1 l_NII122 (nr__) ! [<] Limit flag on NII122
real*4 NII122 (nr__) ! (10-8W/m2/sr) ? [N II] 122 {mu}m intensity averaged over
* the 25 PACS spaxels
real*4 e_NII122 (nr__) ! (10-8W/m2/sr) ? Uncertainty in NII122
real*4 RMS (nr__) ! (10-8) ? RMS variation among the 25 spaxels at
* 122 {mu}m wavelength
integer*4 R1 (nr__) ! [0/3]? Detection status (1)
character*1 l_NII205 (nr__) ! [<] Limit flag on NII205
real*4 NII205 (nr__) ! (10-8W/m2/sr) ? [N II] 205 {mu}m intensity averaged over
* the 25 PACS spaxels
real*4 e_NII205 (nr__) ! (10-8W/m2/sr) ? Uncertainty in NII205
real*4 n_e (nr__) ! (cm-3) ? Electron density
real*4 e_n_e (nr__) ! (cm-3) ? Uncertainty in n(e)
real*4 N_N_ (nr__) ! (10+16cm-2) ? N^+^ column density
real*4 e_N_N_ (nr__) ! (10+16cm-2) ? Uncertainty in N(N+)
*Note (1): Detection status as follows:
* 0 = Not observed;
* 1 = Non-detection;
* 2 = Detection in individual spaxel or spaxels;
* 3 = Only the average of all spaxels provides a detection.
* For the non-detections, we give the 3{sigma} upper limits on the intensities.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! PACS Results
integer*4 nr__1
parameter (nr__1=79) ! Number of records
character*91 ar__1 ! Full-size record
C Position composed of: GLON GLAT
character*10 ID_1 (nr__1) ! Line of sight (LOS) label (GLLL.l+B.b)
real*8 GLON_1 (nr__1) ! (deg) Galactic longitude
real*4 GLAT_1 (nr__1) ! (deg) Galactic latitude
integer*8 OBSID_1 (nr__1) ! Herschel OBSID
integer*4 R_1 (nr__1) ! [0/3]? Detection status (1)
character*1 l_NII122_1 (nr__1) ! [<] Limit flag on NII122
real*4 NII122_1 (nr__1) ! (10-8W/m2/sr) ? [N II] 122 {mu}m intensity averaged over
* the 25 PACS spaxels
real*4 e_NII122_1 (nr__1) ! (10-8W/m2/sr) ? Uncertainty in NII122
real*4 RMS_1 (nr__1) ! (10-8) ? RMS variation among the 25 spaxels at
* 122 {mu}m wavelength
integer*4 R1_1 (nr__1) ! [0/3]? Detection status (1)
character*1 l_NII205_1 (nr__1) ! [<] Limit flag on NII205
real*4 NII205_1 (nr__1) ! (10-8W/m2/sr) ? [N II] 205 {mu}m intensity averaged over
* the 25 PACS spaxels
real*4 e_NII205_1 (nr__1) ! (10-8W/m2/sr) ? Uncertainty in NII205
real*4 n_e_1 (nr__1) ! (cm-3) ? Electron density
real*4 e_n_e_1 (nr__1) ! (cm-3) ? Uncertainty in n(e)
real*4 N_N__1 (nr__1) ! (10+16cm-2) ? N^+^ column density
real*4 e_N_N__1 (nr__1) ! (10+16cm-2) ? Uncertainty in N(N+)
*Note (1): Detection status as follows:
* 0 = Not observed;
* 1 = Non-detection;
* 2 = Detection in individual spaxel or spaxels;
* 3 = Only the average of all spaxels provides a detection.
* For the non-detections, we give the 3{sigma} upper limits on the intensities.
C=============================================================================
C Loading file 'table2.dat' ! Observed PACS Integrated Intensities for
* [N II] 122 and 205 {mu}m Transitions are
* Tabulated Together with the Inferred Electron
* Volume Densities and N^+^ Column Densities
C Format for file interpretation
1 format(
+ A10,1X,F8.4,1X,F3.1,1X,I10,1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.2,
+ 1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,70
read(1,'(A91)')ar__
read(ar__,1)
+ ID(i__),GLON(i__),GLAT(i__),OBSID(i__),R(i__),l_NII122(i__),
+ NII122(i__),e_NII122(i__),RMS(i__),R1(i__),l_NII205(i__),
+ NII205(i__),e_NII205(i__),n_e(i__),e_n_e(i__),N_N_(i__),
+ e_N_N_(i__)
if(ar__(36:36) .EQ. '') R(i__) = iNULL__
if(ar__(39:43) .EQ. '') NII122(i__) = rNULL__
if(ar__(45:49) .EQ. '') e_NII122(i__) = rNULL__
if(ar__(51:55) .EQ. '') RMS(i__) = rNULL__
if(ar__(57:57) .EQ. '') R1(i__) = iNULL__
if(ar__(60:64) .EQ. '') NII205(i__) = rNULL__
if(ar__(66:70) .EQ. '') e_NII205(i__) = rNULL__
if(ar__(72:76) .EQ. '') n_e(i__) = rNULL__
if(ar__(78:81) .EQ. '') e_n_e(i__) = rNULL__
if(ar__(83:86) .EQ. '') N_N_(i__) = rNULL__
if(ar__(88:91) .EQ. '') e_N_N_(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),GLON(i__),GLAT(i__),OBSID(i__),R(i__),l_NII122(i__),
+ NII122(i__),e_NII122(i__),RMS(i__),R1(i__),l_NII205(i__),
+ NII205(i__),e_NII205(i__),n_e(i__),e_n_e(i__),N_N_(i__),
+ e_N_N_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! PACS Results
C Format for file interpretation
2 format(
+ A10,1X,F8.4,1X,F3.1,1X,I10,1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.2,
+ 1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,79
read(1,'(A91)')ar__1
read(ar__1,2)
+ ID_1(i__),GLON_1(i__),GLAT_1(i__),OBSID_1(i__),R_1(i__),
+ l_NII122_1(i__),NII122_1(i__),e_NII122_1(i__),RMS_1(i__),
+ R1_1(i__),l_NII205_1(i__),NII205_1(i__),e_NII205_1(i__),
+ n_e_1(i__),e_n_e_1(i__),N_N__1(i__),e_N_N__1(i__)
if(ar__1(36:36) .EQ. '') R_1(i__) = iNULL__
if(ar__1(39:43) .EQ. '') NII122_1(i__) = rNULL__
if(ar__1(45:49) .EQ. '') e_NII122_1(i__) = rNULL__
if(ar__1(51:55) .EQ. '') RMS_1(i__) = rNULL__
if(ar__1(57:57) .EQ. '') R1_1(i__) = iNULL__
if(ar__1(60:64) .EQ. '') NII205_1(i__) = rNULL__
if(ar__1(66:70) .EQ. '') e_NII205_1(i__) = rNULL__
if(ar__1(72:76) .EQ. '') n_e_1(i__) = rNULL__
if(ar__1(78:81) .EQ. '') e_n_e_1(i__) = rNULL__
if(ar__1(83:86) .EQ. '') N_N__1(i__) = rNULL__
if(ar__1(88:91) .EQ. '') e_N_N__1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1(i__),GLON_1(i__),GLAT_1(i__),OBSID_1(i__),R_1(i__),
+ l_NII122_1(i__),NII122_1(i__),e_NII122_1(i__),RMS_1(i__),
+ R1_1(i__),l_NII205_1(i__),NII205_1(i__),e_NII205_1(i__),
+ n_e_1(i__),e_n_e_1(i__),N_N__1(i__),e_N_N__1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end