anchor.diff:  
       CALL h5dget_space_f(air_did,dspace,hdferr)
       CALL h5sget_simple_extent_dims_f (dspace, air_dim, air_maxdim, hdferr)
  -    CALL h5dget_space_f(tem_did,dspace,hdferr)
  -    CALL h5sget_simple_extent_dims_f (dspace, tem_dim, tem_maxdim, hdferr)
       
       ! ---------------------------------------------------------------
       ! Allocate & initialize variables now that we have the dimensions
  
positive.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)
      CALL h5dread_f(tem_did,     H5T_NATIVE_REAL, lut_tem(1:toz_dim(1),1:srf_dim(1),1:alt_lev_dim(2)), &
           tem_dim, hdferr)
      
      ! --------------
      ! Close datasets
  
negative.code:        IF ( ALLOCATED ( lut_pre_lay ) ) DEALLOCATE ( lut_pre_lay )
         IF ( ALLOCATED ( lut_pre_lay_cld ) ) DEALLOCATE ( lut_pre_lay_cld )
         IF ( ALLOCATED ( lut_pre_lev_cld ) ) DEALLOCATE ( lut_pre_lev_cld )
         
         IF ( ALLOCATED ( lut_I0_clr ) ) DEALLOCATE ( lut_I0_clr )
         IF ( ALLOCATED ( lut_I1_clr ) ) DEALLOCATE ( lut_I1_clr )
  