＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
＊　EGS4 Newsletter in Japan          　　　　　　　　　　９５ー３号　＊
＊　　　　　　　　　　　　　　　　　　　　　　１９９５年　９月　７日　＊
＊　　　　　　　　　　　　　　　　　　　　　高エネルギー物理学研究所  ＊
＊　　　　　　　　　　　　　　　　　　　　　　放射線安全管理センター　＊
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
＜第５回ＥＧＳ４研究会に９３名参加＞
第５回ＥＧＳ４研究会は、９３名が参加し無事終了しました。講習に６２名の参加
があり、特に UNIX 機の場合はワークステーション１台で対応したため応答が悪く
なり、申し訳ありませんでした。それでも、一応皆さんにインストールと実行の経
験をしていただけたのではないかと思っています。ぜひ、各自のコンピュータにイ
ンストールをして、自分の計算に利用できる環境を作って頂きたいと思います。
問題点が出てきましたら、遠慮無く問い合わせて下さい。

＜トラブル例＞
 質問：下記の入力データで pegs4 を実行するとエラーがでる。
MIXT                                                                           
 &INP NE=7,RHO=2.26,RHOZ=49.83,1.71,4.56,31.58,1.92,8.26,1.22 &END        
CONCRETE                                                                       
O  NA AL SI K  CA FE                                                          
ENER                                                                           
 &INP AE=1.5,UE=100000.0,AP=.10,UP=100000.0 &END                              
TEST                                                                           
 &INP &END                                                                     
PWLF                                                                           
 &INP &END                                                                     
DECK                                                                           
 &INP &END                                                                     

 回答：ＰＣ−９８で Lahey Fortran を使用した場合には、理由はわかりませんが、
 　　　エラーとなるようです。ＤＯＳーＶではエラーにならない様です。
 　　　最初の行を、
 &INP NE=7,RHOZ=49.83,1.71,4.56,31.58,1.92,8.26,1.22,         
      RHO=2.26 &END 
      とすると、エラーが出なくなります。
      
 質問：EGS4の起動でエラーが起きたので、報告します。
　　　 我々の研究室では、alphaを使っています。
        OSF1 ver2.1では先日いただいた講習テキストどうりに動いたのですが、
        OSF1 ver3.0では、egs4runを起動すると以下のメッセージが出てエラーと
        なります。

         mortran3 procedure entered
         
         f77 procedure entered
 
         ld:
         Warning: random_ defined as GLOBAL COMMON but is defined in a shared 
         lib as a GLOBAL FUNC
         a.out procedure entered

         time: command terminated abnormally.

 回答：バグの様です。ＤＥＣで検討しているように聞いています。
      
＜SLAC EGS4 Listserver からの情報＞
　Liseserver を介して流されて来た情報で有益と思われたものをお知らせします。

＜質問＞

Hi EGS users,

I am involved in a project which is studying energy deposition and e-h
production in semiconductor memory devices.We wish to get a feel for the
ammount of ionisation and e-h production caused by incident electrons
and positrons of between 0.5-2.0 MeV energy.
Some of the slab thicknesses are of the order of 1 micron but the elemental
compositions are relatively simple.

My question is will EGS4 give a reasonably accurate figure for the ammount
of energy deposition in slabs this thin at this energy range ??

Thanks,

 Shaun Harvey

＜Alex の回答＞

Dear Shaun Harvey:

EGS4 ought to be OK for average energy deposition but if you're after
the differential energy deposited you may have to run with very low
values of AE. Even then the spectrum may be approximate either because
you can not push the cutoff's low enough or because of the fact that
the Moller cross sections assume that the electrons are unbound.
Impact ionisation is not modelled by EGS (yet!).

Here's how to get the most reliable results.

Run using standard EGS (no PRESTA) pushing AE with PEGS as low as you
can go.  Look at the energy deposition spectrum. If it converges to
some shape with lowering AE, then that's the best you can do with EGS.

The multiple scattering theory of EGS breaks down in this regime but it
ought not to matter much as the electron detour factors through foils
of this small thickness ought to be negligible.

Best of luck! Please report your findings to the listserv.  I'm sure
many other are interested in your results.

_______________________________________________
| Alex F Bielajew (aka BLIF)                   |
| Ionising Radiation Standards                 |
| Institute for National Measurement Standards |
| National Research Council of Canada          |
| Ottawa, CANADA                               |
| K1A 0R6                                      |
| tel: 613 993 2197                            |
| fax: 613 952 9865                            |
| e-mail: alex@irs.phy.nrc.ca                  |
|______________________________________________|


