mailr19976 - /branches/relax_disp/specific_analyses/relax_disp/__init__.py


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

Header


Content

Posted by edward on June 08, 2013 - 23:33:
Author: bugman
Date: Sat Jun  8 23:33:16 2013
New Revision: 19976

URL: http://svn.gna.org/viewcvs/relax?rev=19976&view=rev
Log:
Added support for the M61 skew model to the relax_disp.select_model user 
function back end.

This is the Meiboom 1961 on-resonance 2-site model for skewed populations (pA 
pB).

This commit follows step 6 of the relaxation dispersion model addition 
tutorial
(http://thread.gmane.org/gmane.science.nmr.relax.devel/3907).


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

Modified: branches/relax_disp/specific_analyses/relax_disp/__init__.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/specific_analyses/relax_disp/__init__.py?rev=19976&r1=19975&r2=19976&view=diff
==============================================================================
--- branches/relax_disp/specific_analyses/relax_disp/__init__.py (original)
+++ branches/relax_disp/specific_analyses/relax_disp/__init__.py Sat Jun  8 
23:33:16 2013
@@ -858,7 +858,7 @@
     def _select_model(self, model=MODEL_R2EFF):
         """Set up the model for the relaxation dispersion analysis.
 
-        @keyword model: The relaxation dispersion analysis type.  This can 
be one of 'R2eff', 'No Rex', 'LM63', 'CR72', 'M61'.
+        @keyword model: The relaxation dispersion analysis type.  This can 
be one of 'R2eff', 'No Rex', 'LM63', 'CR72', 'M61', 'M61 skew'.
         @type model:    str
         """
 
@@ -917,6 +917,14 @@
             for frq in loop_frq():
                 params.append('r2')
             params += ['phi_ex', 'kex']
+
+        # M61 skew model.
+        elif model == MODEL_M61B:
+            print("The Meiboom (1961) on-resonance 2-site model with skewed 
populations (pA >> pB) for R1rho-type experiments.")
+            params = []
+            for frq in loop_frq():
+                params.append('r2')
+            params += ['pA', 'dw', 'kex']
 
         # Invalid model.
         else:




Related Messages


Powered by MHonArc, Updated Sun Jun 09 00:00:02 2013