anchor.diff:         ! Compute Scattering weights in the look up table grid but
          ! with the correct albedo. amfdiag is used to skip pixel
          ! ---------------------------------------------------------
  -       CALL compute_scatt ( nt, nx, albedo, lat, sza, vza, saa, vaa, l2ctp, l2cfr, terrain_height, amfdiag, &
  +       CALL compute_scatt ( nt, nx, albedo, lat, sza, vza, saa, vaa, l2ctp, l2cfr, terrain_height, amfgeo, amfdiag, &
               scattw)
   
          stop
  
positive.code:     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
  
      ! ------------------
      ! Modified variables
  
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
            ! ----------------------------------------------
  