anchor.diff:      INTEGER (KIND=i4),                                INTENT (IN) :: nt, nx
       INTEGER (KIND=i2), DIMENSION (1:nx,0:nt-1),       INTENT (IN) :: amfdiag
       REAL    (KIND=r4), DIMENSION (1:nx,0:nt-1),       INTENT (IN) :: lat, sza, vza, saa, vaa, terrain_height
  -    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1),       INTENT (IN) :: albedo, l2ctp, l2cfr
  +    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1),       INTENT (IN) :: albedo, l2ctp, l2cfr, amfgeo
   
       ! ------------------
       ! Modified variables
  
positive.code:     RETURN
    END SUBROUTINE amf_diagnostic
  
    SUBROUTINE compute_scatt ( nt, nx, albedo, lat, sza, vza, saa, vaa, l2ctp, l2cfr, terrain_height, amfdiag, &
                               scattw)
  
      USE OMSAO_lininterpolation_module
      USE OMSAO_variables_module, ONLY: verb_thresh_lev
  
negative.code:           ! in function of latitude and ozone column.
            ! -----------------------------------------        
            IF ( ABS(lat(ixtrack,itime)) .GT. 60.0 ) THEN
               toms_idx = MINLOC(ABS(hxxx-omi_ozone_amount(ixtrack,itime)/du/amfgeo(ixtrack,itime)), 1) + 0
            ELSE IF ( ABS(lat(ixtrack,itime)) .GT. 30.0 .AND. ABS(lat(ixtrack,itime)) .LE. 60.0 ) THEN
               toms_idx = MINLOC(ABS(mxxx-omi_ozone_amount(ixtrack,itime)/du/amfgeo(ixtrack,itime)), 1) + 17
            ELSE IF ( ABS(lat(ixtrack,itime)) .LE. 30.0 ) THEN
               toms_idx = MINLOC(ABS(lxxx-omi_ozone_amount(ixtrack,itime)/du)/amfgeo(ixtrack,itime), 1) + 11
            ENDIF
            print*, toms_idx
            ! ----------------------------------------------
            ! Work out relative azimuth angle for this pixel
            ! ----------------------------------------------
  