anchor.diff:           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, & 
  -         alt_lev_did, alt_lay_did, pre_lev_did, pre_lay_did, &
  +         alt_lev_did, alt_lay_did, pre_lay_did, &
            dspace, toms_datatype_id, pre_lev_cld_did, pre_lay_cld_did
   
       INTEGER(SIZE_T)                :: size
  
positive.code:        ALLOCATE (lut_pre_lay(1:ansrf,1:anlay),   STAT=estat)
         ALLOCATE (lut_pre_lev_cld(1:ancld,1:anlev),   STAT=estat)
         ALLOCATE (lut_pre_lay_cld(1:ancld,1:anlay),   STAT=estat)
         ALLOCATE (lut_pre_lev(1:ansrf,1:anlev),   STAT=estat)
  
         ALLOCATE (lut_I0_clr(1:anozo,1:ansrf,1:anvza,1:ansza), STAT=estat)
         ALLOCATE (lut_I1_clr(1:anozo,1:ansrf,1:anvza,1:ansza), STAT=estat)
  
negative.code:     
      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_lev_cld_did, H5T_NATIVE_REAL, lut_pre_lev_cld(1:clp_dim(1),1:alt_lev_dim(2)), &
  