9rnsr / FlexID

内部被ばく線量評価コードFlexID (Flexible code for Internal Dosimetry)
MIT License
1 stars 0 forks source link

排泄コンパートメントの積算放射能として常にゼロが出力される #29

Closed 9rnsr closed 1 week ago

9rnsr commented 1 month ago

本件は、#26にてExcretion周りの計算結果の持ち方を整理した際に発見した。

従来からexcコンパートメントの計算では、他コンパートメントから流入する放射能の総計=積算放射能を求めたうえで、これを出力時間メッシュで割ることで平均放射能を計算し、Retentionファイルに出力していた。しかし求めた積算放射能をCumulativeファイルに出力される数値として正しい位置に格納しておらず、その結果として代わりに常にゼロを出力するという動作となっていた。

26では、出力について従来同等とするため、以下に示すようにexcコンパートメントについてのみ意図的にゼロ設定を行う処理を入れている。(示したのはOIR側のコードだがEIRも同様)

https://github.com/9rnsr/FlexID/blob/ff6601e022d74d076bd4b1e8443e627263fffe7a/FlexID.Calc/MainRoutine_OIR.cs#L239-L249

このコードを削除すれば、他コンパートメントと同様に摂取時からの、excコンパートメントでは体外排出された積算放射能が出力できるようになる。

9rnsr commented 1 month ago

要望がありそうなので @9rnsr にAssignした。