qrdir = /work/quickreduce/from_laptop/merge_master_py3/

prep:
	dfits 2024.04.16/*/*33.fits* | fitsort OBJECT FILTER EXPTIME EXPMEAS > all.list
	grep scp all.list  > cals_scp.list


cals:
	$(qrdir)/podi_makecalibrations.py cals_scp.list cals/ -autorotateflat -nonlinearity -keep 


reduce:
#	$(qrdir)/podi_multicollect.py -fromfile=seq1.list -formatout=reduced/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/
	$(qrdir)/podi_multicollect.py -fromfile=sci.list -formatout=reduced/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/

reduce2:
#	$(qrdir)/podi_multicollect.py -fromfile=seq1.list -formatout=reduced/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/
	$(qrdir)/podi_multicollect.py -fromfile=night2.list -formatout=reduced_0419/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/


reduce_0519:
	$(qrdir)/podi_multicollect.py -fromfile=input_0519.list -formatout=reduced_0519/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/
reduce_0527:
	$(qrdir)/podi_multicollect.py -fromfile=input_0527.list -formatout=reduced_0527/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/
reduce_0601:
	$(qrdir)/podi_multicollect.py -fromfile=input_0601.list -formatout=reduced_0601/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/
reduce_0616:
	$(qrdir)/podi_multicollect.py -fromfile=input_0616.list -formatout=reduced_0616/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/


reduce_20250207:
# ls -1 2025.02.07/o*/*33.fits > input_20250207.list
	$(qrdir)/podi_multicollect.py -fromfile=input_20250207.list -formatout=reduced_20250207/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/
reduce_20250210:
#  ls -1 2025.02.10/o*/*33.fits > input_20250210.list
	$(qrdir)/podi_multicollect.py -fromfile=input_20250210.list -formatout=reduced_20250210/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits -nonlinearity -fixwcs -photcalib -nootalevelplots -crj=3 -cals=may2023cals/



reduce_night3:
	$(qrdir)/podi_multicollect.py \
	-fromfile=night_0428.list \
	-formatout=reduced_0428/%OBSID__%OBJECT__%FILTER__%EXPTIME.fits \
	-nonlinearity \
	-fixwcs \
	-photcalib \
	-nootalevelplots \
	-crj=3 \
	-cals=may2023cals/



illumcorr_g:
	$(qrdir)/podi_illumcorr.py \
	-create \
	illumcorr_g1.fits \
	reduced/*odi_g__50.0.fits

illumcorr_r:
	$(qrdir)/podi_illumcorr.py \
	-create \
	illumcorr_r.fits \
	reduced/*odi_r__180.0.fits \
	reduced/*odi_r__120.0.fits



_matchswarp = $(qrdir)/podi_swarpstack.py \
        -pixelscale=0.4 \
        -combine=weighted,sigmaclipmean,sigmaclipmedian \
        -rmguideota \
	-huge

#	-skipota=45,56 

_swarp = $(_matchswarp) -normsky -targetzp=median 

stack_g:
	$(_swarp) stacks/sequence_g1.fits -bgsub=IC,illumcorr_g1.fits @seq1_reduced.list


# prepare filelist
# ls -1 reduced_0419/*__[A-M]g_exposure__odi_g__40.0.fits > 0419_seq1_g.list

# for stacking use this one as template
stack_template:
		$(_swarp) stacks/output_filename.fits -bgsub @input_files.list

stack_0428:
#	$(_swarp) stacks/sequence_0428.fits -bgsub @night_0428.list
#	$(_swarp) stacks/sequence_0428.fits -bgsub @stack_0428_seq1_g.list
	$(_swarp) stacks/sequence_0428_g1.fits -bgsub @stack_0428_seq1_g.list
	$(_swarp) stacks/sequence_0428_i1.fits -bgsub @stack_0428_seq1_i.list
	$(_swarp) stacks/sequence_0428_g2.fits -bgsub @stack_0428_seq2_g.list

stack_5D:
	$(_swarp) stacks/sequence_0419_5D.fits -bgsub @dither_5D.list

# run stack   <output filename> <additional options>   <input list>



ngc6384:
#	$(_matchswarp) -preswarponly stacks/ngc6384/ngc6384_ref.fits @ngc6384_r.list @ngc6384_ha.list
	$(_swarp) -reference=stacks/ngc6384/ngc6384_ref.fits stacks/ngc6384/ngc6384_ha.fits -bgsub=IC,illumcorr_nb659.fits @ngc6384_ha.list
	$(_swarp) -reference=stacks/ngc6384/ngc6384_ref.fits stacks/ngc6384/ngc6384_r.fits -bgsub=IC,illumcorr_r.fits @ngc6384_r.list

ngc3395:
	$(_matchswarp) -preswarponly stacks/ngc3395/ngc3395_ref.fits @ngc3395_r.list @ngc3395_ha.list
	$(_swarp) -reference=stacks/ngc3395/ngc3395_ref.fits stacks/ngc3395/ngc3395_ha.fits -bgsub=IC,illumcorr_nb659.fits @ngc3395_ha.list
	$(_swarp) -reference=stacks/ngc3395/ngc3395_ref.fits stacks/ngc3395/ngc3395_r.fits -bgsub=IC,illumcorr_r.fits @ngc3395_r.list

ngc4151:
	$(_matchswarp) -preswarponly stacks/ngc4151/ngc4151_ref.fits @ngc4151_r.list @ngc4151_ha.list
	$(_swarp) -reference=stacks/ngc4151/ngc4151_ref.fits stacks/ngc4151/ngc4151_ha.fits -bgsub=IC,illumcorr_nb659.fits @ngc4151_ha.list
	$(_swarp) -reference=stacks/ngc4151/ngc4151_ref.fits stacks/ngc4151/ngc4151_r.fits -bgsub=IC,illumcorr_r.fits @ngc4151_r.list

ngc4214:
	$(_matchswarp) -preswarponly stacks/ngc4214/ngc4214_ref.fits @ngc4214_r.list @ngc4214_ha.list
	$(_swarp) -reference=stacks/ngc4214/ngc4214_ref.fits stacks/ngc4214/ngc4214_ha.fits -bgsub=IC,illumcorr_nb659.fits @ngc4214_ha.list
	$(_swarp) -reference=stacks/ngc4214/ngc4214_ref.fits stacks/ngc4214/ngc4214_r.fits -bgsub=IC,illumcorr_r.fits @ngc4214_r.list


prep2025:
	$(_matchswarp) -preswarponly stacks_2025/ref.fits reduced_202502??/*.fits
g1_20250207:
	$(_swarp) -reference=stacks_2025/ref.fits  stacks_2025/stack_20250207_g1.fits -bgsub @stack_20250207_g1.list


i_20250210:
	$(_swarp) -reference=stacks_2025/ref.fits  stacks_2025/stack_20250210_i.fits -bgsub @stack_20250210_odi_i.list
