提交: 9a26d5e27da004230d63e908404f3333827546e2 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修訂: 9a26d5e27da004230d63e908404f3333827546e2 (tree)
• 時間: 2010-10-28 21:15:35
• 作者: lorenzo
• Commiter: lorenzo

Log Message

I updated the code and it now produces what should be the final version of
the two plots.

Change Summary

• modified: Python-codes/plot_polar_clip.py

差異

diff -r 086a0761c1ef -r 9a26d5e27da0 Python-codes/plot_polar_clip.py
--- a/Python-codes/plot_polar_clip.py	Wed Oct 27 16:25:41 2010 +0000
+++ b/Python-codes/plot_polar_clip.py	Thu Oct 28 12:15:35 2010 +0000
@@ -25,7 +25,17 @@
    return abs(cos(y))
 
 
-theta=0.8*pi/2.
+theta=0.5*pi/2.
+
+theta_1=0.
+
+theta_2=-theta
+
+
+theta_4=-pi/2
+
+theta_5=pi/2
+
 
 # make these smaller to increase the resolution
 dx, dy = 0.0025, 0.0025
@@ -42,6 +52,31 @@
 xy_end=(pi/1.5*cos(theta), pi/1.5*sin(theta))
 
 
+xy_surf_1=(pi/2.*cos(theta_1), pi/2.*sin(theta_1))
+
+xy_end_1=(pi/1.5*cos(theta_1), pi/1.5*sin(theta_1))
+
+
+xy_surf_2=(pi/2.*cos(theta_2), pi/2.*sin(theta_2))
+
+xy_end_2=(pi/1.5*cos(theta_2), pi/1.5*sin(theta_2))
+
+
+xy_end_3=(0.,0.)
+
+xy_surf_3=((pi/1.5-pi/2.),0.)
+
+
+xy_surf_4=(pi/2.*cos(theta_4), pi/2.*sin(theta_4))
+
+xy_end_4=(pi/1.5*cos(theta_4), pi/1.5*sin(theta_4))
+
+xy_surf_5=(pi/2.*cos(theta_5), pi/2.*sin(theta_5))
+
+xy_end_5=(pi/1.5*cos(theta_4), pi/1.5*sin(theta_5))
+
+
+
 print("x is, " )
 
 print (x)
@@ -59,8 +94,6 @@
 
 ax = fig.add_subplot(111)
 
-figtext(.63, .65,r'$-{\bfJ}_1^\perp$', fontdict=None,fontsize=20
-        )
 
 
 
@@ -69,13 +102,14 @@
 
 im.set_clip_path(Circle((0,0),pi/2., transform=ax.transData))
 
-annotate("", xy=(0., ini+1.), xytext=(0, -ini), arrowprops=dict(fc="black"))
+annotate("", xy=(0., ini+1.), xytext=(0, -ini), arrowprops=dict(fc="black",
+                                           arrowstyle="-"))
 
-for y in np.arange(-1.4, 1.5, 0.2):
-   annotate("", xy=(pi/2., y), xytext=(ini, y), arrowprops=dict(fc="g"))
+# for y in np.arange(-1.4, 1.5, 0.2):
+#    annotate("", xy=(pi/2., y), xytext=(ini, y), arrowprops=dict(fc="g"))
 
