Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/599/A85 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/599/A85 Structure of Herbig AeBe disks at mas scale (Lazareff+, 2017)
*================================================================================
*Structure of Herbig AeBe disks at milliarcsecond scale.
*A statistical survey in H band using PIONIER-VLTI.
* Lazareff B., Berger J.-P., Kluska J., Le Bouquin J.-B., Benisty M.,
* Malbet F., Koen C., Pinte C., Thi W.-F., Absil O., Baron F., Delboulbe A.,
* Duvert G., Isella A., Jocou L., Juhasz A., Kraus S., Lachaume R.,
* Menard F., Millan-Gabet R., Monnier J.D., Moulin T., Perraut K., Rochat S.,
* Soulez F., Tallon M., Thiebaut E., Traub W., Zins G.
* <Astron. Astrophys. 599, A85 (2017)>
* =2017A&A...599A..85L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sed_fit.dat' ! SED fits
integer*4 nr__
parameter (nr__=51) ! Number of records
character*99 ar__ ! Full-size record
character*12 Name ! Object name
real*4 chi2r ! reduced chi2 of fit
real*4 lnFnsV ! ([mW/m2]) log10(nuFnu), stellar, V band
real*4 e_lnFnsV ! ([mW/m2]) Standard error for lnFnsV
real*4 AV ! (mag) Extinction in V band
real*4 e_AV ! (mag) Standard error for AV
real*8 lnFndH ! ([mW/m2]) log10(nuFnu), dust, H band
real*4 e_lnFndH ! ([mW/m2]) Standard error for lnFndH
real*4 logTd ! ([K]) log10 dust temperature (photometric)
real*4 e_logTd ! ([K]) Standard error for logTd
real*4 fd ! Fractional flux contribution of dust at
* H band
real*4 e_fd ! Standard error for fd
character*4 SpType ! Adopted spectral type
real*4 logTeff ! ([K]) log10 effective temperature
real*4 BCv ! (mag) Bolometric correction from V magnitude
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ellip.dat' ! Ellipsoid fits
integer*4 nr__1
parameter (nr__1=44) ! Number of records
character*108 ar__1 ! Full-size record
character*12 Name_1 ! Object name
real*4 chi2r_1 ! Reduced chi2 of fit
real*4 kc ! ?=9.99[-6.0,+2.0] Spectral index of
* circumstellar component
real*4 e_kc ! ?=9.99 Standard error for kc
real*4 fs ! [0,1] Fractional flux stellar component
real*4 e_fs ! Standard error for fs
real*4 fh ! [0,1] Fractional flux "halo" component
real*4 e_fh ! Standard error for fh
real*4 cosi ! Cosine inclination (1.0 face-on)
real*4 e_cosi ! Standard error for fh
integer*4 PA ! (deg) Position angle of major axis (East from North)
integer*4 e_PA ! (deg) Standard error for PA
real*4 flor ! [0,1] Weighting for radial profile
real*4 e_flor ! Standard error for flor
real*4 logla ! ([mas]) log10 half light radius
real*4 e_logla ! ([mas]) Standard error for logla
real*4 fc ! [0,1] Fractional flux circumstellar component
real*4 e_fc ! Standard error for fc
real*4 logTdi ! ([K]) log10 dust temperature (interferometric)
real*4 e_logTdi ! ([K]) Standard error for logTdi
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ring_m0.dat' ! Ring fits
integer*4 nr__2
parameter (nr__2=27) ! Number of records
character*117 ar__2 ! Full-size record
character*12 Name_2 ! Object name
real*4 chi2r_2 ! Reduced chi2 of fit
real*4 kc_1 ! ?=9.99[-6.0,+2.0] Spectral index of
* circumstellar component
real*4 e_kc_1 ! ?=9.99 Standard error for kc
real*4 fs_1 ! [0,1] Fractional flux stellar component
real*4 e_fs_1 ! Standard error for fs
real*4 fh_1 ! [0,1] Fractional flux "halo" component
real*4 e_fh_1 ! Standard error for fh
real*4 cosi_1 ! cosine inclination (1.0 face-on)
real*4 e_cosi_1 ! Standard error for fh
integer*4 PA_1 ! (deg) Position angle of major axis (East from North)
integer*4 e_PA_1 ! (deg) Standard error for PA
real*4 flor_1 ! [0,1] Weighting for radial profile
real*4 e_flor_1 ! Standard error for flor
real*4 logla_1 ! ([mas]) log10 half light radius
real*4 e_logla_1 ! ([mas]) Standard error for logla
real*4 fr ! [0,1] Fractional flux ring component
real*4 e_fr ! Standard error for fr
real*4 lrkr ! log10 ratio kernel radius
real*4 e_lrkr ! Standard error for lrkr
real*4 logTdi_1 ! ([K]) log10 dust temperature (interferometric)
real*4 e_logTdi_1 ! ([K]) Standard error for logTdi
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ring_m1.dat' ! Ring fits, m=1 azimuthal modulation
integer*4 nr__3
parameter (nr__3=27) ! Number of records
character*140 ar__3 ! Full-size record
character*12 Name_3 ! Object name
real*4 chi2r_3 ! reduced chi2 of fit
real*4 kc_2 ! ?=9.99[-6.0,+2.0] Spectral index of
* circumstellar component
real*4 e_kc_2 ! ?=9.99 Standard error for kc
real*4 fs_2 ! [0,1] Fractionalflux stellar component
real*4 e_fs_2 ! Standard error for fs
real*4 fh_2 ! [0,1] Fractionalflux "halo" component
real*4 e_fh_2 ! Standard error for fh
real*4 cosi_2 ! cosine inclination (1.0 face-on)
real*4 e_cosi_2 ! Standard error for fh
integer*4 PA_2 ! (deg) Position angle of major axis (East from North)
integer*4 e_PA_2 ! (deg) Standard error for PA
real*4 flor_2 ! [0,1] Weighting for radial profile
real*4 e_flor_2 ! Standard error for flor
real*4 logla_2 ! ([mas]) log10 half light radius
real*4 e_logla_2 ! ([mas]) Standard error for logla
real*4 fr_1 ! [0,1] Fractionalflux ring component
real*4 e_fr_1 ! Standard error for fr
real*4 lrkr_1 ! log10 ratio kernel radius
real*4 e_lrkr_1 ! Standard error for lrkr
real*4 mc1 ! cos modulation order 1
real*4 e_mc1 ! Standard error for mc1
real*4 ms1 ! sin modulation order 1
real*4 e_ms1 ! Standard error for ms1
real*4 logTdi_2 ! ([K]) log10 dust temperature (interferometric)
real*4 e_logTdi_2 ! ([K]) Standard error for logTdi
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ring_m2.dat' ! Ring fits, m=2 azimuthal modulation
integer*4 nr__4
parameter (nr__4=27) ! Number of records
character*162 ar__4 ! Full-size record
character*12 Name_4 ! Object name
real*4 chi2r_4 ! reduced chi2 of fit
real*4 kc_3 ! ?=9.99[-6.0,+2.0] Spectral index of
* circumstellar component
real*4 e_kc_3 ! ?=9.99 standard error for kc
real*4 fs_3 ! [0,1] Fractional flux stellar component
real*4 e_fs_3 ! Standard error for fs
real*4 fh_3 ! [0,1] Fractional flux "halo" component
real*4 e_fh_3 ! Standard error for fh
real*4 cosi_3 ! cosine inclination (1.0 face-on)
real*4 e_cosi_3 ! Standard error for fh
integer*4 PA_3 ! (deg) Position angle of major axis (East from North)
integer*4 e_PA_3 ! (deg) Standard error for PA
real*4 flor_3 ! [0,1] Weighting for radial profile
real*4 e_flor_3 ! Standard error for flor
real*4 logla_3 ! ([mas]) log10 half light radius
real*4 e_logla_3 ! ([mas]) Standard error for logla
real*4 fr_2 ! [0,1] Fractional flux ring component
real*4 e_fr_2 ! Standard error for fr
real*4 lrkr_2 ! log10 ratio kernel radius
real*4 e_lrkr_2 ! Standard error for lrkr
real*4 mc1_1 ! cos modulation order 1
real*4 e_mc1_1 ! Standard error for mc1
real*4 ms1_1 ! sin modulation order 1
real*4 e_ms1_1 ! Standard error for ms1
real*4 mc2 ! cos modulation order 2
real*4 e_mc2 ! Standard error for mc2
real*4 ms2 ! sin modulation order 2
real*4 e_ms2 ! Standard error for ms2
real*4 logTdi_3 ! ([K]) log10 dust temperature (interferometric)
real*4 e_logTdi_3 ! ([K]) Standard error for logTdi
C=============================================================================
C Loading file 'sed_fit.dat' ! SED fits
C Format for file interpretation
1 format(
+ A12,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F7.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,A4,1X,F5.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'sed_fit.dat')
write(6,*) '....Loading file: sed_fit.dat'
do i__=1,51
read(1,'(A99)')ar__
read(ar__,1)
+ Name,chi2r,lnFnsV,e_lnFnsV,AV,e_AV,lnFndH,e_lnFndH,logTd,
+ e_logTd,fd,e_fd,SpType,logTeff,BCv
c ..............Just test output...........
write(6,1)
+ Name,chi2r,lnFnsV,e_lnFnsV,AV,e_AV,lnFndH,e_lnFndH,logTd,
+ e_logTd,fd,e_fd,SpType,logTeff,BCv
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ellip.dat' ! Ellipsoid fits
C Format for file interpretation
2 format(
+ A12,1X,F6.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,I3,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'ellip.dat')
write(6,*) '....Loading file: ellip.dat'
do i__=1,44
read(1,'(A108)')ar__1
read(ar__1,2)
+ Name_1,chi2r_1,kc,e_kc,fs,e_fs,fh,e_fh,cosi,e_cosi,PA,e_PA,
+ flor,e_flor,logla,e_logla,fc,e_fc,logTdi,e_logTdi
c ..............Just test output...........
write(6,2)
+ Name_1,chi2r_1,kc,e_kc,fs,e_fs,fh,e_fh,cosi,e_cosi,PA,e_PA,
+ flor,e_flor,logla,e_logla,fc,e_fc,logTdi,e_logTdi
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ring_m0.dat' ! Ring fits
C Format for file interpretation
3 format(
+ A12,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,I3,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'ring_m0.dat')
write(6,*) '....Loading file: ring_m0.dat'
do i__=1,27
read(1,'(A117)')ar__2
read(ar__2,3)
+ Name_2,chi2r_2,kc_1,e_kc_1,fs_1,e_fs_1,fh_1,e_fh_1,cosi_1,
+ e_cosi_1,PA_1,e_PA_1,flor_1,e_flor_1,logla_1,e_logla_1,fr,
+ e_fr,lrkr,e_lrkr,logTdi_1,e_logTdi_1
c ..............Just test output...........
write(6,3)
+ Name_2,chi2r_2,kc_1,e_kc_1,fs_1,e_fs_1,fh_1,e_fh_1,cosi_1,
+ e_cosi_1,PA_1,e_PA_1,flor_1,e_flor_1,logla_1,e_logla_1,fr,
+ e_fr,lrkr,e_lrkr,logTdi_1,e_logTdi_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ring_m1.dat' ! Ring fits, m=1 azimuthal modulation
C Format for file interpretation
4 format(
+ A12,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,I3,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'ring_m1.dat')
write(6,*) '....Loading file: ring_m1.dat'
do i__=1,27
read(1,'(A140)')ar__3
read(ar__3,4)
+ Name_3,chi2r_3,kc_2,e_kc_2,fs_2,e_fs_2,fh_2,e_fh_2,cosi_2,
+ e_cosi_2,PA_2,e_PA_2,flor_2,e_flor_2,logla_2,e_logla_2,fr_1,
+ e_fr_1,lrkr_1,e_lrkr_1,mc1,e_mc1,ms1,e_ms1,logTdi_2,e_logTdi_2
c ..............Just test output...........
write(6,4)
+ Name_3,chi2r_3,kc_2,e_kc_2,fs_2,e_fs_2,fh_2,e_fh_2,cosi_2,
+ e_cosi_2,PA_2,e_PA_2,flor_2,e_flor_2,logla_2,e_logla_2,fr_1,
+ e_fr_1,lrkr_1,e_lrkr_1,mc1,e_mc1,ms1,e_ms1,logTdi_2,e_logTdi_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ring_m2.dat' ! Ring fits, m=2 azimuthal modulation
C Format for file interpretation
5 format(
+ A12,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,I3,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'ring_m2.dat')
write(6,*) '....Loading file: ring_m2.dat'
do i__=1,27
read(1,'(A162)')ar__4
read(ar__4,5)
+ Name_4,chi2r_4,kc_3,e_kc_3,fs_3,e_fs_3,fh_3,e_fh_3,cosi_3,
+ e_cosi_3,PA_3,e_PA_3,flor_3,e_flor_3,logla_3,e_logla_3,fr_2,
+ e_fr_2,lrkr_2,e_lrkr_2,mc1_1,e_mc1_1,ms1_1,e_ms1_1,mc2,e_mc2,
+ ms2,e_ms2,logTdi_3,e_logTdi_3
c ..............Just test output...........
write(6,5)
+ Name_4,chi2r_4,kc_3,e_kc_3,fs_3,e_fs_3,fh_3,e_fh_3,cosi_3,
+ e_cosi_3,PA_3,e_PA_3,flor_3,e_flor_3,logla_3,e_logla_3,fr_2,
+ e_fr_2,lrkr_2,e_lrkr_2,mc1_1,e_mc1_1,ms1_1,e_ms1_1,mc2,e_mc2,
+ ms2,e_ms2,logTdi_3,e_logTdi_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end