anchor.diff:      INTEGER (KIND=i4)                                :: locerrstat, itt, spixx, epixx, current_month
       INTEGER (KIND=i2), DIMENSION (1:nx,0:nt-1)       :: amfdiag
       REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: amfgeo
  -    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: l2cfr, l2ctp
  +    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: l2cfr, l2ctp, amf_pro_err
       REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: albedo, cli_psurface
       REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1,CmETA) :: climatology, climatology_sd
       REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1,CmETA) :: scattw !, akernels
  
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:     cli_psurface   = r8_missval
      scattw         = r8_missval
      saoamf         = r8_missval
      amf_pro_err    = r8_missval
      amfgeo         = r8_missval
      amfdiag        = i2_missval
  
  