Hi everyone,

I'm calling lammps as a shared lib from OpenFOAM. I would like to use

the "include" command as described here"

http://lammps.sandia.gov/doc/include.html " to include some variables

and fixes.

When I run the input script with lammps as standalone it includes the

files but when I call the script with lammps as a shared lib it doesn't

work. I've attached the input script and include files.

Any thoughts or tips? Thanks in advance.

i cannot reproduce it with the current version of LAMMPS (10 March 2017).

if i comment out the 'read_restart' command (for lack of having your

restart file), i can see "include" working via the library interface. e.g.

when using this little python script:

[[email protected] Downloads]$ cat test.py

from lammps import *

lmp=lammps()

lmp.command('echo both')

lmp.file('in.lammps')

which produces the following screen output indicating that the "include"

directive is properly processed.

axel.

[[email protected] Downloads]$ python test.py

LAMMPS (10 Mar 2017)

using 1 OpenMP thread(s) per MPI task

# input script lammpsFoam simulations

package omp 0

using multi-threaded neighbor list subroutines

atom_style sphere

atom_modify map array

boundary f f f

newton off

processors * * *

units cgs

comm_modify vel yes

# Read restart and import variables

variable path getenv "PWD"

print "path is \{path\}"
path is /home/akohlmey/Downloads
\# read\_restart restart\.grow
include {path}/variables

include /home/akohlmey/Downloads/variables

# Variables used in multiple scripts

# domain variables

variable H0 equal 11.5 # [cm]

variable EF equal 3.0 # [-]

variable hb equal v_H0*v_EF # [cm]

variable height equal 1.1*v_H0*v_EF # [cm], add 10% extra height to domain

variable width equal 1.0 # [cm]

variable depth equal 1.0 # [cm]

variable eps equal 1.-$((1.-0.4)/v_EF) # average vol fraction, needed for

particle insertion eqn

variable eps equal 1.-0.19999999999999998335

variable rhop equal 2.8 # [g/cm3]

# particle insertion variables

variable dslice equal 0.2 # thicknes of a single slice

variable slices equal floor(\{hb\}/{dslice}) # nr of insertion slices

variable slices equal floor(34.5/\{dslice\}\)
variable slices equal floor\(34\.5/0\.2\)
variable vslice equal (v_width)*\(v\_width\)\*v\_dslice \# volume of a single
slice
variable vslice equal 1\*(v_width)*v_dslice

variable vslice equal 1*1*v_dslice

variable dred equal 0.5 # reduction factor for insertion size

variable CG equal 2. # particle coarse graining factor

# Growth variables

variable dred equal 0.1 # reduction factor for insertion size

variable growts equal 25000

variable growevery equal 1

variable relaxts1 equal 10000

variable relaxts2 equal 100000

# granular variables (based on eqns from Landry 2003 )

variable dpav equal 0.025

variable mpav equal $(v_rhop*PI/6.*v_dpav*v_dpav*v_dpav)

variable mpav equal 2.2907446432425579132e-05

variable g equal 981.

variable kn equal \(2\.e5\*v\_mpav\*v\_g/v\_dpav\)
variable kn equal 179777\.63960167608457
variable kt equal (2./7.*v_kn)

variable kt equal 51365.039886193138955

variable gamman equal 50*\(sqrt\(v\_g/v\_dpav\)\)
variable gamman equal 50\*198\.09088823063012796
variable gammat equal {gamman}

variable gammat equal 9904.54441153151

variable mu equal 0.15

print "included variables"

included variables

include \{path\}/fixes
include /home/akohlmey/Downloads/fixes
\# pair\_style none
pair\_style gran/hertz/history/omp {kn} \{kt\} {gamman} \{gammat\} {mu} 1

pair_style gran/hertz/history/omp 179777.639601676 \{kt\} {gamman}

\{gammat\} {mu} 1

pair_style gran/hertz/history/omp 179777.639601676 51365.0398861931

\{gamman\} {gammat} \{mu\} 1
pair\_style gran/hertz/history/omp 179777\.639601676 51365\.0398861931
9904\.54441153151 {gammat} \{mu\} 1
pair\_style gran/hertz/history/omp 179777\.639601676 51365\.0398861931
9904\.54441153151 9904\.54441153151 {mu} 1

pair_style gran/hertz/history/omp 179777.639601676 51365.0398861931

9904.54441153151 9904.54441153151 0.15 1

pair_coeff * *

ERROR: Pair_coeff command before simulation box is defined

(../input.cpp:1762)