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 loading all data files into arrays.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-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        (nr__) ! (d) Barycentric Julian date
      real*8        Flux       (nr__) ! Normalised flux
      real*8        e_Flux     (nr__) ! Normalised flux error
      real*4        Xoffset    (nr__) ! (pix) X centroid position offset
      real*8        Yoffset    (nr__) ! (pix) Y centroid position offset
      real*8        Roll       (nr__) ! Roll angle
      real*8        Contam     (nr__) ! (ct) Contamination
      real*8        Bckg       (nr__) ! (ct) Background
      real*4        Smear      (nr__) ! (ct) Smear
      real*8        DeltaT     (nr__) ! 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      (nr__1) ! (d) Barycentric Julian date
      real*8        Flux_1     (nr__1) ! Normalised flux
      real*8        e_Flux_1   (nr__1) ! Normalised flux error
      real*4        Xoffset_1  (nr__1) ! (pix) X centroid position offset
      real*8        Yoffset_1  (nr__1) ! (pix) Y centroid position offset
      real*8        Roll_1     (nr__1) ! Roll angle
      real*8        Contam_1   (nr__1) ! (ct) Contamination
      real*8        Bckg_1     (nr__1) ! (ct) Background
      real*4        Smear_1    (nr__1) ! (ct) Smear
      real*8        DeltaT_1   (nr__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      (nr__2) ! (d) Barycentric Julian date
      real*8        Flux_2     (nr__2) ! Normalised flux
      real*8        e_Flux_2   (nr__2) ! Normalised flux error
      real*4        Xoffset_2  (nr__2) ! (pix) X centroid position offset
      real*8        Yoffset_2  (nr__2) ! (pix) Y centroid position offset
      real*8        Roll_2     (nr__2) ! Roll angle
      real*8        Contam_2   (nr__2) ! (ct) Contamination
      real*8        Bckg_2     (nr__2) ! (ct) Background
      real*4        Smear_2    (nr__2) ! (ct) Smear
      real*8        DeltaT_2   (nr__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      (nr__3) ! (d) Barycentric Julian date
      real*8        Flux_3     (nr__3) ! Normalised flux
      real*8        e_Flux_3   (nr__3) ! Normalised flux error
      real*4        Xoffset_3  (nr__3) ! (pix) X centroid position offset
      real*8        Yoffset_3  (nr__3) ! (pix) Y centroid position offset
      real*8        Roll_3     (nr__3) ! Roll angle
      real*8        Contam_3   (nr__3) ! (ct) Contamination
      real*8        Bckg_3     (nr__3) ! (ct) Background
      real*4        Smear_3    (nr__3) ! (ct) Smear
      real*8        DeltaT_3   (nr__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      (nr__4) ! (d) Barycentric Julian date
      real*8        Flux_4     (nr__4) ! Normalised flux
      real*8        e_Flux_4   (nr__4) ! Normalised flux error
      real*4        Xoffset_4  (nr__4) ! (pix) X centroid position offset
      real*8        Yoffset_4  (nr__4) ! (pix) Y centroid position offset
      real*8        Roll_4     (nr__4) ! Roll angle
      real*8        Contam_4   (nr__4) ! (ct) Contamination
      real*8        Bckg_4     (nr__4) ! (ct) Background
      real*4        Smear_4    (nr__4) ! (ct) Smear
      real*8        DeltaT_4   (nr__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      (nr__5) ! (d) Barycentric Julian date
      real*8        Flux_5     (nr__5) ! Normalised flux
      real*8        e_Flux_5   (nr__5) ! Normalised flux error
      real*4        Xoffset_5  (nr__5) ! (pix) X centroid position offset
      real*8        Yoffset_5  (nr__5) ! (pix) Y centroid position offset
      real*8        Roll_5     (nr__5) ! Roll angle
      real*8        Contam_5   (nr__5) ! (ct) Contamination
      real*8        Bckg_5     (nr__5) ! (ct) Background
      real*4        Smear_5    (nr__5) ! (ct) Smear
      real*8        DeltaT_5   (nr__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      (nr__6) ! (d) Barycentric Julian date
      real*8        Flux_6     (nr__6) ! Normalised flux
      real*8        e_Flux_6   (nr__6) ! Normalised flux error
      real*8        Fluxd      (nr__6) ! Normalised detrended flux
      real*8        e_Fluxd    (nr__6) ! Normalised detrended flux error
      real*8        Fluxpsfmodel(nr__6) ! PSF SCALPELS PSF shape change model
      real*4        Theta0     (nr__6) ! 1st PSF SCALPELS principal component
      real*4        Theta1     (nr__6) ! 2nd PSF SCALPELS principal component
      real*8        Theta2     (nr__6) ! 3rd PSF SCALPELS principal component
      real*4        Theta3     (nr__6) ! 4th PSF SCALPELS principal component
      real*8        Theta4     (nr__6) ! 5th PSF SCALPELS principal component
      real*4        Theta5     (nr__6) ! 6th PSF SCALPELS principal component
      real*4        Theta6     (nr__6) ! 7th PSF SCALPELS principal component
      real*4        Theta7     (nr__6) ! 8th PSF SCALPELS principal component
      real*4        Theta8     (nr__6) ! 9th PSF SCALPELS principal component
      real*4        Theta9     (nr__6) ! 10th PSF SCALPELS principal component
      real*4        Theta10    (nr__6) ! ? 11th PSF SCALPELS principal component
      real*4        Theta11    (nr__6) ! ? 12th PSF SCALPELS principal component
      real*4        Theta12    (nr__6) ! ? 13th PSF SCALPELS principal component
      real*4        Theta13    (nr__6) ! ? 14th PSF SCALPELS principal component
      real*4        Theta14    (nr__6) ! ? 15th PSF SCALPELS principal component
      real*4        Theta15    (nr__6) ! ? 16th PSF SCALPELS principal component
      real*4        Theta16    (nr__6) ! ? 17th PSF SCALPELS principal component
      real*4        Theta17    (nr__6) ! ? 18th PSF SCALPELS principal component
      real*4        Theta18    (nr__6) ! ? 19th PSF SCALPELS principal component
      real*4        Theta19    (nr__6) ! ? 20th PSF SCALPELS principal component
      real*4        Theta20    (nr__6) ! ? 21st PSF SCALPELS principal component
      real*4        Theta21    (nr__6) ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22    (nr__6) ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23    (nr__6) ! ? 24th PSF SCALPELS principal component
      real*4        Theta24    (nr__6) ! ? 25th PSF SCALPELS principal component
      real*4        Theta25    (nr__6) ! ? 26th PSF SCALPELS principal component
      real*4        Theta26    (nr__6) ! ? 27th PSF SCALPELS principal component
      real*4        Theta27    (nr__6) ! ? 28th PSF SCALPELS principal component
      real*4        Theta28    (nr__6) ! ? 29th PSF SCALPELS principal component
      real*4        Theta29    (nr__6) ! ? 30th PSF SCALPELS principal component
      real*4        Theta30    (nr__6) ! ? 31st PSF SCALPELS principal component
      real*4        Theta31    (nr__6) ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32    (nr__6) ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33    (nr__6) ! ? 34th PSF SCALPELS principal component
      real*8        Theta34    (nr__6) ! ? 35th PSF SCALPELS principal component
      real*4        Theta35    (nr__6) ! ? 36th PSF SCALPELS principal component
      real*8        Theta36    (nr__6) ! ? 37th PSF SCALPELS principal component
      real*4        Theta37    (nr__6) ! ? 38th PSF SCALPELS principal component
      real*4        Theta38    (nr__6) ! ? 39th PSF SCALPELS principal component
      real*4        Theta39    (nr__6) ! ? 40th PSF SCALPELS principal component
      real*4        Theta40    (nr__6) ! ? 41st PSF SCALPELS principal component
      real*8        Theta41    (nr__6) ! ? 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      (nr__7) ! (d) Barycentric Julian date
      real*8        Flux_7     (nr__7) ! Normalised flux
      real*8        e_Flux_7   (nr__7) ! Normalised flux error
      real*8        Fluxd_1    (nr__7) ! Normalised detrended flux
      real*8        e_Fluxd_1  (nr__7) ! Normalised detrended flux error
      real*8        Fluxpsfmodel_1(nr__7) ! PSF SCALPELS PSF shape change model
      real*4        Theta0_1   (nr__7) ! 1st PSF SCALPELS principal component
      real*4        Theta1_1   (nr__7) ! 2nd PSF SCALPELS principal component
      real*8        Theta2_1   (nr__7) ! 3rd PSF SCALPELS principal component
      real*4        Theta3_1   (nr__7) ! 4th PSF SCALPELS principal component
      real*8        Theta4_1   (nr__7) ! 5th PSF SCALPELS principal component
      real*4        Theta5_1   (nr__7) ! 6th PSF SCALPELS principal component
      real*4        Theta6_1   (nr__7) ! 7th PSF SCALPELS principal component
      real*4        Theta7_1   (nr__7) ! 8th PSF SCALPELS principal component
      real*4        Theta8_1   (nr__7) ! 9th PSF SCALPELS principal component
      real*4        Theta9_1   (nr__7) ! 10th PSF SCALPELS principal component
      real*4        Theta10_1  (nr__7) ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_1  (nr__7) ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_1  (nr__7) ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_1  (nr__7) ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_1  (nr__7) ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_1  (nr__7) ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_1  (nr__7) ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_1  (nr__7) ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_1  (nr__7) ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_1  (nr__7) ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_1  (nr__7) ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_1  (nr__7) ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_1  (nr__7) ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_1  (nr__7) ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_1  (nr__7) ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_1  (nr__7) ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_1  (nr__7) ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_1  (nr__7) ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_1  (nr__7) ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_1  (nr__7) ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_1  (nr__7) ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_1  (nr__7) ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_1  (nr__7) ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_1  (nr__7) ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_1  (nr__7) ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_1  (nr__7) ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_1  (nr__7) ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_1  (nr__7) ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_1  (nr__7) ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_1  (nr__7) ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_1  (nr__7) ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_1  (nr__7) ! ? 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      (nr__8) ! (d) Barycentric Julian date
      real*8        Flux_8     (nr__8) ! Normalised flux
      real*8        e_Flux_8   (nr__8) ! Normalised flux error
      real*8        Fluxd_2    (nr__8) ! Normalised detrended flux
      real*8        e_Fluxd_2  (nr__8) ! Normalised detrended flux error
      real*8        Fluxpsfmodel_2(nr__8) ! PSF SCALPELS PSF shape change model
      real*4        Theta0_2   (nr__8) ! 1st PSF SCALPELS principal component
      real*4        Theta1_2   (nr__8) ! 2nd PSF SCALPELS principal component
      real*8        Theta2_2   (nr__8) ! 3rd PSF SCALPELS principal component
      real*4        Theta3_2   (nr__8) ! 4th PSF SCALPELS principal component
      real*8        Theta4_2   (nr__8) ! 5th PSF SCALPELS principal component
      real*4        Theta5_2   (nr__8) ! 6th PSF SCALPELS principal component
      real*4        Theta6_2   (nr__8) ! 7th PSF SCALPELS principal component
      real*4        Theta7_2   (nr__8) ! 8th PSF SCALPELS principal component
      real*4        Theta8_2   (nr__8) ! 9th PSF SCALPELS principal component
      real*4        Theta9_2   (nr__8) ! 10th PSF SCALPELS principal component
      real*4        Theta10_2  (nr__8) ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_2  (nr__8) ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_2  (nr__8) ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_2  (nr__8) ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_2  (nr__8) ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_2  (nr__8) ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_2  (nr__8) ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_2  (nr__8) ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_2  (nr__8) ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_2  (nr__8) ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_2  (nr__8) ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_2  (nr__8) ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_2  (nr__8) ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_2  (nr__8) ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_2  (nr__8) ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_2  (nr__8) ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_2  (nr__8) ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_2  (nr__8) ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_2  (nr__8) ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_2  (nr__8) ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_2  (nr__8) ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_2  (nr__8) ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_2  (nr__8) ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_2  (nr__8) ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_2  (nr__8) ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_2  (nr__8) ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_2  (nr__8) ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_2  (nr__8) ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_2  (nr__8) ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_2  (nr__8) ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_2  (nr__8) ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_2  (nr__8) ! ? 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      (nr__9) ! (d) Barycentric Julian date
      real*8        Flux_9     (nr__9) ! Normalised flux
      real*8        e_Flux_9   (nr__9) ! Normalised flux error
      real*8        Fluxd_3    (nr__9) ! Normalised detrended flux
      real*8        e_Fluxd_3  (nr__9) ! Normalised detrended flux error
      real*8        Fluxpsfmodel_3(nr__9) ! PSF SCALPELS PSF shape change model
      real*4        Theta0_3   (nr__9) ! 1st PSF SCALPELS principal component
      real*4        Theta1_3   (nr__9) ! 2nd PSF SCALPELS principal component
      real*8        Theta2_3   (nr__9) ! 3rd PSF SCALPELS principal component
      real*4        Theta3_3   (nr__9) ! 4th PSF SCALPELS principal component
      real*8        Theta4_3   (nr__9) ! 5th PSF SCALPELS principal component
      real*4        Theta5_3   (nr__9) ! 6th PSF SCALPELS principal component
      real*4        Theta6_3   (nr__9) ! 7th PSF SCALPELS principal component
      real*4        Theta7_3   (nr__9) ! 8th PSF SCALPELS principal component
      real*4        Theta8_3   (nr__9) ! 9th PSF SCALPELS principal component
      real*4        Theta9_3   (nr__9) ! 10th PSF SCALPELS principal component
      real*4        Theta10_3  (nr__9) ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_3  (nr__9) ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_3  (nr__9) ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_3  (nr__9) ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_3  (nr__9) ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_3  (nr__9) ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_3  (nr__9) ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_3  (nr__9) ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_3  (nr__9) ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_3  (nr__9) ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_3  (nr__9) ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_3  (nr__9) ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_3  (nr__9) ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_3  (nr__9) ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_3  (nr__9) ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_3  (nr__9) ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_3  (nr__9) ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_3  (nr__9) ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_3  (nr__9) ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_3  (nr__9) ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_3  (nr__9) ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_3  (nr__9) ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_3  (nr__9) ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_3  (nr__9) ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_3  (nr__9) ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_3  (nr__9) ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_3  (nr__9) ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_3  (nr__9) ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_3  (nr__9) ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_3  (nr__9) ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_3  (nr__9) ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_3  (nr__9) ! ? 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     (nr__10) ! (d) Barycentric Julian date
      real*8        Flux_10    (nr__10) ! Normalised flux
      real*8        e_Flux_10  (nr__10) ! Normalised flux error
      real*8        Fluxd_4    (nr__10) ! Normalised detrended flux
      real*8        e_Fluxd_4  (nr__10) ! Normalised detrended flux error
      real*8        Fluxpsfmodel_4(nr__10) ! PSF SCALPELS PSF shape change model
      real*4        Theta0_4   (nr__10) ! 1st PSF SCALPELS principal component
      real*4        Theta1_4   (nr__10) ! 2nd PSF SCALPELS principal component
      real*8        Theta2_4   (nr__10) ! 3rd PSF SCALPELS principal component
      real*4        Theta3_4   (nr__10) ! 4th PSF SCALPELS principal component
      real*8        Theta4_4   (nr__10) ! 5th PSF SCALPELS principal component
      real*4        Theta5_4   (nr__10) ! 6th PSF SCALPELS principal component
      real*4        Theta6_4   (nr__10) ! 7th PSF SCALPELS principal component
      real*4        Theta7_4   (nr__10) ! 8th PSF SCALPELS principal component
      real*4        Theta8_4   (nr__10) ! 9th PSF SCALPELS principal component
      real*4        Theta9_4   (nr__10) ! 10th PSF SCALPELS principal component
      real*4        Theta10_4  (nr__10) ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_4  (nr__10) ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_4  (nr__10) ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_4  (nr__10) ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_4  (nr__10) ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_4  (nr__10) ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_4  (nr__10) ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_4  (nr__10) ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_4  (nr__10) ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_4  (nr__10) ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_4  (nr__10) ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_4  (nr__10) ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_4  (nr__10) ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_4  (nr__10) ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_4  (nr__10) ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_4  (nr__10) ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_4  (nr__10) ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_4  (nr__10) ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_4  (nr__10) ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_4  (nr__10) ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_4  (nr__10) ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_4  (nr__10) ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_4  (nr__10) ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_4  (nr__10) ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_4  (nr__10) ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_4  (nr__10) ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_4  (nr__10) ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_4  (nr__10) ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_4  (nr__10) ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_4  (nr__10) ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_4  (nr__10) ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_4  (nr__10) ! ? 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     (nr__11) ! (d) Barycentric Julian date
      real*8        Flux_11    (nr__11) ! Normalised flux
      real*8        e_Flux_11  (nr__11) ! Normalised flux error
      real*8        Fluxd_5    (nr__11) ! Normalised detrended flux
      real*8        e_Fluxd_5  (nr__11) ! Normalised detrended flux error
      real*8        Fluxpsfmodel_5(nr__11) ! PSF SCALPELS PSF shape change model
      real*4        Theta0_5   (nr__11) ! 1st PSF SCALPELS principal component
      real*4        Theta1_5   (nr__11) ! 2nd PSF SCALPELS principal component
      real*8        Theta2_5   (nr__11) ! 3rd PSF SCALPELS principal component
      real*4        Theta3_5   (nr__11) ! 4th PSF SCALPELS principal component
      real*8        Theta4_5   (nr__11) ! 5th PSF SCALPELS principal component
      real*4        Theta5_5   (nr__11) ! 6th PSF SCALPELS principal component
      real*4        Theta6_5   (nr__11) ! 7th PSF SCALPELS principal component
      real*4        Theta7_5   (nr__11) ! 8th PSF SCALPELS principal component
      real*4        Theta8_5   (nr__11) ! 9th PSF SCALPELS principal component
      real*4        Theta9_5   (nr__11) ! 10th PSF SCALPELS principal component
      real*4        Theta10_5  (nr__11) ! ? 11th PSF SCALPELS principal component
      real*4        Theta11_5  (nr__11) ! ? 12th PSF SCALPELS principal component
      real*4        Theta12_5  (nr__11) ! ? 13th PSF SCALPELS principal component
      real*4        Theta13_5  (nr__11) ! ? 14th PSF SCALPELS principal component
      real*4        Theta14_5  (nr__11) ! ? 15th PSF SCALPELS principal component
      real*4        Theta15_5  (nr__11) ! ? 16th PSF SCALPELS principal component
      real*4        Theta16_5  (nr__11) ! ? 17th PSF SCALPELS principal component
      real*4        Theta17_5  (nr__11) ! ? 18th PSF SCALPELS principal component
      real*4        Theta18_5  (nr__11) ! ? 19th PSF SCALPELS principal component
      real*4        Theta19_5  (nr__11) ! ? 20th PSF SCALPELS principal component
      real*4        Theta20_5  (nr__11) ! ? 21st PSF SCALPELS principal component
      real*4        Theta21_5  (nr__11) ! ? 22nd PSF SCALPELS principal component
      real*4        Theta22_5  (nr__11) ! ? 23rd PSF SCALPELS principal component
      real*4        Theta23_5  (nr__11) ! ? 24th PSF SCALPELS principal component
      real*4        Theta24_5  (nr__11) ! ? 25th PSF SCALPELS principal component
      real*4        Theta25_5  (nr__11) ! ? 26th PSF SCALPELS principal component
      real*4        Theta26_5  (nr__11) ! ? 27th PSF SCALPELS principal component
      real*4        Theta27_5  (nr__11) ! ? 28th PSF SCALPELS principal component
      real*4        Theta28_5  (nr__11) ! ? 29th PSF SCALPELS principal component
      real*4        Theta29_5  (nr__11) ! ? 30th PSF SCALPELS principal component
      real*4        Theta30_5  (nr__11) ! ? 31st PSF SCALPELS principal component
      real*4        Theta31_5  (nr__11) ! ? 32nd PSF SCALPELS principal component
      real*8        Theta32_5  (nr__11) ! ? 33rd PSF SCALPELS principal component
      real*4        Theta33_5  (nr__11) ! ? 34th PSF SCALPELS principal component
      real*8        Theta34_5  (nr__11) ! ? 35th PSF SCALPELS principal component
      real*4        Theta35_5  (nr__11) ! ? 36th PSF SCALPELS principal component
      real*8        Theta36_5  (nr__11) ! ? 37th PSF SCALPELS principal component
      real*4        Theta37_5  (nr__11) ! ? 38th PSF SCALPELS principal component
      real*4        Theta38_5  (nr__11) ! ? 39th PSF SCALPELS principal component
      real*4        Theta39_5  (nr__11) ! ? 40th PSF SCALPELS principal component
      real*4        Theta40_5  (nr__11) ! ? 41st PSF SCALPELS principal component
      real*8        Theta41_5  (nr__11) ! ? 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     (nr__12) ! (d) Barycentric Julian date
      real*8        Flux_12    (nr__12) ! Normalised detrended flux
      real*8        e_Flux_12  (nr__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     (nr__13) ! (d) Barycentric Julian date
      real*8        Flux_13    (nr__13) ! Normalised detrended flux
      real*8        e_Flux_13  (nr__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     (nr__14) ! (d) Barycentric Julian date
      real*8        Flux_14    (nr__14) ! Normalised detrended flux
      real*8        e_Flux_14  (nr__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     (nr__15) ! (d) Barycentric Julian date
      real*8        Flux_15    (nr__15) ! Normalised detrended flux
      real*8        e_Flux_15  (nr__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     (nr__16) ! (d) Barycentric Julian date
      real*8        Flux_16    (nr__16) ! Normalised detrended flux
      real*8        e_Flux_16  (nr__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     (nr__17) ! (d) Barycentric Julian date
      real*8        Flux_17    (nr__17) ! Normalised detrended flux
      real*8        e_Flux_17  (nr__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     (nr__18) ! (d) Barycentric Julian date
      real*8        Flux_18    (nr__18) ! Normalised flux
      real*8        e_Flux_18  (nr__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     (nr__19) ! (d) Barycentric Julian date
      real*8        Flux_19    (nr__19) ! Normalised flux
      real*8        e_Flux_19  (nr__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     (nr__20) ! (d) Barycentric Julian date
      real*8        Flux_20    (nr__20) ! Normalised detrended flux
      real*8        e_Flux_20  (nr__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     (nr__21) ! (d) Barycentric Julian date
      real*8        Flux_21    (nr__21) ! Normalised detrended flux
      real*8        e_Flux_21  (nr__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     (nr__22) ! (d) Barycentric Julian date
      real*8        RV         (nr__22) ! (m/s) Radial velocity
      real*8        e_RV       (nr__22) ! (m/s) Radial velocity error
      real*4        BIS        (nr__22) ! (m/s) Bisector span of the CCF
      real*4        FWHM       (nr__22) ! (m/s) FWHM of the CCF
      real*8        Contrast   (nr__22) ! Contrast of the CCF
      real*4        Sindex     (nr__22) ! S index
      real*4        Halpha     (nr__22) ! H{alpha} activity index
      real*4        NaD1       (nr__22) ! Na D1 activity index
      real*4        NaD2       (nr__22) ! 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     (nr__23) ! (d) Barycentric Julian date
      real*8        RV_1       (nr__23) ! (m/s) Radial velocity
      real*8        e_RV_1     (nr__23) ! (m/s) Radial velocity error
      real*8        Theta0_6   (nr__23) ! 1st CCF SCALPELS principal component
      real*8        Theta1_6   (nr__23) ! 2nd CCF SCALPELS principal component
      real*8        Theta2_6   (nr__23) ! 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     (nr__24) ! (d) Barycentric Julian date
      real*8        Flux_22    (nr__24) ! Normalised flux
      real*8        e_Flux_22  (nr__24) ! 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     (nr__25) ! (d) Barycentric Julian date
      real*8        Flux_23    (nr__25) ! Normalised flux
      real*8        e_Flux_23  (nr__25) ! 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     (nr__26) ! (d) Barycentric Julian date
      real*8        Flux_24    (nr__26) ! Normalised flux
      real*8        e_Flux_24  (nr__26) ! 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     (nr__27) ! (d) Barycentric Julian date
      real*8        Flux_25    (nr__27) ! Normalised flux
      real*8        e_Flux_25  (nr__27) ! 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     (nr__28) ! (d) Barycentric Julian date
      real*8        Flux_26    (nr__28) ! Normalised flux
      real*8        e_Flux_26  (nr__28) ! 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     (nr__29) ! (d) Barycentric Julian date
      real*8        Flux_27    (nr__29) ! Normalised flux
      real*8        e_Flux_27  (nr__29) ! Normalised flux error
      real*8        Airmass    (nr__29) ! Airmass
      real*8        SkyBkg     (nr__29) ! (ct) Sky background
      real*8        X          (nr__29) ! (pix) X centroid position
      real*8        Y          (nr__29) ! (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     (nr__30) ! (d) Barycentric Julian date
      real*8        Flux_28    (nr__30) ! Normalised detrended flux
      real*8        e_Flux_28  (nr__30) ! 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     (nr__31) ! (d) Barycentric Julian date
      real*8        Flux_29    (nr__31) ! Normalised flux
      real*8        e_Flux_29  (nr__31) ! 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     (nr__32) ! (d) Barycentric Julian date
      real*8        Flux_30    (nr__32) ! Normalised flux
      real*8        e_Flux_30  (nr__32) ! 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     (nr__33) ! (d) Barycentric Julian date
      real*8        Flux_31    (nr__33) ! Normalised flux
      real*8        e_Flux_31  (nr__33) ! 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     (nr__34) ! (d) Barycentric Julian date
      real*8        Flux_32    (nr__34) ! Normalised flux
      real*8        e_Flux_32  (nr__34) ! 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     (nr__35) ! (d) Barycentric Julian date
      real*8        Flux_33    (nr__35) ! Normalised flux
      real*8        e_Flux_33  (nr__35) ! 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     (nr__36) ! (d) Barycentric Julian date
      real*8        Flux_34    (nr__36) ! Normalised flux
      real*8        e_Flux_34  (nr__36) ! 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     (nr__37) ! (nm) Wavelength (562 or 832)
      real*8        Sep        (nr__37) ! (arcsec) Separation
      real*8        Deltamag   (nr__37) ! 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   (nr__38) ! (nm) Wavelength (562 or 832)
      real*8        Sep_1      (nr__38) ! (arcsec) Separation
      real*8        Deltamag_1 (nr__38) ! 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(i__),Flux(i__),e_Flux(i__),Xoffset(i__),Yoffset(i__),
     +  Roll(i__),Contam(i__),Bckg(i__),Smear(i__),DeltaT(i__)
c    ..............Just test output...........
        write(6,1)
     +  BJD(i__),Flux(i__),e_Flux(i__),Xoffset(i__),Yoffset(i__),
     +  Roll(i__),Contam(i__),Bckg(i__),Smear(i__),DeltaT(i__)
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(i__),Flux_1(i__),e_Flux_1(i__),Xoffset_1(i__),
     +  Yoffset_1(i__),Roll_1(i__),Contam_1(i__),Bckg_1(i__),
     +  Smear_1(i__),DeltaT_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  BJD_1(i__),Flux_1(i__),e_Flux_1(i__),Xoffset_1(i__),
     +  Yoffset_1(i__),Roll_1(i__),Contam_1(i__),Bckg_1(i__),
     +  Smear_1(i__),DeltaT_1(i__)
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(i__),Flux_2(i__),e_Flux_2(i__),Xoffset_2(i__),
     +  Yoffset_2(i__),Roll_2(i__),Contam_2(i__),Bckg_2(i__),
     +  Smear_2(i__),DeltaT_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  BJD_2(i__),Flux_2(i__),e_Flux_2(i__),Xoffset_2(i__),
     +  Yoffset_2(i__),Roll_2(i__),Contam_2(i__),Bckg_2(i__),
     +  Smear_2(i__),DeltaT_2(i__)
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(i__),Flux_3(i__),e_Flux_3(i__),Xoffset_3(i__),
     +  Yoffset_3(i__),Roll_3(i__),Contam_3(i__),Bckg_3(i__),
     +  Smear_3(i__),DeltaT_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  BJD_3(i__),Flux_3(i__),e_Flux_3(i__),Xoffset_3(i__),
     +  Yoffset_3(i__),Roll_3(i__),Contam_3(i__),Bckg_3(i__),
     +  Smear_3(i__),DeltaT_3(i__)
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(i__),Flux_4(i__),e_Flux_4(i__),Xoffset_4(i__),
     +  Yoffset_4(i__),Roll_4(i__),Contam_4(i__),Bckg_4(i__),
     +  Smear_4(i__),DeltaT_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  BJD_4(i__),Flux_4(i__),e_Flux_4(i__),Xoffset_4(i__),
     +  Yoffset_4(i__),Roll_4(i__),Contam_4(i__),Bckg_4(i__),
     +  Smear_4(i__),DeltaT_4(i__)
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(i__),Flux_5(i__),e_Flux_5(i__),Xoffset_5(i__),
     +  Yoffset_5(i__),Roll_5(i__),Contam_5(i__),Bckg_5(i__),
     +  Smear_5(i__),DeltaT_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  BJD_5(i__),Flux_5(i__),e_Flux_5(i__),Xoffset_5(i__),
     +  Yoffset_5(i__),Roll_5(i__),Contam_5(i__),Bckg_5(i__),
     +  Smear_5(i__),DeltaT_5(i__)
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(i__),Flux_6(i__),e_Flux_6(i__),Fluxd(i__),e_Fluxd(i__),
     +  Fluxpsfmodel(i__),Theta0(i__),Theta1(i__),Theta2(i__),
     +  Theta3(i__),Theta4(i__),Theta5(i__),Theta6(i__),Theta7(i__),
     +  Theta8(i__),Theta9(i__),Theta10(i__),Theta11(i__),
     +  Theta12(i__),Theta13(i__),Theta14(i__),Theta15(i__),
     +  Theta16(i__),Theta17(i__),Theta18(i__),Theta19(i__),
     +  Theta20(i__),Theta21(i__),Theta22(i__),Theta23(i__),
     +  Theta24(i__),Theta25(i__),Theta26(i__),Theta27(i__),
     +  Theta28(i__),Theta29(i__),Theta30(i__),Theta31(i__),
     +  Theta32(i__),Theta33(i__),Theta34(i__),Theta35(i__),
     +  Theta36(i__),Theta37(i__),Theta38(i__),Theta39(i__),
     +  Theta40(i__),Theta41(i__)
        if(ar__6(320:339) .EQ. '') Theta10(i__) = rNULL__
        if(ar__6(341:361) .EQ. '') Theta11(i__) = rNULL__
        if(ar__6(362:382) .EQ. '') Theta12(i__) = rNULL__
        if(ar__6(384:404) .EQ. '') Theta13(i__) = rNULL__
        if(ar__6(406:426) .EQ. '') Theta14(i__) = rNULL__
        if(ar__6(428:448) .EQ. '') Theta15(i__) = rNULL__
        if(ar__6(449:469) .EQ. '') Theta16(i__) = rNULL__
        if(ar__6(471:491) .EQ. '') Theta17(i__) = rNULL__
        if(ar__6(493:512) .EQ. '') Theta18(i__) = rNULL__
        if(ar__6(514:534) .EQ. '') Theta19(i__) = rNULL__
        if(ar__6(536:556) .EQ. '') Theta20(i__) = rNULL__
        if(ar__6(558:578) .EQ. '') Theta21(i__) = rNULL__
        if(ar__6(580:599) .EQ. '') Theta22(i__) = rNULL__
        if(ar__6(601:621) .EQ. '') Theta23(i__) = rNULL__
        if(ar__6(623:643) .EQ. '') Theta24(i__) = rNULL__
        if(ar__6(645:665) .EQ. '') Theta25(i__) = rNULL__
        if(ar__6(667:686) .EQ. '') Theta26(i__) = rNULL__
        if(ar__6(688:707) .EQ. '') Theta27(i__) = rNULL__
        if(ar__6(709:729) .EQ. '') Theta28(i__) = rNULL__
        if(ar__6(731:751) .EQ. '') Theta29(i__) = rNULL__
        if(ar__6(753:773) .EQ. '') Theta30(i__) = rNULL__
        if(ar__6(775:795) .EQ. '') Theta31(i__) = rNULL__
        if(ar__6(797:814) .EQ. '') Theta32(i__) = rNULL__
        if(ar__6(816:836) .EQ. '') Theta33(i__) = rNULL__
        if(ar__6(838:855) .EQ. '') Theta34(i__) = rNULL__
        if(ar__6(857:876) .EQ. '') Theta35(i__) = rNULL__
        if(ar__6(878:895) .EQ. '') Theta36(i__) = rNULL__
        if(ar__6(897:917) .EQ. '') Theta37(i__) = rNULL__
        if(ar__6(919:939) .EQ. '') Theta38(i__) = rNULL__
        if(ar__6(941:961) .EQ. '') Theta39(i__) = rNULL__
        if(ar__6(963:982) .EQ. '') Theta40(i__) = rNULL__
        if(ar__6(984:1001) .EQ. '') Theta41(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  BJD_6(i__),Flux_6(i__),e_Flux_6(i__),Fluxd(i__),e_Fluxd(i__),
     +  Fluxpsfmodel(i__),Theta0(i__),Theta1(i__),Theta2(i__),
     +  Theta3(i__),Theta4(i__),Theta5(i__),Theta6(i__),Theta7(i__),
     +  Theta8(i__),Theta9(i__),Theta10(i__),Theta11(i__),
     +  Theta12(i__),Theta13(i__),Theta14(i__),Theta15(i__),
     +  Theta16(i__),Theta17(i__),Theta18(i__),Theta19(i__),
     +  Theta20(i__),Theta21(i__),Theta22(i__),Theta23(i__),
     +  Theta24(i__),Theta25(i__),Theta26(i__),Theta27(i__),
     +  Theta28(i__),Theta29(i__),Theta30(i__),Theta31(i__),
     +  Theta32(i__),Theta33(i__),Theta34(i__),Theta35(i__),
     +  Theta36(i__),Theta37(i__),Theta38(i__),Theta39(i__),
     +  Theta40(i__),Theta41(i__)
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(i__),Flux_7(i__),e_Flux_7(i__),Fluxd_1(i__),
     +  e_Fluxd_1(i__),Fluxpsfmodel_1(i__),Theta0_1(i__),
     +  Theta1_1(i__),Theta2_1(i__),Theta3_1(i__),Theta4_1(i__),
     +  Theta5_1(i__),Theta6_1(i__),Theta7_1(i__),Theta8_1(i__),
     +  Theta9_1(i__),Theta10_1(i__),Theta11_1(i__),Theta12_1(i__),
     +  Theta13_1(i__),Theta14_1(i__),Theta15_1(i__),Theta16_1(i__),
     +  Theta17_1(i__),Theta18_1(i__),Theta19_1(i__),Theta20_1(i__),
     +  Theta21_1(i__),Theta22_1(i__),Theta23_1(i__),Theta24_1(i__),
     +  Theta25_1(i__),Theta26_1(i__),Theta27_1(i__),Theta28_1(i__),
     +  Theta29_1(i__),Theta30_1(i__),Theta31_1(i__),Theta32_1(i__),
     +  Theta33_1(i__),Theta34_1(i__),Theta35_1(i__),Theta36_1(i__),
     +  Theta37_1(i__),Theta38_1(i__),Theta39_1(i__),Theta40_1(i__),
     +  Theta41_1(i__)
        if(ar__7(320:339) .EQ. '') Theta10_1(i__) = rNULL__
        if(ar__7(341:361) .EQ. '') Theta11_1(i__) = rNULL__
        if(ar__7(362:382) .EQ. '') Theta12_1(i__) = rNULL__
        if(ar__7(384:404) .EQ. '') Theta13_1(i__) = rNULL__
        if(ar__7(406:426) .EQ. '') Theta14_1(i__) = rNULL__
        if(ar__7(428:448) .EQ. '') Theta15_1(i__) = rNULL__
        if(ar__7(449:469) .EQ. '') Theta16_1(i__) = rNULL__
        if(ar__7(471:491) .EQ. '') Theta17_1(i__) = rNULL__
        if(ar__7(493:512) .EQ. '') Theta18_1(i__) = rNULL__
        if(ar__7(514:534) .EQ. '') Theta19_1(i__) = rNULL__
        if(ar__7(536:556) .EQ. '') Theta20_1(i__) = rNULL__
        if(ar__7(558:578) .EQ. '') Theta21_1(i__) = rNULL__
        if(ar__7(580:599) .EQ. '') Theta22_1(i__) = rNULL__
        if(ar__7(601:621) .EQ. '') Theta23_1(i__) = rNULL__
        if(ar__7(623:643) .EQ. '') Theta24_1(i__) = rNULL__
        if(ar__7(645:665) .EQ. '') Theta25_1(i__) = rNULL__
        if(ar__7(667:686) .EQ. '') Theta26_1(i__) = rNULL__
        if(ar__7(688:707) .EQ. '') Theta27_1(i__) = rNULL__
        if(ar__7(709:729) .EQ. '') Theta28_1(i__) = rNULL__
        if(ar__7(731:751) .EQ. '') Theta29_1(i__) = rNULL__
        if(ar__7(753:773) .EQ. '') Theta30_1(i__) = rNULL__
        if(ar__7(775:795) .EQ. '') Theta31_1(i__) = rNULL__
        if(ar__7(797:814) .EQ. '') Theta32_1(i__) = rNULL__
        if(ar__7(816:836) .EQ. '') Theta33_1(i__) = rNULL__
        if(ar__7(838:855) .EQ. '') Theta34_1(i__) = rNULL__
        if(ar__7(857:876) .EQ. '') Theta35_1(i__) = rNULL__
        if(ar__7(878:895) .EQ. '') Theta36_1(i__) = rNULL__
        if(ar__7(897:917) .EQ. '') Theta37_1(i__) = rNULL__
        if(ar__7(919:939) .EQ. '') Theta38_1(i__) = rNULL__
        if(ar__7(941:961) .EQ. '') Theta39_1(i__) = rNULL__
        if(ar__7(963:982) .EQ. '') Theta40_1(i__) = rNULL__
        if(ar__7(984:1001) .EQ. '') Theta41_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  BJD_7(i__),Flux_7(i__),e_Flux_7(i__),Fluxd_1(i__),
     +  e_Fluxd_1(i__),Fluxpsfmodel_1(i__),Theta0_1(i__),
     +  Theta1_1(i__),Theta2_1(i__),Theta3_1(i__),Theta4_1(i__),
     +  Theta5_1(i__),Theta6_1(i__),Theta7_1(i__),Theta8_1(i__),
     +  Theta9_1(i__),Theta10_1(i__),Theta11_1(i__),Theta12_1(i__),
     +  Theta13_1(i__),Theta14_1(i__),Theta15_1(i__),Theta16_1(i__),
     +  Theta17_1(i__),Theta18_1(i__),Theta19_1(i__),Theta20_1(i__),
     +  Theta21_1(i__),Theta22_1(i__),Theta23_1(i__),Theta24_1(i__),
     +  Theta25_1(i__),Theta26_1(i__),Theta27_1(i__),Theta28_1(i__),
     +  Theta29_1(i__),Theta30_1(i__),Theta31_1(i__),Theta32_1(i__),
     +  Theta33_1(i__),Theta34_1(i__),Theta35_1(i__),Theta36_1(i__),
     +  Theta37_1(i__),Theta38_1(i__),Theta39_1(i__),Theta40_1(i__),
     +  Theta41_1(i__)
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(i__),Flux_8(i__),e_Flux_8(i__),Fluxd_2(i__),
     +  e_Fluxd_2(i__),Fluxpsfmodel_2(i__),Theta0_2(i__),
     +  Theta1_2(i__),Theta2_2(i__),Theta3_2(i__),Theta4_2(i__),
     +  Theta5_2(i__),Theta6_2(i__),Theta7_2(i__),Theta8_2(i__),
     +  Theta9_2(i__),Theta10_2(i__),Theta11_2(i__),Theta12_2(i__),
     +  Theta13_2(i__),Theta14_2(i__),Theta15_2(i__),Theta16_2(i__),
     +  Theta17_2(i__),Theta18_2(i__),Theta19_2(i__),Theta20_2(i__),
     +  Theta21_2(i__),Theta22_2(i__),Theta23_2(i__),Theta24_2(i__),
     +  Theta25_2(i__),Theta26_2(i__),Theta27_2(i__),Theta28_2(i__),
     +  Theta29_2(i__),Theta30_2(i__),Theta31_2(i__),Theta32_2(i__),
     +  Theta33_2(i__),Theta34_2(i__),Theta35_2(i__),Theta36_2(i__),
     +  Theta37_2(i__),Theta38_2(i__),Theta39_2(i__),Theta40_2(i__),
     +  Theta41_2(i__)
        if(ar__8(320:339) .EQ. '') Theta10_2(i__) = rNULL__
        if(ar__8(341:361) .EQ. '') Theta11_2(i__) = rNULL__
        if(ar__8(362:382) .EQ. '') Theta12_2(i__) = rNULL__
        if(ar__8(384:404) .EQ. '') Theta13_2(i__) = rNULL__
        if(ar__8(406:426) .EQ. '') Theta14_2(i__) = rNULL__
        if(ar__8(428:448) .EQ. '') Theta15_2(i__) = rNULL__
        if(ar__8(449:469) .EQ. '') Theta16_2(i__) = rNULL__
        if(ar__8(471:491) .EQ. '') Theta17_2(i__) = rNULL__
        if(ar__8(493:512) .EQ. '') Theta18_2(i__) = rNULL__
        if(ar__8(514:534) .EQ. '') Theta19_2(i__) = rNULL__
        if(ar__8(536:556) .EQ. '') Theta20_2(i__) = rNULL__
        if(ar__8(558:578) .EQ. '') Theta21_2(i__) = rNULL__
        if(ar__8(580:599) .EQ. '') Theta22_2(i__) = rNULL__
        if(ar__8(601:621) .EQ. '') Theta23_2(i__) = rNULL__
        if(ar__8(623:643) .EQ. '') Theta24_2(i__) = rNULL__
        if(ar__8(645:665) .EQ. '') Theta25_2(i__) = rNULL__
        if(ar__8(667:686) .EQ. '') Theta26_2(i__) = rNULL__
        if(ar__8(688:707) .EQ. '') Theta27_2(i__) = rNULL__
        if(ar__8(709:729) .EQ. '') Theta28_2(i__) = rNULL__
        if(ar__8(731:751) .EQ. '') Theta29_2(i__) = rNULL__
        if(ar__8(753:773) .EQ. '') Theta30_2(i__) = rNULL__
        if(ar__8(775:795) .EQ. '') Theta31_2(i__) = rNULL__
        if(ar__8(797:814) .EQ. '') Theta32_2(i__) = rNULL__
        if(ar__8(816:836) .EQ. '') Theta33_2(i__) = rNULL__
        if(ar__8(838:855) .EQ. '') Theta34_2(i__) = rNULL__
        if(ar__8(857:876) .EQ. '') Theta35_2(i__) = rNULL__
        if(ar__8(878:895) .EQ. '') Theta36_2(i__) = rNULL__
        if(ar__8(897:917) .EQ. '') Theta37_2(i__) = rNULL__
        if(ar__8(919:939) .EQ. '') Theta38_2(i__) = rNULL__
        if(ar__8(941:961) .EQ. '') Theta39_2(i__) = rNULL__
        if(ar__8(963:982) .EQ. '') Theta40_2(i__) = rNULL__
        if(ar__8(984:1001) .EQ. '') Theta41_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  BJD_8(i__),Flux_8(i__),e_Flux_8(i__),Fluxd_2(i__),
     +  e_Fluxd_2(i__),Fluxpsfmodel_2(i__),Theta0_2(i__),
     +  Theta1_2(i__),Theta2_2(i__),Theta3_2(i__),Theta4_2(i__),
     +  Theta5_2(i__),Theta6_2(i__),Theta7_2(i__),Theta8_2(i__),
     +  Theta9_2(i__),Theta10_2(i__),Theta11_2(i__),Theta12_2(i__),
     +  Theta13_2(i__),Theta14_2(i__),Theta15_2(i__),Theta16_2(i__),
     +  Theta17_2(i__),Theta18_2(i__),Theta19_2(i__),Theta20_2(i__),
     +  Theta21_2(i__),Theta22_2(i__),Theta23_2(i__),Theta24_2(i__),
     +  Theta25_2(i__),Theta26_2(i__),Theta27_2(i__),Theta28_2(i__),
     +  Theta29_2(i__),Theta30_2(i__),Theta31_2(i__),Theta32_2(i__),
     +  Theta33_2(i__),Theta34_2(i__),Theta35_2(i__),Theta36_2(i__),
     +  Theta37_2(i__),Theta38_2(i__),Theta39_2(i__),Theta40_2(i__),
     +  Theta41_2(i__)
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(i__),Flux_9(i__),e_Flux_9(i__),Fluxd_3(i__),
     +  e_Fluxd_3(i__),Fluxpsfmodel_3(i__),Theta0_3(i__),
     +  Theta1_3(i__),Theta2_3(i__),Theta3_3(i__),Theta4_3(i__),
     +  Theta5_3(i__),Theta6_3(i__),Theta7_3(i__),Theta8_3(i__),
     +  Theta9_3(i__),Theta10_3(i__),Theta11_3(i__),Theta12_3(i__),
     +  Theta13_3(i__),Theta14_3(i__),Theta15_3(i__),Theta16_3(i__),
     +  Theta17_3(i__),Theta18_3(i__),Theta19_3(i__),Theta20_3(i__),
     +  Theta21_3(i__),Theta22_3(i__),Theta23_3(i__),Theta24_3(i__),
     +  Theta25_3(i__),Theta26_3(i__),Theta27_3(i__),Theta28_3(i__),
     +  Theta29_3(i__),Theta30_3(i__),Theta31_3(i__),Theta32_3(i__),
     +  Theta33_3(i__),Theta34_3(i__),Theta35_3(i__),Theta36_3(i__),
     +  Theta37_3(i__),Theta38_3(i__),Theta39_3(i__),Theta40_3(i__),
     +  Theta41_3(i__)
        if(ar__9(320:339) .EQ. '') Theta10_3(i__) = rNULL__
        if(ar__9(341:361) .EQ. '') Theta11_3(i__) = rNULL__
        if(ar__9(362:382) .EQ. '') Theta12_3(i__) = rNULL__
        if(ar__9(384:404) .EQ. '') Theta13_3(i__) = rNULL__
        if(ar__9(406:426) .EQ. '') Theta14_3(i__) = rNULL__
        if(ar__9(428:448) .EQ. '') Theta15_3(i__) = rNULL__
        if(ar__9(449:469) .EQ. '') Theta16_3(i__) = rNULL__
        if(ar__9(471:491) .EQ. '') Theta17_3(i__) = rNULL__
        if(ar__9(493:512) .EQ. '') Theta18_3(i__) = rNULL__
        if(ar__9(514:534) .EQ. '') Theta19_3(i__) = rNULL__
        if(ar__9(536:556) .EQ. '') Theta20_3(i__) = rNULL__
        if(ar__9(558:578) .EQ. '') Theta21_3(i__) = rNULL__
        if(ar__9(580:599) .EQ. '') Theta22_3(i__) = rNULL__
        if(ar__9(601:621) .EQ. '') Theta23_3(i__) = rNULL__
        if(ar__9(623:643) .EQ. '') Theta24_3(i__) = rNULL__
        if(ar__9(645:665) .EQ. '') Theta25_3(i__) = rNULL__
        if(ar__9(667:686) .EQ. '') Theta26_3(i__) = rNULL__
        if(ar__9(688:707) .EQ. '') Theta27_3(i__) = rNULL__
        if(ar__9(709:729) .EQ. '') Theta28_3(i__) = rNULL__
        if(ar__9(731:751) .EQ. '') Theta29_3(i__) = rNULL__
        if(ar__9(753:773) .EQ. '') Theta30_3(i__) = rNULL__
        if(ar__9(775:795) .EQ. '') Theta31_3(i__) = rNULL__
        if(ar__9(797:814) .EQ. '') Theta32_3(i__) = rNULL__
        if(ar__9(816:836) .EQ. '') Theta33_3(i__) = rNULL__
        if(ar__9(838:855) .EQ. '') Theta34_3(i__) = rNULL__
        if(ar__9(857:876) .EQ. '') Theta35_3(i__) = rNULL__
        if(ar__9(878:895) .EQ. '') Theta36_3(i__) = rNULL__
        if(ar__9(897:917) .EQ. '') Theta37_3(i__) = rNULL__
        if(ar__9(919:939) .EQ. '') Theta38_3(i__) = rNULL__
        if(ar__9(941:961) .EQ. '') Theta39_3(i__) = rNULL__
        if(ar__9(963:982) .EQ. '') Theta40_3(i__) = rNULL__
        if(ar__9(984:1001) .EQ. '') Theta41_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,10)
     +  BJD_9(i__),Flux_9(i__),e_Flux_9(i__),Fluxd_3(i__),
     +  e_Fluxd_3(i__),Fluxpsfmodel_3(i__),Theta0_3(i__),
     +  Theta1_3(i__),Theta2_3(i__),Theta3_3(i__),Theta4_3(i__),
     +  Theta5_3(i__),Theta6_3(i__),Theta7_3(i__),Theta8_3(i__),
     +  Theta9_3(i__),Theta10_3(i__),Theta11_3(i__),Theta12_3(i__),
     +  Theta13_3(i__),Theta14_3(i__),Theta15_3(i__),Theta16_3(i__),
     +  Theta17_3(i__),Theta18_3(i__),Theta19_3(i__),Theta20_3(i__),
     +  Theta21_3(i__),Theta22_3(i__),Theta23_3(i__),Theta24_3(i__),
     +  Theta25_3(i__),Theta26_3(i__),Theta27_3(i__),Theta28_3(i__),
     +  Theta29_3(i__),Theta30_3(i__),Theta31_3(i__),Theta32_3(i__),
     +  Theta33_3(i__),Theta34_3(i__),Theta35_3(i__),Theta36_3(i__),
     +  Theta37_3(i__),Theta38_3(i__),Theta39_3(i__),Theta40_3(i__),
     +  Theta41_3(i__)
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(i__),Flux_10(i__),e_Flux_10(i__),Fluxd_4(i__),
     +  e_Fluxd_4(i__),Fluxpsfmodel_4(i__),Theta0_4(i__),
     +  Theta1_4(i__),Theta2_4(i__),Theta3_4(i__),Theta4_4(i__),
     +  Theta5_4(i__),Theta6_4(i__),Theta7_4(i__),Theta8_4(i__),
     +  Theta9_4(i__),Theta10_4(i__),Theta11_4(i__),Theta12_4(i__),
     +  Theta13_4(i__),Theta14_4(i__),Theta15_4(i__),Theta16_4(i__),
     +  Theta17_4(i__),Theta18_4(i__),Theta19_4(i__),Theta20_4(i__),
     +  Theta21_4(i__),Theta22_4(i__),Theta23_4(i__),Theta24_4(i__),
     +  Theta25_4(i__),Theta26_4(i__),Theta27_4(i__),Theta28_4(i__),
     +  Theta29_4(i__),Theta30_4(i__),Theta31_4(i__),Theta32_4(i__),
     +  Theta33_4(i__),Theta34_4(i__),Theta35_4(i__),Theta36_4(i__),
     +  Theta37_4(i__),Theta38_4(i__),Theta39_4(i__),Theta40_4(i__),
     +  Theta41_4(i__)
        if(ar__10(320:339) .EQ. '') Theta10_4(i__) = rNULL__
        if(ar__10(341:361) .EQ. '') Theta11_4(i__) = rNULL__
        if(ar__10(362:382) .EQ. '') Theta12_4(i__) = rNULL__
        if(ar__10(384:404) .EQ. '') Theta13_4(i__) = rNULL__
        if(ar__10(406:426) .EQ. '') Theta14_4(i__) = rNULL__
        if(ar__10(428:448) .EQ. '') Theta15_4(i__) = rNULL__
        if(ar__10(449:469) .EQ. '') Theta16_4(i__) = rNULL__
        if(ar__10(471:491) .EQ. '') Theta17_4(i__) = rNULL__
        if(ar__10(493:512) .EQ. '') Theta18_4(i__) = rNULL__
        if(ar__10(514:534) .EQ. '') Theta19_4(i__) = rNULL__
        if(ar__10(536:556) .EQ. '') Theta20_4(i__) = rNULL__
        if(ar__10(558:578) .EQ. '') Theta21_4(i__) = rNULL__
        if(ar__10(580:599) .EQ. '') Theta22_4(i__) = rNULL__
        if(ar__10(601:621) .EQ. '') Theta23_4(i__) = rNULL__
        if(ar__10(623:643) .EQ. '') Theta24_4(i__) = rNULL__
        if(ar__10(645:665) .EQ. '') Theta25_4(i__) = rNULL__
        if(ar__10(667:686) .EQ. '') Theta26_4(i__) = rNULL__
        if(ar__10(688:707) .EQ. '') Theta27_4(i__) = rNULL__
        if(ar__10(709:729) .EQ. '') Theta28_4(i__) = rNULL__
        if(ar__10(731:751) .EQ. '') Theta29_4(i__) = rNULL__
        if(ar__10(753:773) .EQ. '') Theta30_4(i__) = rNULL__
        if(ar__10(775:795) .EQ. '') Theta31_4(i__) = rNULL__
        if(ar__10(797:814) .EQ. '') Theta32_4(i__) = rNULL__
        if(ar__10(816:836) .EQ. '') Theta33_4(i__) = rNULL__
        if(ar__10(838:855) .EQ. '') Theta34_4(i__) = rNULL__
        if(ar__10(857:876) .EQ. '') Theta35_4(i__) = rNULL__
        if(ar__10(878:895) .EQ. '') Theta36_4(i__) = rNULL__
        if(ar__10(897:917) .EQ. '') Theta37_4(i__) = rNULL__
        if(ar__10(919:939) .EQ. '') Theta38_4(i__) = rNULL__
        if(ar__10(941:961) .EQ. '') Theta39_4(i__) = rNULL__
        if(ar__10(963:982) .EQ. '') Theta40_4(i__) = rNULL__
        if(ar__10(984:1001) .EQ. '') Theta41_4(i__) = rNULL__
c    ..............Just test output...........
        write(6,11)
     +  BJD_10(i__),Flux_10(i__),e_Flux_10(i__),Fluxd_4(i__),
     +  e_Fluxd_4(i__),Fluxpsfmodel_4(i__),Theta0_4(i__),
     +  Theta1_4(i__),Theta2_4(i__),Theta3_4(i__),Theta4_4(i__),
     +  Theta5_4(i__),Theta6_4(i__),Theta7_4(i__),Theta8_4(i__),
     +  Theta9_4(i__),Theta10_4(i__),Theta11_4(i__),Theta12_4(i__),
     +  Theta13_4(i__),Theta14_4(i__),Theta15_4(i__),Theta16_4(i__),
     +  Theta17_4(i__),Theta18_4(i__),Theta19_4(i__),Theta20_4(i__),
     +  Theta21_4(i__),Theta22_4(i__),Theta23_4(i__),Theta24_4(i__),
     +  Theta25_4(i__),Theta26_4(i__),Theta27_4(i__),Theta28_4(i__),
     +  Theta29_4(i__),Theta30_4(i__),Theta31_4(i__),Theta32_4(i__),
     +  Theta33_4(i__),Theta34_4(i__),Theta35_4(i__),Theta36_4(i__),
     +  Theta37_4(i__),Theta38_4(i__),Theta39_4(i__),Theta40_4(i__),
     +  Theta41_4(i__)
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(i__),Flux_11(i__),e_Flux_11(i__),Fluxd_5(i__),
     +  e_Fluxd_5(i__),Fluxpsfmodel_5(i__),Theta0_5(i__),
     +  Theta1_5(i__),Theta2_5(i__),Theta3_5(i__),Theta4_5(i__),
     +  Theta5_5(i__),Theta6_5(i__),Theta7_5(i__),Theta8_5(i__),
     +  Theta9_5(i__),Theta10_5(i__),Theta11_5(i__),Theta12_5(i__),
     +  Theta13_5(i__),Theta14_5(i__),Theta15_5(i__),Theta16_5(i__),
     +  Theta17_5(i__),Theta18_5(i__),Theta19_5(i__),Theta20_5(i__),
     +  Theta21_5(i__),Theta22_5(i__),Theta23_5(i__),Theta24_5(i__),
     +  Theta25_5(i__),Theta26_5(i__),Theta27_5(i__),Theta28_5(i__),
     +  Theta29_5(i__),Theta30_5(i__),Theta31_5(i__),Theta32_5(i__),
     +  Theta33_5(i__),Theta34_5(i__),Theta35_5(i__),Theta36_5(i__),
     +  Theta37_5(i__),Theta38_5(i__),Theta39_5(i__),Theta40_5(i__),
     +  Theta41_5(i__)
        if(ar__11(320:339) .EQ. '') Theta10_5(i__) = rNULL__
        if(ar__11(341:361) .EQ. '') Theta11_5(i__) = rNULL__
        if(ar__11(362:382) .EQ. '') Theta12_5(i__) = rNULL__
        if(ar__11(384:404) .EQ. '') Theta13_5(i__) = rNULL__
        if(ar__11(406:426) .EQ. '') Theta14_5(i__) = rNULL__
        if(ar__11(428:448) .EQ. '') Theta15_5(i__) = rNULL__
        if(ar__11(449:469) .EQ. '') Theta16_5(i__) = rNULL__
        if(ar__11(471:491) .EQ. '') Theta17_5(i__) = rNULL__
        if(ar__11(493:512) .EQ. '') Theta18_5(i__) = rNULL__
        if(ar__11(514:534) .EQ. '') Theta19_5(i__) = rNULL__
        if(ar__11(536:556) .EQ. '') Theta20_5(i__) = rNULL__
        if(ar__11(558:578) .EQ. '') Theta21_5(i__) = rNULL__
        if(ar__11(580:599) .EQ. '') Theta22_5(i__) = rNULL__
        if(ar__11(601:621) .EQ. '') Theta23_5(i__) = rNULL__
        if(ar__11(623:643) .EQ. '') Theta24_5(i__) = rNULL__
        if(ar__11(645:665) .EQ. '') Theta25_5(i__) = rNULL__
        if(ar__11(667:686) .EQ. '') Theta26_5(i__) = rNULL__
        if(ar__11(688:707) .EQ. '') Theta27_5(i__) = rNULL__
        if(ar__11(709:729) .EQ. '') Theta28_5(i__) = rNULL__
        if(ar__11(731:751) .EQ. '') Theta29_5(i__) = rNULL__
        if(ar__11(753:773) .EQ. '') Theta30_5(i__) = rNULL__
        if(ar__11(775:795) .EQ. '') Theta31_5(i__) = rNULL__
        if(ar__11(797:814) .EQ. '') Theta32_5(i__) = rNULL__
        if(ar__11(816:836) .EQ. '') Theta33_5(i__) = rNULL__
        if(ar__11(838:855) .EQ. '') Theta34_5(i__) = rNULL__
        if(ar__11(857:876) .EQ. '') Theta35_5(i__) = rNULL__
        if(ar__11(878:895) .EQ. '') Theta36_5(i__) = rNULL__
        if(ar__11(897:917) .EQ. '') Theta37_5(i__) = rNULL__
        if(ar__11(919:939) .EQ. '') Theta38_5(i__) = rNULL__
        if(ar__11(941:961) .EQ. '') Theta39_5(i__) = rNULL__
        if(ar__11(963:982) .EQ. '') Theta40_5(i__) = rNULL__
        if(ar__11(984:1001) .EQ. '') Theta41_5(i__) = rNULL__
c    ..............Just test output...........
        write(6,12)
     +  BJD_11(i__),Flux_11(i__),e_Flux_11(i__),Fluxd_5(i__),
     +  e_Fluxd_5(i__),Fluxpsfmodel_5(i__),Theta0_5(i__),
     +  Theta1_5(i__),Theta2_5(i__),Theta3_5(i__),Theta4_5(i__),
     +  Theta5_5(i__),Theta6_5(i__),Theta7_5(i__),Theta8_5(i__),
     +  Theta9_5(i__),Theta10_5(i__),Theta11_5(i__),Theta12_5(i__),
     +  Theta13_5(i__),Theta14_5(i__),Theta15_5(i__),Theta16_5(i__),
     +  Theta17_5(i__),Theta18_5(i__),Theta19_5(i__),Theta20_5(i__),
     +  Theta21_5(i__),Theta22_5(i__),Theta23_5(i__),Theta24_5(i__),
     +  Theta25_5(i__),Theta26_5(i__),Theta27_5(i__),Theta28_5(i__),
     +  Theta29_5(i__),Theta30_5(i__),Theta31_5(i__),Theta32_5(i__),
     +  Theta33_5(i__),Theta34_5(i__),Theta35_5(i__),Theta36_5(i__),
     +  Theta37_5(i__),Theta38_5(i__),Theta39_5(i__),Theta40_5(i__),
     +  Theta41_5(i__)
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(i__),Flux_12(i__),e_Flux_12(i__)
c    ..............Just test output...........
        write(6,13)BJD_12(i__),Flux_12(i__),e_Flux_12(i__)
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(i__),Flux_13(i__),e_Flux_13(i__)
c    ..............Just test output...........
        write(6,14)BJD_13(i__),Flux_13(i__),e_Flux_13(i__)
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(i__),Flux_14(i__),e_Flux_14(i__)
c    ..............Just test output...........
        write(6,15)BJD_14(i__),Flux_14(i__),e_Flux_14(i__)
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(i__),Flux_15(i__),e_Flux_15(i__)
c    ..............Just test output...........
        write(6,16)BJD_15(i__),Flux_15(i__),e_Flux_15(i__)
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(i__),Flux_16(i__),e_Flux_16(i__)
c    ..............Just test output...........
        write(6,17)BJD_16(i__),Flux_16(i__),e_Flux_16(i__)
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(i__),Flux_17(i__),e_Flux_17(i__)
c    ..............Just test output...........
        write(6,18)BJD_17(i__),Flux_17(i__),e_Flux_17(i__)
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(i__),Flux_18(i__),e_Flux_18(i__)
c    ..............Just test output...........
        write(6,19)BJD_18(i__),Flux_18(i__),e_Flux_18(i__)
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(i__),Flux_19(i__),e_Flux_19(i__)
c    ..............Just test output...........
        write(6,20)BJD_19(i__),Flux_19(i__),e_Flux_19(i__)
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(i__),Flux_20(i__),e_Flux_20(i__)
c    ..............Just test output...........
        write(6,21)BJD_20(i__),Flux_20(i__),e_Flux_20(i__)
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(i__),Flux_21(i__),e_Flux_21(i__)
c    ..............Just test output...........
        write(6,22)BJD_21(i__),Flux_21(i__),e_Flux_21(i__)
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(i__),RV(i__),e_RV(i__),BIS(i__),FWHM(i__),
     +  Contrast(i__),Sindex(i__),Halpha(i__),NaD1(i__),NaD2(i__)
c    ..............Just test output...........
        write(6,23)
     +  BJD_22(i__),RV(i__),e_RV(i__),BIS(i__),FWHM(i__),
     +  Contrast(i__),Sindex(i__),Halpha(i__),NaD1(i__),NaD2(i__)
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(i__),RV_1(i__),e_RV_1(i__),Theta0_6(i__),Theta1_6(i__),
     +  Theta2_6(i__)
c    ..............Just test output...........
        write(6,24)
     +  BJD_23(i__),RV_1(i__),e_RV_1(i__),Theta0_6(i__),Theta1_6(i__),
     +  Theta2_6(i__)
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(i__),Flux_22(i__),e_Flux_22(i__)
c    ..............Just test output...........
        write(6,25)BJD_24(i__),Flux_22(i__),e_Flux_22(i__)
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(i__),Flux_23(i__),e_Flux_23(i__)
c    ..............Just test output...........
        write(6,26)BJD_25(i__),Flux_23(i__),e_Flux_23(i__)
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(i__),Flux_24(i__),e_Flux_24(i__)
c    ..............Just test output...........
        write(6,27)BJD_26(i__),Flux_24(i__),e_Flux_24(i__)
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(i__),Flux_25(i__),e_Flux_25(i__)
c    ..............Just test output...........
        write(6,28)BJD_27(i__),Flux_25(i__),e_Flux_25(i__)
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(i__),Flux_26(i__),e_Flux_26(i__)
c    ..............Just test output...........
        write(6,29)BJD_28(i__),Flux_26(i__),e_Flux_26(i__)
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(i__),Flux_27(i__),e_Flux_27(i__),Airmass(i__),
     +  SkyBkg(i__),X(i__),Y(i__)
c    ..............Just test output...........
        write(6,30)
     +  BJD_29(i__),Flux_27(i__),e_Flux_27(i__),Airmass(i__),
     +  SkyBkg(i__),X(i__),Y(i__)
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(i__),Flux_28(i__),e_Flux_28(i__)
c    ..............Just test output...........
        write(6,31)BJD_30(i__),Flux_28(i__),e_Flux_28(i__)
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(i__),Flux_29(i__),e_Flux_29(i__)
c    ..............Just test output...........
        write(6,32)BJD_31(i__),Flux_29(i__),e_Flux_29(i__)
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(i__),Flux_30(i__),e_Flux_30(i__)
c    ..............Just test output...........
        write(6,33)BJD_32(i__),Flux_30(i__),e_Flux_30(i__)
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(i__),Flux_31(i__),e_Flux_31(i__)
c    ..............Just test output...........
        write(6,34)BJD_33(i__),Flux_31(i__),e_Flux_31(i__)
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(i__),Flux_32(i__),e_Flux_32(i__)
c    ..............Just test output...........
        write(6,35)BJD_34(i__),Flux_32(i__),e_Flux_32(i__)
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(i__),Flux_33(i__),e_Flux_33(i__)
c    ..............Just test output...........
        write(6,36)BJD_35(i__),Flux_33(i__),e_Flux_33(i__)
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(i__),Flux_34(i__),e_Flux_34(i__)
c    ..............Just test output...........
        write(6,37)BJD_36(i__),Flux_34(i__),e_Flux_34(i__)
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(i__),Sep(i__),Deltamag(i__)
c    ..............Just test output...........
        write(6,38)lambda(i__),Sep(i__),Deltamag(i__)
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(i__),Sep_1(i__),Deltamag_1(i__)
c    ..............Just test output...........
        write(6,39)lambda_1(i__),Sep_1(i__),Deltamag_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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