anchor.diff:      
       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_lev_did, H5T_NATIVE_REAL, lut_pre_lev(1:srf_dim(1),1:alt_lev_dim(2)), &
  -         alt_lev_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)), &
  
positive.code:        IF ( ALLOCATED ( lut_toms ) ) DEALLOCATE ( lut_toms )
         
         IF ( ALLOCATED ( lut_alt_lay ) ) DEALLOCATE ( lut_alt_lay )
         IF ( ALLOCATED ( lut_pre_lev ) ) DEALLOCATE ( lut_pre_lev )
         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 )
  
negative.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_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)
  