anchor.diff:      ! Input variables
       ! ---------------
       INTEGER (KIND=i4),                                INTENT(IN) :: nt, nx, CmETA
  -    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1,CmETA), INTENT(IN) :: climatology, &
  -                                                                    cli_temperature, scattw
  -    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1),       INTENT(IN) :: cli_psurface
  +    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1,CmETA), INTENT(IN) :: climatology, scattw
       INTEGER (KIND=i2), DIMENSION (1:nx,0:nt-1),       INTENT(IN) :: amfdiag
   
       ! -----------------------------
  
positive.code:        ! Work out the AMF using the scattering weights and the climatology
         ! Work out Averaging Kernels
         ! -----------------------------------------------------------------
         CALL compute_amf ( nt, nx, CmETA, climatology, cli_temperature, cli_psurface, &
              scattw, saoamf, amfdiag, locerrstat)
  
         ! -----------------------------------------------------------------
  
negative.code:     
    END SUBROUTINE COMPUTE_SCATT
  
    SUBROUTINE compute_amf ( nt, nx, CmETA, climatology, &
                             scattw, saoamf, amfdiag, errstat)
  
      IMPLICIT NONE
  