From 37f2aabb0048b1de9db32fbc91b427ed648ca5b2 Mon Sep 17 00:00:00 2001
From: ivar <amvdw95@gmail.com>
Date: Fri, 10 Aug 2012 10:52:44 +0000
Subject: [PATCH] removed conversion from mol/mol to ppm

---
 da/tools/MixingratioToPermil.py | 40 +++++++++++++++++++++++++++++++++
 1 file changed, 40 insertions(+)
 create mode 100755 da/tools/MixingratioToPermil.py

diff --git a/da/tools/MixingratioToPermil.py b/da/tools/MixingratioToPermil.py
new file mode 100755
index 00000000..91c3f15f
--- /dev/null
+++ b/da/tools/MixingratioToPermil.py
@@ -0,0 +1,40 @@
+#!/usr/bin/env python
+# MixingratioToPermil.py
+
+"""
+Author : ivar 
+
+Revision History:
+File created on 11 May 2012.
+
+"""
+
+
+def MixingratioToPermil(filename,simulated):
+    """ Converts 13C mixing ratios to permil values"""
+    import da.tools.io4 as io
+    import logging
+    import numpy as np
+
+    pdb         = 0.011112
+
+    trlength=len(simulated[0,:])
+    memlength=trlength/2
+    simulated=simulated #*1.e6 #convert to ppm
+ #   np.set_printoptions(threshold=np.nan)
+ #   msg='simulated shape',simulated.shape;logging.info(msg)
+ #   msg='simulated',simulated[0,0],simulated[0,40];logging.info(msg)
+ #   msg='simulated',simulated ;logging.info(msg)
+    simulated=np.float64(simulated)
+    simulated[:,memlength:trlength]=1000.*((simulated[:,memlength:trlength]/simulated[:,0:memlength]/pdb)-1.)
+ #   msg='simulated',simulated[0,0],simulated[0,40],memlength,trlength,pdb;logging.info(msg)
+ #   msg='simulated',simulated ;logging.info(msg)
+
+
+
+
+    return simulated
+
+
+if __name__ == "__main__":
+    pass
-- 
GitLab