powered by
Calculate vertex depth given ray parameter or vice-versa.
ConvP2Vdepth(p, v, r, h, rp, discons) ConvP2Vdepthinv(rpd, v, r)
Ray vertex radius (km)
Planet velocity structure (km/s)
Radii corresponding to v
Ray parameter (s/deg)
Focal radius (km)
Planet radius (km)
Vector of discontinuity radii (km, from FindDiscon)
ConvP2Vdepth: Radius of ray vertex (km)
ConvP2Vdepthinv: Ray parameter (s/deg)
Note that these functions use radii, not depths, so h would be 6371 (or whatever planet radius you're using) - focal depth.
# NOT RUN { data(model) ConvP2Vdepth(7, model$vp, 6371 - model$z, 6361, 6371, FindDiscon(model)) ConvP2Vdepthinv(4881.467, model$vp, 6371 - model$z) # }
Run the code above in your browser using DataLab