Readme Figure 5.

There are 2 data files: SW-Edep-Case1.dat, the stationary Venus case,
and SW-Edep-Case2.dat, the orbiting Venus case.  Each data set is on
a spherical grid Phi, Theta, Radius of dimensions (73,36,71).  The
radial distance covers 100.0-450.0 km in altitude.  The Theta dimension
goes from near the north pole at 2.5 deg to near the south pole at 177.5
deg.  The spacing between theta points is 5 degrees.  The phi points are
in the azimuthal direction and cover the range of -2.5-362.5 deg. or 365
deg. So the values at phi(1) should be set to the values at phi(72) and
the values at phi(73) should be set to the values at phi(2).

To set up the grid read in the data with:

nr = 71
nth = 36
nphi = 73

do k = 1,nr
  do j = 1,nth
    do i = 1,nph
      read(10,*)swedep(i,j,k)
    enddo
  enddo
enddo

Then set the phi boundary cells:

do k = 1,nr
  do j = 1,nth
    swedep(1,j,k) = swedep(nph-1,j,k)
    swedep(nph,j,k) = swedep(2,j,k)
  enddo
enddo

The data was plotted using matlab.  The contour plot was made using
the slice routine though the altitude of 135 km.  The data for the
curves in the other panels where taken by interpolating between the
phi and theta points along the r dimension as needed.