FORTRAN Generation
(/./ftp/cats/J/MNRAS/511/1043)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/511/1043 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 2026-May-18
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/MNRAS/511/1043    Sub-Neptunes transiting TOI-1064             (Wilson+, 2022)
*================================================================================
*A pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 characterized
*with CHEOPS.
*    Wilson T.G., Goffo E., Alibert Y., Gandolfi D., Bonfanti A., Persson C.M.,
*    Collier Cameron A., Fridlund M., Fossati L., Korth J., Benz W., Deline A.,
*    Floren H.-G., Guterman P., Adibekyan V., Hooton M.J., Hoyer S., Leleu A.,
*    Mustill A.J., Salmon S., Sousa S.G., Suarez O., Abe L., Agabi A.,
*    Alonso R., Anglada G., Asquier J., Barczy T., Barrado Navascues D.,
*    Barros S.C.C., Baumjohann W., Beck M., Beck T., Billot N., Bonfils X.,
*    Brandeker A., Broeg C., Bryant E.M., Burleigh M.R., Buttu M., Cabrera J.,
*    Charnoz S., Ciardi D.R., Cloutier R., Cochran W.D., Collins K.A.,
*    Colon K.D., Crouzet N., Csizmadia S., Davies M.B., Deleuil M., Delrez L.,
*    Demangeon O., Demory B.-O., Dragomir D., Dransfield G., Ehrenreich D.,
*    Erikson A., Fortier A., Gan T., Gill S., Gillon M., Gnilka C.L.,
*    Grieves N., Grziwa S., Gudel M., Guillot T., Haldemann J., Heng K.,
*    Horne K., Howell S.B., Isaak K.G., Jenkins J.M., Jensen E.L.N., Kiss L.,
*    Lacedelli G., Lam K., Laskar J., Latham D.W., Lecavelier des Etangs A.,
*    Lendl M., Lester K.V., Levine A.M., Livingston J., Lovis C., Luque R.,
*    Magrin D., Marie-Sainte W., Maxted P.F.L., Mayo A.W., McLean B., Mecina M.,
*    Mekarnia D., Nascimbeni V., Nielsen L.D., Olofsson G., Osborn H.P.,
*    Osborne H.L.M., Ottensamer R., Pagano I., Palle E., Peter G., Piotto G.,
*    Pollacco D., Queloz D., Ragazzoni R., Rando N., Rauer H., Redfield S.,
*    Ribas I., Ricker G.R., Rieder M., Santos N.C., Scandariato G.,
*    Schmider F.-X., Schwarz R.P., Scott N.J., Seager S., Segransan D.,
*    Serrano L.M., Simon A.E., Smith A.M.S., Steller M., Stockdale C., Szabo G.,
*    Thomas N., Ting E.B., Triaud A.H.M.J., Udry S., Van Eylen V.,
*    Van Grootel V., Vanderspek R.K., Viotto V., Walton N., Winn J.N.
*   <Mon. Not. R. Astron. Soc., 511, 1043-1071 (2022)>
*   =2022MNRAS.511.1043W    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'cheoraw1.dat'	! CHEOPS visit 1 raw data

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

      real*8        BJD         ! (d) Barycentric Julian date
      real*8        Flux        ! Normalised flux
      real*8        e_Flux      ! Normalised flux error
      real*4        Xoffset     ! (pix) X centroid position offset
      real*8        Yoffset     ! (pix) Y centroid position offset
      real*8        Roll        ! Roll angle
      real*8        Contam      ! (ct) Contamination
      real*8        Bckg        ! (ct) Background
      real*4        Smear       ! (ct) Smear
      real*8        DeltaT      ! Delta temperature

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

C  Declarations for 'cheoraw2.dat'	! CHEOPS visit 2 raw data

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

      real*8        BJD_1       ! (d) Barycentric Julian date
      real*8        Flux_1      ! Normalised flux
      real*8        e_Flux_1    ! Normalised flux error
      real*4        Xoffset_1   ! (pix) X centroid position offset
      real*8        Yoffset_1   ! (pix) Y centroid position offset
      real*8        Roll_1      ! Roll angle
      real*8        Contam_1    ! (ct) Contamination
      real*8        Bckg_1      ! (ct) Background
      real*4        Smear_1     ! (ct) Smear
      real*8        DeltaT_1    ! Delta temperature

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

C  Declarations for 'cheoraw3.dat'	! CHEOPS visit 3 raw data

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

      real*8        BJD_2       ! (d) Barycentric Julian date
      real*8        Flux_2      ! Normalised flux
      real*8        e_Flux_2    ! Normalised flux error
      real*4        Xoffset_2   ! (pix) X centroid position offset
      real*8        Yoffset_2   ! (pix) Y centroid position offset
      real*8        Roll_2      ! Roll angle
      real*8        Contam_2    ! (ct) Contamination
      real*8        Bckg_2      ! (ct) Background
      real*4        Smear_2     ! (ct) Smear
      real*8        DeltaT_2    ! Delta temperature

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

C  Declarations for 'cheoraw4.dat'	! CHEOPS visit 4 raw data

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

      real*8        BJD_3       ! (d) Barycentric Julian date
      real*8        Flux_3      ! Normalised flux
      real*8        e_Flux_3    ! Normalised flux error
      real*4        Xoffset_3   ! (pix) X centroid position offset
      real*8        Yoffset_3   ! (pix) Y centroid position offset
      real*8        Roll_3      ! Roll angle
      real*8        Contam_3    ! (ct) Contamination
      real*8        Bckg_3      ! (ct) Background
      real*4        Smear_3     ! (ct) Smear
      real*8        DeltaT_3    ! Delta temperature

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

C  Declarations for 'cheoraw5.dat'	! CHEOPS visit 5 raw data

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

      real*8        BJD_4       ! (d) Barycentric Julian date
      real*8        Flux_4      ! Normalised flux
      real*8        e_Flux_4    ! Normalised flux error
      real*4        Xoffset_4   ! (pix) X centroid position offset
      real*8        Yoffset_4   ! (pix) Y centroid position offset
      real*8        Roll_4      ! Roll angle
      real*8        Contam_4    ! (ct) Contamination
      real*8        Bckg_4      ! (ct) Background
      real*4        Smear_4     ! (ct) Smear
      real*8        DeltaT_4    ! Delta temperature

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

C  Declarations for 'cheoraw6.dat'	! CHEOPS visit 6 raw data

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

      real*8        BJD_5       ! (d) Barycentric Julian date
      real*8        Flux_5      ! Normalised flux
      real*8        e_Flux_5    ! Normalised flux error
      real*4        Xoffset_5   ! (pix) X centroid position offset
      real*8        Yoffset_5   ! (pix) Y centroid position offset
      real*8        Roll_5      ! Roll angle
      real*8        Contam_5    ! (ct) Contamination
      real*8        Bckg_5      ! (ct) Background
      real*4        Smear_5     ! (ct) Smear
      real*8        DeltaT_5    ! Delta temperature

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

C  Declarations for 'cheopsf1.dat'	! CHEOPS visit 1 PSF SCALPELS data

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

      real*8        BJD_6       ! (d) Barycentric Julian date
      real*8        Flux_6      ! Normalised flux
      real*8        e_Flux_6    ! Normalised flux error
      real*8        Fluxd       ! Normalised detrended flux
      real*8        e_Fluxd     ! Normalised detrended flux error
      real*8        Fluxpsfmodel ! PSF SCALPELS PSF shape change model
      real*4        Theta0      ! 1st PSF SCALPELS principal component
      real*4        Theta1      ! 2nd PSF SCALPELS principal component
      real*8        Theta2      ! 3rd PSF SCALPELS principal component
      real*4        Theta3      ! 4th PSF SCALPELS principal component
      real*8        Theta4      ! 5th PSF SCALPELS principal component
      real*4        Theta5      ! 6th PSF SCALPELS principal component
      real*4        Theta6      ! 7th PSF SCALPELS principal component
      real*4        Theta7      ! 8th PSF SCALPELS principal component
      real*4        Theta8      ! 9th PSF SCALPELS principal component
      real*4        Theta9      ! 10th PSF SCALPELS principal component
      real*4        Theta10     ! ? 11th PSF SCALPELS principal component
      real*4        Theta11     ! ? 12th PSF SCALPELS principal component
      real*4        Theta12     ! ? 13th PSF SCALPELS principal component
      real*4        Theta13     ! ? 14th PSF SCALPELS principal component
      real*4        Theta14     ! ? 15th PSF SCALPELS principal component
      real*4        Theta15     ! ? 16th PSF SCALPELS principal component
      real*4        Theta16     ! ? 17th PSF SCALPELS principal component
      real*4        Theta17     ! ? 18th PSF SCALPELS principal component
      real*4        Theta18     ! ? 19th PSF SCALPELS principal component
      real*4        Theta19     ! ? 20th PSF SCALPELS principal component
      real*4        Theta20     ! ? 21st PSF SCALPELS principal component
      real*4        Theta21     ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22     ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23     ! ? 24th PSF SCALPELS principal component
      real*4        Theta24     ! ? 25th PSF SCALPELS principal component
      real*4        Theta25     ! ? 26th PSF SCALPELS principal component
      real*4        Theta26     ! ? 27th PSF SCALPELS principal component
      real*4        Theta27     ! ? 28th PSF SCALPELS principal component
      real*4        Theta28     ! ? 29th PSF SCALPELS principal component
      real*4        Theta29     ! ? 30th PSF SCALPELS principal component
      real*4        Theta30     ! ? 31st PSF SCALPELS principal component
      real*4        Theta31     ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32     ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33     ! ? 34th PSF SCALPELS principal component
      real*8        Theta34     ! ? 35th PSF SCALPELS principal component
      real*4        Theta35     ! ? 36th PSF SCALPELS principal component
      real*8        Theta36     ! ? 37th PSF SCALPELS principal component
      real*4        Theta37     ! ? 38th PSF SCALPELS principal component
      real*4        Theta38     ! ? 39th PSF SCALPELS principal component
      real*4        Theta39     ! ? 40th PSF SCALPELS principal component
      real*4        Theta40     ! ? 41st PSF SCALPELS principal component
      real*8        Theta41     ! ? 42nd PSF SCALPELS principal component

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

