examples/balance/in.balance.neigh.rcb

   1 # 3d Lennard-Jones melt
   2
   3 units           lj
   4 atom_style      atomic
   5 processors      * 1 1
   6
   7 lattice         fcc 0.8442
   8 region          box block 0 10 0 10 0 10
   9 create_box      3 box
  10 create_atoms    1 box
  11 mass            * 1.0
  12
  13 region          long block 3 6 0 10 0 10
  14 set             region long type 2
  15
  16 velocity        all create 1.0 87287
  17
  18 pair_style      lj/cut 2.5
  19 pair_coeff      * * 1.0 1.0 2.5
  20 pair_coeff      * 2 1.0 1.0 5.0
  21
  22 comm_style      tiled
  23
  24 neighbor        0.3 bin
  25 neigh_modify    every 2 delay 4 check yes
  26 fix             p all property/atom d_WEIGHT
  27 fix             0 all balance 50 1.0 rcb weight neigh 0.8 weight store WEIGHT
  28 compute         p all property/atom d_WEIGHT
  29 variable        maximb equal f_0[1]
  30 variable        iter   equal f_0[2]
  31 variable        prev   equal f_0[3]
  32 variable        final  equal f_0
  33
  34 #fix            3 all print 50 "${iter} ${prev} ${final} ${maximb}"
  35
  36 fix             1 all nve
  37
  38 #dump           id all atom 50 dump.melt
  39 #dump           id all custom 50 dump.lammpstrj id type x y z c_p
  40
  41 #dump           2 all image 25 image.*.jpg type type &
  42 #               axes yes 0.8 0.02 view 60 -30
  43 #dump_modify    2 pad 3
  44
  45 #dump           3 all movie 25 movie.mpg type type &
  46 #               axes yes 0.8 0.02 view 60 -30
  47 #dump_modify    3 pad 3
  48
  49 thermo          50
  50
  51 run             250
  52 run             250
  53