mailr22451 - /trunk/specific_analyses/relax_disp/disp_data.py


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

Header


Content

Posted by tlinnet on March 12, 2014 - 09:29:
Author: tlinnet
Date: Wed Mar 12 09:29:18 2014
New Revision: 22451

URL: http://svn.gna.org/viewcvs/relax?rev=22451&view=rev
Log:
Added verbosity flag to calc_rotating_frame_params() to allow switching of 
print information.

Regarding sr #3124, (https://gna.org/support/index.php?3124) - Grace graphs 
production for R1rho analysis with R2_eff as function of Omega_eff.

Modified:
    trunk/specific_analyses/relax_disp/disp_data.py

Modified: trunk/specific_analyses/relax_disp/disp_data.py
URL: 
http://svn.gna.org/viewcvs/relax/trunk/specific_analyses/relax_disp/disp_data.py?rev=22451&r1=22450&r2=22451&view=diff
==============================================================================
--- trunk/specific_analyses/relax_disp/disp_data.py     (original)
+++ trunk/specific_analyses/relax_disp/disp_data.py     Wed Mar 12 09:29:18 
2014
@@ -156,7 +156,7 @@
     return intensity
 
 
-def calc_rotating_frame_params(spin=None, spin_id=None):
+def calc_rotating_frame_params(spin=None, spin_id=None, verbosity=0):
     """Calculates and rotating frame parameters, calculated from:
     - The spectrometer frequency.
     - The spin-lock or hard pulse offset.
@@ -173,6 +173,8 @@
     @type spin:         SpinContainer instance
     @keyword spin_id:   The spin ID string.
     @type spin_id:      None or str
+    @keyword verbosity: A flag specifying to print calculations.
+    @type verbosity:    int
     @return:            List with dict() of theta, Domega, w_eff and list of 
dict() keys.
     @rtype:             List of dict() 
     """
@@ -186,8 +188,9 @@
     chemical_shifts, offsets, tilt_angles, Delta_omega, w_eff = 
return_offset_data(spins=[spin], spin_ids=[spin_id], field_count=field_count, 
fields=spin_lock_nu1)
         
     # Loop over the index of spins, then exp_type, frq, offset
-    print("Printing the following")    
-    print("exp_type spin_id frq offset{ppm} offsets[ei][si][mi][oi]{rad/s} 
ei mi oi si di cur_spin.chemical_shift{ppm} 
chemical_shifts[ei][si][mi]{rad/s} spin_lock_nu1{Hz} 
tilt_angles[ei][si][mi][oi]{rad}")
+    if verbosity:
+        print("Printing the following")    
+        print("exp_type spin_id frq offset{ppm} 
offsets[ei][si][mi][oi]{rad/s} ei mi oi si di cur_spin.chemical_shift{ppm} 
chemical_shifts[ei][si][mi]{rad/s} spin_lock_nu1{Hz} 
tilt_angles[ei][si][mi][oi]{rad}")
 
     si = 0
     theta_spin_dic = dict()
@@ -199,7 +202,8 @@
         # Loop over the dispersion points.
         spin_lock_fields = spin_lock_nu1[ei][mi][oi]
         for di in range(len(spin_lock_fields)):
-            print("%-8s %-10s %11.1f %8.4f %12.5f %i  %i  %i  %i  %i %7.3f 
%12.5f %12.5f %12.5f"%(exp_type, spin_id, frq, offset, 
offsets[ei][si][mi][oi], ei, mi, oi, si, di, spin.chemical_shift, 
chemical_shifts[ei][si][mi], spin_lock_fields[di], 
tilt_angles[ei][si][mi][oi][di]))
+            if verbosity:
+                print("%-8s %-10s %11.1f %8.4f %12.5f %i  %i  %i  %i  %i 
%7.3f %12.5f %12.5f %12.5f"%(exp_type, spin_id, frq, offset, 
offsets[ei][si][mi][oi], ei, mi, oi, si, di, spin.chemical_shift, 
chemical_shifts[ei][si][mi], spin_lock_fields[di], 
tilt_angles[ei][si][mi][oi][di]))
             dic_key = return_param_key_from_data(exp_type=exp_type, frq=frq, 
offset=offset, point=spin_lock_fields[di])
             dic_key_list.append(dic_key) 
             theta_spin_dic["%s"%(dic_key)] = tilt_angles[ei][si][mi][oi][di]




Related Messages


Powered by MHonArc, Updated Wed Mar 12 09:40:03 2014