Home > Legacy archive > Visualization > IDL > IDL tools : rep.pro
The routine 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) :
As 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 :
You should see something like this :
Please note that the cartesian view is not as optimized and polished as the polar view. Nevertheless, you can at least have a better view of the planet, using :
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 :
which center the view on the planet, with a radial extend of 0.1 *2, and an azimuthal extend of 1 *2. :
The default value of center is [0,1], so we commonly omit to specify this keyword.
You can also explicitely specify the range. For example, his command will produces the same result as the last one :
Internally, the values of rrange and 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.
The rep.pro routine internally calls reppol.pro or 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.
RSS 2.0