Skip to content

Commit

Permalink
Fix error in fluence for dosxyznrc sources 1 and 7
Browse files Browse the repository at this point in the history
Change minus sign to plus sign for particles which miss the geometry in
the fluence calculation. Previously, particles that missed the geometry
were subtracted from the total fluence (used to normalize the dose).
This bug has been there for years, indicating, hopefully, how rarely
these sources are used.
  • Loading branch information
blakewalters authored and ftessier committed Jan 24, 2017
1 parent c2cb50e commit ae360a8
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions HEN_HOUSE/user_codes/dosxyznrc/dosxyznrc.mortran
Original file line number Diff line number Diff line change
Expand Up @@ -2908,9 +2908,9 @@ ELSEIF((enflag=1)|(enflag=2)|(enflag=3)|(enflag=4))[ ;OUTPUT61 etot/esrc;

"Calculate incident fluence
IF(isource = 0 | isource = 1 | isource = 3 | isource = 7)[
temp2=dble(IHSTRY+ncaseold-nsmiss-nmissm);
IF(beamarea = 0)[ ainflu=dble(IHSTRY+ncaseold-nsmiss-nmissm);]
ELSE[ ainflu=dble(IHSTRY+ncaseold-nsmiss-nmissm)/beamarea;]
temp2=dble(IHSTRY+ncaseold+nsmiss+nmissm);
IF(beamarea = 0)[ ainflu=dble(IHSTRY+ncaseold+nsmiss+nmissm);]
ELSE[ ainflu=dble(IHSTRY+ncaseold+nsmiss+nmissm)/beamarea;]
]
ELSEIF(isource = 2 | isource = 8 )[
ainflu=dble(IHSTRY+ncaseold+nsmiss+
Expand Down

0 comments on commit ae360a8

Please sign in to comment.