anchor.diff:            idx_lat = MINVAL(MINLOC(ABS(latvals(1:Cmlat) - lat(ixtrack,itimes) )))
             idx_lon = MINVAL(MINLOC(ABS(lonvals(1:Cmlon) - lon(ixtrack,itimes) )))
   
  -          local_temperature(ixtrack,itimes,1:CmETA) = Temperature(idx_lon,idx_lat,1:CmETA)
  -          ltmp(1:CmETA)                             = Temperature(idx_lon,idx_lat,1:CmETA)
  +!!$          local_temperature(ixtrack,itimes,1:CmETA) = REAL(Temperature(idx_lon,idx_lat,1:CmETA), KIND=r8)
  +!!$          ltmp(1:CmETA)                             = REAL(Temperature(idx_lon,idx_lat,1:CmETA), KIND=r8)
   
             ! ----------------------------------------------
             ! Convert pixel terrain height to pressure using
             ! Xiong suggested to use pressure altitude:
             !  Z = -16 alog10 (P / Po) Z in km and P in hPa.
  +          ! From OMI L1B file.
             ! ----------------------------------------------
             local_psurf(ixtrack,itimes) = 1013.0_r8 * &
                  (10.0_r8 ** (REAL(terrain_height(ixtrack,itimes),KIND=r8) / 1000.0_r8 / (-16.0_r8)))
   
  +          clima_psurf = REAL(Psurface(idx_lon, idx_lat), KIND=r8)
             ! lpre(0:CmETA) is pressure at layer boundaries, same unit as
             ! local_surf, convert from hPa to Pa, it is a CmETA+1 = CmEp1 array
             ! Ap(1:CmEp1) & Bp(1:CmEp1) are coeff at layer boundary for pressure
             ! calculation
             DO ilevel = 0, CmETA
  -             lpre(ilevel) = (Ap(ilevel+1) + local_psurf(ixtrack,itimes) * Bp(ilevel+1))*1.E2 ! Pa
  +             lpre(CmETA-ilevel) = (Ap(ilevel+1) + local_psurf(ixtrack,itimes) * Bp(ilevel+1))*1.E2 ! Pa
  +             clima_lpre(CmETA-ilevel) = (Ap(ilevel+1) + clima_psurf * Bp(ilevel+1))*1.E2 ! Pa
             END DO
   
   
  
positive.code:     ! ----------------------
      locerrstat = pge_errstat_ok
      
      ! ------------------------------------------------------------
      ! Find the Climatology corresponding to each lat and lon pixel
      ! ------------------------------------------------------------
  
negative.code:           saoamf(ixtrack,itimes) = SUM(scattw(ixtrack, itimes, 1:CmETA) * &
                                          climatology(ixtrack,itimes,1:CmETA))     / &
                                     SUM(climatology(ixtrack,itimes,1:CmETA))
            print*, saoamf(ixtrack,itimes)
         END DO ! Finish xtrack pixel loop
      END DO ! Finish 
      
  