FORTRAN Generation
(/./ftp/cats/J/A_A/390/717)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/390/717 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-Mar-29
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/390/717  Solar Velocities of the Sun in 1989-99 (IRIS++) (Salabert+, 2002)
*================================================================================
*IRIS++ database: Merging of IRIS + Mark-1 + LOWL
*    Salabert D., Fossat E., Gelly B., Tomczyk S., Palle P.L.,
*    Jimenez-Reyes S.J., Cacciani A., Corbard T., Ehgamberdiev S., Grec G.,
*    Hoeksema J.T., Kholikov S., Lazrek M., Schmider F.X.
*   <Astron. Astrophys. 390, 717 (2002)>
*   =2002A&A...390..717S
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v89f.dat'	! Year 1989 pre-merged

      integer*4 nr__
      parameter (nr__=352142)	! Number of records
      character*35 ar__   	! Full-size record

      real*8        MJD        (nr__) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res        (nr__) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans      (nr__) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind        (nr__) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v90f.dat'	! Year 1990 pre-merged

      integer*4 nr__1
      parameter (nr__1=699609)	! Number of records
      character*35 ar__1  	! Full-size record

      real*8        MJD_1      (nr__1) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_1      (nr__1) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_1    (nr__1) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_1      (nr__1) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v91f.dat'	! Year 1991 pre-merged

      integer*4 nr__2
      parameter (nr__2=699101)	! Number of records
      character*35 ar__2  	! Full-size record

      real*8        MJD_2      (nr__2) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_2      (nr__2) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_2    (nr__2) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_2      (nr__2) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v92f.dat'	! Year 1992 pre-merged

      integer*4 nr__3
      parameter (nr__3=701963)	! Number of records
      character*35 ar__3  	! Full-size record

      real*8        MJD_3      (nr__3) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_3      (nr__3) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_3    (nr__3) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_3      (nr__3) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v93f.dat'	! Year 1993 pre-merged

      integer*4 nr__4
      parameter (nr__4=325813)	! Number of records
      character*35 ar__4  	! Full-size record

      real*8        MJD_4      (nr__4) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_4      (nr__4) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_4    (nr__4) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_4      (nr__4) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v94f.dat'	! Year 1994 pre-merged

      integer*4 nr__5
      parameter (nr__5=696614)	! Number of records
      character*35 ar__5  	! Full-size record

      real*8        MJD_5      (nr__5) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_5      (nr__5) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_5    (nr__5) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_5      (nr__5) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v95f.dat'	! Year 1995 pre-merged

      integer*4 nr__6
      parameter (nr__6=698556)	! Number of records
      character*35 ar__6  	! Full-size record

      real*8        MJD_6      (nr__6) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_6      (nr__6) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_6    (nr__6) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_6      (nr__6) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v96f.dat'	! Year 1996 pre-merged

      integer*4 nr__7
      parameter (nr__7=702161)	! Number of records
      character*35 ar__7  	! Full-size record

      real*8        MJD_7      (nr__7) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_7      (nr__7) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_7    (nr__7) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_7      (nr__7) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v97f.dat'	! Year 1997 pre-merged

      integer*4 nr__8
      parameter (nr__8=698622)	! Number of records
      character*35 ar__8  	! Full-size record

      real*8        MJD_8      (nr__8) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_8      (nr__8) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_8    (nr__8) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_8      (nr__8) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v98f.dat'	! Year 1998 pre-merged

      integer*4 nr__9
      parameter (nr__9=699232)	! Number of records
      character*35 ar__9  	! Full-size record

      real*8        MJD_9      (nr__9) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_9      (nr__9) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_9    (nr__9) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_9      (nr__9) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'v99f.dat'	! Year 1999 pre-merged

      integer*4 nr__10
      parameter (nr__10=668952)	! Number of records
      character*35 ar__10 	! Full-size record

      real*8        MJD_10     (nr__10) ! (d) Modified Julian Date (sampling 45 seconds) (1)
      real*8        Res_10     (nr__10) ! (m/s) ?=0.000000 Radial velocity (2)
      real*8        Trans_10   (nr__10) ! ?=0.00 Transparency coefficient (2)
      integer*4     Ind_10     (nr__10) ! ?=0000 Index of information redundancy (3)
*Note (1): The first point in each file corresponds to the first observing
*          point, and not to the beginning of each year.
*    Low-frequencies are high-pass filtered, the cutoff frequency being at
*    1.1 mHz.
*    The integer part of the value is the number of days elapsed since
*    Nov 17, 1858 at 0 UT, and the fractional part is the UT fraction of
*    day. A 64 bits data format allows to store the MJD with a resolution
*    better than 1 sec which is quite enough for us. MJD is deduced from
*    the classical Julian Date by subtracting 2400000.5 days to the JD.
*Note (2): The missing points have a value of zero.
*Note (3): The last column is a synoptic index indicating how many sites
*    have contributed to the corresponding data points. Each site used in
*    the merging sets a given bit inside a 16 bits integer word.
*    The bits meaning is as follow:
*
*    Bit # Value  Site
*        1     1  MK = Izana [Tenerife, Spain] Mark-1 instrument
*        2     2  IZ = Izana [Tenerife, Spain] IRIS instrument
*        3     4  KU = Kumbel [Uzbekistan] IRIS instrument
*        4     8  LA = La Silla [Chile] IRIS instrument
*        5    16  CU = Culgoora [Australia] IRIS instrument
*        6    32  OU = Oukaimden [Morocco] IRIS instrument
*        7    64  NI = Nice [France] IRIS instrument
*        8   128  ST = Stanford [California, USA] IRIS instrument
*        9   256  JP = Jet Prop Lab [California, USA] Cacciani's MOF instrument
*       10   512  ML = Mauna Loa [Hawaii] LOWL instrument
*       11  1024  PK = Parkent [Uzbekistan] IRIS instrument
*
*    hence, a data point merged from Izana alone would have a value of 2,
*    and a data point merged from the Oukaimden IRIS site (OU) and the LOWL
*    site (ML) would have a value of 544.

