＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
＊　EGS4 Newsletter in Japan          　　　　　　　　　　９４ー４号　＊
＊　　　　　　　　　　　　　　　　　　　　　　１９９５年　２月　２日　＊
＊　　　　　　　　　　　　　　　　　　　　　高エネルギー物理学研究所  ＊
＊　　　　　　　　　　　　　　　　　　　　　　放射線安全管理センター　＊
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊

＜Lahey Fortran 90J について＞
　ＰＣでＥＧＳ４を使用する場合に必要な Lahey Fortran 90J がようやく（株）
日本コンピュータ研究所〔〒220  横浜市西区平沼 1-2-14 横浜 NT ビル、
Tel 045-324-2898, Fax 045-322-2798〕から発売されました。
　PC-9821 Ap2 にインストールして使用してみました。ハードディスクをパーティ
ションして使用している場合、コンパイルするディレクトリーと Lahey Fortran 
をインストールするディレクトリーが異なったパーティション上にある場合には、
うまくリンクされないという問題はありますが、その他は特に問題はないようです。
問題を避けるためには、EGS4 を実行するディレクトリーと同じパーティション内の
ディレクトリーにシステムを組み込んだ方が良いと思います。
　コンパイルの最適化は、デフォルトの -o1 が一番良いようです。
　すでに、DOS-V 用の Lahey Fortran を使っておられる方は、新たに買い直す必要
は無いと思います。
　参考のために、egs4run.bat (PRESTA なし) と egs4runp.bat (PRESTA あり) を
紹介します。



REM 使用方法: EGS4RUN %1 %2
REM %1 : ユーザーコード　例：ucsampl4.mor
REM %2 : 断面積データ  　例：sampl4.dat

DEL MORTJOB.*

REM MORTJOB.MORの作成

COPY LISTING.ON+EGS4MAC.MOR+NRCC4MAC.MOR+PAIRMAC.MOR+KEK4MAC.MOR MORTJOB.TMP
COPY MORTJOB.TMP+%1+KEK4.MOR+EGS4BLOK.MOR+EGS4.MOR MORTJOB.MOR

REM 断面積ファイルのコピー
COPY %2 MORTJOB.XSE

REM MORTJOB.MORからMORTJOB.Fへの変換
MORTRAN3

REM MORTJOB.Fのコンパイル、リンク、実行
lf90  d:\egs4\egs4unix\MORTJOB.for -nw
mortjob



REM 使用方法: EGS4RUNP %1 %2
REM %1 : ユーザーコード　例：ucsampl4.mor
REM %2 : 断面積データ  　例：sampl4.datDEL MORTJOB.*

REM MORTJOB.MORの作成
COPY LISTING.ON+EGS4MAC.MOR+NRC4MACP.MOR+PRESTA.MOR+PAIRMAC.MOR MORTJOB.TMP
COPY MORTJOB.TMP+KEK4MAC.MOR+%1+KEK4.MOR+EGS4BLOK.MOR+EGS4.MOR MORTJOB.MOR

REM 断面積ファイルのコピー
COPY %2 MORTJOB.XSE

REM MORTJOB.MORからMORTJOB.Fへの変換
MORTRAN3

REM MORTJOB.Fのコンパイル、リンク、実行
lf90  d:\egs4\egs4unix\MORTJOB.for -nw
mortjob

＜egs4-l から送られてきた Rogers の MCNP に関する情報＞

Re Use of MCNP4A

MCNP is a very powerful Monte Carlo tool.  However, when electron transport
was patched onto it, an early version of the ITS system was used (V1).
Many years ago we showed that there is a significant error in this version
when electron depth-dose curves are being calculated (about 10% effects or
less).  Thus one should not use this code for careful work with
electron beams.

On the other hand, for brem production cases we have shown that ITS and EGS4
give nearly the same results - or at least at the level that can be
experimentally checked.  However the default versions of EGS4 or ITS can
not be used for accurate simulations.  This presumably applies to MCNP4A
as well.

For more complete discussions, see:

Title: Differences in Electron Depth Dose Curves Calculated with EGS and ETRAN
 and Improved Energy Range Relationships
Authors: D.W.O. Rogers and A.F. Bielajew
Source: Medical Physics , 13
Year: 1986
Pages: 687 -- 694

Title: Comparisons of thick-target bremsstrahlung calculations by EGS4 and ITS
Authors: B.A. Faddegon and D.W.O. Rogers
Source: Nuc. Inst. Meth. A 327
Year: 1993
Pages: 556 -- 565

For electron transport calculations, I also expect that well coded EGS4
calculations are much faster - although MCNP should be much faster for
photon calculations where some of their sophisticated variance reduction
techniques can be applied as part of the package.

Dave Rogers

IRS/NRC/Ottaw
dave@irs.phy.nrc.ca

＜pegs4 の入力データに関する質問と答え＞
〔質問〕
I can not make pegs's data about GAS.
I use this *.inp file.
----------------------------------------------------------------
MIXT
 &INP NE=3,RHO=0.0012,RHOZ=0.76,0.23,0.01,GASP=1.0,IRAYL=1 &END
AIR-GAS                       AIR-GAS
N O AR
ENER
 &INP AE=0.512,UE=500000.0,AP=.10,UP=500000.0 &END
TEST
 &INP &END
PWLF
 &INP &END
DECK
 &INP &END
-----------------------------------------------------------------
And pgs4job.pegs4lst said this.
-----------------------------------------------------------------
1                    PEGS4 LISTING FILE
                    (WITH NRCC MODIFICATIONS, Jan 13,1988)


 ************************************************************
 *                                                          *
 *  OPT = MIXT                                              *
 *                                                          *
 ************************************************************


 ------------------------------------------------------------
 MEDIUM=AIR-GAS                 ,STERNHEIMER ID=AIR-GAS       
 ------------------------------------------------------------

 ATOMIC SYMBOLS ARE:  N   A   
  A NOT AN ATOMIC SYMBOL FOR AN ELEMENT WITH Z LE 100
 BAD ATOMIC SYMBOL....JOB ABORTED.

＜答え＞
　pegs4 の入力データにおいて、原子記号は A2,1X のフォーマートで入力する
様になっています。従って、この例の場合には、
N  O  AR
とする必要があります。
N O AR
として指定すると、 A2,1X フォーマートでは、最初の原子記号は N で２番目の
原子記号は A と解釈されたために、上記の様なエラーメッセージが出され、プロ
グラムは終了となったのです。