Read atoms from data file and generation of a specific amount of mass

Submitted by Matteo on Fri, 08/12/2016 - 14:46

Hello,

I am trying to read atoms stored in a dat file by using the command "read_data" but I get the following error:

ERROR on proc 1: Invalid atom type in Atoms section of data file (../atom_vec_sphere.cpp:1066)

---------------------------------------------
My input file is as follow:

modify_timing on
atom_style granular
atom_modify map array
boundary f f f
newton off
communicate single vel yes
units si
read_data data/Particles.dat
---------------------------------------------------
The dat file is:

2 atoms
1 atoms types

#Simulation box
-0.135 0.135 xlo xhi
-0.05 0.335 ylo yhi
-0.05 0.135 zlo zhi

#Body
#atom number, material type, diameter, x, y, z, density
Atoms

1 1 0.00105 4500 0.030 0.3 0.104
2 1 0.00105 4500 0.000 0.3 0.08

Velocities

1 0 0 0 0 0 0
2 0 0 0 0 0 0
--------------------------------------------

I do not understand what is wrong. Do you have any suggestion?

Also, is there a way to repeatedly introduce these atoms to reach a specific amount of mass? For example as the insert/stream command allows to do?
I would like to do this because I need to insert in my simulation a specific amount of particles with a size distribution which cannot be described by a Gaussian distribution. I would then need to specify in a data file the dimensions of my particles, group them into bins, and insert different quantities of them into the model.