anchor.diff:      CALL h5dget_space_f(alt_lay_did,dspace,hdferr)
       CALL h5sget_simple_extent_dims_f (dspace, alt_lay_dim, alt_lay_maxdim, hdferr)
   
  -    CALL h5dget_space_f(pre_lev_cld_did,dspace,hdferr)
  -    CALL h5sget_simple_extent_dims_f (dspace, alt_lev_cld_dim, alt_lev_cld_maxdim, hdferr)
       CALL h5dget_space_f(pre_lay_cld_did,dspace,hdferr)
       CALL h5sget_simple_extent_dims_f (dspace, alt_lay_cld_dim, alt_lay_cld_maxdim, hdferr)
       
  
positive.code:     CALL h5dopen_f(input_file_id,'/Profiles/Altitude Level', alt_lev_did, hdferr)
      CALL h5dopen_f(input_file_id,'/Profiles/Altitude Layer', alt_lay_did, hdferr)
      CALL h5dopen_f(input_file_id,'/Profiles/Pressure Layer', pre_lay_did, hdferr)
      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,'/Scattering Weights/Clear Sky/dI0', dI0_clr_did, hdferr)
  
negative.code:     CALL h5dclose_f(dI1_clr_did, hdferr); CALL h5dclose_f(dI1_cld_did, hdferr)
      CALL h5dclose_f(dI2_clr_did, hdferr); CALL h5dclose_f(dI2_cld_did, hdferr)
      
      CALL h5dclose_f(alt_lay_did, hdferr); CALL h5dclose_f(pre_lay_did, hdferr)
      CALL h5dclose_f(alt_lev_did, hdferr)
      CALL h5dclose_f(pre_lay_cld_did, hdferr)
      
      ! ----------
      ! Close file
  