anchor.diff:         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
  +       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
  
positive.code:        ALLOCATE (lut_wavelength(1:anwav), STAT=estat)
         ALLOCATE (lut_toms(1:anozo),       STAT=estat)
         
         ALLOCATE (lut_alt_lay(1:ansrf,1:anlay),   STAT=estat)
         ALLOCATE (lut_pre_lay(1:ansrf,1:anlay),   STAT=estat)
         ALLOCATE (lut_pre_lay_cld(1:ancld,1:anlay),   STAT=estat)
  
  
negative.code:     CALL h5dget_space_f(dI0_cld_did,dspace,hdferr)
      CALL h5sget_simple_extent_dims_f (dspace, dI0_cld_dim, dI0_cld_maxdim, hdferr)
  
      CALL h5dget_space_f(pre_lay_did,dspace,hdferr)
      CALL h5sget_simple_extent_dims_f (dspace, alt_lay_dim, alt_lay_maxdim, hdferr)
      
      ! ---------------------------------------------------------------
      ! Allocate & initialize variables now that we have the dimensions
  