Mike Lay at Oxford wrote to me:

> 1) In UPHI, the code assumes that the contents of variables A, B, C
> are around from a previous call. This is not strictly true (as my
> compiler has pointed out!), and is certainly not true on machines
> with dynamic memory allocation. A simple SAVE A, B, C command would
> seem to fix this.....

> 2) ...except when you use PRESTA. Which also uses the variable B.
> Wouldn't be a problem except that the variable B is included in UPHI
> as part of the common block USERPR. This one is not so obvious to get
> around. Any suggestions?

The answer may be of interest to the EGS-afficionados out there so:

This was brought to my attention a long time ago.  Only those who delve
into the depth's of EGS ever discover this!

Yes, UPHI depends on A,B,C being static.  I recommend that EGS4 be
compiled with all local variables held static because it was coded in
the days when this was the default for all computers (although never
part of the ANSI Fortran standard).

SAVE A,B,C

certainly works but then there persists the conflict with the B of
Moliere theory that USERPR puts in a common. I've satisfied my self
that this conflict is never a problem and so I have not supplied a
patch, although the next version of PRESTA renames A,B,C to auphi,
buphi, cuphi to avoid the conflict. Try renaming them yourself and try
a repeat of a calculation. I'll bet you a pint there's no difference.

_______________________________________________
| Alex F Bielajew (aka BLIF)                   |
| Ionising Radiation Standards                 |
| Institute for National Measurement Standards |
| National Research Council of Canada          |
| Ottawa, CANADA                               |
| K1A 0R6                                      |
| tel: 613 993 2197                            |
| fax: 613 952 9865                            |
| e-mail: alex@irs.phy.nrc.ca                  |
|______________________________________________|

★WWW の EGS4 homepage は、
  http://ehssun.lbl.gov/egs/egs.html
  です。WWW を利用できる方は、アクセスしてみて下さい。

☆ Microsoft Fortran の最新版により、ＰＣで EGS4 System が使える事が
　判りましたのでお知らせします。

           Microsoft FORTRAN Powerstation によるEGS4実行

　　　　　　　　　　　　　　　　　　　　　　　7 SEP 1995/波戸　芳仁

1. Mortran3 インストール
   ・mortran3.dat作成、check77.mor: 問題なし
   ・FL32 /W0 mortran3.for でコンパイルする。

2. egs4run.batによるucsampl4.mor 実行
   ・CALL EXIT(0);がない。 -> STOP; に変更。(nrcc4mac.mor)
   ・TIMER がない。 -> コメントアウト

3. PEGS4コンパイル
   ・DIMENSION STDAT1(6,20)---　とその後の3行(Equivalent文)の順番が正
　　しくない。  -> これらを、DATA文の前に移動。
   ・SUBROUTINE HPLT1の引数の型が一致しない。
     -> SUBROUTINE HPLT1内で、INTEGER ICAP(12);を CHARCTER ICAP(12);
　　　　に変更。

4. PEGS4実行
   ・「NAMELIST入力において、宣言していない変数を入力している。」
     -> &END を変数と誤解していたので、/に変更。

5. UCNAI.MOR
   ・CALL EXIT(0);がない。 -> STOP; に変更。(nrc4macp.mor)

6. XYZDOS.MOR
   ・「statement out of order」-> EQUVALENT 文とDATA文の順序入れ替え。
　　　(１箇所)
   ・ SUBROUTINE PRINTER の引数が一致しない。-> SUBROUTINE PRINTER;を
      SUBROUTINE PRINTER(I,J,K,L);と変更。

7. EXAMIN MOR
   ・ SUBROUTINE AUSGAB の引数が一致しない。-> SUBROUTINE AUSGAB;を
      SUBROUTINE AUSGAB(I);と変更。

以上の変更により、全チェック項目について妥当な結果が得られた。
ただし、オプティマイズ(/Ox)は、ucsampl4.mor の最終乱数を変化させたので、
原則として使用しないことにする。（この場合、速度は、Lahey Fortran使用時の
約2/3になる。）

MS-Fortran Powerstationは58,000円と、Lahey Fortran(198,000円）に比べて大
幅に安い。Watcom FortranもProsideから、7万円程度で購入できるので、EGS4に
問題なく使用できるか、今後2カ月以内に調査する。