forked from eutactic/rxnlvl
-
Notifications
You must be signed in to change notification settings - Fork 4
/
test.py
executable file
·43 lines (36 loc) · 1.68 KB
/
test.py
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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
#! /usr/bin/python
# Boilerplate
import sys
sys.path.insert(1,"/home/user/bin/rxnlvl/") # Change this to the full path of rxnlvl
from rxnlvl import *
# Plot
p = plot([25.0,10.0],vbuf=10.0,hbuf=5.0,bgcolour=None, qualified='sortof')
p + level(energy( 00, 'kjmol'), 1, '1', 0x0)
p + level(energy( -85.5, 'kjmol'), 2, 'EC1', 0x0)
p + level(energy( 244, 'kjmol'), 3, 'TS1a', 0xFF4444)
p + level(energy( 51, 'kjmol'), 3, 'TS1b', 0x0)
p + level(energy( -102, 'kjmol'), 4, 'DC1', 0x0)
p + level(energy( -82, 'kjmol'), 5, '2', 0x0)
p + level(energy( -111, 'kjmol'), 6, 'EC2', 0x0)
p + level(energy( -83, 'kjmol'), 7, 'TS2b', 0x0)
p + level(energy( 18, 'kjmol'), 7, 'TS2a', 0x44FF44)
p + level(energy( -103, 'kjmol'), 8, 'DC2', 0x0)
p + level(energy( 85, 'kjmol'), 9, '3', 0x0)
p + level(energy( 185, 'kjmol'), 10, '4', 0x0)
p + level(energy( 25, 'kjmol'), 11, '5', 0x0)
p + edge( '1', 'EC1', 0x0, 0.4, 'normal')
p + edge( 'EC1', 'TS1a', 0x0, 0.2, 'normal')
p + edge( 'EC1', 'TS1b', 0x0, 0.4, 'normal')
p + edge( 'TS1a', 'DC1', 0x0, 0.2, 'normal')
p + edge( 'TS1b', 'DC1', 0x0, 0.4, 'normal')
p + edge( 'DC1', '2', 0x0, 0.4, 'normal')
p + edge( '2', 'EC2', 0x0, 0.4, 'normal')
p + edge( 'EC2', 'TS2b', 0x0, 0.4, 'normal')
p + edge( 'EC2', 'TS2a', 0x0, 0.2, 'normal')
p + edge( 'TS2a', 'DC2', 0x0, 0.2, 'normal')
p + edge( 'TS2b', 'DC2', 0x0, 0.4, 'normal')
p + edge( 'DC2', '3', 0x0, 0.4, 'normal')
p + edge( '3', '4', 0x0, 0.4, 'normal')
p + edge( '4', '5', 0x0, 0.4, 'normal')
p + baseline(energy( 0.0, 'kjmol'),colour=0x0,mode='dashed',opacity=0.1)
p.write()