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:           saoamf(ixtrack,itimes) = SUM(scattw(ixtrack, itimes, 1:CmETA) * &
                                          climatology(ixtrack,itimes,1:CmETA))     / &
                                     SUM(climatology(ixtrack,itimes,1:CmETA))
  
         END DO ! Finish xtrack pixel loop
      END DO ! Finish 
  
negative.code:     cli_psurface   = r8_missval
      scattw         = r8_missval
      saoamf         = r8_missval
      amf_pro_err    = r8_missval
      amfgeo         = r8_missval
      amfdiag        = i2_missval
  
  