1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
#!/usr/bin/env python
# Script: logplot.py
# Purpose: use GnuPlot to plot two columns from a LAMMPS log file
# Syntax: logplot.py log.lammps X Y
# log.lammps = LAMMPS log file
# X,Y = plot Y versus X where X,Y are thermo keywords
# once plot appears, you are in Python interpreter, type C-D to exit
# Author: Steve Plimpton (Sandia), sjplimp at gmail.com
import sys,os,code
path = os.environ["LAMMPS_PYTHON_TOOLS"]
sys.path.insert(1,path)
from log import log
from gnu import gnu
if len(sys.argv) != 4:
raise Exception("Syntax: logplot.py log.lammps X Y")
logfile = sys.argv[1]
xlabel = sys.argv[2]
ylabel = sys.argv[3]
lg = log(logfile)
x,y = lg.get(xlabel,ylabel)
g = gnu()
g.plot(x,y)
code.interact()
|
