anchor.diff:      
       CALL h5dread_f(alt_lay_did, H5T_NATIVE_REAL, lut_alt_lay(1:srf_dim(1),1:alt_lay_dim(2)), &
            alt_lay_dim, hdferr)
  -    CALL h5dread_f(pre_lev_did, H5T_NATIVE_REAL, lut_pre_lev(1:srf_dim(1),1:alt_lev_dim(2)), &
  -         alt_lev_dim, hdferr)
       CALL h5dread_f(pre_lay_did, H5T_NATIVE_REAL, lut_pre_lay(1:srf_dim(1),1:alt_lay_dim(2)), &
            alt_lay_dim, hdferr)
       CALL h5dread_f(pre_lev_cld_did, H5T_NATIVE_REAL, lut_pre_lev_cld(1:clp_dim(1),1:alt_lev_dim(2)), &
  
positive.code:   REAL(KIND=r4),    DIMENSION(:),           ALLOCATABLE :: lut_alb, lut_clp, lut_sza, &
         lut_vza, lut_wavelength, lut_srf
    CHARACTER(LEN=5), DIMENSION(:),           ALLOCATABLE :: lut_toms
    REAL(KIND=r4),    DIMENSION(:,:),         ALLOCATABLE :: lut_alt_lay, lut_pre_lay, lut_pre_lev, &
         lut_pre_lev_cld, lut_pre_lay_cld
    REAL(KIND=r4),    DIMENSION(:,:),         ALLOCATABLE :: lut_Sb_clr, lut_Sb_cld
    REAL(KIND=r4),    DIMENSION(:,:,:,:),     ALLOCATABLE :: lut_I0_clr, lut_I1_clr, lut_I2_clr, lut_Ir_clr, &
  
negative.code:     CALL h5dclose_f(dI2_clr_did, hdferr); CALL h5dclose_f(dI2_cld_did, hdferr)
      
      CALL h5dclose_f(alt_lay_did, hdferr); CALL h5dclose_f(alt_lev_did, hdferr)
      CALL h5dclose_f(pre_lay_cld_did, hdferr); CALL h5dclose_f(pre_lev_cld_did, hdferr)
      
      ! ----------
  