;+
; NAME: ircalstrehl
; PURPOSE:
; 	Compute Strehl for an IRCAL image.
;
; USAGE:
;   ircalstrehl,image,strehl,wavelength=wavelength,_extra=_extra
;
;	
; INPUTS:
; KEYWORDS:
; OUTPUTS:
;
; HISTORY:
; 	Began 2005-07-26 03:27:45 by Marshall Perrin 
;-
PRO ircalstrehl,image,strehl,wavelength=wavelength,_extra=_extra

	sz = size(image)
	if sz[1] ne sz[2] then begin
		findmaxstar,image,x,y
		size = x < (sz[1]-x) < y < (sz[2]-y)
		size -= 1
		ic = imcut(image,2*size,x,y)
	endif else ic = image
	;stop

	fixpix,ic,0,ic2,/nan
	ic = temporary(ic2)

	fitstrehl,ic,strehl,wavelength=wavelength,$
		teldiam=3.0,telsecdiam=0.8,platescale=0.0754,_extra=_extra

end