The master script xwambfr computes the background solar wind. It uses the full-rotation WSA maps.
Boundary values of the solar wind radial velocity (top) and proton number density (bottom) for the Carrington rotation (CR) 2095 (begins on 2010-03-26, ends on 2010-04-22) at the boundary (at 0.1 AU).
[Image produced by plot_bnd.pro]
Distribution of the radial solar wind velocity on three slices passing through the Earth:
[Image produced by plot_tim3.pro]
Evolution of the solar wind radial velocity, number density, and magnetic field strength (from top to bottom) at Earth:
[Image produced by plot_evo5o1.pro]
Execution of the master script xwambfr with default parameters produces results in the following sub-directories:
code.ra1low1p3.16-mp1va2cd-1 with executables (reg2grd, wsa2bc, enlil, out2evo)
case.case.ra1low1p3.fro-a6b1-g53x05.20100326 with grid parameters (grd.nc) and boundary conditions (bnd.nc)
run.ra1low1p3.fro-a6b1-g53x05.16-mp1va2d-1.t03d02.20100326 with values at time levels (tim.<rec>.nc) and evolution at observers (evo.<namobs>_<time-range>.nc and evo/evo.<namobs>_<time-range>.dat)
Numerical results evo/evo_earth_20111127T00-20111225T00.dat should be the same as the archived results doc/evo_earth_20100326T00-20100423T00.dat Values correspond to the evolution plot shown above.
The master script xwambfr can be copied and edited for different coronal maps and for different parameters of the ambient solar wind 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:
Start computation: rundate = 2010-04-09 (left) and rundate = 2010-04-22 (right)
[Images produced by plot_evo5o1.pro]