rep.pro is used to display the data, with a lot of options. Start by typing :
which should display the density field over the whole disc, in a polar (r,phi) representation, at the fifth output (hence, 5 orbits in you ran the
template2.par file provided) :
rep.pro internally calls
open.pro, you can change density to whatever you like (pressure, vortensity, etc.). For a detailed list of the possibilities of this command,
you should consult the help provided by issuing:
The (r,phi) representation is not mandatory, one can prefer (phi,r). This is controlled by the repmode keyword :
Another (phi,r) view is available with repmode=0, with an r axis inverted.
This keyword is sticky, so you don’t need to always specify it.
As for a cartesian (x,y) view, just type :
which center the view onto the planet ([1,0]), with a extend of 0.5 *2 in the x and y direction.
Opposite to the repmode, the cartesian mode is not sticky.
As an example, let us have a look at the azimuthal speed.
The center & zoom behaviour is also available in polar mode. You can try :
You can also explicitely specify the range. For example, his command will produces the same result as the last one :
Internally, the values of
phirange are computed every time (even if you only specify center & zoom), and stored, to be sticky.
Now that we have a closer look, we can try to see finer things. To see the perturbed velocity (that is, relative to the output 0), use the keyword relative:
You should perfectly see the wake, but not much details in the horseshoe region. The keyword zoptimal will take care of this, by adjusting the zrange color scale to saturate the wake. While we’re working on the color scale, we can ask it to be symmetric, with the keyword zsymmetric :
The small velocities perturbation inside the horseshoe region are now visible.
Please note that you can still set the zrange color scale by hand, by adding zrange=[z_low,z_high].
But an other use for the zrange keyword is to recall the previous zrange, by setting zrange=1. This is very useful for comparing two outputs. Here is an example :
By looking at the stride on top of the window, you’ll see that the color scale was kept.
We just spoke about the horseshoe region. You can see it via the separatrice keyword, which overplots the streamline corresponding to the separatrix :
On top of this, you can also gain an overview of the flow by plotting streamlines :
This adds 10 streamlines, regularly spaced. The color of the separatrix was changed to distinguish it :
A note about large/complicated runs.
rep.pro routine internally calls
repcart.pro, which makes use of IDL’s
CONTOUR routine to display the data. For slow computers or large grids, this method can be slow. You can speed it up by adding /light to any command, which lowers the number of levels.