Plotting with o2graph
=====================

One can plot the EOS data with the O₂sclpy (see
https://neutronstars.utk.edu/code/o2sclpy) command ``o2graph``. To
begin, first download an EOS table::

  o2graph -download fid_3_5_22.o2 \
  https://isospin.roam.utk.edu/public_data/eos_tables/du21/fid_3_5_22.o2

Then, to summarize the objects in the file::

  o2graph -filelist fid_3_5_22.o2

Which gives output similar to::

  O2scl object "A" of type tensor_grid.
  O2scl object "A_max" of type tensor_grid.
  O2scl object "A_min" of type tensor_grid.
  ...
  Dataset "n_T" of type size_t with value 160.
  Dataset "n_Ye" of type size_t with value 70.
  Dataset "n_nB" of type size_t with value 301.
  Dataset "n_oth" of type size_t with value 11.
  ...

To create a two-dimensional density plot of a tensor_grid object,
for example::

  o2graph -download fid_3_5_22.o2 \
  https://isospin.roam.utk.edu/public_data/eos_tables/du21/fid_3_5_22.o2 \
  -read fid_3_5_22.o2 A -set logx 1 -set logy 1 -set colbar 1 \
  -to-table3d 0 2 slice 0.4 -den-plot slice pcm=True \
  -xtitle "$ n_B~(\mathrm{fm}^{-3}) $" -ytitle "$ T~(\mathrm{MeV}) $" \
  -save A_Ye04.png

plots the average nuclear mass number as a function of density and
temperature at a fixed electron fraction of 0.4:

.. image:: static/A_Ye04.png
   :width: 60%     
   :alt: alt text

The following::

   o2graph -read fid_3_5_22.o2 T_grid -set logy 1 -plot1 \
   marker=+ -save T_grid.png
   
plots the temperature grid:

.. image:: static/T_grid.png
   :width: 60%     
   :alt: alt text