C  Declarations for 'cheopsf2.dat'	! CHEOPS visit 2 PSF SCALPELS data

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

      real*8        BJD_7       ! (d) Barycentric Julian date
      real*8        Flux_7      ! Normalised flux
      real*8        e_Flux_7    ! Normalised flux error
      real*8        Fluxd_1     ! Normalised detrended flux
      real*8        e_Fluxd_1   ! Normalised detrended flux error
      real*8        Fluxpsfmodel_1 ! PSF SCALPELS PSF shape change model
      real*4        Theta0_1    ! 1st PSF SCALPELS principal component
      real*4        Theta1_1    ! 2nd PSF SCALPELS principal component
      real*8        Theta2_1    ! 3rd PSF SCALPELS principal component
      real*4        Theta3_1    ! 4th PSF SCALPELS principal component
      real*8        Theta4_1    ! 5th PSF SCALPELS principal component
      real*4        Theta5_1    ! 6th PSF SCALPELS principal component
      real*4        Theta6_1    ! 7th PSF SCALPELS principal component
      real*4        Theta7_1    ! 8th PSF SCALPELS principal component
      real*4        Theta8_1    ! 9th PSF SCALPELS principal component
      real*4        Theta9_1    ! 10th PSF SCALPELS principal component
      real*4        Theta10_1   ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_1   ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_1   ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_1   ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_1   ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_1   ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_1   ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_1   ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_1   ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_1   ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_1   ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_1   ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_1   ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_1   ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_1   ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_1   ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_1   ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_1   ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_1   ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_1   ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_1   ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_1   ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_1   ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_1   ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_1   ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_1   ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_1   ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_1   ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_1   ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_1   ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_1   ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_1   ! ? 42nd PSF SCALPELS principal component

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

C  Declarations for 'cheopsf3.dat'	! CHEOPS visit 3 PSF SCALPELS data

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

      real*8        BJD_8       ! (d) Barycentric Julian date
      real*8        Flux_8      ! Normalised flux
      real*8        e_Flux_8    ! Normalised flux error
      real*8        Fluxd_2     ! Normalised detrended flux
      real*8        e_Fluxd_2   ! Normalised detrended flux error
      real*8        Fluxpsfmodel_2 ! PSF SCALPELS PSF shape change model
      real*4        Theta0_2    ! 1st PSF SCALPELS principal component
      real*4        Theta1_2    ! 2nd PSF SCALPELS principal component
      real*8        Theta2_2    ! 3rd PSF SCALPELS principal component
      real*4        Theta3_2    ! 4th PSF SCALPELS principal component
      real*8        Theta4_2    ! 5th PSF SCALPELS principal component
      real*4        Theta5_2    ! 6th PSF SCALPELS principal component
      real*4        Theta6_2    ! 7th PSF SCALPELS principal component
      real*4        Theta7_2    ! 8th PSF SCALPELS principal component
      real*4        Theta8_2    ! 9th PSF SCALPELS principal component
      real*4        Theta9_2    ! 10th PSF SCALPELS principal component
      real*4        Theta10_2   ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_2   ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_2   ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_2   ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_2   ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_2   ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_2   ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_2   ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_2   ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_2   ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_2   ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_2   ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_2   ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_2   ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_2   ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_2   ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_2   ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_2   ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_2   ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_2   ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_2   ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_2   ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_2   ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_2   ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_2   ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_2   ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_2   ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_2   ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_2   ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_2   ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_2   ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_2   ! ? 42nd PSF SCALPELS principal component

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

C  Declarations for 'cheopsf4.dat'	! CHEOPS visit 4 PSF SCALPELS data

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

      real*8        BJD_9       ! (d) Barycentric Julian date
      real*8        Flux_9      ! Normalised flux
      real*8        e_Flux_9    ! Normalised flux error
      real*8        Fluxd_3     ! Normalised detrended flux
      real*8        e_Fluxd_3   ! Normalised detrended flux error
      real*8        Fluxpsfmodel_3 ! PSF SCALPELS PSF shape change model
      real*4        Theta0_3    ! 1st PSF SCALPELS principal component
      real*4        Theta1_3    ! 2nd PSF SCALPELS principal component
      real*8        Theta2_3    ! 3rd PSF SCALPELS principal component
      real*4        Theta3_3    ! 4th PSF SCALPELS principal component
      real*8        Theta4_3    ! 5th PSF SCALPELS principal component
      real*4        Theta5_3    ! 6th PSF SCALPELS principal component
      real*4        Theta6_3    ! 7th PSF SCALPELS principal component
      real*4        Theta7_3    ! 8th PSF SCALPELS principal component
      real*4        Theta8_3    ! 9th PSF SCALPELS principal component
      real*4        Theta9_3    ! 10th PSF SCALPELS principal component
      real*4        Theta10_3   ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_3   ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_3   ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_3   ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_3   ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_3   ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_3   ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_3   ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_3   ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_3   ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_3   ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_3   ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_3   ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_3   ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_3   ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_3   ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_3   ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_3   ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_3   ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_3   ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_3   ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_3   ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_3   ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_3   ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_3   ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_3   ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_3   ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_3   ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_3   ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_3   ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_3   ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_3   ! ? 42nd PSF SCALPELS principal component

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

C  Declarations for 'cheopsf5.dat'	! CHEOPS visit 5 PSF SCALPELS data

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

      real*8        BJD_10      ! (d) Barycentric Julian date
      real*8        Flux_10     ! Normalised flux
      real*8        e_Flux_10   ! Normalised flux error
      real*8        Fluxd_4     ! Normalised detrended flux
      real*8        e_Fluxd_4   ! Normalised detrended flux error
      real*8        Fluxpsfmodel_4 ! PSF SCALPELS PSF shape change model
      real*4        Theta0_4    ! 1st PSF SCALPELS principal component
      real*4        Theta1_4    ! 2nd PSF SCALPELS principal component
      real*8        Theta2_4    ! 3rd PSF SCALPELS principal component
      real*4        Theta3_4    ! 4th PSF SCALPELS principal component
      real*8        Theta4_4    ! 5th PSF SCALPELS principal component
      real*4        Theta5_4    ! 6th PSF SCALPELS principal component
      real*4        Theta6_4    ! 7th PSF SCALPELS principal component
      real*4        Theta7_4    ! 8th PSF SCALPELS principal component
      real*4        Theta8_4    ! 9th PSF SCALPELS principal component
      real*4        Theta9_4    ! 10th PSF SCALPELS principal component
      real*4        Theta10_4   ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_4   ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_4   ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_4   ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_4   ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_4   ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_4   ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_4   ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_4   ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_4   ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_4   ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_4   ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_4   ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_4   ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_4   ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_4   ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_4   ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_4   ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_4   ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_4   ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_4   ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_4   ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_4   ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_4   ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_4   ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_4   ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_4   ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_4   ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_4   ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_4   ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_4   ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_4   ! ? 42nd PSF SCALPELS principal component

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

C  Declarations for 'cheopsf6.dat'	! CHEOPS visit 6 PSF SCALPELS data

      integer*4 nr__11
      parameter (nr__11=762)	! Number of records
      character*1001 ar__11 	! Full-size record

      real*8        BJD_11      ! (d) Barycentric Julian date
      real*8        Flux_11     ! Normalised flux
      real*8        e_Flux_11   ! Normalised flux error
      real*8        Fluxd_5     ! Normalised detrended flux
      real*8        e_Fluxd_5   ! Normalised detrended flux error
      real*8        Fluxpsfmodel_5 ! PSF SCALPELS PSF shape change model
      real*4        Theta0_5    ! 1st PSF SCALPELS principal component
      real*4        Theta1_5    ! 2nd PSF SCALPELS principal component
      real*8        Theta2_5    ! 3rd PSF SCALPELS principal component
      real*4        Theta3_5    ! 4th PSF SCALPELS principal component
      real*8        Theta4_5    ! 5th PSF SCALPELS principal component
      real*4        Theta5_5    ! 6th PSF SCALPELS principal component
      real*4        Theta6_5    ! 7th PSF SCALPELS principal component
      real*4        Theta7_5    ! 8th PSF SCALPELS principal component
      real*4        Theta8_5    ! 9th PSF SCALPELS principal component
      real*4        Theta9_5    ! 10th PSF SCALPELS principal component
      real*4        Theta10_5   ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_5   ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_5   ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_5   ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_5   ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_5   ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_5   ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_5   ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_5   ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_5   ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_5   ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_5   ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_5   ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_5   ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_5   ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_5   ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_5   ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_5   ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_5   ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_5   ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_5   ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_5   ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_5   ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_5   ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_5   ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_5   ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_5   ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_5   ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_5   ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_5   ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_5   ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_5   ! ? 42nd PSF SCALPELS principal component

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

C  Declarations for 'cheodet1.dat'	! CHEOPS visit 1 detrended data

      integer*4 nr__12
      parameter (nr__12=440)	! Number of records
      character*59 ar__12 	! Full-size record

      real*8        BJD_12      ! (d) Barycentric Julian date
      real*8        Flux_12     ! Normalised detrended flux
      real*8        e_Flux_12   ! Normalised flux error

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

C  Declarations for 'cheodet2.dat'	! CHEOPS visit 2 detrended data

      integer*4 nr__13
      parameter (nr__13=278)	! Number of records
      character*59 ar__13 	! Full-size record

      real*8        BJD_13      ! (d) Barycentric Julian date
      real*8        Flux_13     ! Normalised detrended flux
      real*8        e_Flux_13   ! Normalised flux error

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

C  Declarations for 'cheodet3.dat'	! CHEOPS visit 3 detrended data

      integer*4 nr__14
      parameter (nr__14=282)	! Number of records
      character*59 ar__14 	! Full-size record

      real*8        BJD_14      ! (d) Barycentric Julian date
      real*8        Flux_14     ! Normalised detrended flux
      real*8        e_Flux_14   ! Normalised flux error

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

