# 1 "/opt/topspin3.0/exp/stan/nmr/lists/pp/user/pjmatpass.ids"
; pjmatpass.yz
;
; newer version than pjmat2b
;
; MATPASS (f1,f2)=(aniso,iso+aniso) experiment
;  -projection-method magic-angle turning experiment
;   w/ two evolution segments
;  -no shearing necessary to obtain PASS spectrum
;
; ref: ZH Gan, RR Ernst, J Magn Reson A 1996, 123, 140.

;NS : 20 * n
;p1 : f1 90 pulse width
;d1 : relaxation delay
;d3 : =DE
;pl1 : f1 90 pulse power
;in8 : =D7 set in_f=3*in8
;in10 : =D7
;in11 : =D7*2
;l1 : # t1 points, set TD1=L1 during 2D acq
;l2 : # total MAT rotor periods, !=3*n, n=integer>1
;cnst0 : # Run a test expt. to set: cnst0 = the first order phasing in degrees when cnst0=0 (~1.31u/dw*180=588.84)
;cnst1 : # current t1 point, integer>0 2*td
;cnst31 : spinning frequency [Hz]

;$COMMENT=projection method magic-angle turning w/ delayed acquisition
;$CLASS=Solids
;$DIM=1D/2D
;$TYPE=spin-1/2
;$OWNER=ivan
# 33 "/opt/topspin3.0/exp/stan/nmr/lists/pp/user/pjmatpass.ids"
;"d3=de";d3 : =DE
;"l1=td1";l1 : # t1 points, set TD1=L1 during 2D acq
"cnst1=2*td1";cnst1 : # current t1 point, integer>0 2*td

;"cnst0=(d0/dw)*180"
"d0= dw*((cnst0+320)/180)" ; Not actually used, for display only.
"d9=0"

"d2 =p1"
"d7 = 1s/(cnst31*3*td1)"            ;d7 : =t1 increment, Tr/3/TD1
"d8 =(1s/cnst31)*(1-cnst1/(3*td1))-p1"
"d10=(1s/cnst31)*(l2+1/3+cnst1/(3*td1))-p1"
"d11=(1s/cnst31)*(2-cnst1*2/(3*td1))-(p1*2)/3.14159-de -d9 + dw*((cnst0+320)/180)"
;"d8 = 1s*cnst1/(cnst31*3*l1)-p1"
;"d10=(1s/(cnst31*3))*(l2-cnst1/l1)-p1"
;"d11= 2s*cnst1/(cnst31*3*l1)-(p1/2)-dw*2"
# 50 "/opt/topspin3.0/exp/stan/nmr/lists/pp/user/pjmatpass.ids"
;"in8=d7";in8 : =D7 set in_f=in8*3
;"in10=d7";in10 : =D7
;"in11=2*d7";in11 : =D7*2
"inf1=d7*3"


d7
# 1 "mc_line 57 file /opt/topspin3.0/exp/stan/nmr/lists/pp/user/pjmatpass.ids dc-measurement inserted automatically"

    dccorr
# 57 "/opt/topspin3.0/exp/stan/nmr/lists/pp/user/pjmatpass.ids"
1 ze
2 d1
  5u pl1:f1
  (3u ph1):f1
  (p1 ph1):f1
    d8
  (p1 ph2):f1
    d10
  (p1 ph3):f1
    d8
  (p1 ph4):f1
    d10
    d8
    d2
  (p1 ph5):f1
    d11
  go=2 ph31
  50u wr #0 if #0 zd
  5u
  "d8=d8+d7"
  "d10=d10-d7"
  "d11=d11+(d7*2)"
;  5u id8
;  5u dd10
;  5u id11
  lo to 2 times td1
exit
# 85 "/opt/topspin3.0/exp/stan/nmr/lists/pp/user/pjmatpass.ids"
;p=0 +1 0 +1 0 -1
ph1= 0
ph2=(20) 0 12  4 16  8
ph3=(20) 0 10
ph4=(20) 0  2  4  6  8 10 12 14 16 18
ph5=(20) 0  1  2  3  4  5  6  7  8  9
        10 11 12 13 14 15 16 17 18 19
ph31=0 1 2 3

;p=0 -1 0 -1 0 -1
;ph1= 0
;ph2=(20) 0  8 16 4 12
;ph3=(20) 0 14  8  2 16 10  4 18 12  6
;ph4=(20) 0  2  4  6  8 10 12 14 16 18
;ph5=(20) 0  1  2  3  4  5  6  7  8  9
;        10 11 12 13 14 15 16 17 18 19
;ph31=0 1 2 3

