Hi Falgun, As I mentioned, you need to load the script on the command line:
pymol arrows.py or when already running Pymol you load the script with: run arrows.py The @ is for loading pymol scripts, which are different from python scripts. Please read more about the pymol syntax. Also, it doesn't hurt to get to know python. That would be a good investment in your future. Cheers, Tsjerk On Fri, Mar 13, 2009 at 8:14 PM, Falgun Shah <fhs...@olemiss.edu> wrote: > dear sir > > Thanks for the script. Sir, i am not much familiar with python. what i did > is i installed it from plugin/install plugin in python. also in command > input area i gave this command: @arrow.py for installation. > > it showed me following error > /home/Amber9.0/Desktop/desktop/pymol/modules/pymol/parser.py, line 285, in > parse > parsing.run_file(exp_path(args[nest][0]),pymol_names,pymol_names) > File "/home/Amber9.0/Desktop/pymol/modules/pymol/parsing.py", line 407, in > run_file > execfile(file,global_ns,local_ns) > File "arrow.py", line 115 > C[0]/sqrt(C[0]**2+C[1]**2) ) ) > ^ > SyntaxError: invalid syntax > > Show i tried this command: dist arrow > PyMOL>cmd.load_cgo([13.668, -9.582, 5.011, 7.347,-31.499, 16.242], "arrow") > > nothing appeared (no error message too!) > > dist arrow, resi 580 and name ca, id 27014, label=0 > dashed line appeared as before > > I think i am making some mistakes in placing the script or executing the > command. Can you please help me out > > i am sorry for asking too much. > > I will appreciate your help. > On Fri, Mar 13, 2009 at 12:44 PM, Tsjerk Wassenaar <tsje...@gmail.com> > wrote: >> >> # Hi Falgun, >> >> # Here's a mail I sent earlier to the lsit. I don't know why it >> doesn't show up when searching. >> # Anyway, this still is a python proof e-mail. You can save the whole >> text as arrow.py and add it to >> # the pymol command line or use 'run' >> >> # Cheers, >> >> print "Tsjerk :p" >> >> #---------- Forwarded message ---------- >> #From: Tsjerk Wassenaar <tsje...@gmail.com> >> #Date: Thu, Jan 24, 2008 at 10:03 AM >> #Subject: Re: [PyMOL] Arrows >> #To: João Costa <ramosco...@hotmail.com> >> #Cc: pymol-users@lists.sourceforge.net >> >> >> ## (: Python proof E-mail :) ## >> >> # Hi João, >> # >> # It indeed appears that the attachments are not preserved for the >> list (or that one at least) >> # You can save the text of this mail as arrow.py and run it instead. >> # >> # Best, >> # >> # Tsjerk >> >> from pymol import cmd >> from pymol.cgo import * >> from math import * >> >> # >> # Some functions to allow drawing arrows (vectors) in Pymol >> # In need of proper documentation... >> # >> # Please don't distribute (parts of) this file, without credits >> # >> # (c)2006 Tsjerk A. Wassenaar, PhD, University of Utrecht >> # >> # t s j e r k w .at. g m a i l .dot. c o m >> # http://nmr.chem.uu.nl/~tsjerk/ >> # >> >> def t( X ): >> if not X: return X >> Y = [] >> for i in range( len( X[0] ) ): >> Y.append( [] ) >> for j in X: >> Y[i].append( j[i] ) >> return Y >> >> def v_add( a, b ): return ( a[0]+b[0], a[1]+b[1], a[2]+b[2] ) >> def v_sub( a, b ): return ( a[0]-b[0], a[1]-b[1], a[2]-b[2] ) >> def vecprod( a, b ): return ( a[1]*b[2]-a[2]*b[1], >> a[2]*b[0]-a[0]*b[2], a[0]*b[1]-a[1]*b[0] ) >> def inprod( a, b=None ): >> if b: return a[0]*b[0] + a[1]*b[1] + a[2]*b[2] >> else: return a[0]*a[0] + a[1]*a[1] + a[2]*a[2] >> def svmult( s, a ): return ( s*a[0], s*a[1], s*a[2] ) >> def norm( a ): return svmult( 1/sqrt(inprod( a )), a ) >> >> def mvmult( R, x ): >> y = [] >> for i in R: y.append( inprod( i, x ) ) >> return tuple(y) >> >> def mv_add( X, a ): >> Y = [] >> for i in X: >> Y.append( v_add( i, a ) ) >> return Y >> >> def mmmult( R, X ): >> Y = [] >> for i in X: Y.append( mvmult( R, i ) ) >> return Y >> >> def smatrix( v ): return [[ v[0], 0, 0 ], [ 0, v[1], 0 ], [ 0, 0, v[2] ]] >> def rmatrix( v ): >> cosx, sinx = cos( v[0] ), sin( v[0] ) >> cosy, siny = cos( v[1] ), sin( v[1] ) >> cosz, sinz = cos( v[2] ), sin( v[2] ) >> >> return mmmult( mmmult( [[1,0,0],[0,cosx,-sinx],[0,sinx,cosx]], >> [[cosy,0,-siny],[0,1,0],[siny,0,cosy]] ), >> [[cosz,-sinz,0],[sinz,cosz,0],[0,0,1]] ) >> >> def block( i, dphi ): >> ddphi = 0.25*dphi >> phi0 = i*dphi >> phi1 = phi0+ddphi >> phi2 = phi1+ddphi >> phi3 = phi2+ddphi >> phi4 = phi3+ddphi >> sqrt2 = sqrt(2) >> return [ (-0.5*sqrt2,-0.5*sqrt2*cos(phi2),0.5*sqrt2*sin(phi2)), >> (1,0,0), >> (0,cos(phi0),-sin(phi0)), >> (0,cos(phi1),-sin(phi1)), >> (0,cos(phi2),-sin(phi2)), >> (0,cos(phi3),-sin(phi3)), >> (0,cos(phi4),-sin(phi4)) >> ] >> >> def cgo_triangle_fan( X ): >> Y = [] >> while ( X ): >> i = X.pop(0) >> Y.extend( [ NORMAL, i[0], i[1], i[2], ] ) >> for i in range( 6 ): >> i = X.pop(0) >> Y.extend( [ VERTEX, i[0], i[1], i[2], ] ) >> return Y >> >> def cgo_arrow1( S, E, r=0.2, hr=0.4, hl=1.0 ): >> P0 = S >> D = v_sub( E, S ) >> DL = inprod( D, D ) >> P1 = v_add( S, svmult( (DL-hl)/DL, D ) ) >> P2 = E >> >> # Define a vector orthogonal to P1-P0 >> V = v_sub( P1, P0 ) >> V = norm( V ) >> if V[2] != 0: >> A = ( 1, 1, -(V[0]+V[1])/V[2] ) >> elif V[1] != 0: >> A = ( 1, -V[0]/V[1], 0 ) >> else: >> A = ( 0, -V[0], 0 ) >> A = norm( A ) >> >> B = vecprod( V, A ) >> print (inprod(V), inprod(B), inprod(A)) >> R = t([ svmult( hl,V ), svmult( hr,A ), svmult( hr,B ) ]) >> >> # Define the transformation matrix (scale and rotation) >> #C = v_sub( P2, P1 ) >> #scale = ( hl, hr, hr ) >> #rotate = ( 0, acos( C[0]/sqrt(C[0]**2+C[2]**2) ), acos( >> C[0]/sqrt(C[0]**2+C[1]**2) ) ) >> #R = mmmult( smatrix( scale ), rmatrix( rotate ) ) >> >> obj = [ >> CYLINDER, S[0], S[1], S[2], P1[0], P1[1], P1[2], r, 1, 1, 1, 1, 1, 1, >> COLOR, 1, 0, 0, >> BEGIN, TRIANGLE_FAN ] >> >> N = 10 >> dphi = 2*pi/N >> crds = [] >> for i in range(N+1): >> crds.extend( block( i, dphi ) ) >> crds = mv_add( mmmult( R, crds ), P1 ) >> obj.extend( cgo_triangle_fan( crds ) ) >> obj.extend( [ END, ] ) >> >> return obj >> >> def cgo_arrow( S, E, r=0.2, hr=0.4, hl=1.0, name="arrow", state=1 ): >> obj = cgo_arrow1( S, E, r=r, hr=hr, hl=hl ) >> cmd.load_cgo( obj, name, state ) >> >> def cgo_arrows( X, r=0.2, hr=0.4, hl=1.0, name="arrow", state=1 ): >> obj = [] >> for i in X: >> obj.extend( cgo_arrow1( (i[0], i[1], i[2]), (i[3], i[4], >> i[5]), r=r, hr=hr, hl=hl ) ) >> >> cmd.load_cgo( obj, name, state ) >> >> >> #-- >> # Tsjerk A. Wassenaar, Ph.D. >> # Junior UD (post-doc) >> # Biomolecular NMR, Bijvoet Center >> # Utrecht University >> # Padualaan 8 >> # 3584 CH Utrecht >> # The Netherlands >> # P: +31-30-2539931 >> # F: +31-30-2537623 >> >> >> >> -- >> Tsjerk A. Wassenaar, Ph.D. >> Junior UD (post-doc) >> Biomolecular NMR, Bijvoet Center >> Utrecht University >> Padualaan 8 >> 3584 CH Utrecht >> The Netherlands >> P: +31-30-2539931 >> F: +31-30-2537623 > > > > -- > Falgun H shah > PhD candidate (3rd year) > Department of Medicinal Chemistry > 2028, Natural Product Center > University of Mississippi > > Ph No: 6629151286(O) > 662 801 5667(M) > email: fhs...@olemiss.edu > -- Tsjerk A. Wassenaar, Ph.D. Junior UD (post-doc) Biomolecular NMR, Bijvoet Center Utrecht University Padualaan 8 3584 CH Utrecht The Netherlands P: +31-30-2539931 F: +31-30-2537623
