Skip to content

RMG-generated network__.py doesn't run in measure.py #40

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
wants to merge 6 commits into from

Conversation

rwest
Copy link
Member

@rwest rwest commented Jun 16, 2011

When I run RMG-py it generates a pdep/network73_1.py file, which contains the lines

isomer("HCCOH")
reactants("C2H", "OH")

The species labels are missing their numbers.
In order to run in standalone measure.py without giving a KeyError on the line

461         network.isomers = [speciesDict[label] for label in network.isomers]

it should read:

isomer("HCCOH(31)")
reactants("C2H(3)", "OH(32)")

rwest added 5 commits June 16, 2011 15:59
The commit e1f226c changed the method for writing
a measure input file, but not the call to it from rmg/model.py
When reading in from the smiles '[CH]' the carbon ends up with a radicalElectrons
count of 0 instead of 3. Strangely, [CH2] and [CH3] work fine.
This causes problems in data/thermo.py when trying to find the thermo for [CH].
…oner.

As soon as the input file has been read in, we write an output file 
which has 0 species and 0 reactions. The species pictures etc. from the 
previous run have been wiped by this point already, but it can be confusing
to have the previous output file still there when you're running a new job.
Whenever the latest chemXXXX.inp file is written, it is copied (or rather hard-linked)
to chemkin/chem.inp, so that you can always just look at this file and get the latest.
It used to write
    isomer("HCCOH")
    reactants("C2H", "OH")
instead of 
    isomer("HCCOH(31)")
    reactants("C2H(3)", "OH(32)")

This closes #40
@rwest
Copy link
Member Author

rwest commented Jun 16, 2011

Hmm. Only the last of these commits actually deals with this issue, so perhaps I shouldn't have turned it into a pull request, but in any event merging it will close the issue.

This is close to something I've been using in RMG-Java.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant