================================================================================
Title: The precession and inclinational variation across the plane of
    the Milky Way revealed by open clusters
Authors: He Z.
================================================================================
Description of contents: This tar.gz archive contains the data and code 
    necessary to reproduce the results in the accepted version of the
    ApJ Letter listed above. The specific files included in this deposit
    are documented below. 

    ----(1) ocs.csv---------------------------------------------------------
    Open Clusters used in the paper (from Gaia DR3, member star matched 
        radius: 0.1 arc-second). Citation catalogs:
        cg20     = Cantat-Gaudin et al. (2020) [2020A&A...640A...1C]
        castro22 = Castro-Ginard et al. (2022) [2022A&A...661A.118C]
        he21b    = He et al. (2021b) [2021RAA....21...93H]
        he22     = He et al. (2022) [2022ApJS..260....8H]
        he23a    = He et al. (2023a) [2023ApJS..264....8H]
        he23b    = He et al. (2023b) [2023arXiv230510269H]
    ------------------------------------------------------------------------
    Cluster                  (Cluster ID)
    Table                    (citation catalog)
    l                        (weighted mean Galactic longitude)
    b                        (weighted mean Galactic  latitude)
    n                        (number of member stars)
    plx                      (mean parallax)
    e_plx                    (parallax error)
    pmx                      (mean proper motion in R.A.*)
    e_pmx                    (pmRA* dispersion)
    pmy                      (mean proper motion in Dec.) 
    e_pmy                    (pmDEC error)
    age                      (logarithmic age ) 
    rv                       (radial velocity) 
    e_RV                     (RV dispersion when NRV>1, or 
                              radial_velocity_error when NRV=1)
    nrv                      (number of stars with RV values, 
                              and radial_velocity_error<10 km/s)
    pml                      (mean proper motion in Galactic longitude)
    pmb                      (mean proper motion in Galactic latitude)
    xx                       (X-axis of the Cassegrain coordinate system)
    yy                       (Y-axis of the Cassegrain coordinate system)
    zz                       (scale height)
    gcphi                    (\phi angel of the GC coordinate system)
    gcr                      (GC radii)
    gctheta                  (\theta angel of the GC coordinate system)
    vphi                     (tangent velocity in \phi ange)
    e_vphi                   (error in vphi)
    vgcr                     (velocity along GC radius)
    e_vgcr                   (error in vgcr)
    vz                       (vertical velocity)
    e_vz                     (error in vz)
    ------------------------------------------------------------------------
    ----(2) results.txt-----------------------------------------------------
    Full Figure 4 data in the paper: 
    ------------------------------------------------------------------------
    rgc                      (GC radii, in bin of +- 1~kpc)
    inc1                     (initial inclination angle, from anti-GC)
    disp_inc1                (uncertainty(sample) in inc1)
    theta_z0                 (angle of differential scale height)
    disp_theta_z0            (uncertainty(sample) in theta_z0)
    phi_lon                  (LON position)
    disp_phi_lon             (uncertainty(sample) in phi_lon)
    err_phi_lon              (uncertainty(vz_sun) in phi_lon)
    wz_p                     (precession rate)
    disp_wz_p                (uncertainty(sample) in wz_p)
    err_wz_p                 (uncertainty(vz_sun) in wz_p)
    wz_n                     (inclinational variation rate)
    disp_wz_n                (uncertainty(sample) in wz_n)
    err_wz_n                 (uncertainty(vz_sun) in wz_n)
    inc                      (inclination angle, from LON)
    disp_inc                 (uncertainty(sample) in inc)
    dz0                      (differential scale height, see Sec.3 in paper)
    disp_dz0                 (uncertainty(sample) in dz0)
    ------------------------------------------------------------------------
    ----(3) ./files/inc_pre_nut_*kpc.csv------------------------------------
        Intermediate process files generated during fitting:
        4-14 kpc fits from code calculate.py
        And read by read_calculate.py
        See explanatory notes in the code 
    ------------------------------------------------------------------------
    max_age                  (logarithmic age limit of the OC sub-sample)
    n                        (OC number of sub-sample)
    inc1                     (initial inclination angle, from anti-GC)
    theta_z0                 (angle of differential scale height)
    phi_lon                  (LON position)
    disp_phi_lon             (uncertainty(vz_sun) in phi_lon)
    wz_p                     (precession rate)
    disp_wz_p                (uncertainty(vz_sun) in wz_p)
    wz_n                     (inclinational variation rate)
    disp_wz_n                (uncertainty(vz_sun) in wz_n)
    inc                      (inclination angle, from LON)
    dz0                      (differential scale height)
    ------------------------------------------------------------------------
    ------------------------------------------------------------------------

System requirements: Python packages: scipy, numpy, pandas

Additional comments: If the reader has any questions, please contact the code
author: Zhihong He (hezh@mail.ustc.edu.cn)
================================================================================