Figure 4a

Data for figure 4a is available from the PDS.  See the Data
Availability Statement in Brecht and Ledvina [2021]. 


Figure 4b & c

The data file contains the sum of the O+ and O2+ densities in
the ecliptic plane.  This is a ASCII data file that contains
313026 data points.  The data was written with the following
statements:

do j = 1,609
   do i = 1,514
      write(30,*)den(i,j)
   enddo
enddo

The x-axis (i index) covers 5.2776 Venus radii, while the
y-axis covers 6 Venus radii.  Thus dx = 0.0102877 Rv and
dy = 0.00986842 Rv.  The planet is located at xp = 3 Rv and
yp = 2.52943 Rv. So to convert to a planet centered coordinate
system xpc = x-xp and ypc = y-yp.