From ae2a4df69228c1e6ab06c1d7bebfe079b3f8d9fc Mon Sep 17 00:00:00 2001
From: Chris Hodges <code@c-hodges.co.uk>
Date: Mon, 26 Jul 2021 12:31:25 +0100
Subject: [PATCH] Add files via upload

Example/test code for time-resolved files
---
 OPUS_TR_example.py | 57 ++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 57 insertions(+)
 create mode 100644 OPUS_TR_example.py

diff --git a/OPUS_TR_example.py b/OPUS_TR_example.py
new file mode 100644
index 0000000..200b898
--- /dev/null
+++ b/OPUS_TR_example.py
@@ -0,0 +1,57 @@
+#!/usr/bin/python3
+
+#This requires python 3, and brukeropusreader, which can be installed with pip3 install brukeropusreader (or sometimes just pip install...)
+# but to be sure it's up to date, get it from https://github.com/qedsoftware/brukeropusreader
+
+
+from brukeropusreader import read_file
+from brukeropusreader import parse_sm
+
+import matplotlib.pyplot as plt
+import numpy as np
+import glob
+import sys
+import os 
+
+plot=False
+
+
+
+def ProcessFile(filename):
+	opus_data = read_file(filename)
+	SC_X=opus_data.get_range("ScSm")
+	IG_X=opus_data.get_range("IgSm")
+	ScSm=parse_sm(opus_data)
+	IgSm=parse_sm(opus_data,"IgSm")
+	format_Ig=["%d"]
+	format_Ig.extend(["%.8e"]*np.shape(IgSm)[1])
+	format_Sc="%.8e"
+	print (format_Ig)
+	print(format_Sc)
+	print (ScSm)
+
+	np.savetxt(filename+"IG.txt",np.c_[IG_X,IgSm], fmt=format_Ig)
+	np.savetxt(filename+"SC.txt",np.c_[SC_X,ScSm], fmt=format_Sc)
+	
+	if plot:
+		plt.plot(SC_X,ScSm[:,0])
+		plt.show()
+		plt.plot(IG_X,IgSm[:,0])
+		plt.show()
+
+for arg in sys.argv:
+	if arg.lower().find('plot')>=0:
+		plot=True
+	elif len(arg.lower())>=1 :
+		InputFileName=arg
+print (InputFileName)		
+if InputFileName.find("*")>=0 or InputFileName.find("?")>=0:#called with wildcards
+	InputFileName=os.path.join(os.getcwd(),InputFileName)
+	for InputFile in glob.glob(InputFileName):
+		print ("Globbed input: loading file "+InputFile)
+		
+		ProcessFile(InputFile)
+else:
+	InputFile=InputFileName
+	print ("Loading file "+InputFile)	
+	ProcessFile(InputFile)		
\ No newline at end of file