anchor.diff:         ALLOCATE (lut_pre_lay_cld(1:ancld,1:anlay),   STAT=estat)
          ALLOCATE (lut_pre_lev(1:ansrf,1:anlev),   STAT=estat)
   
  -       ALLOCATE (lut_tem(1:anozo,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)
          ALLOCATE (lut_I2_clr(1:anozo,1:ansrf,1:anvza,1:ansza), STAT=estat)
  
positive.code:     CALL h5dopen_f(input_file_id,'/Profiles/Cloud pressure Level', pre_lev_cld_did, hdferr)
      CALL h5dopen_f(input_file_id,'/Profiles/Cloud pressure Layer', pre_lay_cld_did, hdferr)
      CALL h5dopen_f(input_file_id,'/Profiles/Ozone Column Layer', ozo_did, hdferr)
      CALL h5dopen_f(input_file_id,'/Profiles/Temperature Level', tem_did, hdferr)
  
      CALL h5dopen_f(input_file_id,'/Scattering Weights/Clear Sky/dI0', dI0_clr_did, hdferr)
      CALL h5dopen_f(input_file_id,'/Scattering Weights/Clear Sky/dI1', dI1_clr_did, hdferr)
  
negative.code:          alt_lay_cld_dim, hdferr)
      CALL h5dread_f(ozo_did,     H5T_NATIVE_REAL, lut_ozo(1:toz_dim(1),1:srf_dim(1),1:alt_lay_dim(2)), &
           air_dim, hdferr)
      
      ! --------------
      ! Close datasets
  