mailr4992 - /branches/N_state_model/specific_fns/n_state_model.py


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

Header


Content

Posted by edward on February 15, 2008 - 09:19:
Author: bugman
Date: Fri Feb 15 09:19:39 2008
New Revision: 4992

URL: http://svn.gna.org/viewcvs/relax?rev=4992&view=rev
Log:
Renamed the cone angle variables.


Modified:
    branches/N_state_model/specific_fns/n_state_model.py

Modified: branches/N_state_model/specific_fns/n_state_model.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/N_state_model/specific_fns/n_state_model.py?rev=4992&r1=4991&r2=4992&view=diff
==============================================================================
--- branches/N_state_model/specific_fns/n_state_model.py (original)
+++ branches/N_state_model/specific_fns/n_state_model.py Fri Feb 15 09:19:39 
2008
@@ -147,11 +147,11 @@
         cdp.rot_CoM = norm(cdp.CoM) * cdp.red_CoM / len_red
 
         # The cone angle for diffusion on an axially symmetric cone.
-        cdp.diff_on_cone_theta = acos(len_red)
+        cdp.theta_diff_on_cone = acos(len_red)
 
         # The cone angle and order parameter for diffusion in an axially 
symmetric cone.
-        cdp.diff_in_cone_theta = acos(2.*len_red - 1.)
-        cdp.S_diff_in_cone = cos(cdp.diff_in_cone_theta) * (1 + 
cos(cdp.diff_in_cone_theta)) / 2.0
+        cdp.theta_diff_in_cone = acos(2.*len_red - 1.)
+        cdp.S_diff_in_cone = cos(cdp.theta_diff_in_cone) * (1 + 
cos(cdp.theta_diff_in_cone)) / 2.0
 
         # Print out.
         print "\nThe pivot point of the domain motions is:\n" + 
`cdp.pivot_point` + "\n"
@@ -160,8 +160,8 @@
         print "The reduced CoM vector is:\n" + `cdp.red_CoM` + "\n"
         print "The full length rotated CoM vector is:\n" + `cdp.rot_CoM` + 
"\n"
         print "The length reduction is:\n" + `len_red` + "\n"
-        print "The cone angle for diffusion on an axially symmetric cone 
is:\n%.5f rad (%.5f deg)\n" % (cdp.diff_on_cone_theta, cdp.diff_on_cone_theta 
/ (2*pi) *360.)
-        print "The cone angle for diffusion in an axially symmetric cone 
is:\n%.5f rad (%.5f deg)\n" % (cdp.diff_in_cone_theta, cdp.diff_in_cone_theta 
/ (2*pi) *360.)
+        print "The cone angle for diffusion on an axially symmetric cone 
is:\n%.5f rad (%.5f deg)\n" % (cdp.theta_diff_on_cone, cdp.theta_diff_on_cone 
/ (2*pi) *360.)
+        print "The cone angle for diffusion in an axially symmetric cone 
is:\n%.5f rad (%.5f deg)\n" % (cdp.theta_diff_in_cone, cdp.theta_diff_in_cone 
/ (2*pi) *360.)
         print "The order parameter (S_cone) for diffusion in an axially 
symmetric cone is:\nS_cone = %.5f\nS^2 = %.5f\n" % (cdp.S_diff_in_cone, 
cdp.S_diff_in_cone**2)
 
 




Related Messages


Powered by MHonArc, Updated Fri Feb 15 09:40:34 2008