Re: Giant .err files

From: Alberto Fasso' <fasso_at_mail.cern.ch>
Date: Mon, 15 Jul 2013 16:12:35 +0200

You should not turn off those diagnostic statements, but should eliminate
their cause.
Their cause is the fact that the direction cosines read by the source.f
routine are not properly normalized (in double precision!)
If the direction cosines are not such that the sum of their squares = 1
(in double precision), the program cannot track the particle correctly
and needs to adjust them, wasting valuable computer time.

Supposing you read the cosines as UUU, VVV, WWW, do the following:
        ROOT = SQRT(UUU**2 + VVV**2 + WWW**2)
        UUU = UUU / ROOT
        VVV = VVV / ROOT
        WWW = WWW / ROOT

Alberto

On Sun, 14 Jul 2013, Joe Mildenberger wrote:

> When I run the second step of a two-step model and read in starting
> values from a source file, I get huge .err files.
>
> I hasten to add that there are apparently no real errors that I can find
> listed in the files, just thousands and thousands of entries of the form
>
> Jtrack, Etrack 8 0.939590109600000
> GEOFAR, TXYZ: 0.999999999694202
> 0.794822116000000 -0.180331337000000 -0.579429385000000
> Nfrom, Nreg, X, Y, Z 1001 51
> -864.880746000000 47.6401516000000 3048.00000000000
> Jtrack, Etrack 8 0.939565630105940
> GEOFAR, TXYZ: 0.999999999775368
> -0.119170720000000 -5.087247140000000E-03-0.992860745000000
> Nfrom, Nreg, X, Y, Z 1001 51
> -880.159423000000 66.2758347000000 3048.00000000000
> Jtrack, Etrack 8 0.939963188000000
> GEOFAR, TXYZ: 1.00000000021697
> 0.654909806000000 -0.106889125000000 -0.748109525000000
> Nfrom, Nreg, X, Y, Z 1001 51
> -900.314958000000 33.6777587000000 3048.00000000000
> Jtrack, Etrack 8 0.939963188000000
>
> These diagnostic statements are not coming from my source.f file, and I
> would like to know how to turn them off, because the amount of output is
> enormous, about 0.5 GB for a run of 2e6 primaries.
>
> Thanks for any help anyone can offer!
>
> Joe
Received on Mon Jul 15 2013 - 18:38:58 CEST