;6.3) Program: 	2D 19F exchange NMR with fast MAS (>25 kHz)
;Description: 	Determine proximity of various 19F sites
;Requires: 	!!;  
;Reference(s):	Chen & KSR, Macromolecules 2004


;from srF2dexwrT2fil.RFDR
;from srF2dexwrT2fil.fast:  with l14/l4 T2 filters in both dimensions
;from srF2dexwr.fast
;from srF2dexwr.toss
;rotor-synchronized evolution (SW=wr), spin diffusion
;from srFcodex2d.toss
;from srFcodex2d.simp
;from srFcodex.simp
;from srcodex.ksr

define delay cd15mp1
define delay cd15mp2

   "d13=l4*0.001s/l30 - p1/2 - p2/2"  ;exactly tr from center to center of pulses
   "d23=l14*0.001s/l30 - p1/2 - p2/2"  ;exactly tr from center to center of pulses
   "d26=0.001s/l30-p2-2*p1"  ;
   "d27=0.0005s/l30-p2/2-p1"
  
define loopcounter td1half
  "td1half=td1/2"			;for States cosine/sine
 
1 ze
  10u pl1:f1		;preselect pl1 drive power level for F1
2 d1 			;recycle delay
  rpp10
  rpp11

  p1:f1 ph1		;90 excitation pulse
   d23
  p2:f1 ph2
   d23
   d4			;isotropic-shift evolution, in4 = tr
  p1:f1 ph8     
			;mixing time
3   d27                 ;tr/2
  p1:f1 ph10
  p2:f1 ph11^		;RFDR:  One 180 pulse per tr
  p1:f1 ph10^
    d26 		;tr
  p1:f1 ph10
  p2:f1 ph11^		;RFDR:  One 180 pulse per tr
  p1:f1 ph10^
    d27			;tr/2
  lo to 3 times l6	

  p1:f1 ph9
   
   d13
  p2:f1 ph3             ;180 
   d3                   ;
 
  2u:f1 ph0		;dead time delay, reset RF phase to detect
  go=2 ph31		;start ADC with ph31 signal routing

 
;1D
;  wr #0

;2D				;States' detection
  200m wr #0 if #0  zd	;write file to disk, increment, zero data
  10m ip8			;cosine to sine
  lo to 2 times 2		;cosine/sine detection
  10m ip8 id4			;increment d4 by in4 = tr
  10m ip8
  lo to 1 times td1half

exit

ph0= 0				;reference phase for detection
ph1= 1 3 2 0 3 1 0 2  	;90 19F excitation
ph2= 2 2 3 3 0 0 1 1   ;180 before t1
ph3= 0 0 1 1 2 2 3 3   ;180 before detect (best phase for short d13)
     0 0 1 1 2 2 3 3
     1 1 2 2 3 3 0 0   ;exorcycle
     1 1 2 2 3 3 0 0   ;exorcycle
     2 2 3 3 0 0 1 1
     2 2 3 3 0 0 1 1
     3 3 0 0 1 1 2 2 
     3 3 0 0 1 1 2 2 

ph8= 3 3 0 0 1 1 2 2    ;+z cosine store for tm
     1 1 2 2 3 3 0 0    ;-z cosine cancel transverse

ph9=  1 1 2 2 3 3 0 0   ;read-out after tm for cosine

ph10= 1 2 1 2 2 1 2 1   ;90 of composite 180 pulse in RFDR
ph11= 0 1 0 1 1 0 1 0   ;xy-16 for RFDR
      2 3 2 3 3 2 3 2

ph31= 0 2 1 3 2 0 3 1   ;RCPH
      2 0 3 1 0 2 1 3   ;transverse cancellation
      2 0 3 1 0 2 1 3   ;with exorcycle
      0 2 1 3 2 0 3 1   ;with exorcycle