The master script xamb3cone computes propagation of a heliospheric disturbance initiated by the Cone model in the ambient solar wind (AMB).
Boundary values of the solar wind velocity, density, temperature, and magnetic field strength (from left to right) at time when a hydrodynamic ejecta passes the boundary (at 0.1 AU):
[Image produced by plot_bnd.pro]
Distribution of the radial solar wind velocity on three slices passing trough the observing point “P03” (at 1 AU) and evolution at this point:
[Image produced by plot_tim3.pro]
Evolution of the solar wind radial velocity, number density, temperature, and magnetic field strength (from top to bottom) at the observing point “P03”:
[Image produced by plot_evo5o1.pro]
Execution of the master script xamb3cone with default parameters produces results in the following sub-directories:
code.ra2low2.16-mcp1um1cd-1 with executables (ana3d3bc, enlil, out2evo)
case.ra2low2.b30d120t03v30-d4t1v900r30p2 with grid parameters (grd.nc) and boundary conditions (bnd.nc)
run.ra2low2.b30d120t03v30-d4t1v900r30p2.16-mcp1um1cd-1.g53t03b02 with values at time levels (tim.<rec>.nc) and evolution at observers (evo.<namobs>.nc and evo/evo.<namobs>.dat)
Numerical results evo/evo1_p03.dat should be the same as doc/evo1_p03.dat.
The master script xamb3cone can be copied and edited for computing with different parameters of the ambient wind, cone, and numerical model.
Alternatively, these parameters can be passed to the master script as arguments when it is executed. These modifications can be collected into a parametric script which executes the master script, under different conditions, as individual commands. Each run produces results into the “code”, “case”, and “run” directories named according to the used parameters.
Examples of parametric scripts:
Boundary values: nsho=1 (left) and ntail=3 (right)
Heliospheric solution: nsho=1 (left) and ntail=3 (right)
[Images produced by plot_tim1e1.pro]
Examples of parametric scripts:
Boundary values: ntail=0 (left), ntail=2 (middle), and ntail=3 (right)
[Images produced by plot_bndtim1.pro]
Heliospheric solution: ntail=0 (left), ntail=2 (middle), and ntail=3 (right)
[Images produced by plot_tim1e1.pro]