anchor.diff:      he5stat = HE5_SWrdlattr ( swath_id, TRIM(ADJUSTL(gasdatafieldname)),&
                 "ScaleFactor", scale_gas       )
   
  +    ! -----------------------------------
  +    ! Read gas profile standard deviation
  +    ! -----------------------------------
  +    he5stat = HE5_SWrdfld (                                &
  +         swath_id, TRIM(ADJUSTL(gasdatafieldname_sd)),        &
  +         he5_start_3d, he5_stride_3d, he5_edge_3d,         &
  +         dummy_gaspr_sd(1:Cmlon,1:Cmlat,1:CmETA) )
  +
  +    ! ------------------------------------------------------------
  +    ! Read gas standard deviation datafield scale factor attribute
  +    ! ------------------------------------------------------------
  +    he5stat = HE5_SWrdlattr ( swath_id, TRIM(ADJUSTL(gasdatafieldname_sd)),&
  +              "ScaleFactor", scale_gas_sd    )
  +
       ! -------------------------------------------------------------------
       ! Finding out the data field for water vapor (to compute air density)
       ! -------------------------------------------------------------------
  
positive.code:   ! ---------------------------------------
    REAL(KIND=r8), DIMENSION(:),     ALLOCATABLE :: latvals, lonvals
    REAL(KIND=r8), DIMENSION(:,:),   ALLOCATABLE :: Psurface
    REAL(KIND=r8), DIMENSION(:,:,:), ALLOCATABLE :: Temperature, Gas_profiles, H2O_profiles
  
    ! ---------------------------------------
    ! Data obtained from Vlidort lookup table
  
negative.code: 
      REAL (KIND=r4), ALLOCATABLE, DIMENSION(:) :: dummy_lon, dummy_lat
      REAL (KIND=r4), ALLOCATABLE, DIMENSION(:,:) :: dummy_psurf
      REAL (KIND=r4), ALLOCATABLE, DIMENSION(:,:,:) :: dummy_temp, dummy_gaspr, dummy_h2opr, dummy_gaspr_sd
  
      REAL (KIND=r4) :: scale_lat, scale_lon, scale_gas, scale_Psurf, &
           scale_temperature, scale_H2O, scale_gas_sd
  
      ! ------------------------
      ! Error handling variables
  