pro testisol,movin,movout  ;mea

;*************** PARAMETERS **********************************
; GMSMO - smoothing parameter to define general mask
; GMTHR - threshold to define general mask
; smo=7
 gmsmo=37 	; 3" 2.5" to smooth the penumbra
 gmthr=1060	; max. penumbral intensity in smoothed frame
;*************************************************************

;mea=fltarr(140)
si=size(movin)
movout=intarr(504,480,si(3))

for i=0,si(3)-1 do begin

   oima=movin(*,*,i)	; original image
;   ima=SMOOTH(oima,smo)	; smoothed image (for unsharp masking)
;   ima=(oima-ima)>0	; difference image (unsharp masking)

   ois=SMOOTH(oima,gmsmo)
   mge=ois lt gmthr		;outer general mask
;   mpe=ois gt 550		;inner general mask
;   ima=ima*mge*mpe		;resulting masked image
   ima=oima*mge
;   mea(i)=mean(ima(where(ima)))
   movout(0,0,i)=ima(51:554,85:564)

endfor

end