parameter string=gal string=field

! Create directories and setup procedure
$/home/holtz/bin/phot /home/holtz/archival/{gal}/{field}/xvista
cd /home/holtz/archival/{gal}/{field}/xvista
printf 'parameter ichip' >procedure/{field}setup.pro
printf 'call archival {gal} {field}' >>procedure/{field}setup.pro
$cat procedure/setup.pro | sed '1d' >>procedure/{field}setup.pro
$ln -s /home/holtz/raw/archive/{gal}/{field} archive

! Setup up the directories and load with default procedure, create
!   setup procedure for desired field by scanning through image list
call archival {gal} {field}

! Process HSTPHOT output
!call gtrans
!goto ending

call docheck 1 1
!goto ending

!call domakepsf 1 4 2 1
!goto ending

! Compute image shifts. Product is cross.dat file
call doshift 1 1
goto ending

! Products are msk/*.msk files
call doprep 1 4 3

! Mask out cosmic rays. Products are red/z*.hdr,img frames with CRs set
!  to a very low value. Also makes individual reduced frames for each
!  chip.
call piccrs

! Combine the CR-flagged frames into 3D images, add in masks in the process
! Product are 3D _m0f.fits frames. The individual frames from piccrs (both
! reduced and CR-flagged) are deleted.
call doflag

! Make the psf
! use a single model psf
!call domakepsf 1 4 -1
!call domakepsf 1 4 -2  !link to hodge10 old PSF

! make custom psf
call domakepsf 1 1 0  !read jitter file
call domakepsf 1 4 2 -1
!call domakepsf 1 1 1   !interactive
!call domakepsf 1 1 2 
!call domakepsf 1 4 3 

! Do the photometry
!call dophot

! Aperture corrections (interactive) and transformations
!call cor

ending:
end