anchor.diff:      
     END SUBROUTINE amf_calculation
     
  -  SUBROUTINE omi_climatology (climatology, local_psurf, terrain_height, &
  +  SUBROUTINE omi_climatology (climatology, climatology_sd, local_psurf, terrain_height, &
                                 lat, lon, nt, nx, amfdiag, locerrstat)
   
       ! =========================================
  
positive.code:     
    END SUBROUTINE COMPUTE_SCATT
  
    SUBROUTINE compute_amf ( nt, nx, CmETA, climatology, &
                             scattw, saoamf, amfdiag, errstat)
  
      IMPLICIT NONE
  
negative.code:     ! ------------------------------------
      ! Initialize variables that are output
      ! ------------------------------------
      albedo         = r8_missval
      climatology    = r8_missval
      climatology_sd = r8_missval
      cli_psurface   = r8_missval
      scattw         = r8_missval
      saoamf         = r8_missval
      amfgeo         = r8_missval
      amfdiag        = i2_missval
  
      ! -----------------------------------------
      ! If amf_wvl < 0.0 then the slant column is
  