Determining mol/L at STP from volume and mol fractions in a CFD program

I'm using some CFD software and the outputs are a bit unusual so I thought I would double check my understanding on how to do these conversions. 


The data available per cell is: cell volume, particle volume fraction, average fluid density, mol fraction of gas species, temperature, pressure


I am trying to come up with the mol/L at STP as a way of normalizing for density effects to compare concentration of a gas species throughout the height of a reactor. Temperature and pressure vary wildly, partially because this is a fluidized bed and solids distribution plays a big role in pressure. 


So far my process has been convert from mol fraction to mass fraction. 


And then multiply (cell volume-cell volume\*particle volume fraction) by average fluid density to get gaseous mass in the cell. Then I can use the mass fraction to determine the amount of mass for a given species. Obviously at this point I have the mols of each species. 


The problem is then I need the density of that specific species at a given temperature and pressure (currently working at temperatures around 950 C and 110kPa but they very cell by cell obviously) to get back to volume. Is the only way to do this to estimate the density? 


The other problem is that let's say I do get the volume by estimating the density. Then I convert to volume at STP by multiplying by 273 K/T and P/101325 Pa or whatever it is. Are the mol fractions I had at non-STP still usable? Or do I just take these volumes and divide by 22.4 and that's the mols I should be using. 


Can someone please straighten out my thinking on this. Also are there any good sources for density of a gas species as a function of pressure and temperature. Everything I look at online seems really bad. 


