[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: Help in MCNP5

One possibility (not quite what you wanted) is to use a programming language such as Visual Basic to modify the input file and run MCNP through multiple cycles.  I did something like this using a DOS batchfile for the verification runs when I installed MCNP, and this allowed me to do them in the night and weekend "down" time for my desktop.



Dave Neil



-----Original Message-----

From: Pramoth CM [mailto:pchandr@CLEMSON.EDU]

Sent: Tuesday, November 23, 2004 10:07 AM

To: radsafe@list.vanderbilt.edu

Subject: Help in MCNP5





Hi Group,

	I am a student working with MCNP5 and using a very small code to obtain a 

tally of gamma ray spectrum of a particular gamma ray energy for a given 

number of events or live time. I want to modify this code so that the code 

would generate spectrums at increments at10 kev interval for 0 - 1000 kev.

The only option of doing this which I can think of is running the original 

code over and over at 10 kev increments - which would be time consuming.

Could anybody suggest the required change in the code, so that I could run 

the code only once and then obtain a single OUTPUT file with a 100 columns 

( each column representing the spectrum at 10 kev intervals) . I am very 

new to MCNP5 and read the required parts of the MCNP5 manual and could not 

figure out how to do this. I will be really grateful if anybody could help. 

I have provided the original code below.





c ************************** Cell definition cards **************************

c RightCircularCylinder Detector Assembly

c Sensor Cell Definition Cards

c

     1   117  -3.67 (-1 2 -3 ) imp:p 1 $ NaI Sensor

c

c  End Cap

c

    20   213  -2.70   (84 -85 -82)                   imp:p 1 $ Al EndCap Window

    21   444  -8.00   (82 -83 -85 81 ):(-80 81 -82 ) imp:p 1 $ SS EndCap

c

c  Vacuum

c

    32     0 (-84 1 -82)                            imp:p 1 $ Vacuum

c

c Everything else...

c

   988 220 -0.0012 (-900 -901 902 )(-81: 83: 85) &

                    &

                   imp:p 1 $ Detector environment

   999  0  (900: 901: -902)   imp:p 0 $ Outside world



c ************************* Surface definition cards *************************

c

c Sensor Surface Definition Cards

c

     1    pz      0.00 $ Top Face of sensor

     2    pz     -5.08 $ Bottom Face of sensor

     3    cz      2.54 $ Outer Radius of sensor

c

c End Cap and End Cap Window

c

    80    pz    -5.08 $ Inside Lower bound of detector end cap

    81    pz    -5.3975 $ Outside Lower bound of detector end cap

    82    cz     2.54 $ Inside radius of end cap

    83    cz     2.8575 $ Outside radius of end cap

    84    pz     0.5 $ Inside detector Face

    85    pz     0.55 $ Outside detector face

c

c Bounding surfaces

   900    cz  3.8575 $ Bounding cylinder

   901    pz  12 $ Bounding top Surface

   902    pz  -6.3975 $ Bounding bottom Surface



c ******************************** Data cards ********************************

c

Mode P

c

c Material Cards

c

m220     6000.01p    -0.000124 $ C   Air, rho = -0.001205 g/cc

          7000.01p    -0.755267 $ N

          8000.01p    -0.231781 $ O

         18000.01p    -0.012827 $ Ar

c

m117    11000.01p   1    $ NaI,  rho = -3.67 g/cc

         53000.01p   1

c

m444    6000.01p   -0.0024    $ Stainless Steel,  rho = -8 g/cc

         14000.01p   -0.064

         15000.01p   -0.0028

         16000.01p   -0.0016

         22000.01p   -0.032

         24000.01p   -1.44

         25000.01p   -0.16

         26000.01p   -5.5

         27000.01p   -0.002

         28000.01p   -0.8

c

m213    13000.01p   1    $ Aluminum,  rho = -2.7 g/cc

c

c Tally defination card(s)

c

f8:p  (1) $ Pulse Height tally on cell 1

c  Gaussian Energy Broadening Term

ft8 geb 0  0.05086  0.30486  $ for NAI

c

c Source definition cards

c

sdef  pos= 0  0  10  erg=d3

c

c Remember! The calculated efficiency is per photon emitted, so if you have 

more

c than one photon energy, and they are all equally probable, multiply by N

c

  #   si3    sp3

      L       D

    1.0000   1.565E+15

c

c  Energy Range and bin width [in MeV!]

c

e0 0 1e-5 8e-3 512i 4.112 $ 8 KeV bins up to 4096 KeV

c

c Number of particles to follow

c

nps 67900000

c

c .OR. Minutes to run

c

c ctme   1

c

c Print & Dump Cycle card to create MCTAL output file

c

prdmp 2j 1



Thanks,

Pramoth





************************************************************************

You are currently subscribed to the Radsafe mailing list. To

unsubscribe, send an e-mail to Majordomo@list.vanderbilt.edu  Put the

text "unsubscribe radsafe" (no quote marks) in the body of the e-mail,

with no subject line. You can view the Radsafe archives at

http://www.vanderbilt.edu/radsafe/



************************************************************************

You are currently subscribed to the Radsafe mailing list. To

unsubscribe, send an e-mail to Majordomo@list.vanderbilt.edu  Put the

text "unsubscribe radsafe" (no quote marks) in the body of the e-mail,

with no subject line. You can view the Radsafe archives at

http://www.vanderbilt.edu/radsafe/