# Make a PAM from coeff files.
#
# HRC version.
#
# Quick test version, Richard Hook, January 2004
#
#access the appropriate tasks
stsdas
dither

# Delete temp files
delete tmp*
delete hrc_pam.fits

# First create oversize input
ctio
imcreate tmp.fits 2 1100 1100

# Put some values into the header
hedit tmp.fits exptime 1.0 add+ update+

# Drizzle
unlearn drizzle
drizzle.outnx=1200
drizzle.outny=1200
drizzle tmp.fits tmp_drz.fits outweig=tmp_wht.fits coeffs="HRC_coeffs"

# Blot back
unlearn blot
blot.outnx=1024
blot.outny=1024
blot.expout=1.0
blot tmp_wht.fits tmp_blt.fits coeffs="HRC_coeffs"

imarith 1.0 / tmp_blt.fits hrc_pam.fits

delete tmp*
#