---C. Beghein (08/2026)---

In each depth subfolder, you will find two files: dcg and dap. dcg contains the isotropic part of the velocity model. dap contains the 2-psi anisotropy. 
plot_2psi.sh is an example of GMT script to plot the model 

-dcg has the following content:
col1: longitude (deg)
col2: latitude (deg)
col3: dV (m/s)
col4: dV/V_ref (%)  (where V_ref is PREM in this case)

-dap contains:
col1: longitude (deg)
col2: latitude (deg)
col3: fast-axis azimuth (deg, measured from geographic North, clockwise)
col4: bar length for GMT (km) (=psol/psolmax)·ddd/1.8
col5: constant 1.0 — bar width (km)
col6: 2-psi amplitude (m/s)    = sqrt(a^2+b^2)
col7: relative 2-psi amplitude (% of pvref)

Columns 1–5 of dap are GMT psxy -SE arguments (lon, lat, azimuth, major, minor)
