anchor.diff:      RETURN
     END SUBROUTINE amf_diagnostic
   
  -  SUBROUTINE compute_scatt ( nt, nx, albedo, sza, vza, l2ctp, l2cfr, terrain_height, cli_heights, amfdiag, &
  +  SUBROUTINE compute_scatt ( nt, nx, albedo, lat, sza, vza, l2ctp, l2cfr, terrain_height, amfdiag, &
                                scattw)
   
       USE OMSAO_lininterpolation_module
       USE OMSAO_variables_module, ONLY: verb_thresh_lev
  +    USE OMSAO_omidata_module, ONLY: omi_ozone_amount
   
       IMPLICIT NONE
   
  
positive.code:           local_sza    = cos(d2r*REAL(sza(ixtrack,itime), KIND = r8))  ! JED fix
            local_vza    = cos(d2r*REAL(vza(ixtrack,itime), KIND = r8))  ! JED fix
            local_thg    = REAL(terrain_height(ixtrack,itime), KIND = r8)
            local_chg(:) = cli_heights(ixtrack,itime,:)
  
            ! ----------------------------------------------
            ! Convert pixel terrain height to pressure using
  
negative.code:     ! This should be moved to the program that creates the
      ! look up tables
      ! ------------------------------- ------------------------
      DO issza = 1, sza_dim(1)
         re_sza(issza) = cos(d2r*REAL(lut_sza(sza_dim(1)+1-issza), KIND = r8))  ! JED fix
      END DO
      DO isvza = 1, vza_dim(1)
         re_vza(isvza) = cos(d2r*REAL(lut_vza(vza_dim(1)+1-isvza), KIND = r8)) ! JED fix
      END DO
      DO ispre = 1, clp_dim(1)
         re_pre(ispre) = REAL(lut_clp(clp_dim(1)+1-ispre), KIND = r8)
      END DO
      DO isalt = 1, alt_lay_dim(1)
         re_alt(isalt) = REAL(lut_pre_lay(isalt), KIND=8)
      END DO
  
      ! ---------------
  