The master script xwambdu computes the background solar wind. It uses the daily-updated WSA maps.
Boundary values of the solar wind radial velocity (top) and proton number density (bottom) for realization R000 at the boundary (at 0.1 AU) for the realization R0000 and the first map.
[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 xwambdu with default parameters produces results in the following sub-directories:
code.ra1low1p3.16-mcp1um1cd-1 with executables (reg2grd, wsa2bc, enlil, out2evo)
case.ra1low1ip3.t45gR000-a7b1-g53x05.20181101 with grid parameters (grd.nc) and boundary conditions (bnd.nc)
run.ra1low1p3.t45gR000-a7b1-g53x05.16-mp1va2d-1.t03d02.20181101 with values at time levels (tim.<rec>.nc) and evolution at observers (evo.<namobs>.nc and evo/evo.<namobs>.dat)
Numerical results evo/evo_earth.dat should be the same as the archived results doc/evo_earth_20111127T00-20111227T00.dat Values correspond to the evolution plot shown above.
The master script xwambdu 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:
Single-map: gongb-wsa22 (left), gongz-wsa22 (middle), gongz-wsa52 (right)
Multi-maps: gongb-wsa22 (left), gongz-wsa22 (middle), gongz-wsa52 (right)
[Images produced by plot_evo5o1.pro]