Patching Enzo with react_primordialZ

[christopher.jessop]

Hello,

Firstly, apologies if this is a topic that should be posted on the enzo users group, I wasn't certain which group would be best to post in.

I've been trying to patch enzo with krome using the react_primordialZ network after I last enquired about recreating the Omukai 05 plots with this network in the krome collapse test, which I managed to do successfully. 

I created the patch with: ./krome -enzo -compact -cooling ATOMIC,H2,Z -n networks/react_primordialZ -useTabs -noSinkCheck

Following the steps for adding non-primordial species to enzo, I am now trying to do the same recreation for the enzo onezone freefall test, but am not sure how to set the fields so that the initial abundances are the same as Omuakai 05 for oxygen and carbon. In krome, I set CI and OI as:

x(KROME_idx_C) = 0.927d-4 * 1d1**zs(jz)

and 

x(KROME_idx_O) = 3.568d-4 * 1d1**zs(jz) 

In the enzo Grid_OneZoneFreefallTestInitializeGrid.C I try initialise them as:

BaryonField[CINum][i] = TestProblemData.CI_Fraction * BaryonField[0][i], and the same for OINum, where I select a value for CI and OI fraction that will give the same value as the krome abundance. I don't think this is correct at all, as the species are initialised as mass densities in enzo, but initialised in krome as n(i)/n_H. Is there a way to translate these krome initial abundances into the enzo freefall collapse test? 

There is very odd behaviour when I plot any non-primordial densities when I run the test and I hope properly setting the fields will partly help with this, but I think there must be more issues that I have introduced by adding 18 non-primordial species.  Any help would be appreciated.

Cheers,

Chris