anchor.diff:           lut_dI2_cld(1:toz_dim(1),1:alt_lay_dim(2),1:clp_dim(1),1:vza_dim(1),1:sza_dim(1)), &
            dI0_cld_dim, hdferr)
       
  -    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_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_lay_cld_did, H5T_NATIVE_REAL, lut_pre_lay_cld(1:clp_dim(1),1:alt_lay_dim(2)), &
  -         alt_lay_cld_dim, hdferr)
  +         alt_lay_dim, hdferr)
       
       ! --------------
       ! Close datasets
  
positive.code:        wav_dim, wav_maxdim
    INTEGER(HSIZE_T), DIMENSION(2) :: Sb_clr_dim, Sb_clr_maxdim, &
         alt_lay_dim, alt_lay_maxdim, &
         alt_lay_cld_dim, alt_lay_cld_maxdim, Sb_cld_dim, Sb_cld_maxdim
    INTEGER(HSIZE_T), DIMENSION(4) :: I0_clr_dim,  I0_clr_maxdim, I0_cld_dim, I0_cld_maxdim
    INTEGER(HSIZE_T), DIMENSION(5) :: dI0_cld_dim, dI0_cld_maxdim
    INTEGER(HSIZE_T), DIMENSION(6) :: dI0_clr_dim,  dI0_clr_maxdim
  
negative.code:          I0_cld_did, I1_cld_did, I2_cld_did, Ir_cld_did, Sb_cld_did,        &
           dI0_clr_did, dI1_clr_did, dI2_clr_did, & 
           dI0_cld_did, dI1_cld_did, dI2_cld_did, & 
           pre_lay_did, &
           dspace, toms_datatype_id, pre_lay_cld_did
  
      INTEGER(SIZE_T)                :: size
  