RE: USRYIELD binary output

From: Lee, Kerry T. \(JSC-SF\)[LMIT] (kerry.t.lee@nasa.gov)
Date: Mon Dec 11 2006 - 16:00:07 CET

Next message: me@marychin.org: "Low energy neutron xsec not found for some media"

Previous message: me@marychin.org: "Re: miraculous zero photoneutron production?"
In reply to: Alfredo Ferrari: "RE: USRYIELD binary output"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This is a multi-part message in MIME format.

------_=_NextPart_001_01C71D35.0BE304A4
Content-Type: text/plain;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

Hi Alfredo,

I have determined that the algorithm is doing things right and even =
though the what1 of the input file does not follow the manual's =
instructions it turns out to give the right results anyway. The =
difference is that instead of making "ia" a positive number only, if one =
makes it negative the algorithm produces the right result anyway. The =
following code does the work right.

          AAUSYL =3D 0.01D+00 * DBLE (ITUSYL(NY))
          IAUSYL (NY) =3D NINT (AAUSYL)
          IEUSYL (NY) =3D ITUSYL (NY) - IAUSYL (NY) * 100

If you put in ITUSYL =3D -123, then IAUSYL =3D -1 and IEUSYL =3D -23, =
but for ITUSYL =3D 77 you get IAUSYL =3D 1 and IEUSYL =3D -23, but I =
assume since IA represents a single bin the sign is ignored and when =
what1 is packed back up before being written to the USRYIELD output the =
algorithm ignores the sign again by doing

ITUSYL(NY) =3D IEUSYL(NY) + ABS(IAUSYL(NY))*100

Although this is a little confusing on how I was getting the right =
results while not following the explicit instructions of the manual, and =
that the written what1 was different than the supplied what1, the fault =
was mine. Even so, it might be better to not ignore the sign of IA when =
writing to the USRYIELD output so that it will agree with what the user =
input.

Thanks
Kerry

-----Original Message-----
From: Alfredo Ferrari [mailto:alfredo.ferrari@cern.ch]
Sent: Mon 12/11/2006 7:35 AM
To: Lee, Kerry T. (JSC-SF)[LMIT]
Cc: fluka-discuss@fluka.org
Subject: RE: USRYIELD binary output
=20
Hi Kerry

the code tries to unpack from what(1) the infos about which two physical =

quantities you want to deal with. Probably a nonsense number is unpacked
into something strange (I should verify).

Ciao
Alfredo
On Fri, 8 Dec 2006, Lee, Kerry T. (JSC-SF)[LMIT] wrote:

