Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/582/A68 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-19
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/582/A68 Database of circumstellar OH masers (Engels+, 2015)
*================================================================================
*A database of circumstellar OH masers.
* Engels D., Bunzel F.
* <Astron. Astrophys. 582, A68 (2015)>
* =2015A&A...582A..68E (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'oh_maser.dat' ! OH maser observations
integer*4 nr__
parameter (nr__=13655) ! Number of records
character*195 ar__ ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 Obs ! [1/3]? Prioritized observations (1)
integer*4 IDobs ! Observation identification number (1)
integer*4 Pap ! Reference identification number (1)
integer*4 ID ! Object identification number (1)
character*16 Name ! Object name
integer*4 Freq ! (MHz) Frequency
character*1 SpType ! [DSI] Spectral type of maser profile (2)
integer*4 RAh ! (h) Right ascension (J2000.0)
integer*4 RAm ! (min) Right ascension (J2000.0)
real*4 RAs ! (s) Right ascension (J2000.0)
character*1 DE_ ! Declination sign (J2000.0)
integer*4 DEd ! (deg) Declination (J2000.0)
integer*4 DEm ! (arcmin) Declination (J2000.0)
real*4 DEs ! (arcsec) Declination (J2000.0)
real*4 vel_blue ! (km/s) ? Velocity of blue-shifted maser peak
real*4 vel_red ! (km/s) ? Velocity of red-shifted maser peak
real*4 e_vel ! (km/s) ? Velocity resolution
real*4 vel_cent ! (km/s) ? Radial velocity
real*4 vel_shell ! (km/s) ? Expansion velocity of the maser shell
character*1 Flag ! [><=:.&] Flag for detected maser emission (3)
real*8 Sblue ! (Jy) ? Flux density of blue-shifted maser peak (4)
real*8 Sred ! (Jy) ? Flux density of red-shifted maser peak
real*8 Iblue ! (Jy) ? Integral flux of blue-shifted maser peak (4)
real*8 Ired ! (Jy) ? Integral flux of red-shifted maser peak (4)
character*6 Ref ! Acronym of the reference (1)
real*8 RAdeg ! (deg) Right ascension (J2000.0)
real*8 DEdeg ! (deg) Declination (J2000.0)
integer*4 Monitor ! [0/1] Presence in monitor.dat file (5)
integer*4 Maps ! [0/1] Presence in maps.dat file (6)
character*19 BibCode ! Bibliographic Code
*Note (1): Numbers/designations used for administration of the database
*Note (2): Spectral type of maser profile as follows:
* D = Double peaked maser profile
* S = Single peaked maser profile
* I = Irregular maser profile
*Note (3): Marker for detected maser emission ad follows:
* < = non-detections
* > = lower flux density limit
* = = no velocities available
* . = blue flux density peak questionable
* : = red flux density peak questionable
* & = both flux density peaks questionable
*Note (4): 3-sigma upper limit in case of non-detections
*Note (5): 1 if monitor.dat contains entry for this source and frequency.
*Note (6): 1 if maps.dat contains entry for this source and frequency.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'maps.dat' ! Interferometric OH maser observations (oh_maser_maps)
integer*4 nr__1
parameter (nr__1=271) ! Number of records
character*162 ar__1 ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 IDmeas ! Observation identification number (1)
integer*4 Pap_1 ! Reference identification number (1)
integer*4 ID_1 ! Object identification number (1)
character*15 Name_1 ! Object name
integer*4 RAh_1 ! (h) Right ascension (J2000.0)
integer*4 RAm_1 ! (min) Right ascension (J2000.0)
real*4 RAs_1 ! (s) Right ascension (J2000.0)
character*1 DE__1 ! Declination sign (J2000.0)
integer*4 DEd_1 ! (deg) Declination (J2000.0)
integer*4 DEm_1 ! (arcmin) Declination (J2000.0)
real*4 DEs_1 ! (arcsec) Declination (J2000.0)
integer*4 Freq_1 ! (MHz) Frequency
character*6 Inst ! Name of radio interferometer
real*4 Res ! (mas) ? Spatial resolution
character*1 Flag_1 ! [<] Limit flag for detected maser emission (2)
real*4 rms ! (mJy) ? Sensitivity of the observation
character*1 Pol ! [1/0] Marker for polarization observations (3)
integer*4 Nmaps ! Number of interferometric maps
character*10 ObsdateS ! ("date") Observing date of the first map (4)
character*10 ObsdateE ! ("date") Observing date of the last map (4)
character*6 Ref_1 ! Acronym of the reference (1)
real*8 RAdeg_1 ! (deg) Right ascension (J2000.0)
real*8 DEdeg_1 ! (deg) Declination (J2000.0)
character*19 BibCode_1 ! Bibliographic Code
*Note (1): Numbers/designations used for administration of the database
*Note (2): '<' for non-detections
*Note (3): Marker for polarization observations as follows:
* 1 = polarization measured
* 0 = polarization not measured
*Note (4): Format: YYYY-MM-DD, missing information is replaced by blank.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'monitor.dat' ! OH maser monitoring programs (oh_maser_monitor)
integer*4 nr__2
parameter (nr__2=180) ! Number of records
character*151 ar__2 ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 IDmeas_1 ! Observation identification number(1)
integer*4 Pap_2 ! Reference identification number (1)
integer*4 ID_2 ! Object identification number (1)
character*12 Name_2 ! Object name
integer*4 RAh_2 ! (h) Right ascension (J2000.0)
integer*4 RAm_2 ! (min) Right ascension (J2000.0)
real*4 RAs_2 ! (s) Right ascension (J2000.0)
character*1 DE__2 ! Declination sign (J2000.0)
integer*4 DEd_2 ! (deg) Declination (J2000.0)
integer*4 DEm_2 ! (arcmin) Declination (J2000.0)
real*4 DEs_2 ! (arcsec) Declination (J2000.0)
integer*4 Freq_2 ! (MHz) Frequency
character*10 ObsdateS_1 ! ("date") Start of the monitoring program (2)
character*10 ObsdateE_1 ! ("date") End of the monitoring program (2)
character*1 l_Smaxpk ! [<] Marker for detected maser emission (3)
real*8 Smaxpk ! (Jy) ? Maximum flux density detected
character*1 l_Sminpk ! Marker for detected maser emission (3)
character*8 Sminpk ! (Jy) Minimum flux density detected
character*6 Ref_2 ! Acronym of the reference (1)
real*8 RAdeg_2 ! (deg) Right ascension (J2000.0)
real*8 DEdeg_2 ! (deg) Declination (J2000.0)
character*19 BibCode_2 ! Bibliographic Code
*Note (1): Numbers/designations used for administration of the database
*Note (2): Format: YYYY-MM-DD, missing information is replaced by '00'
*Note (3): '<' for non-detections
C=============================================================================
C Loading file 'oh_maser.dat' ! OH maser observations
C Format for file interpretation
1 format(
+ I1,1X,I5,1X,I3,1X,I5,1X,A16,1X,I4,1X,A1,1X,I2,1X,I2,1X,F6.3,
+ 1X,A1,I2,1X,I2,1X,F5.2,2X,F6.1,2X,F6.1,1X,F4.2,2X,F6.1,3X,
+ F5.1,1X,A1,1X,F8.3,1X,F7.3,1X,F8.3,1X,F7.3,1X,A6,1X,F11.7,1X,
+ F11.7,1X,I1,1X,I1,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'oh_maser.dat')
write(6,*) '....Loading file: oh_maser.dat'
do i__=1,13655
read(1,'(A195)')ar__
read(ar__,1)
+ Obs,IDobs,Pap,ID,Name,Freq,SpType,RAh,RAm,RAs,DE_,DEd,DEm,DEs,
+ vel_blue,vel_red,e_vel,vel_cent,vel_shell,Flag,Sblue,Sred,
+ Iblue,Ired,Ref,RAdeg,DEdeg,Monitor,Maps,BibCode
if(ar__(1:1) .EQ. '') Obs = iNULL__
if(ar__(70:75) .EQ. '') vel_blue = rNULL__
if(ar__(78:83) .EQ. '') vel_red = rNULL__
if(ar__(85:88) .EQ. '') e_vel = rNULL__
if(ar__(91:96) .EQ. '') vel_cent = rNULL__
if(ar__(100:104) .EQ. '') vel_shell = rNULL__
if(ar__(108:115) .EQ. '') Sblue = rNULL__
if(ar__(117:123) .EQ. '') Sred = rNULL__
if(ar__(125:132) .EQ. '') Iblue = rNULL__
if(ar__(134:140) .EQ. '') Ired = rNULL__
c ..............Just test output...........
write(6,1)
+ Obs,IDobs,Pap,ID,Name,Freq,SpType,RAh,RAm,RAs,DE_,DEd,DEm,DEs,
+ vel_blue,vel_red,e_vel,vel_cent,vel_shell,Flag,Sblue,Sred,
+ Iblue,Ired,Ref,RAdeg,DEdeg,Monitor,Maps,BibCode
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'maps.dat' ! Interferometric OH maser observations (oh_maser_maps)
C Format for file interpretation
2 format(
+ I5,1X,I3,1X,I5,1X,A15,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,
+ F5.2,1X,I4,1X,A6,1X,F6.1,1X,A1,F5.3,1X,A1,1X,I3,1X,A10,1X,A10,
+ 1X,A6,1X,F11.7,1X,F11.7,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'maps.dat')
write(6,*) '....Loading file: maps.dat'
do i__=1,271
read(1,'(A162)')ar__1
read(ar__1,2)
+ IDmeas,Pap_1,ID_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,
+ DEs_1,Freq_1,Inst,Res,Flag_1,rms,Pol,Nmaps,ObsdateS,ObsdateE,
+ Ref_1,RAdeg_1,DEdeg_1,BibCode_1
if(ar__1(71:76) .EQ. '') Res = rNULL__
if(ar__1(79:83) .EQ. '') rms = rNULL__
c ..............Just test output...........
write(6,2)
+ IDmeas,Pap_1,ID_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,
+ DEs_1,Freq_1,Inst,Res,Flag_1,rms,Pol,Nmaps,ObsdateS,ObsdateE,
+ Ref_1,RAdeg_1,DEdeg_1,BibCode_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'monitor.dat' ! OH maser monitoring programs (oh_maser_monitor)
C Format for file interpretation
3 format(
+ I5,1X,I3,1X,I5,1X,A12,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,
+ F5.2,1X,I4,1X,A10,1X,A10,1X,A1,F7.3,1X,A1,A8,1X,A6,1X,F11.7,
+ 1X,F11.7,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'monitor.dat')
write(6,*) '....Loading file: monitor.dat'
do i__=1,180
read(1,'(A151)')ar__2
read(ar__2,3)
+ IDmeas_1,Pap_2,ID_2,Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,
+ DEm_2,DEs_2,Freq_2,ObsdateS_1,ObsdateE_1,l_Smaxpk,Smaxpk,
+ l_Sminpk,Sminpk,Ref_2,RAdeg_2,DEdeg_2,BibCode_2
if(ar__2(84:90) .EQ. '') Smaxpk = rNULL__
c ..............Just test output...........
write(6,3)
+ IDmeas_1,Pap_2,ID_2,Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,
+ DEm_2,DEs_2,Freq_2,ObsdateS_1,ObsdateE_1,l_Smaxpk,Smaxpk,
+ l_Sminpk,Sminpk,Ref_2,RAdeg_2,DEdeg_2,BibCode_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end