/*UoM_Project

T2 measurement using Periodic Refocusing of J evolution by Coherance Transfer (PROJECT) using perfect spin echo

Developed By NMR Methodology Group
Department of Chemistry
University of Manchester
United Kingdom
2019

J.A. Aguilar, M. Nilsson, G. Bodenhausen, G. Morris, Chem. Commun.,2012,48,811-813.*/


#include <standard.h>
#include <chempack.h>

static int phs1[16] = {0,2,0,2,1,3,1,3,2,0,2,0,3,1,3,1},  							// 90 degree excitation pulse

	phs2[16] = {1,1,1,1,0,0,0,0,3,3,3,3,2,2,2,2},									// 180 degree refocusing pulse 

	phs3[16] = {1,1,3,3,0,0,2,2,3,3,1,1,2,2,0,0},									//90 degree J modulation refocusing pulse

	phs4[16] = {0,2,0,2,1,3,1,3,2,0,2,0,3,1,3,1};									//Receiver Phase
	
	
	
pulsesequence()
{
	
double 			d1 = getval("d1"),
			d20 = getval("d20"),
			cycles = getval("cycles");
			
	settable(t1,16,phs1);
	settable(t2,16,phs2);
	settable(t3,16,phs3);
	settable(t4,16,phs4);


	getelem(t1,ct,v1);
	getelem(t2,ct,v2);
	getelem(t3,ct,v3);
	getelem(t4,ct,oph);
	

	initval(cycles,v10);


/*START PULSE SEQUENCE*/
obspower(tpwr);

delay(d1);



status (B);

rgpulse(pw,v1,rof1,rof1);

	loop(v10,v11);

		delay(d20);

		rgpulse(2.0*pw,v2,rof1,rof1);

		delay(d20);

		rgpulse(pw,v3,rof1,rof1); 

		delay(d20);

		rgpulse(2.0*pw,v2,rof1,rof2);  

		delay(d20);

	endloop(v11);
	
status (C);

}