Ambient solar wind using the GONG-ADAPT-WSA maps can be computed with the master script xwsada.
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 trough the Earth:
[Image produced by plot_tim2e5.pro]
Execution of the master script xwsada with default parameters produces results in the following sub-directories:
code.ra1low1p5.16-mcp1um1cd-1 with executables (reg2grd, wsa2bc, enlil, out2evo)
case.ra1low1ip5.t45gR000-a7b1-g53x05.20181101 with grid parameters (grd.nc) and boundary conditions (bnd.nc)
run.ra1low1p5.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.dat Values correspond to the evolution plot shown above.
The master script xwsada 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:
Slow solar wind: runpar = g53, g15, g53h100000, g53t03 (from left to right)
[Images produced by plot_rad1.pro]