Re-implement Photon Source Normalization Factor in source.F90

svalinn / r2s-act

Rigorous 2 Step Activation Workflow

4 stars 7 forks source link

Re-implement Photon Source Normalization Factor in source.F90 #43

Closed erelson closed 11 years ago

erelson commented 11 years ago

Currently, the write_gammas.py script calculates and applies a normalization factor to the source information from phtn_src (which is first tagged to the mesh, e.g. for visualization). With the implementation of reading the MOAB mesh directly into source.F90, this code needs to be adapted from Python to Fortran, and then verified.

erelson commented 11 years ago

This was overall simple to do. Testing verifies that the same results are seen in the problem summary table for:

Uniform sampling (source_gamma.F90)
Structured voxel sampling (source_gamma.F90)
Hexahedra sampling with source_moab.F90 (on structured hexahedra mesh)
Tet sampling with source_moab.F90 (Starting from unstructured neutron mesh tally of problem)

(The above tests were done with the |concrete|h20|concrete| slab problem).

Simplicity and testing make me confident this was re-implemented correctly, so closing the issue.