anchor.diff:         ! ----------------------------
          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)),INT(alt_lev_dim(2)),errstat)
  +            INT(alt_lay_dim(2)),errstat)
   
          ! -----------------------------------------------------------------
          ! Work out the AMF using the scattering weights and the climatology
  
positive.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(alt_lev_did,dspace,hdferr)
      CALL h5sget_simple_extent_dims_f (dspace, alt_lev_dim, alt_lev_maxdim, hdferr)
      CALL h5dget_space_f(alt_lay_did,dspace,hdferr)
      CALL h5sget_simple_extent_dims_f (dspace, alt_lay_dim, alt_lay_maxdim, hdferr)
  
  
negative.code:     ! Input variables
      ! ---------------
      CHARACTER (LEN=1), INTENT (IN) :: ad
      INTEGER (KIND=i4), INTENT (IN) :: anozo, ansrf, ancld, analb, ansza, anvza, anwav, anlay
  
      ! ------------------
      ! Modified variables
  