parameter string=gal string=field istop single string=version

! Create directories and setup procedure

if single==1
  $/home/holtz/bin/phot /home/holtz/archival/single/{field}/xvista
  cd /home/holtz/archival/single/{field}/xvista
else
  $/home/holtz/bin/phot /home/holtz/archival/{gal}/{field}/xvista
  cd /home/holtz/archival/{gal}/{field}/xvista
end_if
printf 'parameter ichip' >procedure/{field}setup.pro
printf 'call archival {gal} {field} {version}' >>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
!if istop==0|istop>=1|istop==-1
  call archival {gal} {field} {version}
!end_if

! Compute image shifts. Product is cross.dat file
if istop==0|istop>=2|istop==-2
  call docheck 1 1
  call doshift 1 1
end_if

! Products are msk/*.msk files
if istop==0|istop>=3|istop==-3
  call doprep 1 4 3
end_if

! 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.
if istop==0|istop>=4|istop==-4
  call piccrs
end_if

! 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.
if istop==0|istop>=5|istop==-5
  call doflag
end_if

! Make the psf
! use a single model psf
!call domakepsf 1 4 -1
! use the old hodge 10 psf
!call domakepsf 1 4 -2

! make custom psf
if istop==0|istop>=6|istop==-6
  call domakepsf 1 1 0  !jitter files
  !call domakepsf 1 1 1   !interactive
  !call domakepsf 1 1 2 
  !call domakepsf 1 4 3 
end_if
  
! Do the photometry
if istop==0|istop>=7|istop==-7
  call dophot
end_if

! Aperture corrections (interactive) and transformations
if istop==0|istop>=8|istop==-8
  call cor
end_if

! Process HSTPHOT output
if istop==0|istop>=9|istop==-9
  call gtrans
end_if

ending:
end