r4422 - /branches/N_state_model/generic_fns/align_tensor.py -- January 07, 2008 - 14:20

Author: bugman
Date: Mon Jan  7 14:20:42 2008
New Revision: 4422

URL: http://svn.gna.org/viewcvs/relax?rev=4422&view=rev
Log:
Updated the alignment tensor code to use the new gyromagnetic names as per 
the 1.3 line.

The old names used were {gC, gH, gN, gO, gP} whereas the new names are {g13C, 
g1H, g15N, g17O,
g31P}.


Modified:
    branches/N_state_model/generic_fns/align_tensor.py

Modified: branches/N_state_model/generic_fns/align_tensor.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/N_state_model/generic_fns/align_tensor.py?rev=4422&r1=4421&r2=4422&view=diff
==============================================================================
--- branches/N_state_model/generic_fns/align_tensor.py (original)
+++ branches/N_state_model/generic_fns/align_tensor.py Mon Jan  7 14:20:42 
2008
@@ -32,7 +32,7 @@
 from data import Data as relax_data_store
 from data.align_tensor import AlignTensorData
 from data.data_classes import SpecificData
-from physical_constants import gC, gH, gN, gO, gP, h_bar, mu0
+from physical_constants import g13C, g1H, g15N, g17O, g31P, h_bar, mu0
 import pipes
 from relax_errors import RelaxError, RelaxNoPipeError, RelaxNoTensorError, 
RelaxStrError, RelaxTensorError, RelaxUnknownParamCombError, 
RelaxUnknownParamError
 
@@ -544,27 +544,27 @@
 
     # Gyromagnetic ratio of the first nucleus.
     if nuc1 == 'C':
-        gI = gC
+        gI = g13C
     elif nuc1 == 'H':
-        gI = gH
+        gI = g1H
     elif nuc1 == 'N':
-        gI = gN
+        gI = g15N
     elif nuc1 == 'O':
-        gI = gO
+        gI = g17O
     elif nuc1 == 'P':
-        gI = gP
+        gI = g31P
 
     # Gyromagnetic ratio of the second nucleus.
     if nuc2 == 'C':
-        gS = gC
+        gS = g13C
     elif nuc2 == 'H':
-        gS = gH
+        gS = g1H
     elif nuc2 == 'N':
-        gS = gN
+        gS = g15N
     elif nuc2 == 'O':
-        gS = gO
+        gS = g17O
     elif nuc2 == 'P':
-        gS = gP
+        gS = g31P
 
     # Kappa.
     return -3.0/(8.0*pi**2) * gI * gS * mu0 * h_bar