mailr21538 - /branches/relax_disp/specific_analyses/relax_disp/cpmgfit.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on November 20, 2013 - 10:38:
Author: bugman
Date: Wed Nov 20 10:38:38 2013
New Revision: 21538

URL: http://svn.gna.org/viewcvs/relax?rev=21538&view=rev
Log:
Bug fix for the relax_disp.cpmgfit_input user function.

The nu_CPMG values need to be doubled and then divided by 1e3 to obtain the 
1/tau_CPMG values in ms.


Modified:
    branches/relax_disp/specific_analyses/relax_disp/cpmgfit.py

Modified: branches/relax_disp/specific_analyses/relax_disp/cpmgfit.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/specific_analyses/relax_disp/cpmgfit.py?rev=21538&r1=21537&r2=21538&view=diff
==============================================================================
--- branches/relax_disp/specific_analyses/relax_disp/cpmgfit.py (original)
+++ branches/relax_disp/specific_analyses/relax_disp/cpmgfit.py Wed Nov 20 
10:38:38 2013
@@ -270,8 +270,11 @@
         # Tesla units.
         B0 = frq * 2.0 * pi / g1H
 
+        # The X value of 1/tcp (or 1/tau_CPMG) in ms.  This assumes Art's 
usage of the definition that nu_CPMG = 1 / (2 * tau_CPMG).
+        x = 2.0 * point / 1000.0
+
         # Write out the data and error.
-        file.write("%-20f %-20f %-20f %-20f\n" % (point/1000.0, 
spin.r2eff[param_key], spin.r2eff_err[param_key], B0))
+        file.write("%-20f %-20f %-20f %-20f\n" % (x, spin.r2eff[param_key], 
spin.r2eff_err[param_key], B0))
 
     # Close the file and return its name.
     file.close()




Related Messages


Powered by MHonArc, Updated Wed Nov 20 11:40:02 2013