#!/bin/csh -f

set freq1=88559
set freq2=89115
set inname=n159f_mos
set outname=n159
alias deconv.csh ../../../scripts/deconv.csh
alias idl /Applications/exelis/idl82/bin/idl

goto start

# Map the HCN & HCO+ data:
if (-e deconv.log) mv deconv.log deconv.bak.log
rm -f mossdi.log
foreach freq ($freq1 $freq2)
    foreach file (../20*/${inname}.$freq.cal)
        set dir = $file:h:t
        rm -f $dir.$freq.cal
        ln -s $file $dir.$freq.cal
    end
    deconv.csh freq=$freq start=invert doself=n outfile=$outname.$freq \
        vfirst=222 vlast=250 interact=n delv=0.4 msklev=3 usemaj=6.3 usemin=5.7 usepa=50
    deconv.csh freq=$freq start=cleanup doself=n outfile=$outname.$freq
end
goto end

start:
# Generate the moment maps:
idl<<EOF
makemom,'$outname.$freq2.cmmsk.fits',errfile='$outname.$freq2.cmnse.fits',baseroot='${outname}_hcop_dil',thresh=3,edge=2,/PVMOM0
makemom,'$outname.$freq2.cmmsk.fits',errfile='$outname.$freq2.cmnse.fits',baseroot='${outname}_hcop_smo',thresh=3.5,edge=2,smopar=[9,5],/PVMOM0
makemom,'$outname.$freq1.cmmsk.fits',errfile='$outname.$freq1.cmnse.fits',baseroot='${outname}_hcn_dil',thresh=3,edge=2,/PVMOM0
makemom,'$outname.$freq1.cmmsk.fits',errfile='$outname.$freq1.cmnse.fits',baseroot='${outname}_hcn_smo',thresh=3.5,edge=2,smopar=[9,5],/PVMOM0
EOF
#goto end

# Generate pdf files:
foreach type (mom01 mom0pv)
	set pdflist=()
	foreach file (${outname}*${type}.eps)
		if ($type == 'mom01') then
			epspdf $file
			echo "$file" | enscript -r -s 1 -B -f Courier-Bold@12 -o- | \
				ps2pdf - | pdftk $file:r.pdf stamp - output output.pdf
		else
			epspdf -b $file
			echo "$file" | enscript -r -s 1 -B -i 82 -f Courier-Bold@12 -o- | \
				ps2pdf - | pdftk $file:r.pdf stamp - output output.pdf
		endif
		if (-e output.pdf) mv output.pdf $file:r.pdf
		set pdflist=($pdflist $file:r.pdf)
	end
	if ($#pdflist != 0) then
		gs -q -dSAFER -dNOPAUSE -dBATCH -sOutputFile=$type.pdf \
			-sDEVICE=pdfwrite -c .setpdfwrite -f $pdflist
		echo Plotted images output to $type.pdf
		rm -f $pdflist
	endif
end
goto end


end: