Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/625/A69 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 2026-Oct-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/625/A69 Plage area composite series (Chatzistergos+, 2019)
*================================================================================
*Analysis of full disc Ca II K spectroheliograms.
*II. Towards an accurate assessment of long-term variations in plage areas.
* Chatzistergos T., Ermolli I., Krivova N.A., Solanki S.K.
* <Astron. Astrophys. 625, A69 (2019)>
* =2019A&A...625A..69C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'comp_a.dat' ! Plage series, annual values
integer*4 nr__
parameter (nr__=126) ! Number of records
character*20 ar__ ! Full-size record
integer*4 Year (nr__) ! (yr) Year of observation
real*8 plagedf (nr__) ! ?=-99 Projected plage areas in disc fraction (1)
real*8 plagemh (nr__) ! ?=-99 Corrected for foreshortening fractional
* plage areas (1)
*Note (1): missing values are denoted with -99.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'comp_m.dat' ! Plage series, monthly values
integer*4 nr__1
parameter (nr__1=1512) ! Number of records
character*25 ar__1 ! Full-size record
integer*4 Year_1 (nr__1) ! (yr) Year of observation
integer*4 Month (nr__1) ! ("month") Month of observation
real*8 plagedf_1 (nr__1) ! ?=-99 Projected plage areas in disc fraction (1)
real*8 plagemh_1 (nr__1) ! ?=-99 Corrected for foreshortening fractional
* plage areas (1)
*Note (1): missing values are denoted with -99.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'comp_d.dat' ! Plage series, daily values
integer*4 nr__2
parameter (nr__2=46872) ! Number of records
character*30 ar__2 ! Full-size record
integer*4 Year_2 (nr__2) ! (yr) Year of observation
integer*4 Month_1 (nr__2) ! ("month") Month of observation
integer*4 Day (nr__2) ! (d) Day of observation
real*8 plagedf_2 (nr__2) ! ?=-99 Projected plage areas in disc fraction (1)
real*8 plagemh_2 (nr__2) ! ?=-99 Corrected for foreshortening fractional
* plage areas (1)
*Note (1): missing values are denoted with -99.
C=============================================================================
C Loading file 'comp_a.dat' ! Plage series, annual values
C Format for file interpretation
1 format(I4,1X,F7.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'comp_a.dat')
write(6,*) '....Loading file: comp_a.dat'
do i__=1,126
read(1,'(A20)')ar__
read(ar__,1)Year(i__),plagedf(i__),plagemh(i__)
c ..............Just test output...........
write(6,1)Year(i__),plagedf(i__),plagemh(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'comp_m.dat' ! Plage series, monthly values
C Format for file interpretation
2 format(I4,3X,I2,1X,F7.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'comp_m.dat')
write(6,*) '....Loading file: comp_m.dat'
do i__=1,1512
read(1,'(A25)')ar__1
read(ar__1,2)
+ Year_1(i__),Month(i__),plagedf_1(i__),plagemh_1(i__)
c ..............Just test output...........
write(6,2)
+ Year_1(i__),Month(i__),plagedf_1(i__),plagemh_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'comp_d.dat' ! Plage series, daily values
C Format for file interpretation
3 format(I4,3X,I2,3X,I2,1X,F7.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'comp_d.dat')
write(6,*) '....Loading file: comp_d.dat'
do i__=1,46872
read(1,'(A30)')ar__2
read(ar__2,3)
+ Year_2(i__),Month_1(i__),Day(i__),plagedf_2(i__),
+ plagemh_2(i__)
c ..............Just test output...........
write(6,3)
+ Year_2(i__),Month_1(i__),Day(i__),plagedf_2(i__),
+ plagemh_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end