**************************************************************************
Repository of data supporting the results in Daher et al. (JGR, submitted)
**************************************************************************

** CONTENT of this part of the repository:

   README_Figures6to11       This file

   Files containing the results needed to plot Figures 6-11:
   integration_results_use_Schindelegger_116Ma_experiments_ode45_nodeLF_v9.mat
   integration_results_use_Schindelegger_252Ma_experiments_ode45_nodeLF_v9.mat
   integration_results_use_Schindelegger_55Ma_experiments_ode45_nodeLF_v9.mat
   integration_results_use_Schindelegger_PD_experiments_ode45_nodeLF_v9.mat
   ExtraDiagnostics_ode45_Schindelegger_results_nodeLF_v9.mat
   WriteOut_d_dt_terms_ode45_Schindelegger_results_nodeLF_v9.mat
   SchindeleggerResults_MonteCarlo
		This last entry is a directory with results from 1000 Monte Carlo simulations

   The first four "integration_results..." files contain the timesteps and the StateVectors
   of the Earth-Moon evolution trajectories that employ simulations of the ocean tides with
   Schindelegger's model and four different ocean basin geometries--PD (present-day), 55Ma,
   116Ma, and 252Ma.  The ocean tide simulations employ Earth rotation rates ranging from
   6-24 hours.  The timesteps in the "integration_results" files have units of seconds and
   should be divided by 86400*365.25*10^9 to obtain units of Ga.  The StateVector variables
   consist of 
   (1) sidereal Earth rotation rate (s^-1).  To obtain Earth rotation period in hours in a 
        Matlab code, type:  (1/3600)*2*pi./StateVector(:,1) 
   (2) obliquity (radians).  To obtain obliquity in degrees in a Matlab code: 
        (180/pi)*StateVector(:,2)   
   (3) semimajor axis (meters).  To obtain semimajor axis in Earth radii in a Matlab code:
        (1/6378136)*StateVector(:,3)
   (4) eccentricty (dimensionless)
   (5) inclination (radians).  To obtain inclination in degrees in a Matlab code:
        (180/pi)*StateVector(:,5)          

   The "ExtraDiagnostics" file contains "ExtraDiagnostics" structures for the Earth-Moon evolution
   results that employ the four different ocean basin geometries, as described above.  The
   "ExtraDiagnostics" structures include the timesteps (units seconds) as well as several
   other variables:
   n (lunar mean motion, units s^-1):  (1/86400)*2*pi./ExtraDiagnostics.n_ode45 gives the 
      period of the lunar mean motion in days 
   domegabardt (rate of change of longitude of perigee, units s^-1):  
      (1/(86400*365.25))*2*pi./ExtraDiagnostics.domegabar_dt_ode45 gives the period in years
   dOmega_dt (nodal rate of change, units s^-1):  
      -(1/(86400*365.25))*2*pi./ExtraDiagnostics.dOmega_dt_ode45 gives the period in years;
      negative sign added to make the results positive
   xi_Moon_coremantledissipation:  dimensionless parameter xi (see text)
   KoverC_Moon_coremantledissipation:  as with xi, this is another key parameter in the luner core-mantle boundary (CMB) terms (see text) 
   I (lunar equatorial tilt, units radians):  (180/pi)*ExtraDiagnostics.I_lunartilt_ode45
      gives the lunar equatorial tilt in degrees
   dpsi_dt (precession period, units of s^-1):  
      (1/(86400*365.25*1000))*2*pi./ExtraDiagnostics.dpsi_dt_ode45 gives the precession period
      in kiloyears
   ksinchi_diurnal (diurnal ksinchi values, dimensionless):  
      ExtraDiagnostics.ksinchi_diurnal_ode45 gives the diurnal ksinchi values
   ksinchi_semidiurnal (semi-diurnal ksinchi values, dimensionless):  
      ExtraDiagnostics.ksinchi_semidiurnal_ode45 gives the semi-diurnal ksinchi values
   M2_powerdissipation_ode45 (M2 power dissipation, units Watts):  
      (1/(10^12))*ExtraDiagnostics.M2_powerdissipation_ode45 gives M2 dissipation in Terrawatts
   O1_powerdissipation_ode45 (O1 power dissipation, units Watts):  
      (1/(10^12))*ExtraDiagnostics.O1_powerdissipation_ode45 gives O1 dissipation in Terrawatts
   


   The "WriteOut" file contains additional diagnostics for the Earth-Moon trajectories that employ
   ocean tide simulations with the four ocean basin geometries.  These diagnostics are used in 
   Figure 9, which plots d/dt of a, i, and e due to the tides on Earth and to tides and core-mantle
   dissipation within the Moon.  These diagnostics are also used in Figure 10, which plots torques.        
 
   The Matlab files
       plot_Figure6.m
       plot_Figures_7_8_9_11.m
       plot_Figure10.m
       create Figures 6-11 of the paper, e.g., the figures representing the main orbital dynamics results

----------------------------------------------------------------------------

** REFERENCE:
   Daher H., Arbic B.K, Williams J.G., Ansong J.K., Boggs D.H., et al. (2021),
   Long-term Earth-Moon evolution with high-level orbit and ocean tide models,
   Journal of Geophysical Research Planets, submitted.

** CONTACT:
   Brian K. Arbic, arbic@umich.edu (Figures 6-11)
   Michael Schindelegger, schindelegger@igg.uni-bonn.de