Run “ENLIL” with “WSA” data and “Cone” parameters from a list:
calib = Calibration run with different model-free parameters
demo = Demonstration runs for selected CME events
valid = Validation runs
ENLIL heliospheric computations are driven by boundary values calculated WSA and CME data for background solar wind and CME-like transients, respectively.
Numerical computations are performed in the computational region (parameter “reg”) with numerical resolution (parameter “res”) for 1 or more computational blocks (parameter “nblk”) using processors in parallel (parameter “procs”). There are two basic numerical models specified by prameter “numo” - numo=mcp1um1: robust model but with numerical diffusion; - numo=mcp1va2: less-robust model but with reduced diffusion. - wsacode=wsadt2bc: WSA daily-updated sequence of maps. Location of these maps is specified by parameters: - cordir, cordata, corfile, crnum, crlon. and boundary values are calculated with model free-parameters listed below and labeled by “amb”.
Transient disturbances are calculated from geometric and kinematic parameters of CMEs fitted by Cone-model analytic formula (using single-perspective observations) or determined by various iterative tools (when multi-perspective observations are available). These parameters are either explicitly listed (if ncmes>0) or given in cmefile (if ncmes=-1); no distrubances are specified if ncmes=0. Transient disturbances are calculated with model free-parameters listed below and labeled by “cld”.
The master scriot can be executed by a single command:
xwcone [-<parameter> <value>]
or, alternatively, by individual steps:
| command | action |
|---|---|
xwcone case [-<parameter> <value>] |
Set input parameters and execute the case codes |
xwcone plot_case [-<parameter> <value>] |
Plot results of the case computation |
xwcone run [-<parameter> <value>] |
Set input parameters and execute the run code |
xwcone data_run [-<parameter> <value>] |
Calculate data results from the run computation |
xwcone plot_run [-<parameter> <value>] |
Plot results of the run computation |
xwcone clean [-<parameter> <value>] |
Clean the project directory |
The following parameters can be used to modify their default values (see the master script for details):
| argument | default value | comment |
|---|---|---|
-project <char> |
Project name: | |
-projrun <int> |
0 | Project run (0=relax+phys, 1=relax, 2=phys): |
-numo <char> |
mcp1um1d | Numerical model setting: |
-nreg <int> |
1 | Computational grid mode: |
-reg <int> |
2 | Computational region (AU): |
-res <char> |
low | Numerical grid resolution: |
-obsdate <char> |
gongz | Use observations at or up to this date: |
-rundate <char> |
2010-04-03 | Start computation at this date: |
-tstop <float> |
168 | Stop computation at this time (h): |
-tstep <float> |
24 | Output with this step in time (h): |
-istep <int> |
1 | Plot with this step in output files: |
-cordata <char> |
Coronal data sub-directory: | |
-corfile <char> |
null | Coronal data file (null if cordata/corymdh used): |
-cormode <char> |
single | Coronal data mode: |
-corobs <char> |
gongz | Solar photospheric field observatory: |
-correa <char> |
null | Realization of the WSADA coronal map: |
-corymdh <char> |
Coronal data sub-sub-directory structure by date: | |
-crnum <char> |
Carrington Rotation number: | |
-crlon <char> |
Carrington Rotation leading longitude: | |
-amb <char> |
a7b1 | Ambient wind conditions setting: |
-bndrot <char> |
synodic | Rotation of the inner boundary: |
-maxcmes <int> |
30 | Maximum number of simulated CMEs (depends on system memory): |
-ncmes <int> |
60 | Number of simulated CMEs: |
-cld <char> |
d4t1x1 | Hydrodynamic cloud conditions setting: |
-timramp <float> |
0.5 | Transition from background to transient (h): |
-dcld <float> |
4. | Cloud density / fast stream density: |
-dcld0 <float> |
200. | Cloud background density (cm-3): |
-tcld <float> |
1. | Cloud temperature / fast stream temperature: |
-tcld0 <float> |
1. | Cloud background temperature (MK): |
-vcld <float> |
900. | Cloud velocity (km/s): |
-xcld <float> |
1. | Cloud trailing elongation over spherical shape: |
-ncld <int> |
2 | Cloud shape (1=piston, 2=sphesymbnd, 3=sphesymcen, 4=sphekincen): |
-dcav <float> |
0.25 | Cavity density / fast stream density: |
-tcav <float> |
4. | Cavity temperature / fast stream temperature: |
-xcav <float> |
0. | Cavity radius / cloud radius: |
-dtail <float> |
0. | Tail density / fast stream density: |
-ttail <float> |
0. | Tail temperature / fast stream temperature: |
-vtail <float> |
0. | Tail velocity (km/s): |
-xtail <float> |
5. | Tail elongation over spherical shape: |
-ntail <int> |
0 | Tail temporal decline (0=none, 1=back, 2=half, 3=full): |
-runpar <char> |
g53t03 | Run parameters: gamma and heating: |
You have to edit in the ENLIL distribution: - environment variables in env.inc script. You have to edit in this script: - project name and directory (where results will go); - data directory (repository of coronal maps).
Default parameters correspond to a low|resolution computation using the reference numerical model and ambient wind settings for the 2015 August period with CMEs in a heliospheric region between 0.1 and 2.1 AU for -336+168 h (-14+7 days).
It is assumed you are in the
Values of the radial velocity and radial magnetic field at the inner heliospheric boundary are written into the bnd.nc file. This is the NetDCF file and enables computation of the background solar wind by ENLIL code.