anchor.diff:         vza_dim, vza_maxdim, &
          wav_dim, wav_maxdim
     INTEGER(HSIZE_T), DIMENSION(2) :: Sb_clr_dim, Sb_clr_maxdim, &
  -       alt_lev_dim, alt_lev_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
  
positive.code:     ! Input variables
      ! ---------------
      CHARACTER (LEN=1), INTENT (IN) :: ad
      INTEGER (KIND=i4), INTENT (IN) :: anozo, ansrf, ancld, analb, ansza, anvza, anwav, anlay, anlev
  
      ! ------------------
      ! Modified variables
  
negative.code:        ! ----------------------------
         CALL vlidort_allocate ("d", INT(toz_dim(1)),INT(srf_dim(1)),INT(clp_dim(1)),INT(alb_dim(1)), &
              INT(sza_dim(1)),INT(vza_dim(1)),INT(wav_dim(1)), &
              INT(alt_lay_dim(2)),errstat)
  
         ! -----------------------------------------------------------------
         ! Work out the AMF using the scattering weights and the climatology
  