path := $libdir/macro,. # ==== setting output and test-flag ==== output := gNhsqc test_ft := n # test-ft for phasing with XEASY reduce_flag := n # reduce size for test-ft by factor of 2 np := 512 ni := 256 # ==== reading of data ==== read vnmr fid {$np}c {$ni}c byteswap # byteswap is needed on Linux write integer tmp readkay integer tmp {$np}c {$ni}c # gradsort for echo/anti-echo # ==== setting of parameters ==== #DIMENSION: 1 (1H) si(1) = 2048 # frequency domain size aqm(1) := RSHc # acquisition mode (ST/TP, TPPI, RSH) w0(1) = 750.5725085 # frequency delta(1) = 0.000095000 # increment in time domain ppmmax(1) = 11.769 # maximum ppm for calibration #DIMENSION: 2 (15N) si(2) = 1024 # frequency domain size aqm(2) := RSH # acquisition mode (ST/TP, TPPI, RSH) w0(2) = 76.0548503 # frequency delta(2) = 0.000433839 # increment in time domain ppmmax(2) = 133.120 # maximum ppm for calibration set n_dim= 2 window1 := "window sin 90" window2 := "window sin 90" ph0_d(1) = -79.3 ph1_d(1) = 0.0 ph0_d(2) = 0 ph1_d(2) = 0 # ==== processing ==== dim 1 phase 90 #suppress cos 64 # optional solvent suppression dim 2 phase 90 mul 0.5 1 # ==== ft ==== strip_ft 1 1..$si(1)/2 # FT and extract left half strip_ft 2 1..$si(2) if ('$test_ft'!='y') then dim 1 cflatt iterative polynom 5 dim 2 cflatt iterative polynom 5 end if dim 1 2 scale noise if ('$test_ft'=='y') then write easy16 $output.rr r r byteswap write easy16 $output.ri r i byteswap write easy16 $output.ir i r byteswap else write easy16 $output byteswap system echo 'Identifier for dimension w1 ... N15' >> $output.3D.param system echo 'Identifier for dimension w2 ... H1' >> $output.3D.param end if