C  Declarations for 'cheodet4.dat'	! CHEOPS visit 4 detrended data

      integer*4 nr__15
      parameter (nr__15=254)	! Number of records
      character*59 ar__15 	! Full-size record

      real*8        BJD_15      ! (d) Barycentric Julian date
      real*8        Flux_15     ! Normalised detrended flux
      real*8        e_Flux_15   ! Normalised flux error

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

C  Declarations for 'cheodet5.dat'	! CHEOPS visit 5 detrended data

      integer*4 nr__16
      parameter (nr__16=231)	! Number of records
      character*59 ar__16 	! Full-size record

      real*8        BJD_16      ! (d) Barycentric Julian date
      real*8        Flux_16     ! Normalised detrended flux
      real*8        e_Flux_16   ! Normalised flux error

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

C  Declarations for 'cheodet6.dat'	! CHEOPS visit 6 detrended data

      integer*4 nr__17
      parameter (nr__17=762)	! Number of records
      character*59 ar__17 	! Full-size record

      real*8        BJD_17      ! (d) Barycentric Julian date
      real*8        Flux_17     ! Normalised detrended flux
      real*8        e_Flux_17   ! Normalised flux error

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

C  Declarations for 'traws13.dat'	! TESS sector 13 raw data

      integer*4 nr__18
      parameter (nr__18=17939)	! Number of records
      character*59 ar__18 	! Full-size record

      real*8        BJD_18      ! (d) Barycentric Julian date
      real*8        Flux_18     ! Normalised flux
      real*8        e_Flux_18   ! Normalised flux error

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

C  Declarations for 'traws27.dat'	! TESS sector 27 raw data

      integer*4 nr__19
      parameter (nr__19=71703)	! Number of records
      character*59 ar__19 	! Full-size record

      real*8        BJD_19      ! (d) Barycentric Julian date
      real*8        Flux_19     ! Normalised flux
      real*8        e_Flux_19   ! Normalised flux error

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

C  Declarations for 'tdets13.dat'	! TESS sector 13 detrended data

      integer*4 nr__20
      parameter (nr__20=17939)	! Number of records
      character*59 ar__20 	! Full-size record

      real*8        BJD_20      ! (d) Barycentric Julian date
      real*8        Flux_20     ! Normalised detrended flux
      real*8        e_Flux_20   ! Normalised flux error

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

C  Declarations for 'tdets27.dat'	! TESS sector 27 detrended data

      integer*4 nr__21
      parameter (nr__21=71703)	! Number of records
      character*59 ar__21 	! Full-size record

      real*8        BJD_21      ! (d) Barycentric Julian date
      real*8        Flux_21     ! Normalised detrended flux
      real*8        e_Flux_21   ! Normalised flux error

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

C  Declarations for 'harpsraw.dat'	! HARPS raw data

      integer*4 nr__22
      parameter (nr__22=26)	! Number of records
      character*98 ar__22 	! Full-size record

      real*8        BJD_22      ! (d) Barycentric Julian date
      real*8        RV          ! (m/s) Radial velocity
      real*8        e_RV        ! (m/s) Radial velocity error
      real*4        BIS         ! (m/s) Bisector span of the CCF
      real*4        FWHM        ! (m/s) FWHM of the CCF
      real*8        Contrast    ! Contrast of the CCF
      real*4        Sindex      ! S index
      real*4        Halpha      ! H{alpha} activity index
      real*4        NaD1        ! Na D1 activity index
      real*4        NaD2        ! Na D2 activity index

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

C  Declarations for 'harpsccf.dat'	! HARPS CCF SCALPELS data

      integer*4 nr__23
      parameter (nr__23=26)	! Number of records
      character*107 ar__23 	! Full-size record

      real*8        BJD_23      ! (d) Barycentric Julian date
      real*8        RV_1        ! (m/s) Radial velocity
      real*8        e_RV_1      ! (m/s) Radial velocity error
      real*8        Theta0_6    ! 1st CCF SCALPELS principal component
      real*8        Theta1_6    ! 2nd CCF SCALPELS principal component
      real*8        Theta2_6    ! 3rd CCF SCALPELS principal component

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

C  Declarations for 'ctiozs.dat'	! LCOGT CTIO Zs data

      integer*4 nr__24
      parameter (nr__24=255)	! Number of records
      character*59 ar__24 	! Full-size record

      real*8        BJD_24      ! (d) Barycentric Julian date
      real*8        Flux_22     ! Normalised flux
      real*8        e_Flux_22   ! Normalised flux error

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

C  Declarations for 'saaob.dat'	! LCOGT SAAO B data

      integer*4 nr__25
      parameter (nr__25=122)	! Number of records
      character*59 ar__25 	! Full-size record

      real*8        BJD_25      ! (d) Barycentric Julian date
      real*8        Flux_23     ! Normalised flux
      real*8        e_Flux_23   ! Normalised flux error

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

C  Declarations for 'saaozs1.dat'	! LCOGT SAAO Zs visit 1 data

      integer*4 nr__26
      parameter (nr__26=453)	! Number of records
      character*59 ar__26 	! Full-size record

      real*8        BJD_26      ! (d) Barycentric Julian date
      real*8        Flux_24     ! Normalised flux
      real*8        e_Flux_24   ! Normalised flux error

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

C  Declarations for 'saaozs2.dat'	! LCOGT SAAO Zs visit 2 data

      integer*4 nr__27
      parameter (nr__27=368)	! Number of records
      character*59 ar__27 	! Full-size record

      real*8        BJD_27      ! (d) Barycentric Julian date
      real*8        Flux_25     ! Normalised flux
      real*8        e_Flux_25   ! Normalised flux error

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

C  Declarations for 'saaozs3.dat'	! LCOGT SAAO Zs visit 4 data

      integer*4 nr__28
      parameter (nr__28=343)	! Number of records
      character*59 ar__28 	! Full-size record

      real*8        BJD_28      ! (d) Barycentric Julian date
      real*8        Flux_26     ! Normalised flux
      real*8        e_Flux_26   ! Normalised flux error

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

C  Declarations for 'ngtsraw.dat'	! NGTS raw data

      integer*4 nr__29
      parameter (nr__29=1811)	! Number of records
      character*75 ar__29 	! Full-size record

      real*8        BJD_29      ! (d) Barycentric Julian date
      real*8        Flux_27     ! Normalised flux
      real*8        e_Flux_27   ! Normalised flux error
      real*8        Airmass     ! Airmass
      real*8        SkyBkg      ! (ct) Sky background
      real*8        X           ! (pix) X centroid position
      real*8        Y           ! (pix) Y centroid position

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

C  Declarations for 'ngtsdet.dat'	! NGTS detrended data

      integer*4 nr__30
      parameter (nr__30=1811)	! Number of records
      character*59 ar__30 	! Full-size record

      real*8        BJD_30      ! (d) Barycentric Julian date
      real*8        Flux_28     ! Normalised detrended flux
      real*8        e_Flux_28   ! Normalised flux error

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

C  Declarations for 'astep1.dat'	! ASTEP visit 1 data

      integer*4 nr__31
      parameter (nr__31=297)	! Number of records
      character*59 ar__31 	! Full-size record

      real*8        BJD_31      ! (d) Barycentric Julian date
      real*8        Flux_29     ! Normalised flux
      real*8        e_Flux_29   ! Normalised flux error

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

C  Declarations for 'astep2.dat'	! ASTEP visit 2 data

      integer*4 nr__32
      parameter (nr__32=421)	! Number of records
      character*59 ar__32 	! Full-size record

      real*8        BJD_32      ! (d) Barycentric Julian date
      real*8        Flux_30     ! Normalised flux
      real*8        e_Flux_30   ! Normalised flux error

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

C  Declarations for 'astep3.dat'	! ASTEP visit 3 data

      integer*4 nr__33
      parameter (nr__33=267)	! Number of records
      character*59 ar__33 	! Full-size record

      real*8        BJD_33      ! (d) Barycentric Julian date
      real*8        Flux_31     ! Normalised flux
      real*8        e_Flux_31   ! Normalised flux error

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

C  Declarations for 'astep4.dat'	! ASTEP visit 3 data

      integer*4 nr__34
      parameter (nr__34=178)	! Number of records
      character*59 ar__34 	! Full-size record

      real*8        BJD_34      ! (d) Barycentric Julian date
      real*8        Flux_32     ! Normalised flux
      real*8        e_Flux_32   ! Normalised flux error

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

C  Declarations for 'asassn.dat'	! ASAS-SN data

      integer*4 nr__35
      parameter (nr__35=401)	! Number of records
      character*59 ar__35 	! Full-size record

      real*8        BJD_35      ! (d) Barycentric Julian date
      real*8        Flux_33     ! Normalised flux
      real*8        e_Flux_33   ! Normalised flux error

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

C  Declarations for 'wasp.dat'	! WASP data

      integer*4 nr__36
      parameter (nr__36=66832)	! Number of records
      character*59 ar__36 	! Full-size record

      real*8        BJD_36      ! (d) Barycentric Julian date
      real*8        Flux_34     ! Normalised flux
      real*8        e_Flux_34   ! Normalised flux error

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

C  Declarations for 'gems562.dat'	! Gemini-South 562 data

      integer*4 nr__37
      parameter (nr__37=119)	! Number of records
      character*19 ar__37 	! Full-size record

      integer*4     lambda      ! (nm) Wavelength (562 or 832)
      real*8        Sep         ! (arcsec) Separation
      real*8        Deltamag    ! Delta magnitude

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