C=============================================================================

C  Loading file 'v89f.dat'	! Year 1989 pre-merged

C  Format for file interpretation

    1 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v89f.dat')
      write(6,*) '....Loading file: v89f.dat'
      do i__=1,352142
        read(1,'(A35)')ar__
        read(ar__,1)MJD(i__),Res(i__),Trans(i__),Ind(i__)
c    ..............Just test output...........
        write(6,1)MJD(i__),Res(i__),Trans(i__),Ind(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v90f.dat'	! Year 1990 pre-merged

C  Format for file interpretation

    2 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v90f.dat')
      write(6,*) '....Loading file: v90f.dat'
      do i__=1,699609
        read(1,'(A35)')ar__1
        read(ar__1,2)MJD_1(i__),Res_1(i__),Trans_1(i__),Ind_1(i__)
c    ..............Just test output...........
        write(6,2)MJD_1(i__),Res_1(i__),Trans_1(i__),Ind_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v91f.dat'	! Year 1991 pre-merged

C  Format for file interpretation

    3 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v91f.dat')
      write(6,*) '....Loading file: v91f.dat'
      do i__=1,699101
        read(1,'(A35)')ar__2
        read(ar__2,3)MJD_2(i__),Res_2(i__),Trans_2(i__),Ind_2(i__)
c    ..............Just test output...........
        write(6,3)MJD_2(i__),Res_2(i__),Trans_2(i__),Ind_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v92f.dat'	! Year 1992 pre-merged

C  Format for file interpretation

    4 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v92f.dat')
      write(6,*) '....Loading file: v92f.dat'
      do i__=1,701963
        read(1,'(A35)')ar__3
        read(ar__3,4)MJD_3(i__),Res_3(i__),Trans_3(i__),Ind_3(i__)
c    ..............Just test output...........
        write(6,4)MJD_3(i__),Res_3(i__),Trans_3(i__),Ind_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v93f.dat'	! Year 1993 pre-merged

C  Format for file interpretation

    5 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v93f.dat')
      write(6,*) '....Loading file: v93f.dat'
      do i__=1,325813
        read(1,'(A35)')ar__4
        read(ar__4,5)MJD_4(i__),Res_4(i__),Trans_4(i__),Ind_4(i__)
c    ..............Just test output...........
        write(6,5)MJD_4(i__),Res_4(i__),Trans_4(i__),Ind_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v94f.dat'	! Year 1994 pre-merged

C  Format for file interpretation

    6 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v94f.dat')
      write(6,*) '....Loading file: v94f.dat'
      do i__=1,696614
        read(1,'(A35)')ar__5
        read(ar__5,6)MJD_5(i__),Res_5(i__),Trans_5(i__),Ind_5(i__)
c    ..............Just test output...........
        write(6,6)MJD_5(i__),Res_5(i__),Trans_5(i__),Ind_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v95f.dat'	! Year 1995 pre-merged

C  Format for file interpretation

    7 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v95f.dat')
      write(6,*) '....Loading file: v95f.dat'
      do i__=1,698556
        read(1,'(A35)')ar__6
        read(ar__6,7)MJD_6(i__),Res_6(i__),Trans_6(i__),Ind_6(i__)
c    ..............Just test output...........
        write(6,7)MJD_6(i__),Res_6(i__),Trans_6(i__),Ind_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v96f.dat'	! Year 1996 pre-merged

C  Format for file interpretation

    8 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v96f.dat')
      write(6,*) '....Loading file: v96f.dat'
      do i__=1,702161
        read(1,'(A35)')ar__7
        read(ar__7,8)MJD_7(i__),Res_7(i__),Trans_7(i__),Ind_7(i__)
c    ..............Just test output...........
        write(6,8)MJD_7(i__),Res_7(i__),Trans_7(i__),Ind_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v97f.dat'	! Year 1997 pre-merged

C  Format for file interpretation

    9 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v97f.dat')
      write(6,*) '....Loading file: v97f.dat'
      do i__=1,698622
        read(1,'(A35)')ar__8
        read(ar__8,9)MJD_8(i__),Res_8(i__),Trans_8(i__),Ind_8(i__)
c    ..............Just test output...........
        write(6,9)MJD_8(i__),Res_8(i__),Trans_8(i__),Ind_8(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v98f.dat'	! Year 1998 pre-merged

C  Format for file interpretation

   10 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v98f.dat')
      write(6,*) '....Loading file: v98f.dat'
      do i__=1,699232
        read(1,'(A35)')ar__9
        read(ar__9,10)MJD_9(i__),Res_9(i__),Trans_9(i__),Ind_9(i__)
c    ..............Just test output...........
        write(6,10)MJD_9(i__),Res_9(i__),Trans_9(i__),Ind_9(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'v99f.dat'	! Year 1999 pre-merged

C  Format for file interpretation

   11 format(F12.6,1X,F9.6,F7.2,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'v99f.dat')
      write(6,*) '....Loading file: v99f.dat'
      do i__=1,668952
        read(1,'(A35)')ar__10
        read(ar__10,11)MJD_10(i__),Res_10(i__),Trans_10(i__),Ind_10(i__)
c    ..............Just test output...........
        write(6,11)MJD_10(i__),Res_10(i__),Trans_10(i__),Ind_10(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end