C CONVDECONV *******************************************************************
C
C       Program to read in a list of amplitudes and phases and calculate
C       a new list of amplitudes and phases which represent either the
C       convolution or deconvolution of the original.  The necessary 
C       information to describe the relative (x,y) coordinates of the layers
C       required in the convolution or known to be present in the image which
C       is to be deconvoluted, must be given as input data.
C
C       Formula used for mirror-deconvolution is
C           (P+iQ) = [A+iB-(s/r)*Trans(A+iB)*(C1+iC2)] / [r-(s*s/r)*Trans(C1+iC2)*(C1+iC2)]
C
C                       where C1 = Sum of cos(2*PI*(x/a+y/b)))
C                       and   C2 = Sum of sin(2*PI*(x/a+y/b)))
C                       where (x,y) are the coordinate shifts to be applied to the
C                               map whose structure factors are read from the 
C                               input file in convolution, and which describe the 
C                               coordinates of the layers in deconvolution relative
C                               to the final deconvoluted map obtained,
C                       and r and s are factors for the two layer intensities,
C                       and a and b are the unit cell dimensions
C                       WIEN is a WIENer term to avoid divide-by-zero
C       Input data cards
C         1.  x/a,y/b   fractional shifts of second layer (e.g. 0.0,0.25)
C                       the first layer is assumed to be centered after phase 
C                       phase origin shift.
C         2.  XORI,YORI phase origin in degrees
C         3.  ra,rb,rs  factor from-to-stepsize for layer intensities
C
C       Input file is a .aph file, consisting of IH,IK,AMP,PHASE,IQ,BACK,CTF
C
C  Vx 1.00      06.02.2000      RH      Original version
C  Vx 1.01      18.10.2000      RH      change spelling to WIEN (Wiener filter)
C  Vx 2.00      27.11.2001      HS      mirror-deconvolution
C
C           remember to change date in first format statement
C
C*******************************************************************************