C  Declarations for 'gems832.dat'	! Gemini-South 832 data

      integer*4 nr__38
      parameter (nr__38=120)	! Number of records
      character*19 ar__38 	! Full-size record

      integer*4     lambda_1    ! (nm) Wavelength (562 or 832)
      real*8        Sep_1       ! (arcsec) Separation
      real*8        Deltamag_1  ! Delta magnitude

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

C  Loading file 'cheoraw1.dat'	! CHEOPS visit 1 raw data

C  Format for file interpretation

    1 format(
     +  F18.10,1X,F18.16,1X,F22.20,1X,E14.7,1X,F12.9,1X,F19.15,1X,
     +  F20.18,1X,F20.18,1X,E22.16,1X,F13.10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheoraw1.dat')
      write(6,*) '....Loading file: cheoraw1.dat'
      do i__=1,440
        read(1,'(A187)')ar__
        read(ar__,1)
     +  BJD,Flux,e_Flux,Xoffset,Yoffset,Roll,Contam,Bckg,Smear,DeltaT
c    ..............Just test output...........
        write(6,1)
     +  BJD,Flux,e_Flux,Xoffset,Yoffset,Roll,Contam,Bckg,Smear,DeltaT
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheoraw2.dat'	! CHEOPS visit 2 raw data

C  Format for file interpretation

    2 format(
     +  F18.10,1X,F18.16,1X,F22.20,1X,E14.7,1X,F12.9,1X,F19.15,1X,
     +  F20.18,1X,F20.18,1X,E22.16,1X,F13.10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheoraw2.dat')
      write(6,*) '....Loading file: cheoraw2.dat'
      do i__=1,278
        read(1,'(A187)')ar__1
        read(ar__1,2)
     +  BJD_1,Flux_1,e_Flux_1,Xoffset_1,Yoffset_1,Roll_1,Contam_1,
     +  Bckg_1,Smear_1,DeltaT_1
c    ..............Just test output...........
        write(6,2)
     +  BJD_1,Flux_1,e_Flux_1,Xoffset_1,Yoffset_1,Roll_1,Contam_1,
     +  Bckg_1,Smear_1,DeltaT_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheoraw3.dat'	! CHEOPS visit 3 raw data

C  Format for file interpretation

    3 format(
     +  F18.10,1X,F18.16,1X,F22.20,1X,E14.7,1X,F12.9,1X,F19.15,1X,
     +  F20.18,1X,F20.18,1X,E22.16,1X,F13.10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheoraw3.dat')
      write(6,*) '....Loading file: cheoraw3.dat'
      do i__=1,282
        read(1,'(A187)')ar__2
        read(ar__2,3)
     +  BJD_2,Flux_2,e_Flux_2,Xoffset_2,Yoffset_2,Roll_2,Contam_2,
     +  Bckg_2,Smear_2,DeltaT_2
c    ..............Just test output...........
        write(6,3)
     +  BJD_2,Flux_2,e_Flux_2,Xoffset_2,Yoffset_2,Roll_2,Contam_2,
     +  Bckg_2,Smear_2,DeltaT_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheoraw4.dat'	! CHEOPS visit 4 raw data

C  Format for file interpretation

    4 format(
     +  F18.10,1X,F18.16,1X,F22.20,1X,E14.7,1X,F12.9,1X,F19.15,1X,
     +  F20.18,1X,F20.18,1X,E22.16,1X,F13.10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheoraw4.dat')
      write(6,*) '....Loading file: cheoraw4.dat'
      do i__=1,254
        read(1,'(A187)')ar__3
        read(ar__3,4)
     +  BJD_3,Flux_3,e_Flux_3,Xoffset_3,Yoffset_3,Roll_3,Contam_3,
     +  Bckg_3,Smear_3,DeltaT_3
c    ..............Just test output...........
        write(6,4)
     +  BJD_3,Flux_3,e_Flux_3,Xoffset_3,Yoffset_3,Roll_3,Contam_3,
     +  Bckg_3,Smear_3,DeltaT_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheoraw5.dat'	! CHEOPS visit 5 raw data

C  Format for file interpretation

    5 format(
     +  F18.10,1X,F18.16,1X,F22.20,1X,E14.7,1X,F12.9,1X,F19.15,1X,
     +  F20.18,1X,F20.18,1X,E22.16,1X,F13.10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheoraw5.dat')
      write(6,*) '....Loading file: cheoraw5.dat'
      do i__=1,231
        read(1,'(A187)')ar__4
        read(ar__4,5)
     +  BJD_4,Flux_4,e_Flux_4,Xoffset_4,Yoffset_4,Roll_4,Contam_4,
     +  Bckg_4,Smear_4,DeltaT_4
c    ..............Just test output...........
        write(6,5)
     +  BJD_4,Flux_4,e_Flux_4,Xoffset_4,Yoffset_4,Roll_4,Contam_4,
     +  Bckg_4,Smear_4,DeltaT_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheoraw6.dat'	! CHEOPS visit 6 raw data

C  Format for file interpretation

    6 format(
     +  F18.10,1X,F18.16,1X,F22.20,1X,E14.7,1X,F12.9,1X,F19.15,1X,
     +  F20.18,1X,F20.18,1X,E22.16,1X,F13.10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheoraw6.dat')
      write(6,*) '....Loading file: cheoraw6.dat'
      do i__=1,762
        read(1,'(A187)')ar__5
        read(ar__5,6)
     +  BJD_5,Flux_5,e_Flux_5,Xoffset_5,Yoffset_5,Roll_5,Contam_5,
     +  Bckg_5,Smear_5,DeltaT_5
c    ..............Just test output...........
        write(6,6)
     +  BJD_5,Flux_5,e_Flux_5,Xoffset_5,Yoffset_5,Roll_5,Contam_5,
     +  Bckg_5,Smear_5,DeltaT_5
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheopsf1.dat'	! CHEOPS visit 1 PSF SCALPELS data

C  Format for file interpretation

    7 format(
     +  F16.8,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,
     +  E21.13,1X,E21.12,F18.15,1X,E21.13,1X,F18.15,1X,E21.13,1X,
     +  E21.13,1X,E20.12,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.12,
     +  E21.13,1X,E21.13,1X,E21.13,1X,E21.13,E21.13,1X,E21.13,1X,
     +  E20.12,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,E20.12,1X,E20.12,1X,E21.13,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,F18.15,1X,E21.13,1X,F18.15,1X,E20.12,1X,
     +  F18.15,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheopsf1.dat')
      write(6,*) '....Loading file: cheopsf1.dat'
      do i__=1,440
        read(1,'(A1001)')ar__6
        read(ar__6,7)
     +  BJD_6,Flux_6,e_Flux_6,Fluxd,e_Fluxd,Fluxpsfmodel,Theta0,
     +  Theta1,Theta2,Theta3,Theta4,Theta5,Theta6,Theta7,Theta8,
     +  Theta9,Theta10,Theta11,Theta12,Theta13,Theta14,Theta15,
     +  Theta16,Theta17,Theta18,Theta19,Theta20,Theta21,Theta22,
     +  Theta23,Theta24,Theta25,Theta26,Theta27,Theta28,Theta29,
     +  Theta30,Theta31,Theta32,Theta33,Theta34,Theta35,Theta36,
     +  Theta37,Theta38,Theta39,Theta40,Theta41
        if(ar__6(320:339) .EQ. '') Theta10 = rNULL__
        if(ar__6(341:361) .EQ. '') Theta11 = rNULL__
        if(ar__6(362:382) .EQ. '') Theta12 = rNULL__
        if(ar__6(384:404) .EQ. '') Theta13 = rNULL__
        if(ar__6(406:426) .EQ. '') Theta14 = rNULL__
        if(ar__6(428:448) .EQ. '') Theta15 = rNULL__
        if(ar__6(449:469) .EQ. '') Theta16 = rNULL__
        if(ar__6(471:491) .EQ. '') Theta17 = rNULL__
        if(ar__6(493:512) .EQ. '') Theta18 = rNULL__
        if(ar__6(514:534) .EQ. '') Theta19 = rNULL__
        if(ar__6(536:556) .EQ. '') Theta20 = rNULL__
        if(ar__6(558:578) .EQ. '') Theta21 = rNULL__
        if(ar__6(580:599) .EQ. '') Theta22 = rNULL__
        if(ar__6(601:621) .EQ. '') Theta23 = rNULL__
        if(ar__6(623:643) .EQ. '') Theta24 = rNULL__
        if(ar__6(645:665) .EQ. '') Theta25 = rNULL__
        if(ar__6(667:686) .EQ. '') Theta26 = rNULL__
        if(ar__6(688:707) .EQ. '') Theta27 = rNULL__
        if(ar__6(709:729) .EQ. '') Theta28 = rNULL__
        if(ar__6(731:751) .EQ. '') Theta29 = rNULL__
        if(ar__6(753:773) .EQ. '') Theta30 = rNULL__
        if(ar__6(775:795) .EQ. '') Theta31 = rNULL__
        if(ar__6(797:814) .EQ. '') Theta32 = rNULL__
        if(ar__6(816:836) .EQ. '') Theta33 = rNULL__
        if(ar__6(838:855) .EQ. '') Theta34 = rNULL__
        if(ar__6(857:876) .EQ. '') Theta35 = rNULL__
        if(ar__6(878:895) .EQ. '') Theta36 = rNULL__
        if(ar__6(897:917) .EQ. '') Theta37 = rNULL__
        if(ar__6(919:939) .EQ. '') Theta38 = rNULL__
        if(ar__6(941:961) .EQ. '') Theta39 = rNULL__
        if(ar__6(963:982) .EQ. '') Theta40 = rNULL__
        if(ar__6(984:1001) .EQ. '') Theta41 = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  BJD_6,Flux_6,e_Flux_6,Fluxd,e_Fluxd,Fluxpsfmodel,Theta0,
     +  Theta1,Theta2,Theta3,Theta4,Theta5,Theta6,Theta7,Theta8,
     +  Theta9,Theta10,Theta11,Theta12,Theta13,Theta14,Theta15,
     +  Theta16,Theta17,Theta18,Theta19,Theta20,Theta21,Theta22,
     +  Theta23,Theta24,Theta25,Theta26,Theta27,Theta28,Theta29,
     +  Theta30,Theta31,Theta32,Theta33,Theta34,Theta35,Theta36,
     +  Theta37,Theta38,Theta39,Theta40,Theta41
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheopsf2.dat'	! CHEOPS visit 2 PSF SCALPELS data

C  Format for file interpretation

    8 format(
     +  F16.8,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,
     +  E21.13,1X,E21.12,F18.15,1X,E21.13,1X,F18.15,1X,E21.13,1X,
     +  E21.13,1X,E20.12,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.12,
     +  E21.13,1X,E21.13,1X,E21.13,1X,E21.13,E21.13,1X,E21.13,1X,
     +  E20.12,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,E20.12,1X,E20.12,1X,E21.13,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,F18.15,1X,E21.13,1X,F18.15,1X,E20.12,1X,
     +  F18.15,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheopsf2.dat')
      write(6,*) '....Loading file: cheopsf2.dat'
      do i__=1,278
        read(1,'(A1001)')ar__7
        read(ar__7,8)
     +  BJD_7,Flux_7,e_Flux_7,Fluxd_1,e_Fluxd_1,Fluxpsfmodel_1,
     +  Theta0_1,Theta1_1,Theta2_1,Theta3_1,Theta4_1,Theta5_1,
     +  Theta6_1,Theta7_1,Theta8_1,Theta9_1,Theta10_1,Theta11_1,
     +  Theta12_1,Theta13_1,Theta14_1,Theta15_1,Theta16_1,Theta17_1,
     +  Theta18_1,Theta19_1,Theta20_1,Theta21_1,Theta22_1,Theta23_1,
     +  Theta24_1,Theta25_1,Theta26_1,Theta27_1,Theta28_1,Theta29_1,
     +  Theta30_1,Theta31_1,Theta32_1,Theta33_1,Theta34_1,Theta35_1,
     +  Theta36_1,Theta37_1,Theta38_1,Theta39_1,Theta40_1,Theta41_1
        if(ar__7(320:339) .EQ. '') Theta10_1 = rNULL__
        if(ar__7(341:361) .EQ. '') Theta11_1 = rNULL__
        if(ar__7(362:382) .EQ. '') Theta12_1 = rNULL__
        if(ar__7(384:404) .EQ. '') Theta13_1 = rNULL__
        if(ar__7(406:426) .EQ. '') Theta14_1 = rNULL__
        if(ar__7(428:448) .EQ. '') Theta15_1 = rNULL__
        if(ar__7(449:469) .EQ. '') Theta16_1 = rNULL__
        if(ar__7(471:491) .EQ. '') Theta17_1 = rNULL__
        if(ar__7(493:512) .EQ. '') Theta18_1 = rNULL__
        if(ar__7(514:534) .EQ. '') Theta19_1 = rNULL__
        if(ar__7(536:556) .EQ. '') Theta20_1 = rNULL__
        if(ar__7(558:578) .EQ. '') Theta21_1 = rNULL__
        if(ar__7(580:599) .EQ. '') Theta22_1 = rNULL__
        if(ar__7(601:621) .EQ. '') Theta23_1 = rNULL__
        if(ar__7(623:643) .EQ. '') Theta24_1 = rNULL__
        if(ar__7(645:665) .EQ. '') Theta25_1 = rNULL__
        if(ar__7(667:686) .EQ. '') Theta26_1 = rNULL__
        if(ar__7(688:707) .EQ. '') Theta27_1 = rNULL__
        if(ar__7(709:729) .EQ. '') Theta28_1 = rNULL__
        if(ar__7(731:751) .EQ. '') Theta29_1 = rNULL__
        if(ar__7(753:773) .EQ. '') Theta30_1 = rNULL__
        if(ar__7(775:795) .EQ. '') Theta31_1 = rNULL__
        if(ar__7(797:814) .EQ. '') Theta32_1 = rNULL__
        if(ar__7(816:836) .EQ. '') Theta33_1 = rNULL__
        if(ar__7(838:855) .EQ. '') Theta34_1 = rNULL__
        if(ar__7(857:876) .EQ. '') Theta35_1 = rNULL__
        if(ar__7(878:895) .EQ. '') Theta36_1 = rNULL__
        if(ar__7(897:917) .EQ. '') Theta37_1 = rNULL__
        if(ar__7(919:939) .EQ. '') Theta38_1 = rNULL__
        if(ar__7(941:961) .EQ. '') Theta39_1 = rNULL__
        if(ar__7(963:982) .EQ. '') Theta40_1 = rNULL__
        if(ar__7(984:1001) .EQ. '') Theta41_1 = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  BJD_7,Flux_7,e_Flux_7,Fluxd_1,e_Fluxd_1,Fluxpsfmodel_1,
     +  Theta0_1,Theta1_1,Theta2_1,Theta3_1,Theta4_1,Theta5_1,
     +  Theta6_1,Theta7_1,Theta8_1,Theta9_1,Theta10_1,Theta11_1,
     +  Theta12_1,Theta13_1,Theta14_1,Theta15_1,Theta16_1,Theta17_1,
     +  Theta18_1,Theta19_1,Theta20_1,Theta21_1,Theta22_1,Theta23_1,
     +  Theta24_1,Theta25_1,Theta26_1,Theta27_1,Theta28_1,Theta29_1,
     +  Theta30_1,Theta31_1,Theta32_1,Theta33_1,Theta34_1,Theta35_1,
     +  Theta36_1,Theta37_1,Theta38_1,Theta39_1,Theta40_1,Theta41_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheopsf3.dat'	! CHEOPS visit 3 PSF SCALPELS data

C  Format for file interpretation

    9 format(
     +  F16.8,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,
     +  E21.13,1X,E21.12,F18.15,1X,E21.13,1X,F18.15,1X,E21.13,1X,
     +  E21.13,1X,E20.12,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.12,
     +  E21.13,1X,E21.13,1X,E21.13,1X,E21.13,E21.13,1X,E21.13,1X,
     +  E20.12,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,E20.12,1X,E20.12,1X,E21.13,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,F18.15,1X,E21.13,1X,F18.15,1X,E20.12,1X,
     +  F18.15,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheopsf3.dat')
      write(6,*) '....Loading file: cheopsf3.dat'
      do i__=1,282
        read(1,'(A1001)')ar__8
        read(ar__8,9)
     +  BJD_8,Flux_8,e_Flux_8,Fluxd_2,e_Fluxd_2,Fluxpsfmodel_2,
     +  Theta0_2,Theta1_2,Theta2_2,Theta3_2,Theta4_2,Theta5_2,
     +  Theta6_2,Theta7_2,Theta8_2,Theta9_2,Theta10_2,Theta11_2,
     +  Theta12_2,Theta13_2,Theta14_2,Theta15_2,Theta16_2,Theta17_2,
     +  Theta18_2,Theta19_2,Theta20_2,Theta21_2,Theta22_2,Theta23_2,
     +  Theta24_2,Theta25_2,Theta26_2,Theta27_2,Theta28_2,Theta29_2,
     +  Theta30_2,Theta31_2,Theta32_2,Theta33_2,Theta34_2,Theta35_2,
     +  Theta36_2,Theta37_2,Theta38_2,Theta39_2,Theta40_2,Theta41_2
        if(ar__8(320:339) .EQ. '') Theta10_2 = rNULL__
        if(ar__8(341:361) .EQ. '') Theta11_2 = rNULL__
        if(ar__8(362:382) .EQ. '') Theta12_2 = rNULL__
        if(ar__8(384:404) .EQ. '') Theta13_2 = rNULL__
        if(ar__8(406:426) .EQ. '') Theta14_2 = rNULL__
        if(ar__8(428:448) .EQ. '') Theta15_2 = rNULL__
        if(ar__8(449:469) .EQ. '') Theta16_2 = rNULL__
        if(ar__8(471:491) .EQ. '') Theta17_2 = rNULL__
        if(ar__8(493:512) .EQ. '') Theta18_2 = rNULL__
        if(ar__8(514:534) .EQ. '') Theta19_2 = rNULL__
        if(ar__8(536:556) .EQ. '') Theta20_2 = rNULL__
        if(ar__8(558:578) .EQ. '') Theta21_2 = rNULL__
        if(ar__8(580:599) .EQ. '') Theta22_2 = rNULL__
        if(ar__8(601:621) .EQ. '') Theta23_2 = rNULL__
        if(ar__8(623:643) .EQ. '') Theta24_2 = rNULL__
        if(ar__8(645:665) .EQ. '') Theta25_2 = rNULL__
        if(ar__8(667:686) .EQ. '') Theta26_2 = rNULL__
        if(ar__8(688:707) .EQ. '') Theta27_2 = rNULL__
        if(ar__8(709:729) .EQ. '') Theta28_2 = rNULL__
        if(ar__8(731:751) .EQ. '') Theta29_2 = rNULL__
        if(ar__8(753:773) .EQ. '') Theta30_2 = rNULL__
        if(ar__8(775:795) .EQ. '') Theta31_2 = rNULL__
        if(ar__8(797:814) .EQ. '') Theta32_2 = rNULL__
        if(ar__8(816:836) .EQ. '') Theta33_2 = rNULL__
        if(ar__8(838:855) .EQ. '') Theta34_2 = rNULL__
        if(ar__8(857:876) .EQ. '') Theta35_2 = rNULL__
        if(ar__8(878:895) .EQ. '') Theta36_2 = rNULL__
        if(ar__8(897:917) .EQ. '') Theta37_2 = rNULL__
        if(ar__8(919:939) .EQ. '') Theta38_2 = rNULL__
        if(ar__8(941:961) .EQ. '') Theta39_2 = rNULL__
        if(ar__8(963:982) .EQ. '') Theta40_2 = rNULL__
        if(ar__8(984:1001) .EQ. '') Theta41_2 = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  BJD_8,Flux_8,e_Flux_8,Fluxd_2,e_Fluxd_2,Fluxpsfmodel_2,
     +  Theta0_2,Theta1_2,Theta2_2,Theta3_2,Theta4_2,Theta5_2,
     +  Theta6_2,Theta7_2,Theta8_2,Theta9_2,Theta10_2,Theta11_2,
     +  Theta12_2,Theta13_2,Theta14_2,Theta15_2,Theta16_2,Theta17_2,
     +  Theta18_2,Theta19_2,Theta20_2,Theta21_2,Theta22_2,Theta23_2,
     +  Theta24_2,Theta25_2,Theta26_2,Theta27_2,Theta28_2,Theta29_2,
     +  Theta30_2,Theta31_2,Theta32_2,Theta33_2,Theta34_2,Theta35_2,
     +  Theta36_2,Theta37_2,Theta38_2,Theta39_2,Theta40_2,Theta41_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheopsf4.dat'	! CHEOPS visit 4 PSF SCALPELS data

C  Format for file interpretation

   10 format(
     +  F16.8,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,
     +  E21.13,1X,E21.12,F18.15,1X,E21.13,1X,F18.15,1X,E21.13,1X,
     +  E21.13,1X,E20.12,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.12,
     +  E21.13,1X,E21.13,1X,E21.13,1X,E21.13,E21.13,1X,E21.13,1X,
     +  E20.12,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,E20.12,1X,E20.12,1X,E21.13,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,F18.15,1X,E21.13,1X,F18.15,1X,E20.12,1X,
     +  F18.15,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheopsf4.dat')
      write(6,*) '....Loading file: cheopsf4.dat'
      do i__=1,254
        read(1,'(A1001)')ar__9
        read(ar__9,10)
     +  BJD_9,Flux_9,e_Flux_9,Fluxd_3,e_Fluxd_3,Fluxpsfmodel_3,
     +  Theta0_3,Theta1_3,Theta2_3,Theta3_3,Theta4_3,Theta5_3,
     +  Theta6_3,Theta7_3,Theta8_3,Theta9_3,Theta10_3,Theta11_3,
     +  Theta12_3,Theta13_3,Theta14_3,Theta15_3,Theta16_3,Theta17_3,
     +  Theta18_3,Theta19_3,Theta20_3,Theta21_3,Theta22_3,Theta23_3,
     +  Theta24_3,Theta25_3,Theta26_3,Theta27_3,Theta28_3,Theta29_3,
     +  Theta30_3,Theta31_3,Theta32_3,Theta33_3,Theta34_3,Theta35_3,
     +  Theta36_3,Theta37_3,Theta38_3,Theta39_3,Theta40_3,Theta41_3
        if(ar__9(320:339) .EQ. '') Theta10_3 = rNULL__
        if(ar__9(341:361) .EQ. '') Theta11_3 = rNULL__
        if(ar__9(362:382) .EQ. '') Theta12_3 = rNULL__
        if(ar__9(384:404) .EQ. '') Theta13_3 = rNULL__
        if(ar__9(406:426) .EQ. '') Theta14_3 = rNULL__
        if(ar__9(428:448) .EQ. '') Theta15_3 = rNULL__
        if(ar__9(449:469) .EQ. '') Theta16_3 = rNULL__
        if(ar__9(471:491) .EQ. '') Theta17_3 = rNULL__
        if(ar__9(493:512) .EQ. '') Theta18_3 = rNULL__
        if(ar__9(514:534) .EQ. '') Theta19_3 = rNULL__
        if(ar__9(536:556) .EQ. '') Theta20_3 = rNULL__
        if(ar__9(558:578) .EQ. '') Theta21_3 = rNULL__
        if(ar__9(580:599) .EQ. '') Theta22_3 = rNULL__
        if(ar__9(601:621) .EQ. '') Theta23_3 = rNULL__
        if(ar__9(623:643) .EQ. '') Theta24_3 = rNULL__
        if(ar__9(645:665) .EQ. '') Theta25_3 = rNULL__
        if(ar__9(667:686) .EQ. '') Theta26_3 = rNULL__
        if(ar__9(688:707) .EQ. '') Theta27_3 = rNULL__
        if(ar__9(709:729) .EQ. '') Theta28_3 = rNULL__
        if(ar__9(731:751) .EQ. '') Theta29_3 = rNULL__
        if(ar__9(753:773) .EQ. '') Theta30_3 = rNULL__
        if(ar__9(775:795) .EQ. '') Theta31_3 = rNULL__
        if(ar__9(797:814) .EQ. '') Theta32_3 = rNULL__
        if(ar__9(816:836) .EQ. '') Theta33_3 = rNULL__
        if(ar__9(838:855) .EQ. '') Theta34_3 = rNULL__
        if(ar__9(857:876) .EQ. '') Theta35_3 = rNULL__
        if(ar__9(878:895) .EQ. '') Theta36_3 = rNULL__
        if(ar__9(897:917) .EQ. '') Theta37_3 = rNULL__
        if(ar__9(919:939) .EQ. '') Theta38_3 = rNULL__
        if(ar__9(941:961) .EQ. '') Theta39_3 = rNULL__
        if(ar__9(963:982) .EQ. '') Theta40_3 = rNULL__
        if(ar__9(984:1001) .EQ. '') Theta41_3 = rNULL__
c    ..............Just test output...........
        write(6,10)
     +  BJD_9,Flux_9,e_Flux_9,Fluxd_3,e_Fluxd_3,Fluxpsfmodel_3,
     +  Theta0_3,Theta1_3,Theta2_3,Theta3_3,Theta4_3,Theta5_3,
     +  Theta6_3,Theta7_3,Theta8_3,Theta9_3,Theta10_3,Theta11_3,
     +  Theta12_3,Theta13_3,Theta14_3,Theta15_3,Theta16_3,Theta17_3,
     +  Theta18_3,Theta19_3,Theta20_3,Theta21_3,Theta22_3,Theta23_3,
     +  Theta24_3,Theta25_3,Theta26_3,Theta27_3,Theta28_3,Theta29_3,
     +  Theta30_3,Theta31_3,Theta32_3,Theta33_3,Theta34_3,Theta35_3,
     +  Theta36_3,Theta37_3,Theta38_3,Theta39_3,Theta40_3,Theta41_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheopsf5.dat'	! CHEOPS visit 5 PSF SCALPELS data

C  Format for file interpretation

   11 format(
     +  F16.8,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,
     +  E21.13,1X,E21.12,F18.15,1X,E21.13,1X,F18.15,1X,E21.13,1X,
     +  E21.13,1X,E20.12,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.12,
     +  E21.13,1X,E21.13,1X,E21.13,1X,E21.13,E21.13,1X,E21.13,1X,
     +  E20.12,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,E20.12,1X,E20.12,1X,E21.13,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,F18.15,1X,E21.13,1X,F18.15,1X,E20.12,1X,
     +  F18.15,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheopsf5.dat')
      write(6,*) '....Loading file: cheopsf5.dat'
      do i__=1,231
        read(1,'(A1001)')ar__10
        read(ar__10,11)
     +  BJD_10,Flux_10,e_Flux_10,Fluxd_4,e_Fluxd_4,Fluxpsfmodel_4,
     +  Theta0_4,Theta1_4,Theta2_4,Theta3_4,Theta4_4,Theta5_4,
     +  Theta6_4,Theta7_4,Theta8_4,Theta9_4,Theta10_4,Theta11_4,
     +  Theta12_4,Theta13_4,Theta14_4,Theta15_4,Theta16_4,Theta17_4,
     +  Theta18_4,Theta19_4,Theta20_4,Theta21_4,Theta22_4,Theta23_4,
     +  Theta24_4,Theta25_4,Theta26_4,Theta27_4,Theta28_4,Theta29_4,
     +  Theta30_4,Theta31_4,Theta32_4,Theta33_4,Theta34_4,Theta35_4,
     +  Theta36_4,Theta37_4,Theta38_4,Theta39_4,Theta40_4,Theta41_4
        if(ar__10(320:339) .EQ. '') Theta10_4 = rNULL__
        if(ar__10(341:361) .EQ. '') Theta11_4 = rNULL__
        if(ar__10(362:382) .EQ. '') Theta12_4 = rNULL__
        if(ar__10(384:404) .EQ. '') Theta13_4 = rNULL__
        if(ar__10(406:426) .EQ. '') Theta14_4 = rNULL__
        if(ar__10(428:448) .EQ. '') Theta15_4 = rNULL__
        if(ar__10(449:469) .EQ. '') Theta16_4 = rNULL__
        if(ar__10(471:491) .EQ. '') Theta17_4 = rNULL__
        if(ar__10(493:512) .EQ. '') Theta18_4 = rNULL__
        if(ar__10(514:534) .EQ. '') Theta19_4 = rNULL__
        if(ar__10(536:556) .EQ. '') Theta20_4 = rNULL__
        if(ar__10(558:578) .EQ. '') Theta21_4 = rNULL__
        if(ar__10(580:599) .EQ. '') Theta22_4 = rNULL__
        if(ar__10(601:621) .EQ. '') Theta23_4 = rNULL__
        if(ar__10(623:643) .EQ. '') Theta24_4 = rNULL__
        if(ar__10(645:665) .EQ. '') Theta25_4 = rNULL__
        if(ar__10(667:686) .EQ. '') Theta26_4 = rNULL__
        if(ar__10(688:707) .EQ. '') Theta27_4 = rNULL__
        if(ar__10(709:729) .EQ. '') Theta28_4 = rNULL__
        if(ar__10(731:751) .EQ. '') Theta29_4 = rNULL__
        if(ar__10(753:773) .EQ. '') Theta30_4 = rNULL__
        if(ar__10(775:795) .EQ. '') Theta31_4 = rNULL__
        if(ar__10(797:814) .EQ. '') Theta32_4 = rNULL__
        if(ar__10(816:836) .EQ. '') Theta33_4 = rNULL__
        if(ar__10(838:855) .EQ. '') Theta34_4 = rNULL__
        if(ar__10(857:876) .EQ. '') Theta35_4 = rNULL__
        if(ar__10(878:895) .EQ. '') Theta36_4 = rNULL__
        if(ar__10(897:917) .EQ. '') Theta37_4 = rNULL__
        if(ar__10(919:939) .EQ. '') Theta38_4 = rNULL__
        if(ar__10(941:961) .EQ. '') Theta39_4 = rNULL__
        if(ar__10(963:982) .EQ. '') Theta40_4 = rNULL__
        if(ar__10(984:1001) .EQ. '') Theta41_4 = rNULL__
c    ..............Just test output...........
        write(6,11)
     +  BJD_10,Flux_10,e_Flux_10,Fluxd_4,e_Fluxd_4,Fluxpsfmodel_4,
     +  Theta0_4,Theta1_4,Theta2_4,Theta3_4,Theta4_4,Theta5_4,
     +  Theta6_4,Theta7_4,Theta8_4,Theta9_4,Theta10_4,Theta11_4,
     +  Theta12_4,Theta13_4,Theta14_4,Theta15_4,Theta16_4,Theta17_4,
     +  Theta18_4,Theta19_4,Theta20_4,Theta21_4,Theta22_4,Theta23_4,
     +  Theta24_4,Theta25_4,Theta26_4,Theta27_4,Theta28_4,Theta29_4,
     +  Theta30_4,Theta31_4,Theta32_4,Theta33_4,Theta34_4,Theta35_4,
     +  Theta36_4,Theta37_4,Theta38_4,Theta39_4,Theta40_4,Theta41_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheopsf6.dat'	! CHEOPS visit 6 PSF SCALPELS data

C  Format for file interpretation

   12 format(
     +  F16.8,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,F17.15,1X,
     +  E21.13,1X,E21.12,F18.15,1X,E21.13,1X,F18.15,1X,E21.13,1X,
     +  E21.13,1X,E20.12,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.12,
     +  E21.13,1X,E21.13,1X,E21.13,1X,E21.13,E21.13,1X,E21.13,1X,
     +  E20.12,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,E20.12,1X,E20.12,1X,E21.13,1X,E21.13,1X,
     +  E21.13,1X,E21.13,1X,F18.15,1X,E21.13,1X,F18.15,1X,E20.12,1X,
     +  F18.15,1X,E21.13,1X,E21.13,1X,E21.13,1X,E20.12,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheopsf6.dat')
      write(6,*) '....Loading file: cheopsf6.dat'
      do i__=1,762
        read(1,'(A1001)')ar__11
        read(ar__11,12)
     +  BJD_11,Flux_11,e_Flux_11,Fluxd_5,e_Fluxd_5,Fluxpsfmodel_5,
     +  Theta0_5,Theta1_5,Theta2_5,Theta3_5,Theta4_5,Theta5_5,
     +  Theta6_5,Theta7_5,Theta8_5,Theta9_5,Theta10_5,Theta11_5,
     +  Theta12_5,Theta13_5,Theta14_5,Theta15_5,Theta16_5,Theta17_5,
     +  Theta18_5,Theta19_5,Theta20_5,Theta21_5,Theta22_5,Theta23_5,
     +  Theta24_5,Theta25_5,Theta26_5,Theta27_5,Theta28_5,Theta29_5,
     +  Theta30_5,Theta31_5,Theta32_5,Theta33_5,Theta34_5,Theta35_5,
     +  Theta36_5,Theta37_5,Theta38_5,Theta39_5,Theta40_5,Theta41_5
        if(ar__11(320:339) .EQ. '') Theta10_5 = rNULL__
        if(ar__11(341:361) .EQ. '') Theta11_5 = rNULL__
        if(ar__11(362:382) .EQ. '') Theta12_5 = rNULL__
        if(ar__11(384:404) .EQ. '') Theta13_5 = rNULL__
        if(ar__11(406:426) .EQ. '') Theta14_5 = rNULL__
        if(ar__11(428:448) .EQ. '') Theta15_5 = rNULL__
        if(ar__11(449:469) .EQ. '') Theta16_5 = rNULL__
        if(ar__11(471:491) .EQ. '') Theta17_5 = rNULL__
        if(ar__11(493:512) .EQ. '') Theta18_5 = rNULL__
        if(ar__11(514:534) .EQ. '') Theta19_5 = rNULL__
        if(ar__11(536:556) .EQ. '') Theta20_5 = rNULL__
        if(ar__11(558:578) .EQ. '') Theta21_5 = rNULL__
        if(ar__11(580:599) .EQ. '') Theta22_5 = rNULL__
        if(ar__11(601:621) .EQ. '') Theta23_5 = rNULL__
        if(ar__11(623:643) .EQ. '') Theta24_5 = rNULL__
        if(ar__11(645:665) .EQ. '') Theta25_5 = rNULL__
        if(ar__11(667:686) .EQ. '') Theta26_5 = rNULL__
        if(ar__11(688:707) .EQ. '') Theta27_5 = rNULL__
        if(ar__11(709:729) .EQ. '') Theta28_5 = rNULL__
        if(ar__11(731:751) .EQ. '') Theta29_5 = rNULL__
        if(ar__11(753:773) .EQ. '') Theta30_5 = rNULL__
        if(ar__11(775:795) .EQ. '') Theta31_5 = rNULL__
        if(ar__11(797:814) .EQ. '') Theta32_5 = rNULL__
        if(ar__11(816:836) .EQ. '') Theta33_5 = rNULL__
        if(ar__11(838:855) .EQ. '') Theta34_5 = rNULL__
        if(ar__11(857:876) .EQ. '') Theta35_5 = rNULL__
        if(ar__11(878:895) .EQ. '') Theta36_5 = rNULL__
        if(ar__11(897:917) .EQ. '') Theta37_5 = rNULL__
        if(ar__11(919:939) .EQ. '') Theta38_5 = rNULL__
        if(ar__11(941:961) .EQ. '') Theta39_5 = rNULL__
        if(ar__11(963:982) .EQ. '') Theta40_5 = rNULL__
        if(ar__11(984:1001) .EQ. '') Theta41_5 = rNULL__
c    ..............Just test output...........
        write(6,12)
     +  BJD_11,Flux_11,e_Flux_11,Fluxd_5,e_Fluxd_5,Fluxpsfmodel_5,
     +  Theta0_5,Theta1_5,Theta2_5,Theta3_5,Theta4_5,Theta5_5,
     +  Theta6_5,Theta7_5,Theta8_5,Theta9_5,Theta10_5,Theta11_5,
     +  Theta12_5,Theta13_5,Theta14_5,Theta15_5,Theta16_5,Theta17_5,
     +  Theta18_5,Theta19_5,Theta20_5,Theta21_5,Theta22_5,Theta23_5,
     +  Theta24_5,Theta25_5,Theta26_5,Theta27_5,Theta28_5,Theta29_5,
     +  Theta30_5,Theta31_5,Theta32_5,Theta33_5,Theta34_5,Theta35_5,
     +  Theta36_5,Theta37_5,Theta38_5,Theta39_5,Theta40_5,Theta41_5
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheodet1.dat'	! CHEOPS visit 1 detrended data

C  Format for file interpretation

   13 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheodet1.dat')
      write(6,*) '....Loading file: cheodet1.dat'
      do i__=1,440
        read(1,'(A59)')ar__12
        read(ar__12,13)BJD_12,Flux_12,e_Flux_12
c    ..............Just test output...........
        write(6,13)BJD_12,Flux_12,e_Flux_12
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheodet2.dat'	! CHEOPS visit 2 detrended data

C  Format for file interpretation

   14 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheodet2.dat')
      write(6,*) '....Loading file: cheodet2.dat'
      do i__=1,278
        read(1,'(A59)')ar__13
        read(ar__13,14)BJD_13,Flux_13,e_Flux_13
c    ..............Just test output...........
        write(6,14)BJD_13,Flux_13,e_Flux_13
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheodet3.dat'	! CHEOPS visit 3 detrended data

C  Format for file interpretation

   15 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheodet3.dat')
      write(6,*) '....Loading file: cheodet3.dat'
      do i__=1,282
        read(1,'(A59)')ar__14
        read(ar__14,15)BJD_14,Flux_14,e_Flux_14
c    ..............Just test output...........
        write(6,15)BJD_14,Flux_14,e_Flux_14
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheodet4.dat'	! CHEOPS visit 4 detrended data

C  Format for file interpretation

   16 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheodet4.dat')
      write(6,*) '....Loading file: cheodet4.dat'
      do i__=1,254
        read(1,'(A59)')ar__15
        read(ar__15,16)BJD_15,Flux_15,e_Flux_15
c    ..............Just test output...........
        write(6,16)BJD_15,Flux_15,e_Flux_15
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheodet5.dat'	! CHEOPS visit 5 detrended data

C  Format for file interpretation

   17 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheodet5.dat')
      write(6,*) '....Loading file: cheodet5.dat'
      do i__=1,231
        read(1,'(A59)')ar__16
        read(ar__16,17)BJD_16,Flux_16,e_Flux_16
c    ..............Just test output...........
        write(6,17)BJD_16,Flux_16,e_Flux_16
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cheodet6.dat'	! CHEOPS visit 6 detrended data

C  Format for file interpretation

   18 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'cheodet6.dat')
      write(6,*) '....Loading file: cheodet6.dat'
      do i__=1,762
        read(1,'(A59)')ar__17
        read(ar__17,18)BJD_17,Flux_17,e_Flux_17
c    ..............Just test output...........
        write(6,18)BJD_17,Flux_17,e_Flux_17
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'traws13.dat'	! TESS sector 13 raw data

C  Format for file interpretation

   19 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'traws13.dat')
      write(6,*) '....Loading file: traws13.dat'
      do i__=1,17939
        read(1,'(A59)')ar__18
        read(ar__18,19)BJD_18,Flux_18,e_Flux_18
c    ..............Just test output...........
        write(6,19)BJD_18,Flux_18,e_Flux_18
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'traws27.dat'	! TESS sector 27 raw data

C  Format for file interpretation

   20 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'traws27.dat')
      write(6,*) '....Loading file: traws27.dat'
      do i__=1,71703
        read(1,'(A59)')ar__19
        read(ar__19,20)BJD_19,Flux_19,e_Flux_19
c    ..............Just test output...........
        write(6,20)BJD_19,Flux_19,e_Flux_19
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tdets13.dat'	! TESS sector 13 detrended data

C  Format for file interpretation

   21 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tdets13.dat')
      write(6,*) '....Loading file: tdets13.dat'
      do i__=1,17939
        read(1,'(A59)')ar__20
        read(ar__20,21)BJD_20,Flux_20,e_Flux_20
c    ..............Just test output...........
        write(6,21)BJD_20,Flux_20,e_Flux_20
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tdets27.dat'	! TESS sector 27 detrended data

C  Format for file interpretation

   22 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tdets27.dat')
      write(6,*) '....Loading file: tdets27.dat'
      do i__=1,71703
        read(1,'(A59)')ar__21
        read(ar__21,22)BJD_21,Flux_21,e_Flux_21
c    ..............Just test output...........
        write(6,22)BJD_21,Flux_21,e_Flux_21
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'harpsraw.dat'	! HARPS raw data

C  Format for file interpretation

   23 format(
     +  F16.8,1X,F16.10,1X,F16.14,1X,F4.1,1X,F6.1,1X,F7.4,1X,F6.4,1X,
     +  F6.4,1X,F6.4,1X,F6.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'harpsraw.dat')
      write(6,*) '....Loading file: harpsraw.dat'
      do i__=1,26
        read(1,'(A98)')ar__22
        read(ar__22,23)
     +  BJD_22,RV,e_RV,BIS,FWHM,Contrast,Sindex,Halpha,NaD1,NaD2
c    ..............Just test output...........
        write(6,23)
     +  BJD_22,RV,e_RV,BIS,FWHM,Contrast,Sindex,Halpha,NaD1,NaD2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'harpsccf.dat'	! HARPS CCF SCALPELS data

C  Format for file interpretation

   24 format(F16.8,1X,F16.10,1X,F16.14,1X,F18.15,1X,F18.15,1X,F18.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'harpsccf.dat')
      write(6,*) '....Loading file: harpsccf.dat'
      do i__=1,26
        read(1,'(A107)')ar__23
        read(ar__23,24)BJD_23,RV_1,e_RV_1,Theta0_6,Theta1_6,Theta2_6
c    ..............Just test output...........
        write(6,24)BJD_23,RV_1,e_RV_1,Theta0_6,Theta1_6,Theta2_6
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'ctiozs.dat'	! LCOGT CTIO Zs data

C  Format for file interpretation

   25 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'ctiozs.dat')
      write(6,*) '....Loading file: ctiozs.dat'
      do i__=1,255
        read(1,'(A59)')ar__24
        read(ar__24,25)BJD_24,Flux_22,e_Flux_22
c    ..............Just test output...........
        write(6,25)BJD_24,Flux_22,e_Flux_22
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'saaob.dat'	! LCOGT SAAO B data

C  Format for file interpretation

   26 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'saaob.dat')
      write(6,*) '....Loading file: saaob.dat'
      do i__=1,122
        read(1,'(A59)')ar__25
        read(ar__25,26)BJD_25,Flux_23,e_Flux_23
c    ..............Just test output...........
        write(6,26)BJD_25,Flux_23,e_Flux_23
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'saaozs1.dat'	! LCOGT SAAO Zs visit 1 data

C  Format for file interpretation

   27 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'saaozs1.dat')
      write(6,*) '....Loading file: saaozs1.dat'
      do i__=1,453
        read(1,'(A59)')ar__26
        read(ar__26,27)BJD_26,Flux_24,e_Flux_24
c    ..............Just test output...........
        write(6,27)BJD_26,Flux_24,e_Flux_24
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'saaozs2.dat'	! LCOGT SAAO Zs visit 2 data

C  Format for file interpretation

   28 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'saaozs2.dat')
      write(6,*) '....Loading file: saaozs2.dat'
      do i__=1,368
        read(1,'(A59)')ar__27
        read(ar__27,28)BJD_27,Flux_25,e_Flux_25
c    ..............Just test output...........
        write(6,28)BJD_27,Flux_25,e_Flux_25
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'saaozs3.dat'	! LCOGT SAAO Zs visit 4 data

C  Format for file interpretation

   29 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'saaozs3.dat')
      write(6,*) '....Loading file: saaozs3.dat'
      do i__=1,343
        read(1,'(A59)')ar__28
        read(ar__28,29)BJD_28,Flux_26,e_Flux_26
c    ..............Just test output...........
        write(6,29)BJD_28,Flux_26,e_Flux_26
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'ngtsraw.dat'	! NGTS raw data

C  Format for file interpretation

   30 format(F16.8,1X,F10.8,1X,F10.8,1X,F8.6,1X,F9.3,1X,F8.3,1X,F8.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'ngtsraw.dat')
      write(6,*) '....Loading file: ngtsraw.dat'
      do i__=1,1811
        read(1,'(A75)')ar__29
        read(ar__29,30)BJD_29,Flux_27,e_Flux_27,Airmass,SkyBkg,X,Y
c    ..............Just test output...........
        write(6,30)BJD_29,Flux_27,e_Flux_27,Airmass,SkyBkg,X,Y
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'ngtsdet.dat'	! NGTS detrended data

C  Format for file interpretation

   31 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'ngtsdet.dat')
      write(6,*) '....Loading file: ngtsdet.dat'
      do i__=1,1811
        read(1,'(A59)')ar__30
        read(ar__30,31)BJD_30,Flux_28,e_Flux_28
c    ..............Just test output...........
        write(6,31)BJD_30,Flux_28,e_Flux_28
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'astep1.dat'	! ASTEP visit 1 data

C  Format for file interpretation

   32 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'astep1.dat')
      write(6,*) '....Loading file: astep1.dat'
      do i__=1,297
        read(1,'(A59)')ar__31
        read(ar__31,32)BJD_31,Flux_29,e_Flux_29
c    ..............Just test output...........
        write(6,32)BJD_31,Flux_29,e_Flux_29
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'astep2.dat'	! ASTEP visit 2 data

C  Format for file interpretation

   33 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'astep2.dat')
      write(6,*) '....Loading file: astep2.dat'
      do i__=1,421
        read(1,'(A59)')ar__32
        read(ar__32,33)BJD_32,Flux_30,e_Flux_30
c    ..............Just test output...........
        write(6,33)BJD_32,Flux_30,e_Flux_30
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'astep3.dat'	! ASTEP visit 3 data

C  Format for file interpretation

   34 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'astep3.dat')
      write(6,*) '....Loading file: astep3.dat'
      do i__=1,267
        read(1,'(A59)')ar__33
        read(ar__33,34)BJD_33,Flux_31,e_Flux_31
