Page MenuHomec4science
Files F65061134

stable.sh
No OneTemporary
Actions

File Metadata

Created
Fri, May 31, 11:47

stable.sh
View Options

#!/bin/bash -l
#SBATCH --job-name=espic2d
#SBATCH --time=12:00:00
#SBATCH -n 6
#SBATCH -c 6
#SBATCH --mem=144G
#srun amplxe-cl -c hotspots -r results_$nbprocs -- ../src/espic2d < stable.in
module purge
module restore espic2d-intel
export I_MPI_PLATFORM=auto
export espicsrc='/home/lebars/espic2d/src'
export espicwk='/home/lebars/espic2d/wk'
res_folder='/scratch/lebars/gt170_refurbished/Deep_well/'
mkdir -p $res_folder
cat > job.in << EOM
A Skeleton for MPI Time-Dependent program
=========================================
T.M. Tran SPC/EPFL
-
&BASIC
job_time=43000.0, extra_time=200.0,
nrun=100, !# of steps
nlres=f,
newres=f,
femorder=2,2,
ngauss=6,6,
nlppform=.TRUE.
partperiodic=f
nlxg=f, ! Display graphical interface
nlclassical=t, ! Solve classical equations of motion
nz=600, !# of intervals in z
nnr=65,80,25 !# of intervals between radii(1) and radii(2), between radii(2) and radii(3) and between radii(3) and radii(4)
it0d=10,
it2d=100,
ittext=100,
itparts=5000,
ittracer=10,
itgraph=100,
nbcelldiag=1,
itcelldiag=20,
resfile='${res_folder}tes30kv_w.h5'
rstfile='${res_folder}res30kv_w.h5'
nbspecies=3,
partfile='electrons_gauss.in','electron_tracers.in','ion_tracers.in',
radii=0.055, 0.076, 0.080,0.091 !radius of inner metallic wall, plasma boundaries and outer metallic wall
plasmadim=-0.009, 0.0155, 0.074, 0.078, !zmin zmax rmin rmax of initial plasma loading
distribtype=7 ! 1: uniform RZ gaussian in V, 2: stable eq 4.85 from Davidson
lz=-0.1,0.185, !Cylinder length
nplasma=2116800, !132300, ! 100000# of superparticles
nlmaxwellsource=t,
nblock=2000,
nlPhis=t, ! if false deactivate calculation of self electric fieldq
! nplasma=100
potinn=-30000, !potential at inner wall [V]
potout=0, !30000, !potential at outer wall [V]
B0=0.23, !Davidson chapter 4 formula 4.89 [T]
Rcurv=1.5,!1.001, !Davidson chapter 4 formula 4.89
width=0.48,
n0=-3e17, dt=2E-12 !Initial plasma density in plasmadim
! n0=-4E16, dt=8E-12
! n0=-5.e15, dt=8E-12,
! n0=-3.e19, dt=8e-13,
! n0=-5.e8
temp=22000, !Initial temperature of plasma [K] (thermal velocity)
H0=3.2e-14,
P0=8.66e-25,
magnetfile='phiBprofile_refurbasg_well'
bscaling=0,
/
&maxwellsourceparams
frequency=5E8,
temperature=22000,
rlimits=0.064,0.081,
zlimits=-0.010,0.034
time_start=-1.0,
radialtype=2
time_end=-1.0
/
&celldiagparams
specieid=1,
rindex=112,
zindex=117,
/
&geomparams
!x=0.028/2*cos(t)+0.012;y=0.002*sin(t)+0.081
r_a=0.06375
z_0=0.012
r_0=0.081
z_r=0.014
r_r=0.002
r_b=0.081
above2=-1
walltype=9 ! ellipse plus flat bottom and flat top
nlweb=t
/
&spldomain
nbsplines=2,
dist_extent=5e-3,
dr=-0.002,
Lz=0.014,
ra=0.06375,
rb=0.081
epsge=1e-5
epsce=1e-10
h5fname='refurb_geom.h5',
/
&neutcolparams
neutdens=2.5e21,
!Neon parameters
Eion=21.56,
scatter_fac=24.2,
nlcol=t,
io_cross_sec_file='${espicwk}/Ne_io_cross_sec.in',
ela_cross_sec_file='${espicwk}/Ne_ela_cross_sec.in',
/
EOM
export OMP_NUM_THREADS=$SLURM_CPUS_PER_TASK
#srun -n $nbprocs ../src/espic2d "$VAR"
srun $espicsrc/espic2d job.in

Event Timeline

c4science · Help