;1.5) Program: 	Recoupled long-range dipolar dephasing
;Description: 	Selects nonprotonated aromatic carbons far from protons (e.g. charcoal)
;Requires: 	!; DP/TOSS or CP/TOSS at ca. 7 kHz MAS
;Reference(s):	Mao & KSR, J. Magn. Reson. 162, 217-227 (2003). 

;sronepgaref.toss
;reference for sronepgare.toss
;alternative to the simpler sronep.tossgare
;from srtossgare.ksr and sronep.toss
;from srtossT2.ksr
;1H recoupled dephasing controlled by l4 (multiple of tr)
;TOtal Sideband Suppression (TOSS) using four 180 13C pulses
;durations between 180 pulses (d25-d29) are calculated from l31 (=SS in Hz)
;pl2 is power level for contact, pl12 is power level for decoupling on F2

"l20=aq/(p8*2)"

  "d22=(0.5s/l31)"		;gated recoupling..
  "d23=(0.5s/l31)-p2/2-1u"

  "d25=(1s/l31)*0.1226-(p2/2)-2u-3u"   ;TOSS...
  "d26=(1s/l31)*0.0773-p2"   ;too short for power level switch
  "d27=(1s/l31)*0.2236-p2-6u"
  "d28=(1s/l31)*1.0433-p2-6u"
  "d29=(1s/l31)*0.7744-(p2/2)-3u-d13"


1 ze
  10u pl1:f1		;preselect pl1 drive power level for F1
2 d1 do:f2		;recycle delay, turn the decoupler off
  10u pl12:f2		;preselect pl2 drive power for F2
  1u:f1 ph2

  p1:f1 ph2		;90 13C excitation pulse

   1u cpds2:f2
6    d22		; = Ntr 
     d22
  lo to 6 times l4
     
     d22 		
     d23

  p2:f1 ph7             ;180 13C for Hahn echo

     d23 		; = Ntr
     d22
 
7    d22
     d22
  lo to 7 times l4

   1u
  
  p1:f1 ph11		;store for gamma integral

   d6 do:f2		;z filter
   1u cw:f2 pl12:f2 pl11:f1 	;turn on F2 TPPM decoupling, set decoupling drive power level

  p1:f1 ph12		;read-out

   d25                   ;
   3u pl22:f2
  p2:f1 ph3             ;180 13C #1
   d26                   ;
  p2:f1 ph4             ;180 13C #2
   3u pl12:f2
   d27                   ;
   3u pl22:f2
  p2:f1 ph5             ;180 13C #3
   3u pl12:f2
  d28                   ;
   3u pl22:f2
  p2:f1 ph6             ;180 13C #4
   3u pl12:f2
  d29 cpd2:f2                  ;

   2u:f1 ph0             ;dead time delay, reset RF phase to detect
   go=2 ph31		;do NS

;regular cw decoupling
;  aq			;allow time for data sampling
;  1m do:f2		;turn decoupler off

;tppm decoupling
;3 p8:f2 ph20			;d2 < 180 degree, ph10= 25 degree
;  p8:f2 ph21			;ph11=335 degree
;  lo to 3 times l20
;  aq do:f2

 
  1m do:f2
  id6			;increment d6 by in6=tr/4: gamma-average for clean TOSS
  lo to 2 times l5      ;then l5 = 4

  wr #0

exit

ph0= 0			        ;reference phase
ph2= 1 3 2 0 3 1 0 2   ;13C 90 excitation

ph3= 1 1 2 2 3 3 0 0    ;180 13C #1
ph4= 3 3 0 0 1 1 2 2    ;180 13C #2
ph5= 1 1 2 2 3 3 0 0    ;180 13C #3
     3 3 0 0 1 1 2 2
ph6= 2 2 1 1 0 0 3 3    ;180 13C #4
     0 0 3 3 2 2 1 1
ph7= 0	0 1 1 2 2 3 3	;180 in first Hahn echo
     0	0 1 1 2 2 3 3
     0 0 1 1 2 2 3 3
     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
      1 1 2 2 3 3 0 0
      1 1 2 2 3 3 0 0
     2 2 3 3 0 0 1 1
     2 2 3 3 0 0 1 1
     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
      3 3 0 0 1 1 2 2   
      3 3 0 0 1 1 2 2   ;128 steps

ph8= 0 1 0 1		;1H 180 sequence
ph11= 3 3 0 0 1 1 2 2   ;90 store
      3 3 0 0 1 1 2 2
      1 1 2 2 3 3 0 0	;cancel transverse
      1 1 2 2 3 3 0 0
;combined with exorcycle of ph7
      1 1 2 2 3 3 0 0	;cancel transverse
      1 1 2 2 3 3 0 0
      3 3 0 0 1 1 2 2   ;90 store
      3 3 0 0 1 1 2 2
ph12= 1 1 2 2 3 3 0 0   ;90 read out
ph20= (360) 10			;TPPM first pulse
ph21= (360) 0			;TPPM second pulse
ph31= 0 2 1 3 2 0 3 1		;RLS
      0 2 1 3 2 0 3 1
      2 0 3 1 0 2 1 3
      2 0 3 1 0 2 1 3