Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/839/42 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/ApJ/839/42 Spectrosc. of the circumstellar gas of WD1145+017 (Redfield+, 2017)
*================================================================================
*Spectroscopic evolution of disintegrating planetesimals: minute to month
*variability in the circumstellar gas associated with WD 1145+017.
* Redfield S., Farihi J., Cauley P.W., Parsons S.G., Gansicke B.T.,
* Duvvuri G. M
* <Astrophys. J., 839, 42 (2017)>
* =2017ApJ...839...42R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Detected circumstellar absorption features and
fit parameters
integer*4 nr__
parameter (nr__=256) ! Number of records
character*408 ar__ ! Full-size record
character*5 Ion ! Ion with detected circumstellar absorption
real*8 lambda ! (0.1nm) [3056.8/8662.2] Rest wavelength of ion (RWave)
real*4 Elow ! (eV) [0/10.5] Energy of lower level
real*4 Eup ! (eV) [0.4/10.8] Energy of upper level
real*4 loggf ! [-4/0.6] log gf of transition
real*4 DepthK1 ! [0.01/0.3]? Normalized depth of circumstellar
* feature, 2015 Apr Keck data (DepthKeckApr15)
real*4 E_DepthK1 ! [0.002/0.04]? Plus Error on DepthK1
* (E_DepthKeckApr15)
real*4 e_DepthK1_1 ! ? Minus Error on DepthK1 (e_DepthKeckApr15_2)
real*4 VminK1 ! (km/s) [-395/133]? Minimum radial velocity of
* circumstellar feature, 2015 Apr Keck data
* (VminKApr15)
real*4 E_VminK1 ! (km/s) ? Plus Error on VminK1 (E_VminKApr15)
real*4 e_VminK1_1 ! (km/s) [1.4/145]? Minus Error on VminK1
* (e_VminKApr15_2)
real*4 VmaxK1 ! (km/s) [-109/521]? Maximum radial velocity of
* circumstellar feature, 2015 Apr Keck data
* (VmaxKeckApr15)
real*4 E_VmaxK1 ! (km/s) [1/140]? Plus Error on VmaxK1 (E_VmaxKeckApr15)
real*4 e_VmaxK1_1 ! (km/s) ? Minus Error on VmaxK1 (e_VmaxKeckApr15_2)
character*1 l_EWK1 ! Limit flag on EWK1 (l_EWKeckApr15)
real*4 EWK1 ! (0.1nm) [0.01/3]? Equivalent width of circumstellar
* feature, 2015 Apr Keck data (EWKeckApr15)
real*4 E_EWK1 ! (0.1nm) ? Plus Error on EWK1 (E_EWKeckApr15)
real*4 e_EWK1_1 ! (0.1nm) [0.004/0.3]? Minus Error on EWK1
* (e_EWKeckApr15_2)
real*4 DepthK2 ! [0.02/3.1]? Normalized depth of circumstellar
* feature, 2015 Nov Keck data (DepthKeckNov15)
real*4 E_DepthK2 ! [0.005/0.8]? Plus Error on DepthK2
real*4 e_DepthK2_1 ! ? Minus Error on DepthK2
real*4 VminK2 ! (km/s) [-433/153]? Minimum radial velocity of
* circumstellar feature, 2015 Nov Keck data
* (VminKeckNov15)
real*4 E_VminK2 ! (km/s) ? Plus Error on VminK2 (E_VminKeckNov15)
real*4 e_VminK2_1 ! (km/s) [1/146]? Minus Error on VminK2
* (e_VminKeckNov15_2)
real*4 VmaxK2 ! (km/s) [-93/528]? Maximum radial velocity of
* circumstellar feature, 2015 Nov Keck data
* (VmaxKeckNov15)
real*4 E_VmaxK2 ! (km/s) ? Plus Error on VmaxK2 (E_VmaxKeckNov15)
real*4 e_VmaxK2_1 ! (km/s) [0/128]? Minus Error on VmaxK2
* (e_VmaxKeckNov15_2)
character*1 l_EWK2 ! Limit flag on EWK2 (l_EWKeckNov15)
real*4 EWK2 ! (0.1nm) [0.01/2]? Equivalent width of circumstellar
* feature, 2015 Nov Keck data (EWKeckNov15)
real*4 E_EWK2 ! (0.1nm) ? Plus Error on EWK2 (E_EWKeckNov15)
real*4 e_EWK2_1 ! (0.1nm) [0.007/0.2]? Minus Error on on EWK2
* (e_EWKeckNov15_2)
real*4 DepthV1 ! [0.01/0.7]? Normalized depth of circumstellar
* feature, 2016 Feb VLT data (DepthVLTFeb16)
real*4 E_DepthV1 ! [0.002/0.7]? Plus Error on DepthV1
* (E_DepthVLTFeb16)
real*4 e_DepthV1_1 ! ? Minus Error on DepthV1 (e_DepthVLTFeb16_2)
real*4 VminV1 ! (km/s) [-963/620]? Minimum radial velocity of
* circumstellar feature, 2016 Feb VLT data
* (VminVLTFeb16)
real*4 E_VminV1 ! (km/s) ? Plus Error on VminV1 (E_VminVLTFeb16)
real*4 e_VminV1_1 ! (km/s) [0.9/523]? Minus Error on VminV1
* (e_VminVLTFeb16_2)
real*4 VmaxV1 ! (km/s) [-299/885]? Maximum radial velocity of
* circumstellar feature, 2016 Feb VLT data
* (VmaxVLTFeb16)
real*4 E_VmaxV1 ! (km/s) [0.5/174]? Plus Error on VmaxV1 (E_VmaxVLTFeb16)
real*4 e_VmaxV1_1 ! (km/s) ? Minus Error on VmaxV1 (e_VmaxVLTFeb16_2)
character*1 l_EWV1 ! Limit flag on EWV1
real*4 EWV1 ! (0.1nm) [0.02/2.3] Equivalent width of circumstellar
* feature, 2016 Feb VLT data (EWVLTFeb16)
real*4 E_EWV1 ! (0.1nm) [0.007/0.4]? Plus Error on EWV1 (E_EWVLTFeb16)
real*4 e_EWV1_1 ! (0.1nm) ? Minus Error on EWV1 (e_EWVLTFeb16_2)
real*4 DepthV2 ! [0.009/3]? Normalized depth of circumstellar
* feature, 2016 Mar VLT data (DepthVLTMar16)
real*4 E_DepthV2 ! [0.002/1.4]? Plus Error on DepthV2
* (E_DepthVLTMar16)
real*4 e_DepthV2_1 ! ? Minus Error on DepthV2 (e_DepthVLTMar16_2)
real*4 VminV2 ! (km/s) [-988/269]? Minimum radial velocity of
* circumstellar feature, 2016 Mar VLT data
* (VminVLTMar16)
real*4 E_VminV2 ! (km/s) ? Plus Error on VminV2 (E_VminVLTMar16)
real*4 e_VminV2_1 ! (km/s) [1/237]? Minus Error on VminV2
* (e_VminVLTMar16_2)
real*4 VmaxV2 ! (km/s) [-185/593]? Maximum radial velocity of
* circumstellar feature, 2016 Mar VLT data
* (VmaxVLTMar16)
real*4 E_VmaxV2 ! (km/s) [0.4/113]? Plus Error on VmaxV2 (E_VmaxVLTMar16)
real*4 e_VmaxV2_1 ! (km/s) ? Minus Error on VmaxV2 (e_VmaxVLTMar16_2)
character*1 l_EWV2 ! Limit flag on EWV2
real*4 EWV2 ! (0.1nm) [0.01/3.4] Equivalent width of circumstellar
* feature, 2016 Mar VLT data (EWVLTMar16)
real*4 E_EWV2 ! (0.1nm) [0.006/0.4]? Plus Error on EWV2 (E_EWVLTMar16)
real*4 e_EWV2_1 ! (0.1nm) ? Minus Error on EWV2 (e_EWVLTMar16_2)
real*4 DepthV3 ! [0.01/0.5]? Normalized depth of circumstellar
* feature, 2016 Apr VLT data (DepthVLTApr16)
real*4 E_DepthV3 ! [0.001/0.2]? Plus Error on DepthV3
real*4 e_DepthV3_1 ! ? Minus Error on DepthV3
real*4 VminV3 ! (km/s) [-974/240]? Minimum radial velocity of
* circumstellar feature, 2016 Apr VLT data
* (VminVLTApr16)
real*4 E_VminV3 ! (km/s) ? Plus Error on VminV3 (E_VminVLTApr16)
real*4 e_VminV3_1 ! (km/s) [0.2/127]? Minus Error on VminV3
* (e_VminVLTApr16_2)
real*4 VmaxV3 ! (km/s) [-73/580]? Maximum radial velocity of
* circumstellar feature, 2016 Apr VLT data
* (VmaxVLTApr16)
real*4 E_VmaxV3 ! (km/s) [0.2/85]? Plus Error on VmaxV3 (E_VmaxVLTApr16)
real*4 e_VmaxV3_1 ! (km/s) ? Minus Error on VmaxV3 (e_VmaxVLTApr16_2)
character*1 l_EWV3 ! Limit flag on EWV3 (l_EWVLTApr16)
real*4 EWV3 ! (0.1nm) [0.02/3.3] Equivalent width of circumstellar
* feature, 2016 Apr VLT data (EWVLTApr16)
real*4 E_EWV3 ! (0.1nm) [0.004/0.2]? Plus Error on EWV3 (E_EWVLTApr16)
real*4 e_EWV3_1 ! (0.1nm) ? Minus Error on EWV3 (e_EWVLTApr16_2)
character*1 Comm ! Comments regarding the fit:
* i=isolated, b=blended
C=============================================================================
C Loading file 'table1.dat' ! Detected circumstellar absorption features and
* fit parameters
C Format for file interpretation
1 format(
+ A5,1X,F8.3,1X,F5.2,1X,F5.2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F6.1,1X,F5.1,1X,F5.1,1X,F6.1,1X,F5.1,1X,F5.1,1X,A1,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.1,1X,F5.1,1X,F5.1,
+ 1X,F5.1,1X,F5.1,1X,F5.1,1X,A1,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F6.1,1X,F5.1,1X,F5.1,1X,F6.1,1X,F5.1,1X,F5.1,
+ 1X,A1,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.1,1X,
+ F5.1,1X,F5.1,1X,F6.1,1X,F5.1,1X,F5.1,1X,A1,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.1,1X,F5.1,1X,F5.1,1X,F5.1,
+ 1X,F5.1,1X,F5.1,1X,A1,F5.3,1X,F5.3,1X,F5.3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,256
read(1,'(A408)')ar__
read(ar__,1)
+ Ion,lambda,Elow,Eup,loggf,DepthK1,E_DepthK1,e_DepthK1_1,
+ VminK1,E_VminK1,e_VminK1_1,VmaxK1,E_VmaxK1,e_VmaxK1_1,l_EWK1,
+ EWK1,E_EWK1,e_EWK1_1,DepthK2,E_DepthK2,e_DepthK2_1,VminK2,
+ E_VminK2,e_VminK2_1,VmaxK2,E_VmaxK2,e_VmaxK2_1,l_EWK2,EWK2,
+ E_EWK2,e_EWK2_1,DepthV1,E_DepthV1,e_DepthV1_1,VminV1,E_VminV1,
+ e_VminV1_1,VmaxV1,E_VmaxV1,e_VmaxV1_1,l_EWV1,EWV1,E_EWV1,
+ e_EWV1_1,DepthV2,E_DepthV2,e_DepthV2_1,VminV2,E_VminV2,
+ e_VminV2_1,VmaxV2,E_VmaxV2,e_VmaxV2_1,l_EWV2,EWV2,E_EWV2,
+ e_EWV2_1,DepthV3,E_DepthV3,e_DepthV3_1,VminV3,E_VminV3,
+ e_VminV3_1,VmaxV3,E_VmaxV3,e_VmaxV3_1,l_EWV3,EWV3,E_EWV3,
+ e_EWV3_1,Comm
if(ar__(35:39) .EQ. '') DepthK1 = rNULL__
if(ar__(41:45) .EQ. '') E_DepthK1 = rNULL__
if(ar__(47:51) .EQ. '') e_DepthK1_1 = rNULL__
if(ar__(53:58) .EQ. '') VminK1 = rNULL__
if(ar__(60:64) .EQ. '') E_VminK1 = rNULL__
if(ar__(66:70) .EQ. '') e_VminK1_1 = rNULL__
if(ar__(72:77) .EQ. '') VmaxK1 = rNULL__
if(ar__(79:83) .EQ. '') E_VmaxK1 = rNULL__
if(ar__(85:89) .EQ. '') e_VmaxK1_1 = rNULL__
if(ar__(92:96) .EQ. '') EWK1 = rNULL__
if(ar__(98:102) .EQ. '') E_EWK1 = rNULL__
if(ar__(104:108) .EQ. '') e_EWK1_1 = rNULL__
if(ar__(110:114) .EQ. '') DepthK2 = rNULL__
if(ar__(116:120) .EQ. '') E_DepthK2 = rNULL__
if(ar__(122:126) .EQ. '') e_DepthK2_1 = rNULL__
if(ar__(128:133) .EQ. '') VminK2 = rNULL__
if(ar__(135:139) .EQ. '') E_VminK2 = rNULL__
if(ar__(141:145) .EQ. '') e_VminK2_1 = rNULL__
if(ar__(147:151) .EQ. '') VmaxK2 = rNULL__
if(ar__(153:157) .EQ. '') E_VmaxK2 = rNULL__
if(ar__(159:163) .EQ. '') e_VmaxK2_1 = rNULL__
if(ar__(166:170) .EQ. '') EWK2 = rNULL__
if(ar__(172:176) .EQ. '') E_EWK2 = rNULL__
if(ar__(178:182) .EQ. '') e_EWK2_1 = rNULL__
if(ar__(184:188) .EQ. '') DepthV1 = rNULL__
if(ar__(190:194) .EQ. '') E_DepthV1 = rNULL__
if(ar__(196:200) .EQ. '') e_DepthV1_1 = rNULL__
if(ar__(202:207) .EQ. '') VminV1 = rNULL__
if(ar__(209:213) .EQ. '') E_VminV1 = rNULL__
if(ar__(215:219) .EQ. '') e_VminV1_1 = rNULL__
if(ar__(221:226) .EQ. '') VmaxV1 = rNULL__
if(ar__(228:232) .EQ. '') E_VmaxV1 = rNULL__
if(ar__(234:238) .EQ. '') e_VmaxV1_1 = rNULL__
if(ar__(247:251) .EQ. '') E_EWV1 = rNULL__
if(ar__(253:257) .EQ. '') e_EWV1_1 = rNULL__
if(ar__(259:263) .EQ. '') DepthV2 = rNULL__
if(ar__(265:269) .EQ. '') E_DepthV2 = rNULL__
if(ar__(271:275) .EQ. '') e_DepthV2_1 = rNULL__
if(ar__(277:282) .EQ. '') VminV2 = rNULL__
if(ar__(284:288) .EQ. '') E_VminV2 = rNULL__
if(ar__(290:294) .EQ. '') e_VminV2_1 = rNULL__
if(ar__(296:301) .EQ. '') VmaxV2 = rNULL__
if(ar__(303:307) .EQ. '') E_VmaxV2 = rNULL__
if(ar__(309:313) .EQ. '') e_VmaxV2_1 = rNULL__
if(ar__(322:326) .EQ. '') E_EWV2 = rNULL__
if(ar__(328:332) .EQ. '') e_EWV2_1 = rNULL__
if(ar__(334:338) .EQ. '') DepthV3 = rNULL__
if(ar__(340:344) .EQ. '') E_DepthV3 = rNULL__
if(ar__(346:350) .EQ. '') e_DepthV3_1 = rNULL__
if(ar__(352:357) .EQ. '') VminV3 = rNULL__
if(ar__(359:363) .EQ. '') E_VminV3 = rNULL__
if(ar__(365:369) .EQ. '') e_VminV3_1 = rNULL__
if(ar__(371:375) .EQ. '') VmaxV3 = rNULL__
if(ar__(377:381) .EQ. '') E_VmaxV3 = rNULL__
if(ar__(383:387) .EQ. '') e_VmaxV3_1 = rNULL__
if(ar__(396:400) .EQ. '') E_EWV3 = rNULL__
if(ar__(402:406) .EQ. '') e_EWV3_1 = rNULL__
c ..............Just test output...........
write(6,1)
+ Ion,lambda,Elow,Eup,loggf,DepthK1,E_DepthK1,e_DepthK1_1,
+ VminK1,E_VminK1,e_VminK1_1,VmaxK1,E_VmaxK1,e_VmaxK1_1,l_EWK1,
+ EWK1,E_EWK1,e_EWK1_1,DepthK2,E_DepthK2,e_DepthK2_1,VminK2,
+ E_VminK2,e_VminK2_1,VmaxK2,E_VmaxK2,e_VmaxK2_1,l_EWK2,EWK2,
+ E_EWK2,e_EWK2_1,DepthV1,E_DepthV1,e_DepthV1_1,VminV1,E_VminV1,
+ e_VminV1_1,VmaxV1,E_VmaxV1,e_VmaxV1_1,l_EWV1,EWV1,E_EWV1,
+ e_EWV1_1,DepthV2,E_DepthV2,e_DepthV2_1,VminV2,E_VminV2,
+ e_VminV2_1,VmaxV2,E_VmaxV2,e_VmaxV2_1,l_EWV2,EWV2,E_EWV2,
+ e_EWV2_1,DepthV3,E_DepthV3,e_DepthV3_1,VminV3,E_VminV3,
+ e_VminV3_1,VmaxV3,E_VmaxV3,e_VmaxV3_1,l_EWV3,EWV3,E_EWV3,
+ e_EWV3_1,Comm
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end