Follow the steps below to get EMSES up and running for the first time.
mkdir /LARGE0/gr20001/$USER
ln -s /LARGE0/gr20001/$USER ~/large0
.bashrc
, then log out and log back ingrep -qxF 'module load intel-python' ~/.bashrc || echo 'module load intel-python' >> ~/.bashrc
grep -qxF 'export PATH="$PATH:$HOME/.local/bin"' ~/.bashrc || echo 'export PATH="$PATH:$HOME/.local/bin"' >> ~/.bashrc
~/large0
and re-enter your passwordmkdir -p ~/large0/Github
cd ~/large0/Github
git clone https://github.com/CS12-Laboratory/MPIEMSES3D.git
cd MPIEMSES3D
make
cd ~/large0/Github
git clone https://github.com/CS12-Laboratory/EMSES-tutorials.git
cd EMSES-tutorials
pip install -r requirements.txt
dshield*
directorycp ~/large0/Github/MPIEMSES3D/bin/mpiemses3d dshield0/
cp ~/large0/Github/MPIEMSES3D/bin/mpiemses3d dshield1/
cp ~/large0/Github/MPIEMSES3D/bin/mpiemses3d dshield2/
dshield0
cd ~/large0/Github/EMSES-tutorials/dshield0
mysbatch job.sh
mysbatch: Custom command (camptools: https://github.com/Nkzono99/camptools)
Reads the “nodes(:)” entry in plasma.inp, sets the number of processes in job.sh, and then runs sbatch job.sh.
qs # Check job status (custom camptools command)
squeue # Check job status (Slurm)
qgroup # Check resource availability
latestjob # Custom command (camptools: https://github.com/Nkzono99/camptools)
# Shows the latest job’s log file (= tail -n 5 stdout.*****.log)
NOTE: To cancel a submitted job, run
scancel <job-id>
NOTE: To view the stdout & stderr of a running job:
stdout.****.log
: standard outputstderr.****.log
: standard error
squeue
stdout.*****.log
and stderr.****.log
for errors..mypython/plot.py
)These should be saved as data/*.png
and data/gif/*.gif
. For example:
phisp_2d_xy.png
dshield0/plot_example.ipynb
See the following references for how to visualize:
In dshield0/plasma.inp
, nstep
is set to 100
.
This only captures the very early interaction between plasma and the object.
TODO: Increase
nstep
to simulate over a longer time, then run again.
For dshield1
and dshield2
, open their respective plasma.inp
files and run each one.
Compare what changes between the cases.
dshield0
and discuss the results.An isolated negative charge in various plasma environments:
NOTE: If you set
wp = 0.0d0
in ds0, the visualization library emout will not work properly. Therefore, when visualizing ds0, setwp
to something like1.0d0
inplasma.inp
before running the visualization.
NOTE: If you run into issues following the simulation steps, first explain your situation to an M1 student and get advice. If you still can’t resolve the problem and need to ask your advisor or instructors, include all of the following information:
- What you asked the senior student
- What answer you received
- What you actually did based on that answer
- What still remains unresolved
After running the simulations, at minimum check:
nd1p
)nd2p
)phisp
)Discuss the physical meaning of these results with your B4 group members and be prepared to share your thoughts in the next tutorial.