anchor.diff:    END SUBROUTINE COMPUTE_SCATT
   
     SUBROUTINE compute_amf ( nt, nx, CmETA, climatology, climatology_sd, &
  -                           scattw, saoamf, amfdiag, errstat)
  +                           scattw, saoamf, amf_pro_err, amfdiag, errstat)
   
       IMPLICIT NONE
   
  
positive.code:        ! Work out Averaging Kernels
         ! -----------------------------------------------------------------
         CALL compute_amf ( nt, nx, CmETA, climatology, climatology_sd, &
              scattw, saoamf, amfdiag, locerrstat)
  
         ! -----------------------------------------------------------------
         ! Write out scattering weights, altitude grid and averaging kernels
  
negative.code:     ! -----------------------------
      ! Output and modified variables
      ! -----------------------------
      REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1),       INTENT (INOUT) :: saoamf, amf_pro_err
      INTEGER (KIND=i4),                                INTENT (INOUT) :: errstat
  
      ! ---------------
  