pulsesequence() { char shname1[MAXSTR], f1180[MAXSTR], SE_flg[MAXSTR]; int icosel = 0, t1_counter, phase; double d2_init=0.0, pwS4,pwS5,pwS6,pwS7, lambda = getval("lambda"), gzlvl1 = getval("gzlvl1"), gzlvl2 = getval("gzlvl2"), gzlvl3 = getval("gzlvl3"), gzlvl5 = getval("gzlvl5"), gt1 = getval("gt1"), gt3 = getval("gt3"), gt5 = getval("gt5"), gstab = getval("gstab"), shpw1,shlvl1=getval("shlvl1"), pwClvl = getval("pwClvl"), pwNlvl = getval("pwNlvl"), pwN = getval("pwN"), ni = getval("ni"), d2 = getval("d2"), tau1 = getval("tau1"); getstr("shname1", shname1); getstr("SE_flg",SE_flg); getstr("f1180",f1180); phase = (int) (getval("phase") + 0.5); settable(t3,2,phi3); settable(t10,1,phi10); settable(t12,2,phi12); settable(t13,2,phi13); /* INITIALIZE VARIABLES */ shpw1 = getval("shpw1"); pwS4 = h_shapedpw("eburp2",4.0,3.5,zero, 0.0, 0.0); pwS5 = h_shapedpw("pc9f",4.0,3.5,zero, 2.0e-6, 0.0); pwS6 = h_shapedpw("reburp",4.0,3.5,zero, 0.0, 0.0); pwS7 = c_shapedpw2("isnob5",40.0,-125.0,"isnob5",40.0,0.0 , two, 0.0, 0.0); if (SE_flg[0] =='y') { if (phase == 2) {tsadd(t10,2,4); icosel = +1;} else icosel = -1; } else { if (phase == 2) {tsadd(t3,1,4); icosel=1;} } tau1 = d2; if((f1180[A] == 'y') ) { tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; } tau1 = tau1; if( ix == 1) d2_init = d2; t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 ); if(t1_counter % 2) { tsadd(t3,2,4);tsadd(t12,2,4); tsadd(t13,2,4); } status(A); decpower(pwClvl); dec2power(pwNlvl); set_c13offset("co"); zgradpulse(gzlvl3, gt3); delay(d1-gt3); lk_hold(); rcvroff(); h_shapedpulse("pc9f",4.0,3.5,zero, 2.0e-6, 0.0); delay(lambda-pwS5*0.5-pwS6*0.4); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0); delay(lambda-pwS5*0.5-pwS6*0.4); h_shapedpulse("pc9f_",4.0,3.5,one, 0.0, 0.0); obspower(shlvl1); /**************************************************************************/ /** Sensitivity enhanced version **********************/ /**************************************************************************/ if (SE_flg[0] == 'y') { shaped_pulse(shname1,shpw1,zero,0.0,0.0); dec2rgpulse(pwN,t3,0.0,0.0); delay(tau1*0.5); shaped_pulse(shname1,shpw1,two,0.0,0.0); c_shapedpulse2("isnob5",40.0,-125.0,"isnob5",40.0,0.0 , two, 0.0, 0.0); delay(tau1*0.5); zgradpulse(-gzlvl1, gt1); delay(pwS4-shpw1-pwS7-gt1); dec2rgpulse(pwN*2.0,zero,0.0,0.0); h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0); dec2rgpulse(pwN, t10, 0.0, 0.0); } /**************************************************************************/ /** standard INEPT-based version **********************/ /**************************************************************************/ else { if (ni < 1.0) { dec2rgpulse(pwN,t3,0.0,0.0); dec2rgpulse(pwN*2.0, zero, 0.0, 0.0); dec2rgpulse(pwN, zero, 0.0, 0.0); } else { if (tau1 < shpw1) { shaped_pulse(shname1,shpw1,two,0.0,0.0); shaped_pulse(shname1,shpw1,zero,0.0,0.0); dec2rgpulse(pwN,t3,0.0,0.0); delay(tau1); dec2rgpulse(pwN*3.0, zero, 0.0, 0.0); } else { if (tau1*0.5 < (pwS7+shpw1*0.5)) { shaped_pulse(shname1,shpw1,two,0.0,0.0); dec2rgpulse(pwN,t3,0.0,0.0); delay(tau1*0.5-shpw1*0.5); shaped_pulse(shname1,shpw1,zero,0.0,0.0); delay(tau1*0.5-shpw1*0.5); dec2rgpulse(pwN*3.0, zero, 0.0, 0.0); } else { shaped_pulse(shname1,shpw1,two,0.0,0.0); dec2rgpulse(pwN,t3,0.0,0.0); delay(tau1*0.5-shpw1*0.5); shaped_pulse(shname1,shpw1,zero,0.0,0.0); c_shapedpulse2("isnob5",40.0,-125.0,"isnob5",40.0,0.0 , two, 0.0, 0.0); delay(tau1*0.5-shpw1*0.5-pwS7); dec2rgpulse(pwN*3.0, zero, 0.0, 0.0); } } } } /**************************************************************************/ if (SE_flg[0] == 'y') { zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); dec2rgpulse(pwN, one, 0.0, 0.0); h_shapedpulse("eburp2_",4.0,3.5,one, 0.0, 0.0); txphase(zero); dec2phase(zero); delay(lambda-pwS4*0.5-pwS6*0.4); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); dec2phase(t10); delay(lambda-pwS4*0.5-pwS6*0.4); h_shapedpulse("eburp2",4.0,3.5,zero, 0.0, 0.0); delay((gt1/10.0) + 1.0e-4 +gstab + 2.0*GRADIENT_DELAY + POWER_DELAY); h_shapedpulse("reburp",4.0,3.5,zero, 0.0, 0.0); zgradpulse(icosel*gzlvl2, gt1/10.0); /* 2.0*GRADIENT_DELAY */ delay(gstab); } else { h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5-POWER_DELAY-1.0e-4); } dec2power(dpwr2); /* POWER_DELAY */ lk_sample(); if (SE_flg[0] == 'y') setreceiver(t13); else setreceiver(t12); statusdelay(C,1.0e-4 ); }
pulsesequence() { char shname1[MAXSTR], f1180[MAXSTR], f2180[MAXSTR], SE_flg[MAXSTR], dec_flg[MAXSTR], CT_flg[MAXSTR]; int icosel=1, t1_counter, t2_counter, ni2 = getval("ni2"), phase; double d2_init=0.0, d3_init=0.0, pwS1,pwS2,pwS3,pwS4,pwS5,pwS6, lambda = getval("lambda"), CTdelay = getval("CTdelay"), gzlvl1 = getval("gzlvl1"), gzlvl2 = getval("gzlvl2"), gzlvl3 = getval("gzlvl3"), gzlvl4 = getval("gzlvl4"), gzlvl5 = getval("gzlvl5"), gzlvl6 = getval("gzlvl6"), gt1 = getval("gt1"), gt2 = getval("gt2"), gt3 = getval("gt3"), gt4 = getval("gt4"), gt5 = getval("gt5"), gt6 = getval("gt6"), gstab = getval("gstab"), tpwrsf = getval("tpwrsf"), shlvl1, shpw1, pwClvl = getval("pwClvl"), pwC = getval("pwC"), pwNlvl = getval("pwNlvl"), pwN = getval("pwN"), dpwr2 = getval("dpwr2"), d2 = getval("d2"), shbw = getval("shbw"), shofs = getval("shofs")-4.77, timeTN = getval("timeTN"), tauC = getval("tauC"), tau1 = getval("tau1"), tau2 = getval("tau2"), taunh = getval("taunh"); getstr("shname1", shname1); getstr("SE_flg",SE_flg); getstr("dec_flg",dec_flg); getstr("CT_flg",CT_flg); getstr("f1180",f1180); getstr("f2180",f2180); phase = (int) (getval("phase") + 0.5); settable(t1,2,phi1); settable(t3,4,phi3); settable(t10,1,phi10); settable(t12,4,phi12); settable(t13,4,phi13); /* INITIALIZE VARIABLES */ shpw1 = pw*8.0; shlvl1 = tpwr; pwS1 = c13pulsepw("co", "ca", "sinc", 90.0); pwS2 = c13pulsepw("co", "ca", "sinc", 180.0); pwS3 = c13pulsepw("ca", "co", "square", 180.0); pwS4 = h_shapedpw("eburp2",shbw,shofs,zero, 0.0, 0.0); pwS6 = h_shapedpw("reburp",shbw,shofs,zero, 0.0, 0.0); pwS5 = h_shapedpw("pc9f",shbw,shofs,zero, 2.0e-6, 0.0); if (SE_flg[0] == 'y') { if ( ni2*1/(sw2)/2.0 > (timeTN-pwS2*0.5-pwS4)) { printf(" ni2 is too big. Make ni2 equal to %d or less.\n", ((int)((timeTN-pwS2*0.5-pwS4)*2.0*sw2))); psg_abort(1);} } else { if ( ni2*1/(sw2)/2.0 > (timeTN-pwS2*0.5)) { printf(" ni2 is too big. Make ni2 equal to %d or less.\n", ((int)((timeTN-pwS2*0.5)*2.0*sw2))); psg_abort(1);} } if (CT_flg[0] == 'y') { if ( ni*1/(sw1) > (CTdelay)) { printf(" ni is too big. Make ni equal to %d or less.\n", ((int)(CTdelay*sw1))); psg_abort(1);} } if (phase == 1) ; if (phase == 2) tsadd(t1,1,4); if (SE_flg[0] =='y') { if (phase2 == 2) {tsadd(t10,2,4); icosel = +1;} else icosel = -1; } else { if (phase2 == 2) tsadd(t3,1,4); } tau1 = d2; if((f1180[A] == 'y') ) { tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; } tau1 = tau1; tau2 = d3; if((f2180[A] == 'y') ) { tau2 += ( 1.0 / (2.0*sw2) ); if(tau2 < 0.2e-6) tau2 = 0.0; } tau2 = tau2; if( ix == 1) d2_init = d2; t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 ); if(t1_counter % 2) { tsadd(t1,2,4); tsadd(t12,2,4); tsadd(t13,2,4); } if( ix == 1) d3_init = d3; t2_counter = (int) ( (d3-d3_init)*sw2 + 0.5 ); if(t2_counter % 2) { tsadd(t3,2,4); tsadd(t12,2,4); tsadd(t13,2,4); } status(A); rcvroff(); decpower(pwClvl); decoffset(dof); dec2power(pwNlvl); dec2offset(dof2); obspwrf(tpwrsf); decpwrf(4095.0); obsoffset(tof); set_c13offset("co"); zgradpulse(gzlvl2, gt2); delay(1.0e-4); delay(d1-gt2); lk_hold(); h_shapedpulse("pc9f",shbw,shofs,zero, 2.0e-6, 0.0); delay(lambda-pwS5*0.5-pwS6*0.4); h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0); delay(lambda-pwS5*0.5-pwS6*0.4); h_shapedpulse("pc9f_",shbw,shofs,one, 0.0, 0.0); obspower(shlvl1); /**************************************************************************/ dec2rgpulse(pwN,zero,0.0,0.0); zgradpulse(gzlvl6,gt6); delay(timeTN-pwS2*0.5-gt6); sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6); delay(timeTN-pwS2*0.5-taunh*0.5-shpw1); shaped_pulse(shname1,shpw1,two,0.0,0.0); zgradpulse(gzlvl6, gt6); delay(taunh*0.5-gt6); dec2rgpulse(pwN,zero,0.0,0.0); shaped_pulse(shname1,shpw1,zero,0.0,0.0); /**************************************************************************/ /*** CO -> CA transfer *********************************/ /**************************************************************************/ c13pulse("co", "ca", "sinc", 90.0, zero, 2.0e-6, 0.0); zgradpulse(-gzlvl4, gt4); decphase(zero); delay(tauC - gt4 - 0.5*10.933*pwC); decrgpulse(pwC*158.0/90.0, zero, 0.0, 0.0); decrgpulse(pwC*171.2/90.0, two, 0.0, 0.0); decrgpulse(pwC*342.8/90.0, zero, 0.0, 0.0); /* Shaka 6 composite */ decrgpulse(pwC*145.5/90.0, two, 0.0, 0.0); decrgpulse(pwC*81.2/90.0, zero, 0.0, 0.0); decrgpulse(pwC*85.3/90.0, two, 0.0, 0.0); zgradpulse(-gzlvl4, gt4); decphase(one); delay(tauC - gt4 - 0.5*10.933*pwC - WFG_START_DELAY - 2.0e-6); c13pulse("co", "ca", "sinc", 90.0, one, 2.0e-6, 0.0); /**************************************************************************/ /* xxxxxxxxxxxxxxxxxxxx 13Ca EVOLUTION xxxxxxxxxxxxxxxxxx */ /**************************************************************************/ set_c13offset("ca"); if (CT_flg[0] == 'y') { c13pulse("ca", "co", "square", 90.0, t1, 2.0e-6, 0.0); delay(tau1*0.5); sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0); delay(CTdelay*0.5); c13pulse("cab", "co", "square", 180.0, zero, 0.0, 0.0); delay((CTdelay-tau1)*0.5); c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0); c13pulse("ca", "co", "square", 90.0, zero, 0.0, 0.0); } else { if (dec_flg[0] == 'y') { c13pulse("ca", "co", "square", 90.0, t1, 2.0e-6, 2.0e-6); delay(tau1*0.5); c_shapedpulse2("isnob5",30.0,-31.0,"isnob5",30.0,120.0 , zero, 0.0, 0.0); dec2rgpulse(pwN*2.0,zero,0.0,0.0); delay(tau1*0.5); c13pulse("ca", "co", "square", 180.0, zero, 2.0e-6, 2.0e-6); dec2rgpulse(pwN*2.0,two,0.0,0.0); c_shapedpulse2("isnob5",30.0,-31.0,"isnob5",30.0,120.0 , two, 0.0, 0.0); c13pulse("ca", "co", "square", 90.0, zero, 2.0e-6, 2.0e-6); } else { c13pulse("ca", "co", "square", 90.0, t1, 2.0e-6, 2.0e-6); delay(tau1*0.5); sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6); delay(tau1*0.5); c13pulse("ca", "co", "square", 180.0, zero, 2.0e-6, 2.0e-6); c13pulse("co", "ca", "sinc", 180.0, zero, 2.0e-6, 2.0e-6); if (pwN*2.0 > pwS2) delay(pwN*2.0 - pwS2); c13pulse("ca", "co", "square", 90.0, zero, 2.0e-6, 2.0e-6); } } set_c13offset("co"); zgradpulse(gzlvl3, gt3); delay(1.0e-4); /**************************************************************************/ /*** CO -> CA transfer *********************************/ /**************************************************************************/ c13pulse("co", "ca", "sinc", 90.0, one, 2.0e-6, 0.0); zgradpulse(gzlvl4, gt4*0.7); delay(tauC - gt4*0.7 - 0.5*10.933*pwC); decrgpulse(pwC*158.0/90.0, zero, 0.0, 0.0); decrgpulse(pwC*171.2/90.0, two, 0.0, 0.0); decrgpulse(pwC*342.8/90.0, zero, 0.0, 0.0); /* Shaka 6 composite */ decrgpulse(pwC*145.5/90.0, two, 0.0, 0.0); decrgpulse(pwC*81.2/90.0, zero, 0.0, 0.0); decrgpulse(pwC*85.3/90.0, two, 0.0, 0.0); zgradpulse(gzlvl4, gt4*0.7); delay(tauC - gt4*0.7 - 0.5*10.933*pwC - WFG_START_DELAY - 2.0e-6); /* WFG_START_DELAY */ c13pulse("co", "ca", "sinc", 90.0, zero, 2.0e-6, 0.0); /**************************************************************************/ obspower(shlvl1); shaped_pulse(shname1,shpw1,zero,0.0,0.0); dec2rgpulse(pwN,t3,0.0,0.0); if (SE_flg[0] == 'y') { delay(tau2*0.5); c13pulse("ca", "co", "square", 180.0, zero, 0.0, 0.0); delay(taunh*0.5-pwS3); shaped_pulse(shname1,shpw1,two,0.0,0.0); delay(timeTN-pwS2*0.5-shpw1-taunh*0.5-gt1-1.0e-4); zgradpulse(-gzlvl1, gt1); delay(1.0e-4); sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6); delay(timeTN-pwS2*0.5-tau2*0.5-pwS4); h_shapedpulse("eburp2",shbw,shofs,zero, 2.0e-6, 0.0); dec2rgpulse(pwN, t10, 0.0, 0.0); } else { delay(tau2*0.5); c13pulse("ca", "co", "square", 180.0, zero, 0.0, 0.0); delay(taunh*0.5-pwS3); shaped_pulse(shname1,shpw1,two,0.0,0.0); delay(timeTN-pwS2*0.5-taunh*0.5-shpw1); sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6); delay(timeTN-pwS2*0.5-tau2*0.5); dec2rgpulse(pwN, zero, 0.0, 0.0); } /**************************************************************************/ if (SE_flg[0] == 'y') { zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); dec2rgpulse(pwN, one, 0.0, 0.0); h_shapedpulse("eburp2_",shbw,shofs,one, 0.0, 0.0); txphase(zero); dec2phase(zero); delay(lambda-pwS4*0.5-pwS6*0.4); h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); dec2phase(t10); delay(lambda-pwS4*0.5-pwS6*0.4); h_shapedpulse("eburp2",shbw,shofs,zero, 0.0, 0.0); delay((gt1/10.0) + 1.0e-4 +gstab + 2.0*GRADIENT_DELAY + POWER_DELAY); h_shapedpulse("reburp",shbw,shofs,zero, 0.0, rof2); zgradpulse(icosel*gzlvl2, gt1/10.0); /* 2.0*GRADIENT_DELAY */ delay(gstab); } else { h_shapedpulse("eburp2",shbw,shofs,zero, 2.0e-6, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5-POWER_DELAY-1.0e-4); } dec2power(dpwr2); /* POWER_DELAY */ lk_sample(); if (SE_flg[0] == 'y') setreceiver(t13); else setreceiver(t12); rcvron(); statusdelay(C,1.0e-4 ); }
pulsesequence() { char shname1[MAXSTR], f1180[MAXSTR], f2180[MAXSTR], SE_flg[MAXSTR]; int icosel = 0, t1_counter, t2_counter, ni2 = getval("ni2"), phase; double d2_init=0.0, d3_init=0.0, pwS1,pwS2,pwS3,pwS4,pwS5,pwS6,pwS7,pwS8, lambda = getval("lambda"), gzlvl1 = getval("gzlvl1"), gzlvl2 = getval("gzlvl2"), gzlvl3 = getval("gzlvl3"), gzlvl4 = getval("gzlvl4"), gzlvl5 = getval("gzlvl5"), gzlvl6 = getval("gzlvl6"), gt1 = getval("gt1"), gt3 = getval("gt3"), gt4 = getval("gt4"), gt5 = getval("gt5"), gt6 = getval("gt6"), gstab = getval("gstab"), shlvl1 = getval("shlvl1"), shpw1 = getval("shpw1"), pwClvl = getval("pwClvl"), pwNlvl = getval("pwNlvl"), pwN = getval("pwN"), d2 = getval("d2"), timeTN = getval("timeTN"), Delta, tauNCO = getval("tauNCO"), tauC = getval("tauC"), tau1 = getval("tau1"), tau2 = getval("tau2"), taunh = getval("taunh"); getstr("shname1", shname1); getstr("SE_flg",SE_flg); getstr("f1180",f1180); getstr("f2180",f2180); phase = (int) (getval("phase") + 0.5); settable(t1,2,phi1); settable(t2,1,phi2); settable(t3,4,phi3); settable(t4,8,phi4); settable(t10,1,phi10); settable(t12,8,phi12); settable(t13,8,phi13); /* INITIALIZE VARIABLES */ Delta = timeTN-tauNCO; pwS1 = c13pulsepw("ca", "co", "square", 90.0); pwS2 = c13pulsepw("ca", "co", "square", 180.0); pwS3 = c13pulsepw("co", "ca", "sinc", 180.0); pwS7 = c13pulsepw("co", "ca", "sinc", 90.0); pwS8 = c_shapedpw("reburp",60.0 ,-135.0,zero, 0.0, 0.0); pwS4 = h_shapedpw("eburp2",4.0,3.5,zero, 0.0, 0.0); pwS6 = h_shapedpw("reburp",4.0,3.5,zero, 0.0, 0.0); pwS5 = h_shapedpw("pc9f",4.0,3.5,zero, 2.0e-6, 0.0); if (SE_flg[0] == 'y') { if ( ni2*1/(sw2)/2.0 > (timeTN-Delta-pwS3-pwS4)) { printf(" ni2 is too big. Make ni2 equal to %d or less.\n", ((int)((timeTN-Delta-pwS3-pwS4)*2.0*sw2))); psg_abort(1);} } else { if ( ni2*1/(sw2)/2.0 > (timeTN-Delta-pwS3)) { printf(" ni2 is too big. Make ni2 equal to %d or less.\n", ((int)((timeTN-Delta-pwS3)*2.0*sw2))); psg_abort(1);} } if (phase == 1) ; if (phase == 2) {tsadd(t1,1,4);} if (SE_flg[0] =='y') { if (phase2 == 2) {tsadd(t10,2,4); icosel = +1;} else icosel = -1; } else { if (phase2 == 2) {tsadd(t3,1,4); icosel = 1;} } tau1 = d2; if((f1180[A] == 'y') ) { tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; } tau1 = tau1; tau2 = d3; if((f2180[A] == 'y') ) { tau2 += ( 1.0 / (2.0*sw2) ); if(tau2 < 0.2e-6) tau2 = 0.0; } tau2 = tau2; if( ix == 1) d2_init = d2; t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 ); if(t1_counter % 2) { tsadd(t1,2,4); tsadd(t12,2,4); tsadd(t13,2,4); } if( ix == 1) d3_init = d3; t2_counter = (int) ( (d3-d3_init)*sw2 + 0.5 ); if(t2_counter % 2) { tsadd(t3,2,4); tsadd(t12,2,4); tsadd(t13,2,4); } status(A); decpower(pwClvl); dec2power(pwNlvl); set_c13offset("co"); zgradpulse(gzlvl6, gt6); delay(1.0e-4); delay(d1-gt6); lk_hold(); rcvroff(); h_shapedpulse("pc9f",4.0,3.5,zero, 2.0e-6, 0.0); delay(lambda-pwS5*0.5-pwS6*0.4); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0); delay(lambda-pwS5*0.5-pwS6*0.4); h_shapedpulse("pc9f_",4.0,3.5,one, 2.0e-6, 0.0); obspower(shlvl1); /**************************************************************************/ /* xxxxxxxxxxxxxxxxxxxxxx N-> CA transfer xxxxxxxxxxxxxxxxxx */ /**************************************************************************/ set_c13offset("ca"); dec2rgpulse(pwN,zero,0.0,0.0); delay(timeTN); sim3_c13pulse("", "ca", "co", "square", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6); delay(Delta); c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0); delay(timeTN-Delta-pwS3-taunh*0.5-shpw1); shaped_pulse(shname1,shpw1,two,0.0,0.0); delay(taunh*0.5); dec2rgpulse(pwN,zero,0.0,0.0); shaped_pulse(shname1,shpw1,zero,0.0,0.0); /**************************************************************************/ /* xxxxxxxxxxxxxxxxxxxxxx CA->CO TRANSFER xxxxxxxxxxxxxxxxxx */ /**************************************************************************/ set_c13offset("ca"); c13pulse("ca", "co", "square", 90.0, t2, 2.0e-6, 0.0); /* initval(0.0, v2); decstepsize(1.0); dcplrphase(v2); */ zgradpulse(gzlvl4, gt4); delay(tauC*0.5-gt4); c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,zero,0.0,0.0); zgradpulse(gzlvl4, gt4); delay(tauC-gt4); c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,two,0.0,0.0); delay(tauC*0.5); c13pulse("ca", "co", "square", 90.0, one, 0.0, 0.0); zgradpulse(gzlvl3, gt3*3.5); delay(1.0e-4); /**************************************************************************/ /* xxxxxxxxxxxxxxxxxxxxxx 13CO EVOLUTION xxxxxxxxxxxxxxxxxx */ /**************************************************************************/ set_c13offset("co"); c13pulse("co", "ca", "sinc", 90.0, t1, 2.0e-6, 0.0); delay(tau1*0.5); sim3_c13pulse(shname1, "ca", "co", "square", "",shpw1, 180.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0); delay(tau1*0.5); c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0); sim3_c13pulse(shname1, "ca", "co", "square", "",shpw1, 180.0, 0.0, two, zero, zero, 0.0, 0.0); if (pwN*2.0 > pwS2) delay(pwN*2.0-pwS2); c13pulse("co", "ca", "sinc", 90.0, t4, 0.0, 0.0); /***************************************************************************/ /* CA->CO transfer */ /***************************************************************************/ set_c13offset("ca"); c13pulse("ca", "co", "square", 90.0, zero, 2.0e-6, 0.0); initval(0.0, v2); decstepsize(1.0); dcplrphase(v2); zgradpulse(gzlvl3, gt3*2.0); delay(tauC*0.5-gt3*2.0); c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,two,0.0,0.0); zgradpulse(gzlvl3, gt3*2.0); delay(tauC-gt3*2.0); c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,zero,0.0,0.0); delay(tauC*0.5); c13pulse("ca", "co", "square", 90.0, one, 0.0, 0.0); /* dcplrphase(v2); */ /**************************************************************************/ obspower(shlvl1); shaped_pulse(shname1,shpw1,zero,0.0,0.0); dec2rgpulse(pwN,t3,0.0,0.0); delay(tau2*0.5+taunh*0.5); shaped_pulse(shname1,shpw1,two,0.0,0.0); delay(timeTN-shpw1-taunh*0.5-gt1-1.0e-4); zgradpulse(-gzlvl1, gt1); delay(1.0e-4); sim3_c13pulse("", "ca", "co", "square", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6); if (SE_flg[0] == 'y') { delay(Delta); c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0); delay(timeTN-tau2*0.5-pwS4-Delta-pwS3); h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0); dec2rgpulse(pwN, t10, 0.0, 0.0); } else { delay(Delta); c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0); delay(timeTN-tau2*0.5-Delta-pwS3); dec2rgpulse(pwN, zero, 0.0, 0.0); } /**************************************************************************/ if (SE_flg[0] == 'y') { zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); dec2rgpulse(pwN, one, 0.0, 0.0); h_shapedpulse("eburp2_",4.0,3.5,one, 0.0, 0.0); txphase(zero); dec2phase(zero); delay(lambda-pwS4*0.5-pwS6*0.4); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); dec2phase(t10); delay(lambda-pwS4*0.5-pwS6*0.4); h_shapedpulse("eburp2",4.0,3.5,zero, 0.0, 0.0); delay((gt1/10.0) + 1.0e-4 +gstab + 2.0*GRADIENT_DELAY + POWER_DELAY); h_shapedpulse("reburp",4.0,3.5,zero, 0.0, 0.0); zgradpulse(icosel*gzlvl2, gt1/10.0); /* 2.0*GRADIENT_DELAY */ delay(gstab); } else { h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5); h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0); zgradpulse(gzlvl5, gt5); delay(lambda-pwS6*0.4 - gt5-POWER_DELAY-1.0e-4); } dec2power(dpwr2); /* POWER_DELAY */ lk_sample(); if (SE_flg[0] == 'y') setreceiver(t13); else setreceiver(t12); rcvron(); statusdelay(C,1.0e-4 ); }