I’m exploring your axisymmetric code to model a simple spherically-focused disc transducer we’re building, but the p_min and p_max outputs don’t seem to align with what I see in p_final.
Here's a link to a Matlab figure showing p_final alongside p_max and p_min. Apparently p_max and p_min pressures are shifted left, and therefore do not display the focal spot pressures. Near the transducer at the top of the image, you can see where the pressures should lie on top of the transducer, but also appear shifted left:
https://fusf.sharepoint.com/:i:/s/FUSF-Development/Ebn2GZi9SX5Khd7gh80_L4YBRqmwubNlQBw1ZVhP7aEi_Q?e=cgvlb1
Simply altering the kspaceFirstOrderAS to kspaceFirstOrder2D yields pressure results that align as expected. Happy to provide code or 2D screen-captures if those are helpful.