Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Beginners quide, FAQ, examples? #159

Open
erosenst opened this issue Mar 15, 2022 · 4 comments
Open

Beginners quide, FAQ, examples? #159

erosenst opened this issue Mar 15, 2022 · 4 comments

Comments

@erosenst
Copy link

Hello is it possible to add some runable examples in sources directory? (simple .sh with input file?)

I'm currently stuck point where tests gone 100% ok. But how inport/define my endf library to njoy21? Simple commandline example?

I just try simple njoy21 -i input.txt but got:
[error] Encountered unexpected argument
[info] Error while parsing line 8

’PENDF tape for ENDFB-6.8 92235(0)’


input.txt
moder
1 -21
*PENDF tape for ENDFB-6.8 92235(0)*/
20 9228
0/
reconr
-21 -22
’PENDF tape for ENDFB-6.8 92235(0)’
9228 1 0 /
.001/
’92235(0) from ENDFB-6.8’/
0/
broadr
-21 -22 -23
9228 1/
.001/
600/
0/
stop

So I just try dithering resolution.. Is there any hidden commandline switches or enviroment variables to define where to look input file?
Or simple explanation how to use that program..

This is possibly explained somewhere in manual but not yet find it.. 
Thank you.  
@nathangibson14
Copy link
Contributor

nathangibson14 commented Mar 15, 2022 via email

@erosenst
Copy link
Author

Thank you for reply. I found by strace that it looks input file in same directory and ln -s library file to working dir.'tape20'
Then it run test run.. (removed comments, they seem be unrealible feature)
Your link helps lot.. Examples dir don't helped for my question how link input files but figured it out..

Then try some more complicated and stuck to thermr:

[info] Mat Number: 9228
[warning] Material specific default value not found in table
default value of 25 eV used.

[error] Encountered invalid value for mtref
[info] Error while parsing line 23

0 9228 8 1 1 0 1 221 2


[info] 
mt values are labels for endf sections, where the sections are usually
used to hold different reactions.

The mtref argument is the mt for the desired inelastic option. These
mt values must be between 221 and 250, inclusively.

[info] Trouble while validating THERMR Card 2
[info] Trouble while validating THERMR input
[info] Error while running NJOY21

I used input:
thermr
 0 25 42
0 9228 8 1 1 0 1 221 2
600
0.001 1.0

So it says last element 9th=2 is wrong MT number, range but for documentation it says 8th should be MT and 9th is print option.
Is input syntax changed in njoy21 or what is wrong?

     *  card 2                                                        *
     *     matde      material desired on endf tape                   *
     *     matdp      material desired on pendf tape                  *
     *     nbin       number of equi-probable angles                  *
     *     ntemp      number of temperatures                          *
     *     iinc       inelastic options                               *
     *                   0     none                                   *
     *                   1     compute as free gas                    *
     *                   2     reserved                               *
     *                   3     reserved                               *
     *                   4     read s(a,b) and compute matrix         *
     *     icoh       elastic options                                 *
     *                   0     none                                   *
     *                   1     graphite                               *
     *                   2     beryllium                              *
     *                   3     beryllium oxide                        *
     *                  11     polyethylene                           *
     *                  12     h(zrh)                                 *
     *                  13     zr(zrh)                                *
     *     natom      number of principal atoms                       *
     *     mtref      mt for inelastic reaction (201-250 only)        *
     *     iprint     print option (0=minimum, 1=maximum,             *
     *                2=max. normal + intermediate results)           *
     *                (default=0)                                     *

@erosenst
Copy link
Author

Is input syntax changed in njoy21 or what is wrong?

That was just case. njoy21 have ten variables in thermr card 2:

@whaeck
Copy link
Member

whaeck commented Mar 16, 2022

Hi.

NJOY2012 actually modified the NJOY input for card 2 in THERMR - which is causing your issue. It added the iform argument between icoh and natom:

   !  card 2
   !     matde      material desired on endf tape
   !     matdp      material desired on pendf tape
   !     nbin       number of equi-probable angles
   !     ntemp      number of temperatures (default = 1)
   !     iinc       inelastic options
   !                   0     none
   !                   1     compute as free gas
   !                   2     read s(a,b) and compute matrix
   !     icoh       elastic options
   !                   0     none
   !                   1     compute using ENDF6 format data
   !                   --------or for earlier formats
   !                   1     graphite
   !                   2     beryllium
   !                   3     beryllium oxide
   !                  11     polyethylene
   !                  12     h(zrh)
   !                  13     zr(zrh)
   !     iform      output format for inelastic distributions
   !                  0      E-E'-mu ordering (MF6 special)
   !                  1      E-mu-E' ordering (MF6/Law7)
   !     natom      number of principal atoms
   !     mtref      mt for inelastic reaction (221-250 only)
   !     iprint     print option (0=minimum, 1=maximum,
   !                2=max. normal + intermediate results)
   !                (default=0)

This is what is currently in the NJOY2016 manual and in the NJOY2016 source code.

Best regards,
Wim

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants