procedure tmulti (imgFilQ, extA, extB, extC, extD, skyFilQ, a1Q, a2Q, a3Q, a4Q, SNliQ, gaiQ, rnQ, XcQ, YcQ, radQ, XmaxQ, YmaxQ, maskZerQ, displFiQ, fQ, outA, outRMA, outB, outRMB, outC, outRMC, outD, outRMD, fldQ)

string imgFilQ	{"list3",prompt="List of images (abbreviated names)"}
string extA	{"gs", prompt="    Set A: Extension for imgs (- for none)"}
string extB	{"gc", prompt="    Set B: Extension for imgs (- for none)"}
string extC	{"us", prompt="    Set C: Extension for imgs (- for none)"}
string extD	{"uc", prompt="    Set D: Extension for imgs (- for none)"}
string skyFilQ  {"regions1.01a8", prompt="File of sky regions."}
real a1Q	{5, prompt="Aperture 1 (diameter)"}
real a2Q	{8, prompt="Aperture 2 (diameter)"}
real a3Q	{12, prompt="Aperture 3 (diameter)"}
real a4Q	{20, prompt="Aperture 4 (diameter)"}
real SNliQ	{1, prompt="Photometry limit (as a S/N ratio)"}
real gaiQ	{2.74, prompt="CCD gain (e-/ADU)"}
real rnQ	{4.8, prompt="Read noise of CCD (rms e-)"}
real XcQ	{470.0, prompt="X-coord of circular aperture centre"}
real YcQ	{515.0, prompt="Y-coord of circular aperture centre"}
real radQ	{440.0, prompt="Radius of circular aperture"}
int XmaxQ	{987, prompt="Image dimension in X"}
int YmaxQ	{988, prompt="Image dimension in Y"}
string maskZerQ	{"zero01a8",prompt="Zero mask image (full name)"}
bool displFiQ	{yes,prompt="Display mask-applied images, one-by-one?"}
int fQ		{1, prompt="    Ximtool frame to display into"}
string outA	{"apsGS", prompt="   Set A: Output pref for flux file 1"}
string outRMA   {"skyGS", prompt="   Set A: Output pref for flux file 2"}
string outB	{"apsGC", prompt="   Set B: Output pref for flux file 1"}
string outRMB   {"skyGC", prompt="   Set B: Output pref for flux file 2"}
string outC	{"apsUS", prompt="   Set C: Output pref for flux file 1"}
string outRMC   {"skyUS", prompt="   Set C: Output pref for flux file 2"}
string outD	{"apsUC", prompt="   Set D: Output pref for flux file 1"}
string outRMD   {"skyUC", prompt="   Set D: Output pref for flux file 2"}
string fldQ  {"01a8", prompt="Output field extn for ALL files"}




begin

	# Variable Declarations:


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

	if (! defpac ("utilities")) {
		print("UTILITIES package not loaded")
		bye
	}

       	if (! defpac ("artdata")) {
		print("ARTDATA package not loaded")
		bye
	}


	print (" ")
	print (" ")
	print ("  TMULTI Start time:")
	time
	print (" ")



	print (" ")
	print ("========================================================")
	print (" ")
	print (" ")



	if (extA!="-") {
	tedge (imgFilX=imgFilQ, abbNamX+, extX=extA, skyFilX=skyFilQ, a1X=a1Q, a2X=a2Q, a3X=a3Q, a4X=a4Q, SNliX=SNliQ, gaiX=gaiQ, rnY=rnQ, outTX=outA//"."//fldQ, outRMX=outRMA//"."//fldQ, XcY=XcQ, YcY=YcQ, radY=radQ, XmaxY=XmaxQ, YmaxY=YmaxQ, maskZerZ=maskZerQ, displFin=displFiQ, fm=fQ)
	}

	if (extB!="-") {
	tedge (imgFilX=imgFilQ, abbNamX+, extX=extB, skyFilX=skyFilQ, a1X=a1Q, a2X=a2Q, a3X=a3Q, a4X=a4Q, SNliX=SNliQ, gaiX=gaiQ, rnY=rnQ, outTX=outB//"."//fldQ, outRMX=outRMB//"."//fldQ, XcY=XcQ, YcY=YcQ, radY=radQ, XmaxY=XmaxQ, YmaxY=YmaxQ, maskZerZ=maskZerQ, displFin=displFiQ, fm=fQ)
	}

	if (extC!="-") {
	tedge (imgFilX=imgFilQ, abbNamX+, extX=extC, skyFilX=skyFilQ, a1X=a1Q, a2X=a2Q, a3X=a3Q, a4X=a4Q, SNliX=SNliQ, gaiX=gaiQ, rnY=rnQ, outTX=outC//"."//fldQ, outRMX=outRMC//"."//fldQ, XcY=XcQ, YcY=YcQ, radY=radQ, XmaxY=XmaxQ, YmaxY=YmaxQ, maskZerZ=maskZerQ, displFin=displFiQ, fm=fQ)
	}

	if (extD!="-") {
	tedge (imgFilX=imgFilQ, abbNamX+, extX=extD, skyFilX=skyFilQ, a1X=a1Q, a2X=a2Q, a3X=a3Q, a4X=a4Q, SNliX=SNliQ, gaiX=gaiQ, rnY=rnQ, outTX=outD//"."//fldQ, outRMX=outRMD//"."//fldQ, XcY=XcQ, YcY=YcQ, radY=radQ, XmaxY=XmaxQ, YmaxY=YmaxQ, maskZerZ=maskZerQ, displFin=displFiQ, fm=fQ)
	}




	print (" ")
	print (" ")
	print ("========================================================")
	print (" ")

	print (" ")
	print (" ")
	print ("  TMULTI End time:")
	time
	print (" ")
	print (" ")

	print (" ")
	print (" ")
	print (" Done with TMULTI.")




end