Conversion of standardized ReadMe file for
catalog 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-26
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/483/661 Refractive index for silicon carbide (SiC) (Pitman+, 2008)
*================================================================================
*Optical properties of silicon carbide for astrophysical applications.
*I. New laboratory infrared reflectance spectra and optical constants.
* Pitman K.M., Hofmeister A. M., Corman A.B., Speck A.K.
* <Astron. Astrophys. 483, 661 (2008)>
* =2008A&A...483..661P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'betasic.dat' ! Optical constants for 3C-SiC chemical vapor
deposited (CVD) wafer (table 3)
integer*4 nr__
parameter (nr__=105) ! Number of records
character*37 ar__ ! Full-size record
real*4 lambda (nr__) ! (um) Wavelength
real*4 n (nr__) ! Real part of complex refractive index
real*4 k (nr__) ! Imaginary part of complex refractive index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'grysicpa.dat' ! Optical constants for gray 6H-SiC powder,
extraordinary ray (table 4)
integer*4 nr__1
parameter (nr__1=105) ! Number of records
character*37 ar__1 ! Full-size record
real*4 lambda_1 (nr__1) ! (um) Wavelength
real*4 n_1 (nr__1) ! Real part of complex refractive index
real*4 k_1 (nr__1) ! Imaginary part of complex refractive index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'grysicpe.dat' ! Optical constants for gray 6H-SiC powder,
ordinary ray (table 5)
integer*4 nr__2
parameter (nr__2=105) ! Number of records
character*37 ar__2 ! Full-size record
real*4 lambda_2 (nr__2) ! (um) Wavelength
real*4 n_2 (nr__2) ! Real part of complex refractive index
real*4 k_2 (nr__2) ! Imaginary part of complex refractive index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'moisicpa.dat' ! Optical constants for moissanite gem,
extraordinary ray (table 6 )
integer*4 nr__3
parameter (nr__3=105) ! Number of records
character*37 ar__3 ! Full-size record
real*4 lambda_3 (nr__3) ! (um) Wavelength
real*4 n_3 (nr__3) ! Real part of complex refractive index
real*4 k_3 (nr__3) ! Imaginary part of complex refractive index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'moisicpe.dat' ! Optical constants for moissanite gem,
ordinary ray (table 7)
integer*4 nr__4
parameter (nr__4=105) ! Number of records
character*37 ar__4 ! Full-size record
real*4 lambda_4 (nr__4) ! (um) Wavelength
real*4 n_4 (nr__4) ! Real part of complex refractive index
real*4 k_4 (nr__4) ! Imaginary part of complex refractive index
C=============================================================================
C Loading file 'betasic.dat' ! Optical constants for 3C-SiC chemical vapor
* deposited (CVD) wafer (table 3)
C Format for file interpretation
1 format(E11.6,1X,E11.7,1X,E13.6)
C Effective file loading
open(unit=1,status='old',file=
+'betasic.dat')
write(6,*) '....Loading file: betasic.dat'
do i__=1,105
read(1,'(A37)')ar__
read(ar__,1)lambda(i__),n(i__),k(i__)
c ..............Just test output...........
write(6,1)lambda(i__),n(i__),k(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'grysicpa.dat' ! Optical constants for gray 6H-SiC powder,
* extraordinary ray (table 4)
C Format for file interpretation
2 format(E11.6,1X,E11.7,1X,E13.6)
C Effective file loading
open(unit=1,status='old',file=
+'grysicpa.dat')
write(6,*) '....Loading file: grysicpa.dat'
do i__=1,105
read(1,'(A37)')ar__1
read(ar__1,2)lambda_1(i__),n_1(i__),k_1(i__)
c ..............Just test output...........
write(6,2)lambda_1(i__),n_1(i__),k_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'grysicpe.dat' ! Optical constants for gray 6H-SiC powder,
* ordinary ray (table 5)
C Format for file interpretation
3 format(E11.6,1X,E11.7,1X,E13.6)
C Effective file loading
open(unit=1,status='old',file=
+'grysicpe.dat')
write(6,*) '....Loading file: grysicpe.dat'
do i__=1,105
read(1,'(A37)')ar__2
read(ar__2,3)lambda_2(i__),n_2(i__),k_2(i__)
c ..............Just test output...........
write(6,3)lambda_2(i__),n_2(i__),k_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'moisicpa.dat' ! Optical constants for moissanite gem,
* extraordinary ray (table 6 )
C Format for file interpretation
4 format(E11.6,1X,E11.7,1X,E13.6)
C Effective file loading
open(unit=1,status='old',file=
+'moisicpa.dat')
write(6,*) '....Loading file: moisicpa.dat'
do i__=1,105
read(1,'(A37)')ar__3
read(ar__3,4)lambda_3(i__),n_3(i__),k_3(i__)
c ..............Just test output...........
write(6,4)lambda_3(i__),n_3(i__),k_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'moisicpe.dat' ! Optical constants for moissanite gem,
* ordinary ray (table 7)
C Format for file interpretation
5 format(E11.6,1X,E11.7,1X,E13.6)
C Effective file loading
open(unit=1,status='old',file=
+'moisicpe.dat')
write(6,*) '....Loading file: moisicpe.dat'
do i__=1,105
read(1,'(A37)')ar__4
read(ar__4,5)lambda_4(i__),n_4(i__),k_4(i__)
c ..............Just test output...........
write(6,5)lambda_4(i__),n_4(i__),k_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end