From 440dfe2d64fda9904a055aa72ac941a56d6bd52f Mon Sep 17 00:00:00 2001
From: gac-x04sa <gac-x04sa@x04sa-cons-2.psi.ch>
Date: Wed, 3 Jul 2019 13:46:48 +0200
Subject: [PATCH] Script execution

---
 script/users/bill/nno-fourcv.py | 75 +++++++++++++++++++++++++++++++++
 1 file changed, 75 insertions(+)
 create mode 100644 script/users/bill/nno-fourcv.py

diff --git a/script/users/bill/nno-fourcv.py b/script/users/bill/nno-fourcv.py
new file mode 100644
index 0000000..89ca13c
--- /dev/null
+++ b/script/users/bill/nno-fourcv.py
@@ -0,0 +1,75 @@
+###################################################################################################\
+#Geometry
+###################################################################################################
+set_geometry("fourcv", True)
+
+###################################################################################################\
+#Constraints
+###################################################################################################
+#hkl.con('a_eq_b')
+
+#equivalent to 'setmode 0' + 'freeze 5'
+hkl.con( 'alpha', 5) #alpha constant
+
+
+#equivalent to 'setmode 1' + 'freeze 5'
+#hkl.con( 'beta', 5) #alpha constant
+
+#equivalent to 'setmode 2''
+#hkl.con( 'a_eq_b', 0) #alpha constant
+
+#hkl.con( 'eta', 0) 
+
+
+setup_axis(alpha, alpha.getMinValue(), alpha.getMaxValue())
+setup_axis(delta, delta.getMinValue(), 90) #delta.getMaxValue())
+setup_axis(gamma, gamma.getMinValue(), gamma.getMaxValue())
+setup_axis(omegaV, omegaV.getMinValue(), omegaV.getMaxValue())
+    
+
+###################################################################################################\
+#Orientation
+###################################################################################################
+
+"""
+Enter real space lattice parameters:
+ Lattice a (1.305)? 7.723
+ Lattice b (1.305)? 7.707
+ Lattice c (1.305)? 7.723
+ Lattice alpha (90)? 90
+ Lattice beta (90)? 89.265
+ Lattice gamma (90)? 90
+"""
+
+
+newub('test1')
+
+setlat('cubic', 7.723, 7.707, 7.723, 90.0, 89.265, 90.0)
+#en = 9.5 
+#en = 12.4 #!!!!!
+
+
+
+#Only available in newer Diffcalc
+#ub.setmiscut(0)
+
+
+#Calculated with newer Diffcalc
+#   UB matrix:  0.81363   0.00000  -0.01044
+#               0.00000   0.81526   0.00000
+#               0.00000   0.00000   0.81357
+#ub.setub([[ 0.81363,   0.00000,  0.00000], [0.00000,  0.81526,   0.00000], [-0.01044,   0.00000,   0.81357]])
+
+
+addref([0, 1, 2], [5.0, 9.379, 19.7895, 102.6162], energy.take())
+addref([2, 0, 2], [5.0, 19.0754, 20.1865, 11.9693], energy.take())
+addref([2, 2, 2], [5.0, 27.1234, 21.2368, 60.0354], energy.take())
+showub()
+
+
+checkub()
+print getub()
+
+
+print  hklci((44.3400,    0.0000,  123.7322 , 0.0))[0]
+print hklca((0, 1, 2))[0]