./paper_data_code/ contains:
------------------------------------------------------------------------
----(1) ocs.csv---------------------------------------------------------
OCs used in the paper (from Gaia DR3, member star matched 
radius: 0.1 arc-second).
citation catalogs:
cg20:Cantat-Gaudin et al. (2020)
castro22:Castro-Ginard et al. (2022)
he21b:He et al. (2021b)
he22:He et al. (2022)
he23a:He et al. (2023a)
he23b:He et al. (2023b)
------------------------------------------------------------------------
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 inc1)
dz0                      (differential scale height, see Sec.3 in paper)
disp_dz0                 (uncertainty(sample) in inc1)
------------------------------------------------------------------------
----(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)
------------------------------------------------------------------------
------------------------------------------------------------------------