Suddenly, after multiple succesful including of “in.ion” lammps can not open it. What is the reason?
script 1:C_HHG_film0.in
units lj
atom_style charge
dimension 2
boundary m m p
lattice sq2 0.00874 origin 0.0 0.0 0.0
atom_modify map array
variable omega equal 0.057
variable Field equal 0.2
variable I1 equal 9.225/27.21
variable I2 equal 20.521/27.21
variable I3 equal 37.37/27.21
variable I4 equal 54.18/27.21
variable I5 equal 70.99/27.21
variable I6 equal 87.81/27.21
variable I7 equal 104.6/27.21
variable I8 equal 123.2/27.21
region 1 block -1 1 -1 1 -1 1
region 3 block -10 10 -1 1 -1 1
region ghostatoms block 0.0 0.1 0.0 0.0 0.0 0.0 units box
create_box 5 3
create_atoms 4 region 1
mass 1 10
mass 2 10
mass 3 1600
mass 4 10
mass 5 10
group ionsA type 4
group electron type 2
group ions union ionsA
#variable numel0 equal count(ionsA)*8
#create_atoms 1 random ${numel0} 1242313 ghostatoms
mass 1 1
#group electron0 type 1
#set group electron0 charge 0
set group ions charge 0.0
variable CoulCutoff equal 8
pair_style coul/debye 0.01 10
pair_coeff * *
timestep 0.005
variable time equal (step)*0.005
variable ksi atom ${omega}*v_time
variable EFx atom ${Field}*cos(v_ksi)
fix 2 all efield v_EFx 0.0 0.0
fix 1 all nve
thermo 10
dump 4 all custom 10 dump_HHGA.dat id type x y z mass q
fix for2d all enforce2d
variable j equal 0
label rerun0
variable j equal $j+1
run 20
include ionizationAuCD.in
label out1
if “${j} < 100000” then “jump C_HHG_film0.in rerun0”
script 2: ionizationAuCD.in
variable allN equal count(all)
label ionizationCycle
variable i loop ${allN}
include in.ion
next i
jump ionizationAuCD.in ionizationCycle
jump C_HHG_film0.in out1
script 3: in.ion
#HERE is in.ion!!
Sorry for large 1 script, i tried to delete all which is not necessary.