REM ***********************************************************************
REM *                                                                     *
REM *                           COSYPH - GMB 2000                         *
REM *                                                                     *
REM * Homonuclear shift correlated spectrosopy (phase cycled). This is    *
REM * a slightly modified version of the cosy pulse program provided with *
REM * the instrument. The pulse sequence is:                              *
REM *                                                                     *
REM *                     rd - p1 - d0 - p2 - acquire                     *
REM *                                                                     *
REM * TD=SI, number of t1 experiments, x=TD/4, and zero fill in the       *
REM * second dimension. p0 is 90 pw in usec, p1 is either 90 or 45, rd=1  *
REM * to 5 times T1, d is the dwell time and, the time increment in the   *
REM * t1 dimension. The phase cycle is adecuate to give phase-sensitive   *
REM * data instead of magnitude mode.                                     *
REM *                                                                     *
REM ***********************************************************************
REM read in cosyph conditions (8 scans...)
re "c:\eft\h1\cosyph.ini"
showwin
cls
print "COSY acquisition program (phase cycled)."
print 
print "The relaxation delay needs to be set in the range 1 to 5"
print "times the expected T1 for the sample. 2.0 sec is a good"
print "general purpose value."
print
input "Enter a relaxation delay.",d0
cls
print "COSY acquisition program (phase cycled)."
print 
print "Enter a file name for the data. The default file"
print "my_cosyph can be used by responding with an <Enter>."
print "Any file with the same name will be overwritten"
print "without warning."
input "Enter output file name, e.g. data\my_cosyph<Enter>. ",f$
if f$="" then f$="data\my_cosyph"
cls
print "COSY acquisition program (phase cycled)."
print 
print "The preparation pulse is 90 degrees."
define p0 h90 pulses.txt
if p0>0 then print "The 90 degree pulsewidth chosen is ",p0
if p0=0 then input "Enter 90 degree pulsewidth: ",p0
REM Here ask if this is a COSY45 or COSY90...
print
print "Select the pulsewidth of the mixing pulse. Enter"
print "1 for a 90 degree pulse, or 2 for a 45 degree pulse."
input "Pulswidth of the mixing pulse (1) or (2): ",m
REM Set the default p1 to 90 if what the user entered is garbage...
p1=p0
if m=1 then p1=p0
if m=2 then p1=p0/2
SLEEP 1
REM first pulse phase
a$="00112233"
REM second pulse phase
b$="02132031"
REM receiver phase
c$="00112233"
REM will use 3*dw/16 delay before turning receiver on
e=3*d/16
REM will use dw/16 delay before turning digitizer on
f=d/16
REM set up the program for GS
PROGRAM "clear"
PROGRAM 0 d0 S
PROGRAM + 2.0 uS    
PROGRAM + p0  uS TX a$
PROGRAM + e uS
PROGRAM + f uS RX
PROGRAM + 1.0 uS ADC RX
PROGRAM + a S RX
PROGRAM "coadd" c$
PROGRAM + 2.0 uS
ZERO
cls
print "Running GS..."
print
print "If the FID is RED, reduce the gain. <Ctrl+G>"
print
print "<Ctrl+Q> to exit GS and continue with COSY."
STARTGS
REM now run the COSYPH experiment
REM initial t1 delay
d1=1.0 'microseconds
REM Number of t1 slices SI/4=TD/4, TD=t
x=t/4
REM calculate the run time with the revised RD.
w=n*x*(a+d0+x*d/2e6)/60
REM print decent significant figures.
w=(10*w - (10*w)%1)/10
REM define nucleus for 2nd dimension
h$="H1"
REM Write the 2D file header
WH2D f$ x s h$
cls
print "Running COSYPH..."
print
print "Time, in minutes, to complete the experiment. ",w
print
PROGRAM "clear"
PROGRAM 0 d0 S
PROGRAM + 2.0 uS    
PROGRAM + p0 uS TX a$
PROGRAM + d1 uS    
PROGRAM + p1 uS TX b$
PROGRAM + e uS
PROGRAM + f uS RX
PROGRAM + 1.0 uS ADC RX
PROGRAM + a S RX
PROGRAM "coadd" c$
PROGRAM + 2.0 uS
FOR i=1 to x
ZERO
START
REM Write the data to the 2D file
WF f$ i 0
REM Increment the t1 delay by d and update pulse programer...
d1=d1+d
PROGRAM 3 d1 uS
cls
print "Running COSYPH..."
print
REM calculate the remaining time.
w=n*(x-i)*(a+d0+x*d/2e6)/60
REM print decent significant figures.
w=(10*w - (10*w)%1)/10
print "Slice",i,"of",x,"completed - Time remaining",w,"minutes."
NEXT
cls
print "COSYPH Experiment Completed."
print
print "Data written to file ",f$
print
print "Use NutsMacro aii_cosy.mac, <Ctrl+F5>, for processing."
input "<Enter> to exit the pulse program.",g$
END