Skip to content
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

NLTE fixes - discussion needed before merge - relates to #109 #112

Merged
merged 18 commits into from
Apr 16, 2014

Conversation

ssim
Copy link
Contributor

@ssim ssim commented Mar 26, 2014

summary:

  • bug fix in plasma_array, wrong indentation on line 587
    (calculate_nlte_level_populations)
  • simplification of nlte rates - treatment of stimulated emission no
    longer handled as negative absorption (correction factor removed
    from calculate_nlte_level_populations / reformulated terms)

LTE check made using this version He I and He II, but reproducible test case needed - this needs to be discussed/explained how to insert set up.

@wkerzendorf look at the change proposed in the formulation of the rate equations. Does this look problematic? This relates to the suggestion in #109

  • suggestion is to replace use of the stimulated correction factor in the rate equations (that is treating stimulated emission as negative recombination) with direct use of stimulated emission as term. Seems simpler to follow and behaves better in the LTE tests
  • we need to decide if this is ok, test results compared to previous NLTE calculations (rerun one or more from paper) and put in place some mechanism to make sure we have a solution to these equations as one of our tests (I can propose what that test would be but need guidance on how to implement it)

Do not recommend merging this until it has been discussed fully.

@@ -92,7 +92,9 @@ class BasePlasmaArray(object):
"""

@classmethod
def from_abundance(cls, abundance_dict, density, atom_data, time_explosion, nlte_config=None, saha_treatment='lte'):
def from_abundance(cls, abundance_dict, density, atom_data, time_explosion,
nlte_config=None, ionization_mode='lte',
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

shouldn't this be 1 instead of None. This is like in #53

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

for i in range(20):
   print i

@ssim
Copy link
Contributor Author

ssim commented Apr 15, 2014

Tests are now in place. Bug fixes have been attempted. This version should now be considered possible for merging.

@ssim
Copy link
Contributor Author

ssim commented Apr 16, 2014

passed: Travis says "yes"

ssim added a commit that referenced this pull request Apr 16, 2014
NLTE fixes - discussion needed before merge - relates to #109
@ssim ssim merged commit 78f64f4 into tardis-sn:master Apr 16, 2014
@unoebauer unoebauer mentioned this pull request Nov 30, 2015
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants