Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/644/A64 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 2026-Aug-11
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/644/A64 CO observations of the IC443 SNR (Dell'Ova+, 2020)
*================================================================================
*Interstellar anatomy of the TeV gamma-ray peak in the IC443 supernova remnant.
* Dell'Ova P., Gusdorf A., Gerin M., Riquelme D., Gusten R.,
* Noriega-Crespo A., Tram L.N., Houde M., Guillard P., Lehmann A.,
* Lesaffre P., Louvet F., Marcowith A., Padovani M.
* <Astron. Astrophys. 644, A64 (2020)>
* =2020A&A...644A..64D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! List of datacubes
integer*4 nr__
parameter (nr__=9) ! Number of records
character*134 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg ! (deg) Right Ascension of center (J2000)
real*8 DEdeg ! (deg) Declination of center (J2000)
integer*4 Nx ! Number of pixels along X-axis
integer*4 Ny ! Number of pixels along Y-axis
integer*4 Nz ! Number of pixels along Z-axis
integer*4 bVRAD ! (m/s) Lower value of VRAD interval
integer*4 BVRAD_1 ! (m/s) Upper value of VRAD interval
real*8 dVRAD ! (m/s) VRAD resolution
integer*4 size ! (Kibyte) Size of FITS file
character*23 FileName ! Name of FITS file, in subdirectory fits
character*49 Title ! Title of the FITS file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableg1.dat' ! Infrared photometric magnitudes of YSO candidates
detected by WISE in the extended G region
integer*4 nr__1
parameter (nr__1=65) ! Number of records
character*90 ar__1 ! Full-size record
integer*4 ID ! [1/65] Sequential number
character*19 Name ! Catalog name (JHHMMSS.ss+DDMMSS.s)
real*4 W1mag ! (mag) Wise W1 (3.4um) magnitude
real*4 e_W1mag ! (mag) Wise W1 (3.4um) magnitude error
real*4 W2mag ! (mag) Wise W2 (4.6um) magnitude
real*4 e_W2mag ! (mag) Wise W2 (4.6um) magnitude error
real*4 W3mag ! (mag) Wise W3 (12um) magnitude
real*4 e_W3mag ! (mag) Wise W3 (12um) magnitude error
real*4 W4mag ! (mag) Wise W4 (22um) magnitude
real*4 e_W4mag ! (mag) ? Wise W4( 22um) magnitude error
character*8 Label ! Figure 13 label (Class I or Class II) (1)
real*4 alpha ! Spectral index alpha,
* {alpha}=dlog({lambda}F_{lambda}_)dlog{lambda}
integer*4 Flag ! [0/5] Flag (G1)
*Note (1): Label as follows:
* Class I = {alpha}>=3
* Class II = -1.6<={alpha}<=-0.3
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableg2.dat' ! Infrared photometric magnitudes of YSO candidates
detected by 2MASS in the extended G region
integer*4 nr__2
parameter (nr__2=79) ! Number of records
character*74 ar__2 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 ID_1 ! [1/79] Sequential number
integer*4 RAh ! (h) Right ascension (J2000)
integer*4 RAm ! (min) Right ascension (J2000)
real*4 RAs ! (s) Right ascension (J2000)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) Declination (J2000)
integer*4 DEm ! (arcmin) Declination (J2000)
real*4 DEs ! (arcsec) Declination (J2000)
real*4 Jmag ! (mag) 2MASS J (1.25um) magnitude
real*4 e_Jmag ! (mag) 2MASS J (1.25um) magnitude error
real*4 Hmag ! (mag) 2MASS H (1.65um) magnitude
real*4 e_Hmag ! (mag) 2MASS H (1.65um) magnitude error
real*4 Kmag ! (mag) 2MASS K (2.17um) magnitude
real*4 e_Kmag ! (mag) 2MASS K (2.17um) magnitude error
character*5 Label_1 ! Figure 13 label (YSO, CTTS or HAeBe)
integer*4 Flag_1 ! [0/2] Flag (G1)
C=============================================================================
C Loading file 'list.dat' ! List of datacubes
C Format for file interpretation
1 format(
+ F9.5,F9.5,1X,I3,1X,I3,1X,I4,1X,I7,1X,I6,1X,F7.3,1X,I5,1X,A23,
+ 1X,A49)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,9
read(1,'(A134)')ar__
read(ar__,1)
+ RAdeg,DEdeg,Nx,Ny,Nz,bVRAD,BVRAD_1,dVRAD,size,FileName,Title
c ..............Just test output...........
write(6,1)
+ RAdeg,DEdeg,Nx,Ny,Nz,bVRAD,BVRAD_1,dVRAD,size,FileName,Title
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableg1.dat' ! Infrared photometric magnitudes of YSO candidates
* detected by WISE in the extended G region
C Format for file interpretation
2 format(
+ I2,1X,A19,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,A8,1X,F5.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'tableg1.dat')
write(6,*) '....Loading file: tableg1.dat'
do i__=1,65
read(1,'(A90)')ar__1
read(ar__1,2)
+ ID,Name,W1mag,e_W1mag,W2mag,e_W2mag,W3mag,e_W3mag,W4mag,
+ e_W4mag,Label,alpha,Flag
if(ar__1(69:73) .EQ. '') e_W4mag = rNULL__
c ..............Just test output...........
write(6,2)
+ ID,Name,W1mag,e_W1mag,W2mag,e_W2mag,W3mag,e_W3mag,W4mag,
+ e_W4mag,Label,alpha,Flag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableg2.dat' ! Infrared photometric magnitudes of YSO candidates
* detected by 2MASS in the extended G region
C Format for file interpretation
3 format(
+ I2,1X,I1,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F6.3,1X,F5.3,
+ 1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,A5,2X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'tableg2.dat')
write(6,*) '....Loading file: tableg2.dat'
do i__=1,79
read(1,'(A74)')ar__2
read(ar__2,3)
+ ID_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,
+ e_Kmag,Label_1,Flag_1
RAdeg_1 = rNULL__
DEdeg_1 = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg_1=RAh*15.
if(RAm .GT. -180) RAdeg_1=RAdeg_1+RAm/4.
if(RAs .GT. -180) RAdeg_1=RAdeg_1+RAs/240.
if(DEd .GE. 0) DEdeg_1=DEd
if(DEm .GE. 0) DEdeg_1=DEdeg_1+DEm/60.
if(DEs .GE. 0) DEdeg_1=DEdeg_1+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,3)
+ ID_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,
+ e_Kmag,Label_1,Flag_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end