C*****************************************************************************
C
C       CTFFIND2 - determines defocus and astigmatism for images of
C       arbitrary size (MRC format). Astigmatic angle is measured form
C       x axis (same convetiones as in the MRC 2D image processing 
C       programs).
C
C       CARD 1: Input file name for image
C       CARD 2: Output file name to check result
C       CARD 3: CS[mm], HT[kV], AmpCnst, XMAG, DStep[um]
C       CARD 4: Box, ResMin[A], ResMax[A], dFMin[A], dFMax[A], FStep
CHENN>
C       card 5: iNoAstig,dFStartX,dFStartY,dFStartAst
CHENN<
C
C               The output image file to check the result of the fitting
C               shows the filtered average power spectrum of the input
C               image in one half, and the fitted CTF (squared) in the
C               other half. The two halfs should agree very well for a
C               sucessfull fit.
C
C               CS: Spherical aberration coefficient of the objective in mm
C               HT: Electron beam voltage in kV
C               AmpCnst: Amount of amplitude contrast (fraction). For ice
C                        images 0.07, for negative stain about 0.15.
C               XMAG: Magnification of original image
C               DStep: Pixel size on scanner in microns
C               Box: Tile size. The program devides the image into square
C                    tiles and calculates the average power spectrum. Tiles
C                    with a significatly higher or lower variance are 
C                    excluded; these are parts of the image which are unlikely
C                    to contain useful information (beam edge, film number 
C                    etc). IMPORTANT: Box must have a value of power of 2.
C               ResMin: Low resolution end of data to be fitted.
C               ResMaX: High resolution end of data to be fitted.
C               dFMin: Starting defocus value for grid search in Angstrom. 
C                      Positive values represent an underfocus. The program
C                      performs a systematic grid search of defocus values 
C                      and astigmatism before fitting a CTF to machine 
C                      precision.
C               dFMax: End defocus value for grid search in Angstrom.
C               FStep: Step width for grid search in Angstrom.
CHENN>
C               iNoAstig: 0=usual operation
C                         1=searches only defocus, does not vary astigmatism.
C               dfStart*: Initial defocus values for iNoAstig=1 situation
CHENN<
C
C*****************************************************************************
C       example command file (UNIX):
C
C       #!/bin/csh -x
C       #
C       #   ctffind2
C       #
C       time /public/image/bin/ctffind2.exe << eof
C       image.mrc
C       power.mrc
C       2.6,200.0,0.07,60000.0,28.0
C       128,100.0,15.0,30000.0,90000.0,5000.0
CHENN>
C       0,1,1,1
CHENN<
C       eof
C       #
C
C       or
C
C       #!/bin/csh -x
C       #
C       #   ctffind3
C       #
C       time /public/image/bin/ctffind2.exe << eof
C       image.mrc
C       power.mrc
C       2.6,200.0,0.07,60000.0,28.0
C       128,100.0,15.0,30000.0,90000.0,5000.0
C       1,5780.0,4650.0,27.0
C       eof
C       #
C
C*****************************************************************************