anchor.diff:         ! Now it is only needed to interpolate to the pixels of the granule.
          ! It was read there to obtain the dimensions of the number of levels.
          ! ---------------------------------------------------------------------
  -       CALL omi_climatology (climatology, cli_psurface, terrain_height, &
  +       CALL omi_climatology (climatology, climatology_sd, cli_psurface, terrain_height, &
               lat, lon, nt, nx, amfdiag, locerrstat)
          
          ! -------------------------------------
  
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:     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)       :: 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
  
  
  