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) = 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:
  
positive.code:           DO ilevel = 1, CmETA
               CALL ezspline_1d_interpolation (INT(CmETA,KIND=i4), clima_local_heights, &
                    re_tmp, &
                    one, local_heights(ixtrack,itimes,ilevel), &
                    ltmp(ilevel), status)
               CALL ezspline_1d_interpolation (INT(CmETA,KIND=i4), clima_local_heights, &
                    re_gas, &
                    one, local_heights(ixtrack,itimes,ilevel), &
                    lgas(ilevel), status)
               CALL ezspline_1d_interpolation (INT(CmETA,KIND=i4), clima_local_heights, &
                    re_h2o, &
                    one, local_heights(ixtrack,itimes,ilevel), &
                    lh2o(ilevel), status)
            END DO         
           
  
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 
      
  