> This is a multi-part message in MIME format.
>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0
> Content-Type: text/plain;
> charset=3D"iso-8859-1"
> Content-Transfer-Encoding: quoted-printable
>
> I seem to have found that I was in error when explaining the problem. =
I =3D
> missed the fact that what12 is stored as the next variable after =
what1, =3D
> but the value stored where what1 should be is 77 and the value from =
the =3D
> input file is -123. How is this not being transferred correctly to =
the =3D
> output file?
>
> Thanks
> Kerry
>
>
> -----Original Message-----
> From: Lee, Kerry T. (JSC-SF)[LMIT]
> Sent: Thu 12/7/2006 10:27 AM
> To: fluka-discuss@fluka.org
> Subject: USRYIELD binary output
> =3D20
>
> Dear FLUKA authors,
>
> I am reading a USRYIELD output that was written according to the =3D
> following USRYIELD command
>
> USRYIELD -123.00000 201.00000-50.000000 4.0000000 5.0000000 =3D
> 1.0000000LET+0.1 =3D20
> USRYIELD 1000.0000 0.0000000 100.00000 1.5000000 0.0000000 2603.0000 =
& =3D
> =3D20
> =3D20
> I am trying to parse the file using a C++ routine, so it is imperitive =
I =3D
> understand the details of what is written where. I followed the code =
in =3D
> the manual on reading the USRYIELD output =3D
> (http://www.fluka.org/manual/sect/s090/text.html) and see that the =3D
> variable ITUSYL(NY) is supposed to correspond to the What1. When I =
read =3D
> this file the variable reads a value of 2603, which is my What12 and =
not =3D
> What1. I have taken into account the extra bytes that fortran writes =
=3D
> and all the other values are read in correctly, so is this a bug or am =
I =3D
> missing something?
>
> Thanks
> Kerry
>
>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0
> Content-Type: text/html;
> charset=3D"iso-8859-1"
> Content-Transfer-Encoding: quoted-printable
>
> <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
> <HTML>
> <HEAD>
> <META HTTP-EQUIV=3D3D"Content-Type" CONTENT=3D3D"text/html; =3D
> charset=3D3Diso-8859-1">
> <META NAME=3D3D"Generator" CONTENT=3D3D"MS Exchange Server version =3D
> 6.5.7651.34">
> <TITLE>RE: USRYIELD binary output</TITLE>
> </HEAD>
> <BODY>
> 
>
> I seem to have found that I was in error when =3D
> explaining the problem.  I missed the fact that what12 is stored =
as =3D
> the next variable after what1, but the value stored where what1 should =
=3D
> be is 77 and the value from the input file is -123.  How is this =
=3D
> not being transferred correctly to the output file? 
> 
> Thanks 
> Kerry 
> 
> 
> -----Original Message----- 
> From: Lee, Kerry T. (JSC-SF)[LMIT] 
> Sent: Thu 12/7/2006 10:27 AM 
> To: fluka-discuss@fluka.org 
> Subject: USRYIELD binary output 
> 
> 
> Dear FLUKA authors, 
> 
> I am reading a USRYIELD output that was written according to the =3D
> following USRYIELD command 
> 
> USRYIELD  -123.00000 201.00000-50.000000 4.0000000 5.0000000 =3D
> 1.0000000LET+0.1   
> USRYIELD   1000.0000 0.0000000 100.00000 1.5000000 0.0000000 =
=3D
> 2603.0000 &        
> 
> I am trying to parse the file using a C++ routine, so it is imperitive =
I =3D
> understand the details of what is written where.  I followed the =
=3D
> code in the manual on reading the USRYIELD output (<A =3D
> =
HREF=3D3D"http://www.fluka.org/manual/sect/s090/text.html">http://www.flu=
ka=3D
> .org/manual/sect/s090/text.html</A>) and see that the variable =3D
> ITUSYL(NY) is supposed to correspond to the What1.  When I read =
=3D
> this file the variable reads a value of 2603, which is my What12 and =
not =3D
> What1.  I have taken into account the extra bytes that fortran =
=3D
> writes and all the other values are read in correctly, so is this a =
bug =3D
> or am I missing something? 
> 
> Thanks 
> Kerry 
> 
> 
> 
>
> </BODY>
> </HTML>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0--
>

--=20

+------------------------------------------------------------------------=
----+
| Alfredo Ferrari || Tel.: +41.22.767.6119 =
    |
| CERN-AB || Fax.: +41.22.767.7555 =
    |
| 1211 Geneva 23 || e-mail: Alfredo.Ferrari@cern.ch =
    |
| Switzerland || Alfredo.Ferrari@mi.infn.it =
    |
+------------------------------------------------------------------------=
----+

------_=_NextPart_001_01C71D35.0BE304A4
Content-Type: text/html;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
charset=3Diso-8859-1">
<META NAME=3D"Generator" CONTENT=3D"MS Exchange Server version =
6.5.7651.34">
<TITLE>RE: USRYIELD binary output</TITLE>
</HEAD>
<BODY>

Hi Alfredo, 
 
I have determined that the algorithm is doing things right and even =
though the what1 of the input file does not follow the manual's =
instructions it turns out to give the right results anyway.  The =
difference is that instead of making "ia" a positive number =
only, if one makes it negative the algorithm produces the right result =
anyway.  The following code does the work right. 
 
          AAUSYL =3D =
0.01D+00 * DBLE (ITUSYL(NY)) 
          IAUSYL (NY) =3D =
NINT (AAUSYL) 
          IEUSYL (NY) =3D =
ITUSYL (NY) - IAUSYL (NY) * 100 
 
If you put in ITUSYL =3D -123, then IAUSYL =3D -1 and IEUSYL =3D -23, =
but for ITUSYL =3D 77 you get IAUSYL =3D 1 and IEUSYL =3D -23, but I =
assume since IA represents a single bin the sign is ignored and when =
what1 is packed back up before being written to the USRYIELD output the =
algorithm ignores the sign again by doing 
 
ITUSYL(NY) =3D IEUSYL(NY) + ABS(IAUSYL(NY))*100 
 
Although this is a little confusing on how I was getting the right =
results while not following the explicit instructions of the manual, and =
that the written what1 was different than the supplied what1, the fault =
was mine.  Even so, it might be better to not ignore the sign of IA =
when writing to the USRYIELD output so that it will agree with what the =
user input. 
 
Thanks 
Kerry 
 
-----Original Message----- 
From: Alfredo Ferrari [<A =
HREF=3D"mailto:alfredo.ferrari@cern.ch">mailto:alfredo.ferrari@cern.ch</A=
>] 
Sent: Mon 12/11/2006 7:35 AM 
To: Lee, Kerry T. (JSC-SF)[LMIT] 
Cc: fluka-discuss@fluka.org 
Subject: RE: USRYIELD binary output 
 
Hi Kerry 
 
the code tries to unpack from what(1) the infos about which two =
physical 
quantities you want to deal with. Probably a nonsense number is =
unpacked 
into something strange (I should verify). 
 
            &=
nbsp;    Ciao 
            &=
nbsp;   Alfredo 
On Fri, 8 Dec 2006, Lee, Kerry T. (JSC-SF)[LMIT] wrote: 
 
> This is a multi-part message in MIME format. 
> 
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0 
> Content-Type: text/plain; 
>       charset=3D"iso-8859-1" 
> Content-Transfer-Encoding: quoted-printable 
> 
> I seem to have found that I was in error when explaining the =
problem.  I =3D 
> missed the fact that what12 is stored as the next variable after =
what1, =3D 
> but the value stored where what1 should be is 77 and the value from =
the =3D 
> input file is -123.  How is this not being transferred =
correctly to the =3D 
> output file? 
> 
> Thanks 
> Kerry 
> 
> 
> -----Original Message----- 
> From: Lee, Kerry T. (JSC-SF)[LMIT] 
> Sent: Thu 12/7/2006 10:27 AM 
> To: fluka-discuss@fluka.org 
> Subject: USRYIELD binary output 
> =3D20 
> 
> Dear FLUKA authors, 
> 
> I am reading a USRYIELD output that was written according to the =
=3D 
> following USRYIELD command 
> 
> USRYIELD  -123.00000 201.00000-50.000000 4.0000000 5.0000000 =
=3D 
> 1.0000000LET+0.1  =3D20 
> USRYIELD   1000.0000 0.0000000 100.00000 1.5000000 =
0.0000000 2603.0000 & =3D 
>      =3D20 
> =3D20 
> I am trying to parse the file using a C++ routine, so it is =
imperitive I =3D 
> understand the details of what is written where.  I followed =
the code in =3D 
> the manual on reading the USRYIELD output =3D 
> (<A =
HREF=3D"http://www.fluka.org/manual/sect/s090/text.html">http://www.fluka=
.org/manual/sect/s090/text.html</A>) and see that the =3D 
> variable ITUSYL(NY) is supposed to correspond to the What1.  =
When I read =3D 
> this file the variable reads a value of 2603, which is my What12 =
and not =3D 
> What1.  I have taken into account the extra bytes that fortran =
writes =3D 
> and all the other values are read in correctly, so is this a bug or =
am I =3D 
> missing something? 
> 
> Thanks 
> Kerry 
> 
> 
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0 
> Content-Type: text/html; 
>       charset=3D"iso-8859-1" 
> Content-Transfer-Encoding: quoted-printable 
> 
> <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML =
3.2//EN"> 
> <HTML> 
> <HEAD> 
> <META HTTP-EQUIV=3D3D"Content-Type" =
CONTENT=3D3D"text/html; =3D 
> charset=3D3Diso-8859-1"> 
> <META NAME=3D3D"Generator" CONTENT=3D3D"MS =
Exchange Server version =3D 
> 6.5.7651.34"> 
> <TITLE>RE: USRYIELD binary output</TITLE> 
> </HEAD> 
> <BODY> 
>  
> 
> I seem to have found that I was in =
error when =3D 
> explaining the problem.&nbsp; I missed the fact that what12 is =
stored as =3D 
> the next variable after what1, but the value stored where what1 =
should =3D 
> be is 77 and the value from the input file is -123.&nbsp; How =
is this =3D 
> not being transferred correctly to the output file? 
> 
> Thanks 
> Kerry 
> 
> 
> -----Original Message----- 
> From: Lee, Kerry T. (JSC-SF)[LMIT] 
> Sent: Thu 12/7/2006 10:27 AM 
> To: fluka-discuss@fluka.org 
> Subject: USRYIELD binary output 
> 
> 
> Dear FLUKA authors, 
> 
> I am reading a USRYIELD output that was written according to the =
=3D 
> following USRYIELD command 
> 
> USRYIELD&nbsp; -123.00000 201.00000-50.000000 4.0000000 =
5.0000000 =3D 
> 1.0000000LET+0.1&nbsp;&nbsp; 
> USRYIELD&nbsp;&nbsp; 1000.0000 0.0000000 100.00000 =
1.5000000 0.0000000 =3D 
> 2603.0000 =
&amp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&amp=
;nbsp; 
> 
> I am trying to parse the file using a C++ routine, so it is =
imperitive I =3D 
> understand the details of what is written where.&nbsp; I =
followed the =3D 
> code in the manual on reading the USRYIELD output (<A =3D 
> HREF=3D3D"<A =
HREF=3D"http://www.fluka.org/manual/sect/s090/text.html">http://www.fluka=
.org/manual/sect/s090/text.html</A>"><A =
HREF=3D"http://www.fluka=3D">http://www.fluka=3D> 
> .org/manual/sect/s090/text.html</A>) and see that the =
variable =3D 
> ITUSYL(NY) is supposed to correspond to the What1.&nbsp; When I =
read =3D 
> this file the variable reads a value of 2603, which is my What12 =
and not =3D 
> What1.&nbsp; I have taken into account the extra bytes that =
fortran =3D 
> writes and all the other values are read in correctly, so is this a =
bug =3D 
> or am I missing something? 
> 
> Thanks 
> Kerry 
> 
> 
> 
> 
> </BODY> 
> </HTML> 
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0-- 
> 
 
-- 
 
+------------------------------------------------------------------------=
----+ 
|  Alfredo =
Ferrari           =
     ||  Tel.: =
+41.22.767.6119         &nbs=
p;        | 
|  =
CERN-AB           =
             =
||  Fax.: =
+41.22.767.7555         &nbs=
p;        | 
|  1211 Geneva =
23           &nbsp=
;     ||  e-mail: =
Alfredo.Ferrari@cern.ch        | 
|  =
Switzerland          &n=
bsp;         =
||          =
Alfredo.Ferrari@mi.infn.it     | 
+------------------------------------------------------------------------=
----+

</BODY>
</HTML>
------_=_NextPart_001_01C71D35.0BE304A4--

Next message: me@marychin.org: "Low energy neutron xsec not found for some media"

Previous message: me@marychin.org: "Re: miraculous zero photoneutron production?"
In reply to: Alfredo Ferrari: "RE: USRYIELD binary output"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This archive was generated by hypermail 2.1.6 : Wed Dec 13 2006 - 13:59:05 CET