anchor.diff:                  ! Determine phenomenological rate coefficients
                   call network_applyApproximateMethod(net, nIsom, nReac, nProd, Elist, nGrains, T, P, method, K(u,v,:,:))
   
  +            end do
  +
  +            do i = 1, nIsom
  +                deallocate(net%isomers(i)%densStates)
  +                deallocate(net%isomers(i)%eqDist)
  +            end do
  +            do w = 1, size(net%reactions)
  +                deallocate(net%reactions(w)%kf)
  +                deallocate(net%reactions(w)%kb)
               end do
   
           end do
  
positive.code:         ! Return the chosen energy grains
          call network_getEnergyGrains(Emin, Emax, grainSize, nGrains, Elist)
  
  
      end subroutine
  
  
negative.code: 
              ! Calculate microcanonical rate coefficients at the current conditions
              do w = 1, size(net%reactions)
                  allocate(net%reactions(w)%kf(1:nGrains))
                  allocate(net%reactions(w)%kb(1:nGrains))
                  call reaction_calculateMicrocanonicalRates(net%reactions(w), T, Elist, nGrains, net%isomers, net%species, &
                      nIsom, nReac, nProd)
              end do
  