anchor.diff:      he5stat = HE5_SWrdfld (                                &
            swath_id, TRIM(ADJUSTL(gasdatafieldname)),        &
            he5_start_3d, he5_stride_3d, he5_edge_3d,         &
  -         H2O_profiles(1:Cmlon,1:Cmlat,1:CmETA) )
  +         dummy_h2opr(1:Cmlon,1:Cmlat,1:CmETA) )
   
       ! -----------------------------------------
       ! Read gas datafield scale factor attribute
  
positive.code:     CHARACTER (LEN=15), PARAMETER :: cli_Psurf_field       = 'SurfacePressure'
      CHARACTER (LEN=18), PARAMETER :: cli_Temperature_field = 'TemperatureProfile'
  
      REAL      (KIND=r4) :: scale_lat, scale_lon, scale_gas, scale_Psurf, &
                             scale_temperature, scale_H2O
  
      ! ------------------------
      ! Error handling variables
  
negative.code:           ! ----------------------------------------------------
            ! Just selecting the closest location to lat lon pixel
            ! ----------------------------------------------------
            idx_lat = MINVAL(MINLOC(ABS(latvals(1:Cmlat) - REAL(lat(ixtrack,itimes),KIND=8) )))
            idx_lon = MINVAL(MINLOC(ABS(lonvals(1:Cmlon) - REAL(lon(ixtrack,itimes),KIND=8) )))
  
            ! ----------------------------------------------
            ! Convert pixel terrain height to pressure using
            ! USSA (good below 50 km)
            ! ----------------------------------------------
  !!$          CALL ussa_press(REAL(terrain_height(ixtrack,itimes),KIND=r8) / 1000.0_r8, &
  !!$               local_psurf(ixtrack,itimes) )
            
            ! ---------------------------------------------------------
            ! For OMH2O I should use terrain heigh from the climatology
            ! ---------------------------------------------------------
            local_psurf(ixtrack,itimes) = linInterpol(Cmlon, Cmlat, &
                 lonvals(1:Cmlat), latvals(1:Cmlon), &
                 Psurface(1:Cmlon,1:Cmlat), &
                 REAL(lon(ixtrack,itimes),KIND=8), REAL(lat(ixtrack,itimes),KIND=8), &
                 status=status)
            clima_psurf = Psurface(idx_lon, idx_lat)
            
            ! 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
  