-ax.annotate('', xy=(-.3, 2.4),  xycoords='data',
-                xytext=(-.4, 2.2), # textcoords='offset points',
+ax.annotate('', xy=(-.3, 2.6),  xycoords='data',
+                xytext=(-.4, 2.4), # textcoords='offset points',
                 arrowprops=dict(arrowstyle="->",
                                 connectionstyle="angle3,angleA=0,angleB=-10"),
                 )
@@ -83,8 +117,50 @@
 annotate("", xy=xy_surf,
          xytext=xy_end, arrowprops=dict(fc="g"))
 
+annotate("", xy=xy_surf_1,
+         xytext=xy_end_1, arrowprops=dict(fc="g"))
 
-figtext(.5, 0.7,r'$\hat{\bf s}$', fontdict=None,fontsize=20
+annotate("", xy=xy_surf_2,
+         xytext=xy_end_2, arrowprops=dict(fc="g"))
+
+annotate("", xy=xy_surf_3,
+         xytext=xy_end_3, arrowprops=dict(fc="g"))
+
+
+annotate("", xy=xy_surf_4,
+         xytext=xy_end_4, arrowprops=dict(fc="g"))
+
+annotate("", xy=xy_surf_5,
+         xytext=xy_end_5, arrowprops=dict(fc="g"))
+
+
+
+
+# figtext(.5, 0.7,r'$\hat{\bf J}_1$', fontdict=None,fontsize=20
+#         )
+
+
+
+figtext(.73, .66,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+figtext(.73, 0.22,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+figtext(.79, 0.43,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+
+figtext(.53, 0.43,r'$\hat{\bf s}_{\perp}$', fontdict=None,fontsize=20
+        )
+
+
+
+figtext(.53, .1,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+
+figtext(.53, .7,r'${\bf J}_1$', fontdict=None,fontsize=20
         )
 
 
@@ -96,14 +172,17 @@
 
 # colorbar(ticks=[] )
 
-# ax = plt.axes() 
-cb = plt.colorbar()
-cb.ax.set_yticks((0.0, 1))
-cb.ax.set_yticklabels(( r"$\min(-{\bf J}_1^{\perp}\cdot \hat{\bf s})$",
-                        r"$\max(-{\bf J}_1^{\perp}\cdot \hat{\bf s})$"))
+
+# cb = plt.colorbar()
+# cb.ax.set_yticks((0.0, 1))
+# cb.ax.set_yticklabels(( r"$\min(-{\bf J}_1^{\perp}\cdot \hat{\bf s})$",
+#                         r"$\max(-{\bf J}_1^{\perp}\cdot \hat{\bf s})$"))
 
 
 savefig("first-plot.pdf", bbox_inches='tight')
+
+
+
  
 clf()    
 
@@ -111,12 +190,13 @@
 
 
 
-theta=0.4*pi/2.
 
 
-xy_surf=(pi/2.*cos(theta), pi/2.*sin(theta))
 
-xy_end=(pi/1.5*cos(theta), pi/1.5*sin(theta))
+
+xy_end_3=(0.,0.)
+
+xy_surf_3=(0.,(pi/1.5-pi/2.))
 
 
 Z = func5(X, Y)
@@ -127,47 +207,98 @@
 
 ax = fig.add_subplot(111)
 
-figtext(.65, .67,r'$-{\bf J}_1^{\parallel}$', fontdict=None,fontsize=20)
+
 
 im = imshow(Z,cmap=cm.jet, extent=(-pi, pi, -pi, pi))
 im.set_interpolation('bilinear')
 
 im.set_clip_path(Circle((0,0),pi/2., transform=ax.transData))
 
-for x in np.arange(-1.4, 1.5, 0.2):
-   if (abs(x)>1e-4):
-      annotate("", xy=(x,-pi/2.), xytext=(x, -ini), arrowprops=dict(fc="g"))
+# for x in np.arange(-1.4, 1.5, 0.2):
+#    if (abs(x)>1e-4):
+#       annotate("", xy=(x,-pi/2.), xytext=(x, -ini), arrowprops=dict(fc="g"))
 
-for x in np.arange(-1.4, 1.5, 0.2):
-   if (abs(x)>=1e-4):
-      annotate("", xy=(x,pi/2.), xytext=(x, ini), arrowprops=dict(fc="g"))
+# for x in np.arange(-1.4, 1.5, 0.2):
+#    if (abs(x)>=1e-4):
+#       annotate("", xy=(x,pi/2.), xytext=(x, ini), arrowprops=dict(fc="g"))
 
-annotate("", xy=(0., ini+1.), xytext=(0, -ini), arrowprops=dict(fc="black"))
+annotate("", xy=(0., ini+1.), xytext=(0, -ini), arrowprops=dict(fc="black",
+                                                                arrowstyle="-"))
 
-ax.annotate('', xy=(-.3, 2.4),  xycoords='data',
-                xytext=(-.4, 2.2), # textcoords='offset points',
+ax.annotate('', xy=(-.3, 2.6),  xycoords='data',
+                xytext=(-.4, 2.4), # textcoords='offset points',
                 arrowprops=dict(arrowstyle="->",
                                 connectionstyle="angle3,angleA=0,angleB=-10"),
                 )
 
+
 annotate("", xy=xy_surf,
          xytext=xy_end, arrowprops=dict(fc="g"))
 
+annotate("", xy=xy_surf_1,
+         xytext=xy_end_1, arrowprops=dict(fc="g"))
 
-figtext(.65, 0.61,r'$\hat{\bf s}$', fontdict=None,fontsize=20
+annotate("", xy=xy_surf_2,
+         xytext=xy_end_2, arrowprops=dict(fc="g"))
+
+annotate("", xy=xy_surf_3,
+         xytext=xy_end_3, arrowprops=dict(fc="g"))
+
+
+annotate("", xy=xy_surf_4,
+         xytext=xy_end_4, arrowprops=dict(fc="g"))
+
+
+annotate("", xy=xy_surf_5,
+         xytext=xy_end_5, arrowprops=dict(fc="g"))
+
+
+figtext(.73, .66,r'${\bf J}_1$', fontdict=None,fontsize=20
         )
 
+figtext(.73, 0.22,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+figtext(.77, 0.43,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+
+
+
+
+figtext(.53, .1,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+
+figtext(.53, .7,r'${\bf J}_1$', fontdict=None,fontsize=20
+        )
+
+
+
+
+
 
 ax.set_xlim(-ini, ini)
 ax.set_ylim(-ini, ini+1)
 ax.set_xticks([]) # this is needed for bbox_inches
 ax.set_yticks([])
 
-cb = plt.colorbar()
-cb.ax.set_yticks((0.0, 1))
-cb.ax.set_yticklabels(( "min","max"))
-cb.ax.set_yticklabels(( r"$\min(-{\bf J}_1^{\parallel}\cdot \hat{\bf s})$",
-                        r"$\max(-{\bf J}_1^{\parallel}\cdot \hat{\bf s})$"))
+# cb = plt.colorbar()
+# cb.ax.set_yticks((0.0, 1))
+# cb.ax.set_yticklabels(( "min","max"))
+# cb.ax.set_yticklabels(( r"$\min(-{\bf J}_1^{\parallel}\cdot \hat{\bf s})$",
+#                         r"$\max(-{\bf J}_1^{\parallel}\cdot \hat{\bf s})$"))
+
+
+
+
+
+
+
+figtext(.52, 0.43,r'$\hat{\bf s}_{\parallel}$', fontdict=None,fontsize=20
+        )
+
+
 
 
 savefig("second-plot.pdf", bbox_inches='tight')