An Example
The output of a 1D hydro simulation is given by a variety of
scalars as a function of radius, the examples here are
from the file wrbub.1067.
This file can be read into IDL using the READCOL routine
and the various scalars can be plotted in desired units, e.g.,
as done by:
To get a sense of the 3D nature of these 1D scalars, their values
are used to weight the intensity of a spherically symmetric model
created as a 3D cube of cells. These 3D values can be viewed as
a "solid"-ish model in 3D ("thick"), cut in half to better show the
radial structure.
They can also be viewed as a 2D projection on the sky ("thin") -
this is the way they'd be seen
if viewed by a telescope (and assuming they are optically thin.)
Such images are presented below for the sample model file
wrbub.1067 and were produced using the IDL procedure:
The several scalar values from the data file are plotted and
imaged below; note that "pressure" looks to be
proportional to "density * temperature"
in this case.
The "X-ray flux" value is, as a very crude approximation,
taken to be proportional to the density squared and
is zero in regions where T is below 1.e6 K.
| Scalar w/1D Plot
| 3D view, cut in half
| 2D projected intensity (blue to red)
|
Velocity(r)
| 
| 
|
Density(r)
| 
| 
|
Density(r)
version 2, 4/12/06
| 
| 
|
Temp(r)
| 
| 
|
Pressure(r)
| 
| 
|
X-ray_flux(r)
~ density^2 for T > 1.e+6
| 
| 
|
X-ray_flux(r)
version 2, 4/12/06
~ SQRT( density^2 / sqrt(Temp) ) for T > 1.e+6
| 
| 
|
dd@space.mit.edu