anchor.diff:            saoamf(ixtrack,itimes) = SUM(scattw(ixtrack, itimes, 1:CmETA) * &
                                           climatology(ixtrack,itimes,1:CmETA))     / &
                                      SUM(climatology(ixtrack,itimes,1:CmETA))
  +          print*, saoamf(ixtrack,itimes)
          END DO ! Finish xtrack pixel loop
       END DO ! Finish 
       
  
positive.code:     ! ------------------------------------
      albedo       = r8_missval
      climatology  = r8_missval
      cli_heights  = r8_missval
      cli_psurface = r8_missval
      scattw       = r8_missval
      saoamf       = r8_missval
  
negative.code:     
    END SUBROUTINE amf_calculation
    
    SUBROUTINE omi_climatology (climatology, local_psurf, local_temperature, terrain_height, &
                                lat, lon, nt, nx, xtrange, locerrstat)
  
      ! =========================================
  