Using a normalizing flow model as a likelihood

I have some multivariate data which doesn’t match well with any of the usual parametric distributions because of nonlinear correlations between variables. I’m thinking of trying to fit a normalizing flow model to the raw data samples and then using this as a likelihood for inferring the parameters of a latent variable model. I see that pymc3 has some NF capabilities but is it possible to use a fitted NF model as a likelihood in another model?