procedure tsingle (imgScns,b1,b2,b3,s1,e1,s2,e2, s3,e3,inA,outA,clnA,inB,outB,clnB,inC, outC,clnC,inD,outD,clnD,clipV, gnV, rnV, oot)

int imgScns	{3,prompt="Number of image scans (<=3)"}
string b1 	{"u", prompt="    Scan 1: Image name prefix"}
string b2 	{"u", prompt="    Scan 2: Image name prefix"}
string b3 	{"u", prompt="    Scan 3: Image name prefix"}
int s1		{343, prompt="    Scan 1: Image start number"}
int e1		{351, prompt="    Scan 1: Image finish number"}
int s2		{352, prompt="    Scan 2: Image start number"}
int e2		{360, prompt="    Scan 2: Image finish number"}
int s3		{361, prompt="    Scan 3: Image start number"}
int e3		{369, prompt="    Scan 3: Image finish number"}
string inA  {"r", prompt="set A: Extn for input imgs (- for none)"}
string outA {"summed", prompt="set A: Output (cleaned|summed)"}
string clnA {"us", prompt="set A: Output image extension?"}
string inB  {"r", prompt="   set B: Extn for input imgs (- for none)"}
string outB	{"cleaned", prompt="   set B: Output (cleaned|summed)"}
string clnB {"uc", prompt="   set B: Output image extension?"}
string inC  {"g", prompt="set C: Extn for input imgs (- for none)"}
string outC {"summed", prompt="set C: Output (cleaned|summed)"}
string clnC {"gs", prompt="set C: Output image extension?"}
string inD  {"g", prompt="   set D: Extn for input imgs (- for none)"}
string outD	{"cleaned", prompt="   set D: Output (cleaned|summed)"}
string clnD {"gc", prompt="   set D: Output image extension?"}
string clipV	{"avsigclip", prompt="Cleaning method (avsigclip|pclip|crreject)"}
real gnV	{2.74, prompt="  -> crreject: CCD gain (e-/ADU)"}
real rnV	{4.8, prompt="  -> crreject: CCD read noise (e-)"}
string oot 	{"u", prompt="Prefix for all output images?"}





begin


	# Variable Declarations:


	if (! defpac ("images")) {
		print("Package IMAGES needs to be loaded.")
		bye
	}


	print (" ")
	print (" ")
	print ("  TSINGLE Start time:")
	time
	print (" ")



	if ((inA!="-")&&(outA=="summed")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inA,outType="summed",clnNum=1,clnOver-,clnExtn="", clipmode=clipV, gnW=gnV, rnW=rnV, sumNum=1,sumOver-,sumExtn=clnA,outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}

	if ((inA!="-")&&(outA=="cleaned")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inA,outType="cleaned",clnNum=1,clnOver-,clnExtn=clnA, clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn="",outpref=outpref, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}




	if ((inB!="-")&&(outB=="summed")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inB,outType="summed",clnNum=1,clnOver-,clnExtn="", clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn=clnB,outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}

	if ((inB!="-")&&(outB=="cleaned")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inB,outType="cleaned",clnNum=1,clnOver-,clnExtn=clnB, clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn="",outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}




	if ((inC!="-")&&(outC=="summed")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3,  inExtn=inC,outType="summed",clnNum=1,clnOver-,clnExtn="", clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn=clnC,outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}

	if ((inC!="-")&&(outC=="cleaned")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inC,outType="cleaned",clnNum=1,clnOver-,clnExtn=clnC, clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn="",outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}





	if ((inD!="-")&&(outD=="summed")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inD,outType="summed",clnNum=1,clnOver-,clnExtn="", clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn=clnD,outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}

	if ((inD!="-")&&(outD=="cleaned")) {
	tcombine (nScans=imgScns,impref1=b1,impref2=b2,impref3=b3, start1=s1,end1=e1,start2=s2,end2=e2,start3=s3,end3=e3, inExtn=inD,outType="cleaned",clnNum=1,clnOver-,clnExtn=clnD, clipmode=clipV, gnW=gnV, rnW=rnV,sumNum=1,sumOver-,sumExtn="",outpref=oot, maskType="none",maskZer="-",clnmaRMS="--",summaRMS="---",logname="STDOUT", convFits-)
	}


	print (" ")
	print (" ")
	print ("  TSINGLE End time:")
	time
	print (" ")




	print (" ")
	print ("============================================================")
	print (" ")
	print (" ")
	print ("Done with tsingle.")
	print (" ")

end