The routines in pixon_IDL_source have been 
modified from the original distribution in the following
ways:

010326 dd
  - files of the form "<some name>1.pro"
have been copied into "<some name>.pro"
so that IDL will find the file as needed.

  - an extra "Control Z" was removed from the 
ends of (new) files ndshift.pro and wrap.pro

  - the routine "minf_parabol_D.pro"
was copied into "minf_parabol_d.pro"
so IDL could find it.

  - In the (new) routine pxn.pro the line (273):

 F_PSF = FFT( WRAP(PSF/(TOTAL(PSF))), 1 )

was replaced with:

 ddtemp = WRAP(PSF/(TOTAL(PSF)))
 F_PSF = FFT( ddtemp, 1 )

because IDL kept saying WRAP was not defined variable
even though the function WRAP was compiled...

  - Edit pxncal.pro to change its name from pxncal3
to pxncal in the code.

010327 dd
  - Copy radii_lg.pro into pixon_IDL_source to keep all
needed routines in one place.

  - Copy conv.pro into pixon_IDL_source to keep all
needed routines in one place.


  - Comment out two lines in (the new) pxncal.pro :

;;fmt1  = "(/,'      Number of Pixons = ',f10.2)"
fmt2  = "(/,'        Smallest Pixon = ', f7.3, ' pixels')"
fmt3  =   "('        Largest  Pixon = ', f7.3, ' pixels',/)"


;;print, npxns(pxn), format=fmt1
print, min(pxn),   format=fmt2
print, max(pxn),   format=fmt3

Not sure what was intended here!?!  Maybe uniq values of pxn?


  - Add keyword "maxdels" to pxncal.pro and modify
code so that the range of the dels can be input:

  if n_elements(maxdels) EQ 0 then maxdels=32.0
  dels  = radii_lg(1.0, maxdels, 20)