c    ..............Just test output...........
        write(6,34)BJD_33,Flux_31,e_Flux_31
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'astep4.dat'	! ASTEP visit 3 data

C  Format for file interpretation

   35 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'astep4.dat')
      write(6,*) '....Loading file: astep4.dat'
      do i__=1,178
        read(1,'(A59)')ar__34
        read(ar__34,35)BJD_34,Flux_32,e_Flux_32
c    ..............Just test output...........
        write(6,35)BJD_34,Flux_32,e_Flux_32
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'asassn.dat'	! ASAS-SN data

C  Format for file interpretation

   36 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'asassn.dat')
      write(6,*) '....Loading file: asassn.dat'
      do i__=1,401
        read(1,'(A59)')ar__35
        read(ar__35,36)BJD_35,Flux_33,e_Flux_33
c    ..............Just test output...........
        write(6,36)BJD_35,Flux_33,e_Flux_33
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'wasp.dat'	! WASP data

C  Format for file interpretation

   37 format(F18.10,1X,F18.16,1X,F21.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'wasp.dat')
      write(6,*) '....Loading file: wasp.dat'
      do i__=1,66832
        read(1,'(A59)')ar__36
        read(ar__36,37)BJD_36,Flux_34,e_Flux_34
c    ..............Just test output...........
        write(6,37)BJD_36,Flux_34,e_Flux_34
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gems562.dat'	! Gemini-South 562 data

C  Format for file interpretation

   38 format(I3,1X,F7.5,1X,F7.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gems562.dat')
      write(6,*) '....Loading file: gems562.dat'
      do i__=1,119
        read(1,'(A19)')ar__37
        read(ar__37,38)lambda,Sep,Deltamag
c    ..............Just test output...........
        write(6,38)lambda,Sep,Deltamag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'gems832.dat'	! Gemini-South 832 data

C  Format for file interpretation

   39 format(I3,1X,F7.5,1X,F7.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gems832.dat')
      write(6,*) '....Loading file: gems832.dat'
      do i__=1,120
        read(1,'(A19)')ar__38
        read(ar__38,39)lambda_1,Sep_1,Deltamag_1
c    ..............Just test output...........
        write(6,39)lambda_1,Sep_1,Deltamag_1
c    .......End.of.Just test output...........
      end do
      close(1)

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