units        real
dimension    3
boundary    p p f

atom_style          atomic
neighbor          5.0 bin
neigh_modify    delay 5



#######box################
##region                box block 0 100.0 0 100.0 0 4500.0 units box
##create_box          2 box 
read_restart 6.restart
##########################

#####################region#######################
region                vapor block  5    95   5    95 5 2998 units box
region                liquid block 5    95   5    95 3002 4495 units box
region                vapornvt block EDGE EDGE EDGE EDGE EDGE 3000 units box
region                liquidnvt block EDGE EDGE EDGE EDGE 3000 EDGE units box
##################################################

##################create atoms####################
######masses##############
mass        1 39.948
mass        2 39.948
##########################

##########create atoms####
##create_atoms 1 random 104940 6000000 vapor
##create_atoms 2 random 314142 4000000 liquid
##########################

##########LJ potentials###
pair_style      lj/cut 12.0
pair_coeff      * * 0.2379 3.405 12.0
##########################

###################################################



#####################group########################
group    vapor1    type 1
group    liquid1    type 2

group vapor dynamic all region vapornvt every 1
group liquid dynamic all region liquidnvt every 1
##################################################

#####################Wall#########################
fix zwalls all wall/reflect zlo EDGE zhi EDGE
fix wallhi vapor wall/lj126 zhi 3000 1.0 1.0 1 units box
fix wallhi2 liquid wall/lj126 zlo 3000 1.0 1.0 1 units box
##################################################


#####################GCMC#########################
fix vapor1 vapor1 gcmc 100 1000 1000 1 123456543 164.65 0 1.0 region vapor pressure 106.81  tfac_insert 1
fix liquid1 liquid1 gcmc 100 1000 1000 2 123456543 86.0 0 1.0 region liquid pressure 111.66 tfac_insert 1
##################################################

###################Thermostat#####################
#######vapornvt region##########
fix md1 vapor nve
compute mdtemp1 vapor temp/region vapornvt
compute_modify mdtemp1 dynamic/dof yes
fix mdnvt1 vapor langevin 164.65 164.65 100.0 699483
fix_modify mdnvt1 dynamic/dof yes
fix_modify mdnvt1 temp mdtemp1

#######liquidnvt region##########
fix md2 liquid nve
compute mdtemp3 liquid temp/region liquidnvt
compute_modify mdtemp3 dynamic/dof yes
fix mdnvt3 liquid langevin 86 86 100.0 699483
fix_modify mdnvt3 dynamic/dof yes
fix_modify mdnvt3 temp mdtemp3

###################################################

##################Compute##########################
compute cc2 all chunk/atom bin/1d z lower 20.0
fix temppro2 all ave/chunk 10 100 1000 cc2 temp norm sample file temp20.profile

compute den2 all chunk/atom bin/1d z 0.0 20.0
fix denspro all ave/chunk 10 100 1000 den2 density/mass norm sample file density20.profile


###################################################

timestep       5
thermo        10
thermo_style      custom step atoms cpu temp etotal pe ke press c_mdtemp1 c_mdtemp3 



restart 10000 5.restart 6.restart




run        1000000