Custom astropy unit for PSI radial coordinate grids.
One PSI_rsun equals RSUN = 6.96 × 10¹⁰ cm. Radial coordinate arrays
in MAS and POT3D HDF files are stored in unit of solar radii; this unit allows
astropy to convert them to cm, km, or AU automatically.