Home

ATOMIC ENERGY & S & L`ENERGIE ATOMIQUE OF

image

Contents

1. 22 that k is Boltzmann s constant k 8 6275 x 107 gt eV K therefore 5kT 0 1264 eV The Maxwellian spectrum is placed in the region below xkT 3 THETA fission spectrum temperature 1 27 10 ev 4 FCUT point at which a fission spectrum is joined to spectrum 4 selected by IW 67 4 x 10 eV The fission spectrum is placed in the region above FCUT 5 Not used enter zero 0 6 Not used enter zero 0 7 Not used enter zero 0 8 Not used enter zero 0 9 Not used enter zero 0 10 Not used enter zero 0 T Terminate Block 3 Note that if none of the arrays in Block 3 is entered a T must still be entered The data for the NNUC ENDF materials to be processed are stacked one set after the other A set of data for a nuclide consists of a title card and Blocks 4 and 5 as follows Title Card 20A4 format This card is a title card for the nuclide The title will be carried throughout AMPX with the neutron cross sections The first 48 characters will be used as a title for ANISN cross section sets subse quently produced by the NITAWL module Block 4 70 Nuclide General Information 9 3 ID19 MATNO NTEMP LORDER NFY MATID 23 The identification number for this nuclide This number is carried as the identifier of the data on the master neutron cross section library produced by XLACS2 The ENDF B material number MAT number of the set of d
2. GETAP GXWIMS LAMDA NPTRD Error WB N mm OR DWN DA w bt 765 Description Problem with 14 burnup chain edit input cannot find referenced NIN 1 array read error ITYPE invalid must be 1 or 2 Available core exhausted t Invalid WIMSTAR tape unit no Can t find WIMSTAR record NT unit no IDN nuclide identifier IT record type T temperature ISF resonance table type Invalid data source for new nuclide Invalid NIN detected Invalid RIN detected Can t find a resonable value for APA 1 array read error 2 array read error Temperature or a limit exceeded NTE lt 10 and NTE NSP lt 20 3 4 5 array read error 6 array read error Program error WIMSTAR tape structure invalid Energy structure does not match existing WIMSTAR tape 7 8 9 10 Total Fast Resonance Thermal Group structure doesn t match point data Interpolation limit of 25 exceeded on NXS tape Available core exhausted t Temperature on NXS tape does not match 4 array entry for nuclide ID19 t See footnote on page 79 Routine PPGBW P1SCAT RESCAL RESON RPLCE SCAN SIGSET SMOOTH Error FF whe Ee N Le 2 3 4 8 9 10 2 6 LT 2 Description NXS tape problem invalid p detected No NAM tape and energy structure not card input NAM tape has wrong no of groups NWO and NEP are inconsistent 6 array read error 10 array read error NEP
3. where GRGN defines the core region in K bytes GTIME sets the task time limit in minutes and any of the scratch units 14 through 19 that are required are activated as shown Each load module required is concate nated to the steplib as G0 STEPLIB DD DD NPTXS load module DD WIMSTAR load module 11 DD XLACS2 load module Each data file required is also defined to the procedure as GO FTnnF001 DD file description where nn is the two digit FORTRAN logical unit number Finally the card data is entered as GO SYSIN DD data cards The file requirements and data card input for each module are listed in Section 3 2 for the AMPX II modules and Section 3 3 for WIMSTAR The core region parameter is discussed in the next section All modules write printed output to unit 6 and if necessary punched output to unit 7 3 5 DYNAMIC CORE ALLOCATION WIMSTAR uses the AMPX system of dynamic core allocation whereby all available core in the user s region is allocated to the program The use of this core is dynamically determined at execution time by the 51 requirements of the specific task The size required is wholly dependent on the particular nuclide being processed and the processing required If the size requested through the GRGN parameter of the AMPX2 procedure is insufficient a WIMSTAR error message is printed if the system detects this problem an IHO240I message is printed and the user simply inc
4. J 1 NO Records 5 JC JB AA k y JJ k k 1 JA to LAS where JA JC 2 1 and may differ for each record 57 Energy group boundaries in order of decreasing energy Fission source spectrum sums to 1 0 A record for each nuclide describing all nuclides produced by burnup of this nuclide This record contains Length of record JC NIN JB Yield normally 1 0 and NIN of nuclide formed by neutron capture in nuclide JB Decay constant of nuclide JB and NIN of nuclide formed by decay Fission energy yield and indi cator NFA see below Pairs of numbers giving yield and NIN of fission products produced by fission of JB Pairs of numbers 0 0 and NIN of any nuclides formed only indirectly 238 240 from JB e g U ntn Pu s Bus de gt Sur requiring specification of Pu 0 Sm respectively The indicator NFA is equal to NF for nuclide JB unless nuclide JB is a fission product in this case NFA is 1 if nuclide JB does not have a resonance tabulation associated with it and 2 if it does If yield and NIN of nuclide formed by capture in nuclide JB are both set to zero no capture product will be assumed in any subsequent calculation of burnup File Mark 58 Nuclide Files For each of the nuclides specified in record 2 of the general index file there is a file of data containing the following informa f tion 7 Record 1 J AW TAN NF NT NZZ NIN at
5. 65 Record 5 9 zeros Groups of records 1 5 are repeated for each MAT on the tape System EOF mark The variables have the same meaning as for the ENDF B tapes see Reference 2 4 3 4 PXS FD Tape Format This tape contains finitely dilute point cross sections for the unresolved resonance region Record 1 1 1D19 nuclide identification no 2 T temperature 3 SIGP oy value 4 10 zero Record 2 N E 1 0 1 1 1 N Record 2 is repeated for Fay Ff Oy Records 1 and 2 are repeated for each Ip temperature and nuclide on the tape System EOF mark 2 cross sections a el fp Gaa are total elastic scattering ps ssion ar n y respectively 66 5 SUMMARY WIMSTAR Version 4 is a FORTRAN IV computer program developed for inclusion as a module in the AMPX II ENDF B processing system In conjunction with existing AMPX modules WIMSTAR provides the capability of generating library data files for the WIMS lattice code and updating the WIMS library tape WIMSTAR adheres to the conventions of the AMPX system and can easily be implemented on a computer currently supporting the AMPX system Several minor modifications to AMPX II module NPTXS were necessary AMPX II currently accesses ENDF B IV tapes modifications to AMPX modules XLACS2 and NPTXS and to WIMSTAR will be required when ENDF B V is released for general use WIMSTAR is divided into six separately executable segments to accomplish
6. Transport cross section Ter N1 N2 Absorption cross section RE N1 N2 Chi x not used in WIMSC or later versions GXWIMS places zeros in this array which can be overwritten using cards if desired N2 Goldstein Cohen parameters N2 Fission yield cross section Vag Required only if NF gt 2 N1 N2 Fission cross section Ops Required only if NF gt 2 N1 4N2 PO condensed scattering matrix from fast and resonance groups NP4 Temperatures for thermal data tabulations The number and values of the temperatures may be different from the old WIMS tape However if this is the case all thermal data must be replaced In other words if ID4 lt 0 NT must equal NT of old tape and 23 is not required if ID4 gt 0 23 is necessary NT Length of thermal PO scattering matrix for each temperature Required only if ID4 0 NT Thermal edit flags NT entries corresponding to each thermal temperature 0 no edit or input 1 edit or input Block 7 is input for each 1 in 25 array NT Terminate Block 6 Block 7 Repeat for each 1 in 25 array and NES gt 0 18 19 20 21 Transport cross section oe N3 Absorption cross section CAE N3 Fission yield cross section vog Required only if NF gt 2 N3 Fission cross section ee Required only if NF gt 2 N3 42 26 PO condensed scattering matrix from thermal groups Entries in 24 array T Terminate Block 7 Sets of Bl
7. absorption or fission yield depending on whether this is the first or second tabulation Ordered oy within temperature NTExNSP Terminate Block 10 Sets of Blocks 8 9 and 10 are repeated NER times Block 11 NEP gt 0 6 T P1 Scattering Matrix Edit Control 1 1 ID6 Data source O old WIMS tape O0 card input gt 0 WIMSTAR tape ID6 is identification no of data A Terminate Block 11 Block 12 NEP gt 0 10 T Pl scattering matrix ordered to groups within from groups NGxNG Terminate Block 12 Sets of Blocks 11 and 12 are repeated NEP times e Note There are four Pl matrices on a WIMS tape one each for hydrogen deuterium oxygen and carbon in that order If NEP 0 they are copied from old to new tape directly and Blocks 11 and 12 are not required If NEP 4 each is edited in turn End of Data Notes 1 File requirements WIMS tape input WIMS tape output WIMSTAR data tape input scratch units 14 15 3 3 3 SCAN Block 1 1 Run Parameters 3 1 NNUC No of nuclides to be printed 2 Nwo WIMS data tape unit no 3 INEX Inclusive exclusive print option 0 print section if either nuclide is included in NLIST or section is included in IOUT 1 print section only if both nuclide is included in NLIST and section is included in IOUT T Terminate Block 1 Block 2 2 IOUT set of 10 flags 0 1 no print print one for each of the fol
8. the resonance group cross sections are infinitely dilute i e no resonance shielding Note that these cross sections can in general be consistent only with one resonance tabulation K AA I I 1 K Condensed scattering matrix for scattering from fast and resonance groups The vector AA I may be split into N1 N2 blocks of the form AS AL VECTOR I I 1 L where the numbers in VECTOR are the non zero scattering cross sections from a group L AL is the number of these cross sec tions and AS is the position of the self scattter term in the array VECTOR Note that AS 1 if there is no upscatter TEMP J J 1 NT Temperatures K at which thermal daca is tabulated in ascending order If 60 NT 1 there is one thermal data tabulation which is then used for all tempera tures in WIMS case TEMP 1 is usually in this arbitrarily set equal to 300 K Following this record there are three records for each temperature as follows Record 6 Record 7 Record 8 TR J J N1 N2 1 N ABS J J N1 N2 1 N XNUFIS J J N1 N2 1 N FIS J J N1 N2 1 N o nd transport and tr a P R absorption cross sections for thermal groups for the appropriate temperature vo and Tes fission yield f and fission cross sections for thermal groups Record 7 is present only for fissile nuclides i e NF gt 2 KA AAT I 1 KA File Mark I 1 Condensed scattering m
9. 1961 70 TABLE 1 69 GROUP WEIGHTING SPECTRUM REGION A Groups 1 to 14 inclusive Ener Dry 69 Group Wet 69 Group By Weighting Spectrum Weighting Spectrum 11 3315 Mev 0 000077 0 0003 8 82497 0 000288 0 0012 6 87289 0 001082 0 0043 5 35261 0 003415 0 0125 4 16862 0 007253 0 0250 3 24652 0 015299 0 0477 2 52839 0 032235 0 0821 1 96912 0 040766 0 0925 1 53355 0 068804 0 1270 1 19433 0 092740 0 1446 0 930145 0 151457 0 1831 0 724397 0 293853 0 2321 0 564161 0 359617 0 2432 0 439639 0 273385 0 2287 0 342181 0 579860 0 3326 0 266491 0 797033 0 3665 0 207543 1 020714 0 4096 0 161635 1 108708 0 4687 0 125881 1 392871 0 5421 98 0635 Kev 1 538518 0 6346 76 3509 2 025529 0 7520 59 4621 2 256529 0 8839 46 3092 2 717019 1 0847 36 0656 2 929716 1 2975 28 0879 2 913247 1 6002 21 8749 5 728346 2 0012 17 0362 5 345382 2 4685 13 2678 6 613002 3 1004 10 3330 8 260331 3 8690 8 04733 8 314061 4 8358 The first spectrum given above was used as the weighting spectrum for GALAXY in producing the 14 fast groups for every element except hydrogen and oxygen The latter elements used the second spectrum REGION B In Groups 15 to 56 inclusive the weighting is proportional to ED REGION C In Groups 57 to 69 inclusive the weighting is proportional to E kT 1 a 2 Eo where kT 2 522 x 10 8 Mev T This table was taken from Reference 3 JL TABLE 2 69 GROUP ENERGY BOUNDARIES FOR WIMS Lethargy Energy
10. 3 2 1 NPTXS Block 1 0 Broadening Method 1 SIGMAL method default PSI CHI method 1 ID O 1 1 Case Description 1 1 NNUC number of ENDF nuclides to be processed T Terminate Block 1 Repeat Blocks 2 and 3 for NNUC nuclides Block 2 2 Nuclide Selection Array 6 1 MATNO ENDF MAT number for nuclide 2 NDFB unit number for ENDF tape 11 3 IDTAP O 1 no check check ENDF library label 1 4 MODE ENDF library mode 0 2 binary BCD 0 5 NSIGP b of o values 1 6 1D19 nuclide identification number on point cross section tapes 15 3 Nuclide Options 4 1 RFACT the resolved resonance region energy mesh is based directly on this parameter Points are chosen such that the ratio of total cross section from point to point is roughly RFACT 0 9 2 SFACT number of practical widths over which the RFACT scheme is used for a particular resonance 10 3 SIGP o for unresolved calculation Enter 0 if NSIGP gt 1 and use 4 instead 10 4 TDEGIC temperature in Kelvin for Doppler broadening 0 T Terminate Block 2 Block 3 NSIGP gt 1 4 array for the unresolved calculation if the number of a S is greater than 1 Enter 25 values high to low maximum number is 11 NSIGP T Terminate Block 3 End of Data Notes 1 The Le in this report and ge in the AMPX II user s manual both refer to the potential scattering cross section See the defini tion
11. Lethargy Energy Width Sroup Energy Width Width Width MeV 1 10 0 6 0655 3 9345 0 49997 28 4 00 3 30 0 700 0 19237 2 6 0655 3 679 2 3865 0 49998 29 3 30 2 60 0 700 0 23841 3 3 679 2 231 1 448 0 50019 30 2 60 2 10 0 500 0 21357 4 2 231 1 353 0 878 0 50013 31 2 10 1 50 0 600 0 33647 5 1 353 0 821 0 532 0 49956 32 1 50 1 30 0 200 0 14310 6 0 821 0 509 0 321 0 49592 33 1 30 1 15 0 150 0 12260 7 0 500 0 3025 0 1975 0 50253 34 1 15 1 123 0 027 0 02376 8 0 3025 0 183 0 1195 0 50260 35 1 123 1 097 0 026 0 02342 9 0 183 0 1110 0 072 0 49996 36 1 097 1 071 0 026 0 02399 10 0 1110 0 06734 0 04366 0 49978 37 1 071 1 045 0 026 0 02458 il 0 06734 0 04085 0 02649 0 49985 38 1 045 1 020 0 025 0 02421 12 0 04085 0 02478 0 01607 0 49987 39 1 020 0 996 0 024 0 02381 13 0 02478 0 01503 0 00975 0 49999 40 0 996 0 972 0 024 0 02439 14 0 01503 0 009118 0 005912 0 49980 41 0 972 0 950 0 022 0 02289 42 0 950 0 910 0 040 0 04302 43 0 910 0 850 0 060 0 06821 ev 44 0 850 0 780 0 070 0 08594 45 0 780 0 625 0 155 0 22154 15 9118 0 5530 0 3588 0 0 50006 46 0 625 0 500 0 125 0 22314 16 5530 0 3519 1 2010 9 0 45198 47 0 500 0 400 0 100 0 22314 17 3519 1 2239 45 1279 65 0 45198 48 0 400 0 350 0 050 0 13353 18 2239 45 1425 1 814 35 0 45199 49 0 350 0 320 0 030 0 08961 19 1425 1 906 898 518 202 0 45197 50 0 320 0 300 0 020 0 06454 20 906 898 367 262 539 636 0 90395 51 0 300 0 280 0 020 0 06899 21 367 262 148 728
12. NR approximation and a value of 0 yields the Infinite Mass IM approximation Using the p int cross sections for infinite dilution GXWIMS attempts to find a value of that yields the same resonance integral for absorption as that obtained by solving the slowing down equation In some cases the IR approximation is not suffi ciently accurate to yield a value of between O and 1 and the user must supply a value manually The second phase of GXWIMS obtains the weighted multigroup cross sections from XLACS2 and calculates the following quantities CE total P scattering group i to j elastic inelastic 2 n 2n 1 elastic P scattering group i to j Sij 1 GT total of E 6 gnr2n i i aij i J a ST o 0 absorption o E O i i z ij 3 of fission yield v of i i o potential scattering 54 slowing down power amp o u u i j ij i l1 i There are several options for calculating the transport cross section 7 29 The recommended procedure is the following below 4 eV a row sum correction is applied tr _ a o 1 1 oi di f dij 3 z Si j j above 4 eV a weighted column sum correction is applied ot of 209 z2 Ir VW ot i i j ij 3 j ji ji where between 4 eV and 9 2 KeV a 1 E current is assumed i e Wij hu Au Au lethargy width of group i and above 9 2 KeV the weights are obtained as Wij J J where J neutron current in group i GXWIMS has no facility to p
13. System for one resonance absorber and one moderator hydrogen INTGRC performs the integration to calcu late the resonance integral LAMDA calculates the Goldstein Cohen A values for the first temperature and smallest 2 For each successive o5 the unresolved resonance region cross sections are replaced provided a PXS FD tape is available The following service subroutines are used to manipulate the point cross sections ADDXS COMPRS EXPAND LINEAR MULT RESINT RPLCE and THNFIT Results are printed plotted and written to the WIMSTAR tape by RESOUT The multigroup phase of GXWIMS is handled by WIMSXS TREAD reads the AMPX master library and VECTOR and MATRIX accumulate the required one and two dimensional cross sections respectively WIMS cross sections are calculated as per Section 2 3 1 condensed scattering matrices are generated by WIMFMT and the results are printed plotted and written to the WIMSTAR tape UPDATE is the main WIMS library update routine of segment WIMLIB Each of the four sections burnup chains multigroup data resonance data and P1 scattering matrices are updated separately by subroutines BURNUP SMOOTH RESON and PISCAT respectively In each case the program compares the NIN and RIN lists input by the user with those on the old tape to determine whether the data is to be added edited or deleted Before each record is written to the new WIMS tape FIDO input routine FIDAS is called to allow the user to ov
14. data generation post generation data manipulation and WIMS tape update and format conversion Multigroup cross sections are calculated by AMPX II module XLACS2 and placed in the WIMS format by WIMSTAR The major development effort was devoted to resonance processing The slowing down equation is solved numerically using point cross sections tabulated on a very fine energy mesh to obtain resonance integrals for the absorber mixed with pure hydrogen WIMSTAR can generate all the WIMS data required with the exception of fission source spectra burnup chain data and x Future improvements could include providing group structure condensation an upgraded Goldstein Cohen A calculation and better facility for display of cross sections 67 6 ACKNOWLEDGEMENTS The author would like to express his thanks to Dr D Hamel for his guidance in the completion of this work The author would also like to acknowledge the following contributions the AMPX II code system obtained from the Radiation Shielding Information Centre Oak Ridge module RESPU of the RSYST code system also obtained from RSIC and parts of the program LINEAR obtained from Lawrence Livermore Laboratory 10 11 12 68 REFERENCES J R Askew F J Fayers P B Kemshell A General Description of the Lattice Code WIMS J Brit Nucl Energy Soc 5 564 1966 M K Drake ed Data Formats and Procedures for the ENDF Neutron Cross Section Lib
15. device 20 The 6th 7th and 9th entries of the 70 array are used to identify the thermal data when the fast and thermal data reside on separate devices 10 LCSM If IW 4 enter 0 If IW 4 enter the number of regions used to specify the arbitrary weighting function in the 5 array 0 11 MCSM If IW 4 enter 0 If IW 4 enter the number of points used to specify the arbitrary weighting function in the 6 array 0 12 NDNP Not used Enter 0 0 38 Six triggers I PT i to specify amount of printed output desired 6 IGPT i Correspondence Resolved and unresolved resonance processing Averaged cross sections by energy group and process Elastic scattering matrices Thermal scattering matrices O if applicable as processed Inelastic scattering matrices gt gt 0 D g g H w o u ui v ua tal uw wo o A a 4 4 ke est R n 2n scattering matrices Edit Po data in APMX formats Edit in convenient formats x The Po data if applicable are edited with an IMPT i 1 a value of 2 includes the Po data x The terminology convenient is used to imply a convenient form from a programming viewpoint k The default value for I PT i is zero T Terminate Block 2 21 Block 3 4 5 6 7 8 ENDF MT identifiers of the reactions for which cross sections are to be punched IDTAP If IW 4 enter the interpolation s
16. for WIMS TApe Revisor The program can calculate all the data required by WIMS with the exception of the fission source spectrum burnup chains and x used in calculating f p in WIMSB x is not used in WIMSC or later versions and can be set to zero These data can be taken from the old WIMS tape or created elsewhere and added to the new WIMS tape via the update routine Several built in weighting spectra for producing the multi group cross sections from point data are provided However if he prefers the user can substitute his own or use the spectra provided by Winfrith Calculation of the transport cross section supports the methods suggested by Winfrith 8910 The current version of AMPX II accesses ENDF B IV tapes When Version yl becomes more widely used AMPX II and WIMSTAR 4 will be modified accordingly WIMSTAR 4 and other modules of the AMPX II system are cur rently being run at an IBM installation Use on another system would be difficult due to the complications in implementing the AMPX system Section 2 describes the structure of WIMSTAR 4 and the AMPX II modules required Section 3 provides instructions for using the current versions of the programs Section 4 contains information to aid in making modifications to WIMSTAR 4 Section 5 summarizes the current state of the program and discusses future improvements Appendix A describes the error messages that may be received when something goes wrong while Appendix
17. must equal 0 or 4 Can t find correct point data on NXS tape Available core exhausted AWA invalid Incorrect temperature on data from NSD tape WIMS tape problem in resonance tables Problem with resonance table edit data check RIN s or program error 10 11 12 array read error Available core exhausted 13 array read error Problem with resonance table edit data check RIN s 5 array read error Energy range of NSD tape does not match that of NXS tape 1 array read error 2 3 array read error Can t find material MAT on ENDF tape 4 array read error Problem with smooth edit control data check NIN lists for order and content See footnote on page 79 Routine SUMARY TAPMAN TERP1 TRANS TRANS2 TRANS3 UPDATE Error 78 Description Error in reading one of 10 11 12 13 14 15 16 17 23 248 25 26 arrays Attempt to change no of temperatures invalid Error in reading one of 18 19 20 21 26 arrays Program error WIMSTAR tape structure invalid 1 srray read error 2 3 array read error Interpolation code out of range Zero or negative value can t be interpolated by logs 1 2 array read error Program error WIMSTAR tape structure invalid Energy structure on WIMS tape does not match existing WIMSTAR tape 3 4 5 6 Total Fast Resonance Thermal Group boundaries on WIMS tape don t match existing WIMSTAR tape Maximum no
18. resonance region These point cross sections are generated by module NPTXS Module XLACS2 then calcu lates multigroup cross sections for use in creating the second set above while module WIMSTAR generates the third set Intermediate results are stored in data files external to each module Through the organization of the AMPX II system any number of modules or segments of WIMSTAR 4 can be executed in one job step The steps in generating WIMS data are summarized below Module NPTXS evaluates the ENDF B resonance parameters to produce point cross sections File PXS ID contains resolved resonance region cross sections and infinitely dilute unre solved resonance region cross sections while file PXS FD contains unresolved cross sections at several user selected finite dilutions Module XLACS2 calculates weighted multigroup cross sections in the WIMS energy group structure and writes these to a standard AMPX master library tape Module RADE can be used to check and partially list the data and module AJAX can be used to edit the data on these master tapes Module WIMSTAR 4 performs the remaining steps 4 2 2 Segment GXWIMS collects the required data from the ENDF B tape the AMPX master library and the two point cross section tapes and generates the WIMS library data which are placed on a WIMSTAR data tape Segment WIMLIB then uses the WIMSTAR data tape to update the WIMS library by creating a new WIMS tape Se
19. 0I error occurs accompanied by a Run Com pletion code of SOC1 Increase REGION and rerun the step 80 APPENDIX B FIDO INPUT SYSTEM The FIDO input method is especially devised to allow the entry or modification of large data arrays with minimum effort Special advantage is taken of patterns of repetition or symmetry wherever pos sible The FIDO system was patterned after the input method used with the FLOCO coding system at Los Alamos and was first applied to the DTF II code Since that time numerous features requested by users have been added a free field option has been developed and the application of FIDO has spread to innumerable codes This description was taken from Reference 4 of the main text The data are entered in units called arrays An array com prises a group of contiguous storage locations which are to be filled with data at one time These arrays usually correspond on a one to one basis with FORTRAN arrays used in the program A group of one or more arrays read with a single call to the FIDO package forms a block and a special delimiter is required to signify the end of each block Arrays within a block may be read in any order with respect to each other but an array belonging to one block must not be shifted to another block The same array can be entered repeatedly within the same block For example an array could be filled with O using a special option and then a few scattered locations could be ch
20. 2 8 LNIN NG NO N1 NZ N3 NNF NNFP 8 LNIN nuclide identification nos 9 LRIN resonance table identification nos 10 NG 1 energy group boundaries 11 NO fission source spectrum 3 3 NIN3 ID3 LENB 10 LENB burnup chain 14 LENB x3 burnup chain edit control Repeat Blocks 3 and 4 NEB times TABLE D 2 continued Segment Block 10 No 4x 10 11 12 13 14 15 16 17 26 23 24 25 18 19x 20 21 26 Repeat Repeat 5 10 11 12 13 Repeat Repeat Array Length N2 N2 N1 N2 N1 N2 N2 N2 N1 N2 N14N2 NP4 NT NT NT N3 N3 N3 N3 24 95 Data NIN4 ID4 TFR AW IAN NF NT NZZ NP4 P Eo fT FTR condensed fast and res Po scat matrix thermal temperatures lengths of thermal Po scat matrices Block 7 input flags one per temp ITR g a Vae Tf condensed thermal Po scat matrices Block 7 for each I in 25 Blocks 5 7 NES times 5 N2 NTE NSP NTExNSP RIN5S ID5 ISF NTE NSP Block 10 input flags one per res group temperatures increasing only if ID5 0 Sp increasing only if ID5 0 resonance table LS within temperature Block 10 for each 1 in 10 Blocks 8 10 NER times 96 TABLE D 2 concluded Array Segment Block No Length Data 11 6 1 ID6 12 10 NGxNG P scat matrix to within from groups Repeat Blocks 11 and 12 NEP times SCAN 1 1 3 NNUC NWO INEX 2 2 10 data section print flags 3 NNUC nuclide id
21. 218 534 0 90396 52 0 280 0 250 0 030 0 11333 22 148 728 75 5014 73 2266 0 67797 53 0 250 0 220 0 030 0 12783 23 75 5014 48 052 27 4494 0 45187 54 0 220 0 180 0 040 0 20067 24 48 052 27 700 20 352 0 55085 55 0 180 0 140 0 040 0 25131 25 27 700 15 968 11 732 0 55085 56 0 140 0 100 0 040 0 33647 26 15 968 9 877 6 091 0 48038 57 0 100 0 080 0 020 0 22314 27 9 877 4 00 5 877 0 90391 58 0 080 0 067 0 013 0 17733 59 0 067 0 058 0 009 0 14425 60 0 058 0 050 0 008 0 14842 61 0 050 0 042 0 008 0 17435 62 0 042 0 035 0 007 0 18232 63 0 035 0 030 0 005 0 15415 64 0 030 0 025 0 005 0 18232 65 0 025 0 020 0 005 0 22314 66 0 020 0 015 0 005 0 28768 67 0 015 0 010 0 005 0 40547 0 010 0 69315 0 005 72 TABLE 3 OVERLAY MODULE CONTENTS Module 1 Module 2 Module 3 MAIN HDPRNT READ BURNUP GXWIMS ALOCAT HEADER REVS BURN1 PPGBW CORE ICONV SAVE FIND PLOT CORE2 LIMITS SETCOR FLAGS WXCAL DATE MESAGE SUMARY PISCAT CORE1 FIDAS NWDID WRITE RESON GETAP GMR PRINT SMOOTH UPDATE WIMLIB Module 4 Module 5 INDEX SIGP ADDXS EXPAND MULT SIALUP MATGEN TREAD ALPHA INTEGR NPTRD SIGSET MATRIX VECTOR ANS INTERP RECH2 SIGUP1 MWORD WIMSXS ANUSET INTGRC RESCAL SUAI REAL WIMFMT COMPRS ITPOS RESINT TERP1 DIRINT LAMDA RESOUT THNFIT EMESH LINEAR RPLCE ZWEIG1 Module 6 Module 7 Module 8 Module 9 CHECK PRINT1 BCDBIN 10 PRINT2 BINBCD SCAN TAPMAN CONVRT SCAN1 WIAPE1 CONV1 WEP WTAPE2 WIMSCN Common Block Calls other s
22. 8 1 and 3 3238 2 323802 T 3238 2 resonance tables Note that transfer of burnup chain and smooth data is suppressed for the second transfer by using the negative value for NIN 3238 4 File requirements WIMS tape input WIMSTAR tape output 3 3 5 TAPMAN Block 1 1 Run Parameters 5 1 IDTAPE Identification no of the new WIMSTAR tape 2 NWD Unit no of the new WIMSTAR tape Enter 0 to suppress transfer 3 NTR No of entries in the record transfer list Enter 0 to suppress transfer 4 NTAPE No of input WIMSTAR tapes 5 IOUT Data print flag 1 no transfer or listings produce only a summary of the contents of WIMSTAR tapes 0 no listings list only data written to IDTAPE 2 list all data read Terminate Block 1 48 Block 2 2 List of unit nos for the input WIMSTAR tapes NTAPE 3 List of record iransfer flags 0 1 no transfer transfer data record Required only if NTR gt 0 NTR T Terminate Block 2 End of Data Notes 1 The energy structure NG Nl N2 N3 of all WIMSTAR tapes processed together must be identical Unless IOUT 1 in which case no processing is done 24 The length of the 3 array should reflect the total number of records from all input WIMSTAR tapes The list should consist of a set of O s and 1 s 0 No transfer 1 transfer to IDTAPE which are applied successively to the input records If the 3 list is exhausted transfer stops unused e
23. B describes the FIDO input system used by all AMPX modules Appendix C contains the input for a sample case and Appendix D provides a convenient summary of the input data and the required tape units 2 PROGRAM DESCRIPTIONS 2 1 GENERAL REMARKS The generation of WIMS library data is a multi step process requiring the execution of three AMPX II modules NPTXS XLACS2 and WIMSTAR 4 Two utility modules AJAX and RADE provide useful support and have been included in this discussion The procedure is illustrated in Figure 1 and briefly explained here Detailed discussions of each module follow Each nuclide on the WIMS data tape has three sets of informa tion The first is a record describing all the nuclides produced by burnup of this nuclide This data must be provided to WIMSTAR by the user The second set contains the basic cross sections at one tempera ture for the fast and resonance regions and at one or several tempera tures for the thermal region A resonance nuclide requires a third set of data containing absorption and fission yield if fissile cross sections for the resonance region These cross sections are temperature and concentration dependent T os and replace the resonance region cross section of the second set in the WIMS calculations following evaluation of the case data The creation of the resonance tables requires a point cross section curve in the resolved resonance region and a curve for each in the unresolved
24. C 1 Sample Case Input 90 WIMLIB UPDATE TH232 0 0 1 1 18 17 25 1 1 1 0 2 94 69 27 14 13 42 13 37 T 8 2001 6001 2002 8002 3 4 6 7 9 10 1010 11 12 1012 2012 2212 14 16 19 23 27 29 52 55 56 1056 58 63 2063 91 93 112 83 95 99 101 1103 103 105 1105 108 109 113 115 127 131 133 134 135 1135 143 145 147 1147 2147 148 1148 149 150 151 152 153 154 155 1155 157 164 902 176 178 181 207 232 1232 2232 9233 1233 234 235 1235 236 237 3238 939 3239 1240 241 242 941 942 943 1000 2000 1999 9 232 1 1232 1 1232 2 1232 3 2232 1 9233 9233 1233 234 235 2 235 2 235 3 235 3 235 4 235 4 1235 2 1235 2 1235 3 1235 3 1235 4 1235 4 236 3238 5 3239 1 3239 1 T 3 2232 23201 OT T 4 2232 23202 300 232 0381 90 2110T 23 300 T 5 2232 1 23202 100T T SCAN WIMS TAPE FOR ERRORS 1 0 18 OT 2 FO T CONVRT BINARY TO BCD 1 1 18 50 T FIGURE C 1 concluded 91 APPENDIX D INPUT SUMMARIES This appendix contains summaries of the input requirements discussed in Section 3 It is intended as a quick reference after the user becomes familiar with the complete instructions 92 TABLE D 1 AMPX II INPUT SUMMARY Module Block ST Data 1 o 1 broadening method 1 1 NNUC 2 2 6 MATNO NDFB IDTAP MODE NSIGP ID19 3 4 RFACT SFACT SIGP TDEGIC 3 4k NSIGP 5 values only if NSIGP gt 1 Repeat Blocks 2 and 3 NNUC times 1 Five title cards 2 1 5 ID NNUC MAXG NEG IW 2 12 LSLAB LCYL LUNR MSN NPE N
25. ENDF B by G B Wilkin ABSTRACT WIMSTAR Version 4 is a FORTRAN IV computer program developed to generate data files for the WIMS lattice code library from the ENDF B data base The program must be used in conjunction with the AMPX II system and has been designed for implementation as a module of that system This report describes the structure implementation and use of the AMPX WIMSTAR system Atomic Energy of Canada Limited Whiteshell Nuclear Research Establishment Pirawa Manitoba ROE 1L0 1981 August AECL 6809 CONTENTS INTRODUCTION PROGRAM DESCRIPTIONS 2 1 GENERAL REMARKS 2 2 AMPX MODULES 2 3 PROGRAM WIMSTAR INSTRUCTIONS FOR USING THE PROGRAMS 3 1 INPUT DECK SETUP 3 2 AMPX MODULE INPUT 3 3 WIMSTAR INPUT 3 4 JOB CONTROL CONSIDERATIONS 3 5 DYNAMIC CORE ALLOCATION 3 6 ERROR PROCESSING PROGRAM MAINTENANCE INFORMATION 4 1 OVERLAY STRUCTURE 4 2 PROGRAMMING DETAILS 4 3 DATA FILE FORMATS SUMMARY ACKNOWLEDGEMENTS REFERENCES TABLES FIGURES Page NP amp 12 13 14 31 49 50 52 52 r 55 66 67 68 70 73 cont APPENDIX A APPENDIX B APPENDIX C APPENDIX D CONTENTS concluded WIMSTAR ERROR MESSAGES FIDO INPUT SYSTEM SAMPLE CASE INPUT SUMMARTES 80 88 91 1 INTRODUCTION WIMSTAR Version 4 is a FORTRAN IV compucer program developed to generate data files for the WIMS lattice code library tape from the ENDF B 2 nuclear data
26. IMSTAR data tape these arrays are required only if changes are necessary Standard FIDO commands are used to access selected values In some cases the length of these arrays is entered on cards These values are only used if the complete array is input via cards If input is via old WIMS tape or WIMSTAR tape the length cannot be changed via card input 10 Energy group boundaries NG 1 11 Fission source spectrum Sums to 1 0 NO T Terminate Block 2 Block 3 NEB gt 0 3 Burnup Chain Edit Control 3 1 NIN3 NIN of nuclide to be edited Order must match NIN list of 8 array 2 ID3 Data source lt 0 old WIMS tape T 39 0 card input gt O WIMSTAR tape ID3 is identification no of the data 3 LENB Length of burnup chain Required only if ID3 0 If ID3 0 LENB is used to indicate the negative of the no of chain edits entered in 14 For example if ID3 lt 0 and 3 edits are required enter LENB 3 Terminate Block 3 Block 4 NEB gt 0 10 14 Burnup chain for NIN3 The first two words of the chain are replaced by LENB and NIN3 respectively and should be entered as zero LENB Burnup chain edit control A set of 3 entries is input for each edit operation as follows 134 1 2 2 2135 NIN2135 is added for the chain after NIN134 with yield 1 2E 2 134 0 0 0 NIN134 is deleted The length of the chain is adjusted internally The first eight entries in the chain canno
27. L length 4 N3 Data Record Tir Sar Wes Of each of length N3 a PO Scattering Matrix Fast and Resonance Groups Control Record de IDR 5 2 IDN nuclide identification no 36 T temperature 4 9 zero 10 L length N 1 Data Record N D I I 1 N condensed scattering matrix 6 PO Scattering Matrix Thermal Groups The format is identical to record type 5 except that IDR 6 7 Resonance Table Control Record 1 IDR 7 2 IDN nuclide identification no 3 zero 4 ISE type 1 2 absorption fission yiela 5 NTE no of temperatures 6 NSP no of values 64 7 9 zero 10 L length NTE NSP NTE NSP N2 Data Record T L I 1 NTE S I 1 1 NSP D 1 J K I 1 NSP J 1 NTE K 1 N2 resonance table data 8 Pl Scattering Matrix Control Record 1 JDR 8 2 IDN nuclide identification no 3 T temperature 4 9 zero 10 L length NG NG Data Record D J 1 J 1 NG I 1 NG Pl scattering data sequenced to within from groups System EOF mark ends the tape data 4 3 3 PXS ID Tape Format This tape contains infinitely dilute point cross sections for the resolved and unresolved resonance regions Record 1 MAT MF MT ZA AWR zero LFS zero zero Record 2 MAT MF MI T SIGP zero zero N1 N2 NBT I INT I I 1 N1 X I Y I 1 N2 r Record 3 MAT MF 7 zeros Records 1 2 and 3 are repeated for each MT of MAT Record 4 MAT 8 zeros
28. LAS 20500 FH AECL 6809 ATOMIC ENERGY L ENERGIE ATOMIQUE OF CANADA LIMITED OU CANADA LIMITEE WIMSTAR 4 A COMPUTER PROGRAM FOR GENERATING WIMS LIBRARY DATA FROM ENDF B WIMSTAR 4 UN PROGRAMME D ORDINATEUR DESTINE A PRODUIRE DES DONNEES DE BIBLIOTHEQUE WIMS A PARTIR DE ENDF B G B Wilkin Whiteshell Nuclear Research Etablissement de Recherches Establishment Nucl aires de Whiteshell Pinawa Manitoba ROE 1L0 August 1981 ao t ATOMIC ENERGY OF CANADA LIMITED WIMSTAR 4 A COMPUTER PROGRAM FOR GENERATING WIMS LIBRARY DATA FROM ENDF B by G B Wilkin Whiteshell Nuclear Research Establishment Pinawa Manitoba ROE 1L0 1981 August AECL 6809 WIMSTAR 4 UN PROGRAMME D ORDINATEUR DESTINE A PRODUIRE DES DONNEES DE BIBLIOTHEQUE WIMS A PARTIR DE ENDF B par G B Vilkin RESUME WIMSTAR Version 4 est un programme d ordinateur FORTRAN IV tabli pour produire des fichiers de donn es pour la biblioth que de codes de r seaux WIMS partir de la base de donn es ENDF B On doit utiliser le programme concurrement avec le syst me AMPX II Il a t con u pour tre mis en application comme module de ce syst me Ce rapport d crit la structure la mise en application et l emploi du syst me AMPX WIMSTAR L Energie Atomique du Canada Limit e Etablissement de Recherches Nucl aires de Whiteshell Pinawa Manitoba ROE 1L0 1981 ao t AECL 6809 WIMSTAR 4 A COMPUTER PROGRAM FOR GENERATING WIMS LIBRARY DAA FROM
29. M has the same effect as N except that the sign of each entry in the sequence is reversed each time the sequence is entered For example the entries l 2 3 2M2 would be equivalent to 123 3 2 2 3 This option is also useful in entering quadrature coefficients 2 causes N N locations to be set to 0 The pointer is 1 3 advanced by Ny N gt C causes the position of the last array item entered to be printed This is the position of the pointer less 1 The pointer is not moved o causes the print trigger to be changed The trigger is originally off Successive 0 fields turn it on and off alternately When the trigger is on each card image is listed as it is read gt indicates that the pointer is to skip Ny positions leaving those array positions unchanged If the third subfield is blank the pointer is advanced by N If the third subfield is nonblank that data entry is entered following the skip and the pointer is advanced by Ny 1 A moves the pointer to the position Na specified in the third subfield 85 P fills the remainder of the array with the datum entered in the third subfield E skips over the remainder of the array The array length criterion is always satisfled by an E no matter how many entries have been specified No more entries to an array may be given following an E except that data entry may be restarted with an A The reading of data to an array is t
30. MSTAR INPUT Module WIMSTAR is requested via the WIMSTAR Each segment of WIMSTAR is initiated via a segment request card which has the segment name punched starting in column 1 The rest of the card is read and printed as a title for the segment listing For example 2 39 GXWIMS CALCULATE TH 232 WIMS DATA AT 300 K initiates segment GXWIMS The FIDO input data cards follow each segment request card See the sample case in Appendix C 3 3 1 GXWIMS Block 1 1 Run Counter 1 1 NNUK No of nuclides to be processed T Terminate Block 1 Sets of Blocks 2 though 6 are repeated NNUK times Block 2 2 Run Parameters 16 1 NDF ENDF B tape unit no 2 NWD WIMSTAR tape unit no 3 IOUT Output print flag 0 partial print including all WIMS data 1 full print of NAM contents 4 IDTAPE Identification no of the new WIMSTAR tape If IDTAPE 0 an existing tape is assumed and the data is added to the end 5 MAT ENDF B material no 6 IDN WIMSTAR nuclide identification no 7 NTE No of temperatures to be processed maximum of 10 8 NSP No of op values in the resonance calculation NIE NSP must be lt 20 33 9 NG No of groups 69 10 N1 No of fast groups 14 11 N2 No of resonance groups 13 12 N3 No of thermal groups 42 13 IFL Smooth data calculation flag 0 no yes but don t perform the calculation yes including calc
31. PEP IDTAP MODE NGMA LCSM MCSM NDNP 3 6 output flags 3 4 IDTAP ENDF MT identifiers for punched cross sections 5 2xLCSM interpolation scheme for weighting spectrum 6 2xMCSM weighting spectrum 7 MAXG 1 energy group boundaries decreasing 8 10 T X THETA FCUT 6 entries not used One title card 4 70 9 1D19 MATNO NTEMP LORDER NL NFY MATID KMXB KMXA 71 9 SIGP AJIN RFACT SFACT MATPT NUNIT MME MMI MATEL 5 73 NTEMP temperatures increasing Repeat title card and Blocks 4 and 5 NNUC times AJAX 1 o 2 MMT NMAX 1 1 NFILE 93 TABLE D 1 concluded 2 2 NF IOPT 3 1oPT list of nuclide identifiers added or deleted 4 IOPT list of new nuclide identifiers only if change required Enter Block 3 only if IOPT 0 Repeat Blocks 2 and 3 NFILE times 1 4 MMT 3 entries not used 2 IOPT1 IOPT2 18 entries not used 94 TABLE D 2 WIMSTAR INPUT SUMMARY NNUK 2 16 NDF NWD IOUT IDTAPE MAT IDN NTE NSP NG N1 N2 N3 IFL IFR IFP NWK 3 5 DUMAX SIGPA AWA APA ERR 4 NTE temperatures increasing 5 NSP o increasing 6 8 ID19 NAM NXS NSD ITFP ITFT ITFO ITF1 Repeat Block 4 NTE times 7 N2 8 NC 1 energy group boundaries only if no AMPX tapes Input Block 5 only if IFP 1 x 9 NGxN1 TR weights Repeat Block 6 for each ITFT 3 or 5 Repeat Blocks 2 6 NNUK times o 1 IOUT 1 8 NWO NWN NWD LRIN NEB NES NER NEP
32. T Note Sets with duplicate identifiers will not be entered on MMT The first occurrence of an identifier selects that set for the new library T Terminate Block 2 Block 3 Enter only if IOPT 0 3 List of nuclide identifiers to be added or deleted from NF IOPT 4 List of new identifiers This array allows changing of the iden tifiers given in 3 array when selecting nuclides for new library Enter only if identifiers are to be changed IOPT T Terminate Block 3 End of Data Note 1 File requirements AMPX master library tape output AMPX master library tape s input scratch units 15 16 18 19 3 2 4 RADE Block 1 1 Checking Commands 4 1 MMT AMPX master library unit no 2 MWT 3 MAN not required for WIMSTAR applications enter zeroes 4 IFM 31 2 Options 20 1 IOPT1 no of angles at which a display of differential cross sections is desired These angles will be equally spaced in the cosine range l to 1 2 IOPT2 the accuracy in 1 1000 s of a percent to which checks are made e g 1 0 001 1 3 IOPT3 not used enter zeroes 20 IOPT20 T Terminate Block 1 End of Data Notes 1 RADE can also check ANISN libraries Although this option is not required for the WIMSTAR application the user can refer to the AMPX II user s manual 9 for details 2 File requirements AMPX master library tape input 3 3 module request card scratch units 18 19 WI
33. anged by reading in a new set of data for that array If no entries to the arrays in a block are required the delimiter alone satisfies the input requirement Three major types of input are available fixed field input free field input and user field input 81 Fixed Field Input Each card is divided into six 12 column data fields each of which is divided into three subfields The following sketch illustrates a typical data field The three subfields always comprise 2 1 and 9 columns respectively Subfield 1 Subfield 2 Subfield 3 To begin the first array of a block an array originator field is placed in any field on a card Subfield 1 An integer array identifier lt 100 specifying the data array to be read in Subfield 2 An array type indicator S if the array is integer data sx Gf the array is real data Subfield 3 Blank Data are then placed in succ ssive fields until the required number of entries has been accounted for In entering data it is convenient to think of an index or pointer which is under control of the user and which specifies the position in the array into which the next data entry is to go The 82 pointer is always positioned at array location 1 by entering the array originator field The pointer subsequently moves according to the data operator chosen Blank fields are a special case in that they do not cause any data modification and do not move the poi
34. ata to be processed fast data For multithermal groups NEG gt 1 NTEMP is the number of temperatures at which a thermal scattering kernel will be calculated Doppler broadening is included in the resonance treatments but only at the first temperature specified in the 73 array To Doppler broaden a resonance material in a one thermal group calculation at a temperature other than the default value of 300K set NTEMP 1 and input the temperature in the 73 array Caution setting NTEMP O in a multithermal group calculation results in a temperature of 0 K being used and provides no downscatter data in the thermal range 0 Order of expansion for the scattering matrices above thermal Elastic scattering and all inelastic levels MT 51 90 will have matrices of this order 0 If NEG gt 1 enter the order of expansion for the thermal scattering matrices 0 If NEG gt 1 and if thermal ENDF data are to be mounted on KMXA identify the format of the data as follows 0 Binary formatted data 2 BCD formatted data The default value is zero The meaning of this parameter depends on the value of NEG and the presence of thermal data on KMXA as follows SUD Thermal Data NEG on KMXA MATID lt 1 N A Enter zero gt 1 Yes Enter identification number MAT number of thermal data S a 8 data if different from MATNO gt 1 No Parameter serves as a trigger for an analytic free gas calculation The free gas kern
35. atrices Enter 0 for each matrix not required 4 Terminate Block 1 Block 2 Repeat NNUC times 3 Nuclide Transfer Control 3 1 NIT NIN from WIMS tape to be transferred If NIT lt 0 transfer is suppressed 2 IRN Resonance table indicator to be transferred RIN is generated as RIN NIT desl If IRN lt 0 transfer is suppressed 3 IDN WIMSTAR identification no of data from NIT T Terminate Block 2 End of Data Notes 1 If IDTAPE 0 the energy group structure of WIMSTAR tape is checked against that of the WIMS tape then the tape is spaced to the end and the new data is added If IDTAPE gt 0 a new tape is started by placing the energy structure from the WIMS tape on the WIMSTAR tape When checking the energy structure the program ignores the upper and lower values However if a lower value of zero is placed on the WIMSTAR tape segment GXWIMS will fail when attempting to calculate the lethargy mesh for RESPU To avoid this problem do not start a new WIMSTAR tape using TRANS if the WIMS tape has a zero as the lower energy limit T 47 2 If no Pl matrices are required the 2 array is omitted 3 The WIMSTAR identification nos IDN must be unique if one is duplicated the data of the second occurrence cannot be used by the update segment Thus if two sets of resonance tables are required for the same NIN the 3 array is repeated changing IDN For example 3 3238 1 323801 T transfer 323
36. atrix for scattering from thermal groups The data are stored as described for record 4 After the thermal data for all temperatures specified in record 5 data for this nuclide are terminated This concludes the description of the L nuclide files 2 to 61 Resonance Data Files The resonance data are contained in N2 files one for each resonance group Each file contains one absorption or two absorption and fission yield records for each resonance tabulation depending on the trigger NF in the nuclide data of the nuclide to which the tabulation refers The records are in the order of the NIN s in the index Record 2 General Index File The records contain RIN M1 M2 T JB JB 1 M1 Resonance identification number SIGP JD JD 1 M2 RSIG JD JB JD 1 M2 JB 1 M1 temperatures Ces and resonance number of temperatures and oS cross sections ordered a within temperature Note that the upper value of Le in the tabulation is ignored and replaced by infinity when WIMS interpolates in the RSIG table Following the records for each resonance identification number the data for each group are terminated by One record 0 0 1 1 0 0 0 0 0 0 File Mark Pl Scattering Matrices The Pi scattering matrix data consist of a single Py matrix for hydrogen deuterium oxygen and carbon in that order The matrices are written to tape one row at a time thus the Py scattering matrix file consists of 4N r
37. but if no entry for a field is required no space for it need be left Only columns 1 72 may be used as with fixed field input The array originator field can begin in any position The array identifiers and type indicators are used as in fixed field input The type indicator is entered twice to designate free field input i e or The blank third subfield required in fixed field input is not required For example 31L indicates that array 31 a real data array will follow in free field format Data fields may follow the array origin field immediately The data field entries are identical to the fixed field entries with the following restrictions 1 Any number of blanks may separate fields but at least one blank must follow a third subfield entry if one is used 2 If both first and second subfield entries are used no blanks may separate them i e 24S but not 24 S 3 Numbers written with exponents must not have imbedded blanks i e 1 0E 4 1 0E4 1 0 4 or even 1 4 but not 1 0 E4 4 In third subfield data entries only nine digits including the decimal but not including the exponent field can be used i e 123456 89E07 but not 123456 789E07 87 5 The Z entry must be of the form 738Z not Z738 or 738 Z 6 The or data operators are not needed and are not available 7 The Q N and M entries are restricted 3Q4 1N4 or M4 but not 4Q 4N ox 4M User Field Input If
38. ch causes XLACS2 to calculate an MME based on the mass of the nuclide is generally adequate Number of angles for the Lobatto quadrature used in the calculation of the inelastic scattering transfer array The default value 0 which causes XLACS2 to calculate an MMI based on the mass of the nuclide is generally adequate Elastic scattering processing trigger for resonance nuclides T 2 28 0 Prepare elastic scattering data from the ENDF B data on logical unit NoMA MAT Use the MT 2 data for data set MAT on the point cross section library logical unit 31 to prepare the elastic scattering cross sections Terminate Block 4 Block 5 Required only if NTEMP gt 0 73 Temperature s at which to evaluate the thermal scattering kernel NIEMP Temperature s in K is are input low to high Only one temperature should be input for the WIMSTAR application T Terminate Block 5 End of Data Notes 1 XLACS2 is run once for each temperature desired changing the value of ID19 70 array each time Use the same values of ID19 that were used in NPTXS NTEMP 70 array should be set to 1 and the one temperature value entered in the 73 array XLACS2 can be rerun several times by setting NNUC 1 array to the number desired and changing the temperature for each rerun For a resonance nuclide the XLACS2 resolved resonance processing produces average elastic scattering data in 81 equal lethargy panels per energy
39. chemes for the arbitrary weighting spectrum i e NBT i INT i i 1 LCSM 2LCSM Standard ENDF B conventions for the meaning of NBT and INT must be followed see Reference 2 If IW 4 enter the arbitrary weighting spectrum i e E i W i i 1 MCSM 2 MCSM E i units are eV s W i units are per unit energy Entries should be low to high in energy The relation between 5 and 6 arrays is described in Reference 2 The spec trum used by Winfrith is shown in Table l Energy Group Boundaries MAXG 1 The group boundaries are input high to low in eV Many standard group boundaries are available in a library that is available with the AMPX package the built in group structure library When a standard group structure is desired the 7 array can be omitted The group structure library received with the AMPX II package has been modified for the present implementation to include the WIMS 69 group structure illustrated in Table 2 Weighting Option Specifications 10 The 8 array is obtained by default only if IW 1 2 4 or 6 The array is used to specify the energies where the Maxwellian and or fission spectra are joined to the intermediate energy range weighting spectrum and to specify the temperatures of the Maxwellian and fission spectra 1 T temperature in K for Maxwellian spectrum 300 2 x a multiplier on kT where a join to next portion of weighting spectrum and the Maxwellian is made 5 Note
40. ctrum can be specified for the entire problem energy range by using the appropri ate entries in the 8 array to place the Maxwellian spectrum below the problem energy range and or to place the fission spectrum above the problem energy range 18 section library and the identification number for the point data set to be used are specified in the 71 array 2 Secondary Problem Information 12 This array is required only if one or more of its elements is different from the default value shown in parentheses 1 2 LSLAB LCYL 3 Not used 0 The maximum number of points contained in any of the ENDF B File 3 arrays 5 250 If the MATPT or MATEL entries in the 71 array specify that NPTXS produced point cross section data are to be used arrays dimensioned by LCYL may need to be increased to accommodate the point data Consequently LCYL must be the maximum of the following 1 5 250 2 number of points in the weighting function specified by IW 250 3 maximum number of points from MATPT or MATEL data set for MT 1 2 18 or 102 data 250 The number of points in the weighting function specified by IW is as follows IW Number of Points 0 1 or 4 5 or 6 LUNR 2 Assume 2 000 3 MCSM 2 array 7 Number of points in MATPT data set for MI 1 Maximum number of points in any of the arbitrary weighting functions on logical unit 46 The dimension of an array used in conjunction with cr
41. data entry is to be repeated N times The pointer advances by N I indicates linear interpolation The data numerator Ny indicates the number of int erpolated points to be supplied The data entry in the third subfield is entered followed by Ny interpolated entries equally spaced between that value and the data entry found in the third subfield of the next nonblank field The pointer is advanced by Ny 1 The field following an I field is then processed normally according to i 3 own data operator The I entry is especially valuable for specifying a spatial mesh In arrays interpolated values will be rounded to the nearest integer L indicates logarithmic interpolation The effect is the same as that of I except that the resulting data are evenly separated in log space This is especially convenient for specifying an energy mesh Q is used to repeat sequences of numbers The length of the sequence is given by the third subfield N The sequence of Ny entries is to be repeated Ny times The pointer advances by NN If either Ny or Ny is 0 then a sequence of Ny N is repeated one time only and the pointer advances by Ny Ng This feature is especially valuable for geometry specification The N option has the same effect as Q except that the order of the sequence is reversed each time it is entered This is 84 valuable for che type of symmetry possessed by Sa quadrature coeffi cients
42. de only the values of the parameters that change must be input ID19 must change for each temperature 3 IOUT only controls listing the data read from the AMPX master library tape and produces considerable output It should be set to 1 only if problems develop or the data is suspected to be inaccurate 4 For ITFT 4 or 5 only weights read previously for this nuclide can be accessed as the resonance calculation uses the same work tape that the weights are stored on and are thus overwritten Weights are input Gi i 1 N1 j 1 NG where i ranges over the fast to groups and j over all from groups See Section 2 3 1 and Reference 9 5 File requirements AMPX master library tape s input NPTXS point cross section tape s input ENDF B tape input WIMSTAR data tape output scratch units 14 15 16 6 If this nuclide is hydrogen deuterium oxygen or carbon both ITFO and ITF1 should be set to 1 set to O otherwise 3 3 2 WIMLIB Block 1 OS Output Print Flag 1 1 IOUT Controls printing of edited WIMS data 0 1 No Yes The 0 array may also be entered in Blocks 3 5 8 and 11 thus allowing the print to be switched on and off as required 1 Run Parameters 1 NWO LRIN NEB NES NER NEP 8 Old WIMS tape unit no Not required if all input 37 from NWD or cards New WIMS tape wit WIMSTAR Data tape unit Not required if all input from NWO or cards No No No N
43. e checks for inconsistencies in the structure and unreasonable data values and prints user selected files SCAN should be used after every update operation to look for errors and to ensure that the new tape is readable Only the newly added data need be printed although the entire tape is checked internally SCAN performs the following checks duplicate NIN s nuclide identification numbers energy boundaries decrease and all unique Sap fission spectrum sums to 1 0 burnup chains in correct order and all NIN s referenced are present order of nuclide data files check NF NFA and NZZ for consistency count fissiles and fission products and check NNF and NNFP thermal data temperatures increase order of RIN s resonance identification numbers number of tables for each nuclide NZZ NF correct T and g_ should increase and all tables same size for each nuclide structure of condensed scattering matrices SCAN attempts to read the entire tape despite errors encountered See Section 4 3 1 to interpret the printed data and for a description of the WIMS tape format and meaning of the above variables 2 3 4 TRANS Segment TRANS provides a mechanism for transferring data from a WIMS iibrary tape to a WIMSTAR data tape This operation is useful if it is desired to add another set of data for a nuclide already contained on the WIMS tape and the user wishes to use the same burnup chain for both entries An e
44. ecords System EOF Mark Some early WIMS tapes do not contain NNF NNFP or ALAMDA data items WIMSTAR 4 cannot process these tapes 62 4 3 2 WIMSTAR Data Tape Format The WIMSTAR tape consists of pairs of records The first of each pair is a control vector of length 10 that describes the data record that follows it The format of the burnup chain and scattering matrices is the same as on the WIMS library tape 1 Tape Control appears only once as the first pair of records Control Record 1 IDR 1 2 IDTAPE identification no of tape 3 NG no of energy groups 4 Nl no of fast groups 5 N2 no of resonance groups 6 N3 no of thermal groups 7 79 zero 10 L length of data record NG 1 Data Record D I I 1 L energy boundaries 2 Burnup Chain Control Record 1 IDR 2 2 IDN nuclide identification no 3 9 zero 10 L length of data record N Data Record N D I 1 2 N burnup chain 3 Multigroup Data Fast and Resonance Groups Control Record 1 IDR 3 2 IDN nuclide identification no 3 T temperature 4 9 zero 10 L length 4 N1 8 N2 63 Data Record op N2 Eo t N2 9 pr NIAN2 a N1 N2 x N2 ACN2 vo N1 N2 o N1 N2 Note the values in brackets indicate the length of each cross section 4 Multigroup Data Thermal Groups Control Record 1 IDR 4 2 IDN nuclide identification no 3 T temperature 4 9 zero 10
45. el is generated using routines from the tHERMos 19 code The value of MATID depends on the atomic weight of the material being processed as follows Atomic Recommended ile Weight MATID tems lt 19 0 Causes S a 8 free gas data to be generated from which Po order NL matrices are produced gt 19 with no resonances I A Po analytic free gas in thermal range kernel is calculated gt 19 with resonances in 2 A Po analytic free gas thermal range kernel is calculated and normalized to the ENDF B File 3 elastic data 8 KMXB If NEG gt 1 enter the number of atoms per molecule of principal scatterer for which the S a 8 data in File 7 of the ENDF data set apply For example hydrogen is the principal scatterer in water and has two atoms per molecule therefore KMXB 2 e Examples of thermal ENDF data sets are H bound in H20 MAT 1002 D in D20 1004 C in graphite 1065 C in CH 1011 Be 1064 Cele 1095 H in ZrH 1096 and Zr in ZrH 1097 25 9 KMXA If NEG gt 1 enter the logical unit number of the device which contains the thermal data MATID if not the same device as NGMA 71 Resonance Information 9 Needed for a resonance nuclide only 1 SIG o_ the potential scattering cross section in barns per atom of the resonance nuclide for a nuclide or mixture of nuclides which is admixed with the resonance nuclide in a homogeneous system is defined as follows M E No o i j si p N
46. entification no print list TRANS 1 1 6 NWO NWD IDTAPE NNUC NTE NSP 2 4 P scat matrix WIMSTAR identification list 2 3 3 NIT IRN IDN Repeat Block 2 NNUC times e TAPMAN 1 1 5 IDTAPE NWD NTR NTAPE IOUT 2 NTAPE input WIMSTAR tape unit nos 3 NTR record transfer flags CONVRT 1 1 3 ITYPE NWO NWN 97 TABLE D 3 MODULE INPUT OUTPUT REQUIREMENTS Pare WIMSTAR Input Variable Name a 4 se a et et 10 I input O output AMPX master library WIMS library tape card reader printer ENDF B fast PXS FD finite dilution group boundaries tape PXS ID infinite dilution scratch tape ENDF B thermal WIMSTAR data tape m hit oO D gt Enr ISSN 0067 0367 To identify individual documents in the series we have assigned an AECL number to each Please refer to the AECL number when requesting additional copies of this document from Scientific Document Distribution Office Atomic Energy of Canada Limited Chalk River Ontario Canada KOJ 130 Price 5 00 per copy ISSN 0067 0367 Pour identifier les rapports individuels faisant partie de cette s rie nous avons assign un num ro AECL chacun Veuillez faire mention du num ro AECL si vous demandez d autres exemplaires de ce rapport au Service de Distribution des Documents Officiels L Energie Atomique du Canada Limit e Chalk River Ontario Canada KOJ 1J0 prix 5 00 par exemplaire
47. erform these current calculations thus they must be done elsewhere and the Wijs input via cards An alternative option is to assume a 1 E current for all energies above 4 eV Other options include using the row sum correction for all energies or simply setting a a See Reference 9 for a discussion of the merits of these various transport cross section options Finally for all nuclides on the WIMS tape except those for which a Py scattering matrix is supplied i e hydrogen deuterium oxygen and carbon the self scatter ing term of the Po matrix is adjusted as 10 a 0 z ol gt j ii ii wle 1 The user must direct GXWIMS as to whether or not this adjustment is required GXWIMS then prints and plots the one dimensional multigroup data for inspection and writes all results to a WIMSTAR data tape 2 3 2 WIMLIB Segment WIMLIB reads an old WIMS library tape and under the direction of user supplied update and edit commands uses the results on a WIMSTAR tape to write a new WIMS library tape New data can be input entirely from the WIMSTAR tape or entirely from cards or from a combi nation of both A new WIMS tape can be built entirely from scratch if desired Any entry on the input WIMS tape can be edited although changing the energy group structure currently 69 groups is not per mitted Energy group condensation is a facility that may be added in the future 2 3 3 SCAN Segment SCAN reads a WIMS tap
48. erminated when a new array origin field is supplied or when the block is terminated If an incor rect number of positions have been filled an error edit is given and a flag is set which will later abort execution of the problem FIDO then continues with the next array if an array origin was read Other wise it returns zontrol to the calling program A block termination consists of a field having T in the second subfield All entries following T on a card are ignored and control is returned from FIDO to the calling program Comment cards can be entered within a block by placing an apostrophe in column 1 Then columns 2 80 will be listed with column 2 being used for printer carriage control Such cards have no effect on the data array or pointer Free Field Input With free field input data are written without fixed restric tions as to field and subfield size and positioning on the card The options used with fixed field input are available although some are slightly restricted in form In general fewer data cards are required for a problem the interpreting card is easier to read a card listing is more intelligible the cards are easier to keypunch and certain common keypunch errors are tolerated without affecting the problem 86 Data arrays using fixed and free field input can be intermingled at will within a given block The concept of three subfields per field is still applicable to free field input
49. erride any item via card input On option the edited data files are printed 55 WIMSCN is the main WIMS tape scan routine which calls sub routine CHECK entry points OCHK SCHK VCHK and PCHK to check NIN order data sequencing data values and condensed scattering matrix structure respectively Logical function IO controls printing of selected nuclide files Subroutines TRANS3 and TRANS2 handle transfer of data from WIMS to WIMSTAR tapes for resonance tables and all other data respec tively in segment TRANS Subroutine WIAPE2 performs similar functions for segment TAPMAN SUMARY prints a list of the control records for WIMSTAR tapes Segment CONVRT uses subroutine BINBCD to perform binary to BCD conversion and BCDBIN for BCD to binary conversion 4 3 DATA FILE FORMATS 2 for the The user is directed to the ENDF B user s guide format of the ENDF B library and to the AMPX II user s manual 4 for the format of the AMPX master library and the energy group tape The WIMS library format taken from Reference 7 is included here to aid in interpreting the listings produced by the module SCAN 4 3 1 WIMS Library Format Introduction The library tape is composed of a series of files written in the FORTRAN binary mode There is a general index file followed by a file for each nuclide giving the basic cross sections followed by a file for each resonance group giving the temperature and Le effective potential scatter
50. for SIG in XLACS2 array 71 and the definitions of SIGPA and S in GXWIMS arrays 3 and 5 respectively 2 A reasonable range for RFACT is 0 7 to 1 0 If the input value is outside this range NPTXS resets it to the default value 16 3 If NSIGP 1 the desired ds is input in the 3 array and Block 3 array 4 is not required 4 NPTXS is run once for each temperature changing the value of ID19 2 array each time 5 File requirements ENDF library type input PXS ID tape unit 31 output PXS FD tape unit 41 output scratch units 14 15 16 17 18 19 Input instructions for the remaining AMPX modules XLACS2 AJAX and RADE are taken directly from the AMPX II user s manual 4 3 2 2 XLACS2 Block 1 Five title cards each in 20A4 format The five cards of Hollerith information can be used to describe the neutron library being produced The block terminator is not input Block 2 1 General Problem Information 5 1 ID Identification number for the neutron library Any number will suffice This number is carried as a user identifier on the library but is not actually used by any module in AMPX 2 NNUC Number of ENDF materials to be processed 3 MAXG Total number of neutron energy groups 4 NEG Number of thermal neutron energy groups 17 5 IW Weighting option trigger 1 Fission constant Maxwellian 2 Fission 1 E Maxwellian 3 1 04E each material has a
51. gments TRANS TAPMAN and SCAN provide listing and editing services for WIMS and WIMSTAR tapes Segment CONVRT converts WIMS tapes between binary and BCD formats to allow their transfer to other computer centres AMPX MODULES A brief description of the AMPX modules required is given here For a detailed discussion please refer to the AMPX II user s manual t 2 2 1 NPTXS Module NPTXS preprocesses resonance or nonresonance ENDF nuclides to make point files for total fission elastic scattering and n y cross sections Several modifications have been made to NPTXS for use with WIMSTAR thus the user should refer to Section 3 2 of this report and not the AMPX II user s manual when preparing input for the modified version of NPTXS These changes involve writing the unresolved cross sections for the qo S to tape 41 instead of punching them on cards Also the maximum number of s has been increased from 7 to li The resolved resonance region and infinitely dilute unresolved resonance region point cross sections written to tape l PXS ID are used by both XLACS2 the resonance processing of XLACS2 should be skipped completely and segment GXWIMS of module WIMSTAR 4 If the nuclide contains unresolved resonance parameters the data written to tape 41 PXS FD containing unresolved cross sections at several finite dilutions are used only by GXWIMS NPTXS can process both single and multi level Breit Wigner resolved resonance
52. group If MATEL 0 these data are used to prepare the elastic scattering matrices However an alternative more accurate procedure for calculating elastic scattering trans fer matrices is available The alternative procedure is triggered by setting MATEL equal to the identification number of a data set in an NPTXS generated point cross section library on logical unit NUNIT XLACS2 will bypass the use of internally generated elastic scattering data and will use the MT 2 data from the point cross section library 29 2 File requirements ENDF B fast library unit 11 input ENDF B thermal library unit i2 input AMPX master library unit 23 output AMPX group structure library unit 47 input scratch units 14 15 16 17 18 3 2 3 AJAX Block 1 0 Logical Assignments 2 1 MMT Unit no of new library 1 2 NMAX Unit no of the input file which has the largest buffer requirements 1 Number of Files 1 1 NFILE No of file requests to be made When reordering operations are performed which require the same file to be accessed several times each access is counted to determine NFILE T Terminate Block 1 Blocks 2 and 3 are repeated NFILE times Block 2 2 File and Option Selection 2 1 NF Unit no of input library 2 IOPT option N delete N nuclides from NF to create new library on MMT 30 add all nuclides on library MMT io n N add N nuclides from NF to create new library on MM
53. hain and the new multigroup and resonance data The new WIMS tape was then scanned for errors and finally converted to BCD for shipment 89 EXEC AMPX GRGN 640K GTIME 30 U14 amp amp A U15 amp amp 5 U16 amp amp C U17 amp amp D U18 amp amp E U19 amp amp F STEPLIB DD DD NPTXS LOAD MODULE DD XLACS2 LOAD MODULE DD WIMSTAR LOAD MODULE FTO1F001 DD OLD WIMS LIBRARY FT11F001 DD ENDF B TAPE FT23F001 DD AMPX MASTER LIBRARY T31F001 DD POINT X SECTIONS INFINITE DILUTION FTH1F001 DD POINT X SECTIONS FINITE DILUTIONS FT47F001 DD ENERGY GROUP LIBRARY FT50F001 DD NEW WIMS LIBRARY BCD FORMAT SYSIN DD NPTXS 1 1 T 2 1296 A4 2 9 129601 3 A4 300 T 4 146 10500 3208 1662 842 392 162 102 52 T XLACS2 AMPX2 IBM STANDARD VERSION XLACS2 TH232 MAT 1296 USING NPTXS DATA 300K 69 GROUFS FISS 1 E MAX WEIGHTING 1 1 1 69 42 2 2 A2 11000 26000 AB 2 E 3 FIT 8 300 5 1 27 6 67 444 E T TH232 MAT 1296 USING NPTXS DATA 70 129601 1296 1 3 3 A7 2 E 71 146 A5 129601 31 0 0 129601 T 73 300 T WIMSTAR GXWIMS CAL TH232 DATA 1 1T 2 11 17 0 801111 1296 23202 1 9 69 14 13 42 2 3 1 500 T 3 0 11 2 232 0381 0 0 01 4 300 5 52 102 162 392 842 1662 3208 10500 1 6 T 6 129601 23 31411200T 2 FOL2 T TRANS TH232 WIMS TO WIMSTAR 1 1 170 13 10T 3 1232 3 23201 T TAPMAN CHECK LIST TH232 1 A4 1 1 T 2 17 T FIGURE
54. ighting functions The multigroup data produced are written to an AMPX master library tape XLACS2 should be directed through input items MATPT and MATEL of the 71 array see Section 3 2 to skip all resonance processing and to generate multigroup total fission elastic scattering and n y cross sections from data on the point cross section library PXS ID produced by NPTXS XLACS2 contains the facility of substituting ENDF B thermal scattering law S a 8 data in place of normal ENDF B processing in the thermal energy region sze item MATID of the 70 array Section 3 2 This facility should be used for very light moderators XLACS2 should be run once for each temperature desired with the same nuclide identification numbers item ID19 of the 70 array as were used for the NPTXS data 2 2 3 AJAX Module AJAX is used to combine data from AMPX master libraries Options are provided to allow merging adding deleting or reordering for any number of input files AJAX can be used to change the nuclide identification numbers to match the NPTXS data sets if this was not done at the time of generation 2 2 4 RADE Module RADE is provided to check the multigroup libraries for consistent and reasonable data On option the user can request a display of differential cross sections 2 3 PROGRAM WIMSTAR WIMSTAR is organized into six separately executable segments GXWIMS WIMLIB TRANS TAPMAN SCAN and CONVRT Each modu
55. ing dependent cross sections for each resonance element 56 Each nuclide is identified by a nuclide identification number NIN this NIN is in general chosen so that the NIN modulo 1000 is the atomic mass number of the nuclide in question Each resonance tabulation is identified by the NIN of the nuclide to which it refers together with a single decimal point number thus any nuclide can have up to 10 0 1 9 resonance tabulations associated with it these combined give the resonance identification number RIN which is a floating point value as opposed to the integer value NIN The energy groups are counted from high to low energy similarly ail information is tabulated in order of decreasing energy A positive distinction is made between fast resonance and thermal energy groups the cuts usually being taken as 10 keV and 4 eV although these figures are quite arbitrary General Index File D The first file on the tape contains the following general information Record 1 L N NO N1 N2 N3 NNF NNFP Number of nuclides on this library tape Total number of groups N1 N2 N3 Number of groups into which there is a fission source Numbers of fast resonance and thermal groups Numbers of fissile and fission product nuclides Record 2 IN I I 1 L Nuclide identification numbers for the L nuclides on the tape these form an index to the nuclide files Record 3 GB J J 1 N 1 Record 4 FS J
56. le is now described 2 3 1 GXWIMS There are two phases in the execution of GXWIMS calculation of resonance tables including Goldstein Cohen values and calculation of multigroup cross sections The resonance integral tables in the WIMS library result from calculations performed with the resonance nuclide under consideration in a homogeneous mixture with hydrogen as the moderator The table values are a function of both temperature and dilution o gt Module RESPU of the RSYST Code System 17 is used to solve iia ati the slowing down equation for the situation of a resonance nuclide mixed homoge neously with hydrogen using the total and scattering point cross sections provided by NPTXS and yielding the relative flux pu This flux is used to collapse the absorption and fission yield if the nuclide is fissile point cross sections u9 f o u u du u 1 RI u f u du ai into the WIMS multigroup structure This process is repeated for each using the appropriate unresolved resonance region values and each temperature using appropriately broadened values The Intermediate Resonance IR theory of Goldstein and Cohen 18 shows that the resonance integral can be approximated as Aroa paan TR a a 0 04 A6 a 6G 08 te ou du u where Cia is the resonance scattering cross section and is the Goldstein Cohen parameter taking a value between O0 and 1 A value of 1 produces the Narrow Resonance
57. lements are ignored 3 To produce tape listings summaries without writing a new tape a To summarize all input tapes with no lists set IOUT 1 NWD 0 NTR O and omit 3 array b To list all of one or more input tapes set IOUT 2 NWD 0 NTR 0 and omit 3 array c To list selected records from one or more tapes set IOUT 1 NWD 0 and use NTR and 3 array to select required records 4 File requirements WIMSTAR data tape s input WIMSTAR data tape output 49 3 3 6 CONVRT Block 1 1 Run Parameters 3 1 ITYPE Conversion direction 1 Binary to BCD 2 BCD to Binary 2 NWO Input WIMS tape unit no 3 NWN Output WIMS tape unit no T Terminate Block 1 End of Data Note 1 File requirements WIMS tape input WIMS tape output 3 4 JOB CONTROL CONSIDERATIONS Each of the five AMPX modules required NPTXS XLACS2 AJAX RADE and WIMSTAR 4 has been saved in a separate load module The AMPX II driver program has also been saved as a load module This driver program reads the module request cards and module data cards initiates execution of each requested module and prints a summary of the run The driver program and requested modules are executed via a cataloged procedure called AMPX2 See the sample case in Appendix C for a listing and use of this procedure The procedure is executed via EXEC AMPX2 GRGN 400K GTIME 10 ff U14 amp amp A UL5 amp amp B 50
58. lowing sections 10 45 ds Burnup chains Eo s 2 Fast and resonance groups ds Te e fir Ta Xe 3 Fast and resonance groups VO es Og 4 Fast and resonance groups P scattering matrix 5 Thermal temperatures T 6 hermal groups Oey a 7 Thermal groups VO gs Og 8 Thermal groups Po scattering matrix 9 Resonance tables 10 Pi scattering matrices 3 NLIST List of nuclides to be printed This list is used together with the values of INEX and IOUT to determine whether a given section of data is to be printed NNUC T Terminate Block 2 End of Data Notes 1 SCAN also checks the WIMS tape for errors in structure and unreason able data values This checking is performed on the entire tape regardless of whether the data is printed 2 After an update operation only the nuclides added or edited need be printed Use NNUC and 3 array to select these nuclides omit the 2 array and set INEX 0 3 File requirements 3 3 4 Block 1 WIMS tape input TRANS 1 Run Parameters 6 l 2 NWO Input WIMS tape unit no NWD Output WIMSTAR tape unit no 2 SAGE Identification no of new WIMSTAR tape enter O if 3 IDTAPE an existing tape is to be extended 4 NNUC No of nuclides to be transferred 5 NTE Max no of temperatures in resonance tables 6 NSP Max no of a S in resonance tables List of WIMSTAR identification nos for the four P1 scattering e m
59. n tape PXS ID created by NPTXS module of AMPX Only if IFR gt 0 Tape unit no of tape created by NPTXS module of AMPX containing unresolved point cross sections for each op in 5 array PXS FD Only if IFR gt 2 Data plot flags 1 Q no plotting 1 plot WIMS data 2 plot WIMS data and point curves for On Cu a u for resonance calculation T 7 Bx 35 3 same as 2 including u 6 ITFT Transport cross section calculation 2 0 use Sy 1 use row sum method 2 3 4 5 use row column sum method with weight factors lethargy widths read from 9 array previous weights read previous weights read with changes from 9 array 7 ITFO Po self scattering term adjustment 0 yes i no 8 ITF1 P output to NWD 0 no 1 yes Terminate Block 4 Block 5 Aux Input Only if IFP 1 List of values to be used if A calculation is not performed or fails to generate a usable value N2 Energy structure high to low eV Used only if IFL 0 and no AMPX master libraries are used NG 1 Terminate Block 5 Block 6 Repeat for each occurrence of ITFT 3 or 5 gx List of weighting factors for calculating er for the fast groups ordered to 1 gt N1 within from 1 gt NG NGX N1 Terminate Block 6 End of Data 36 Notes 1 Sets of Blocks 2 through 6 are repeated NNUK times 2 For each occurrence of Block 4 after the first for any given nucli
60. nter A data field has the following form Subfield 1 The data numerator an integer lt 100 We refer to this entry as Ny in the following discussion Subfield 2 One of the special data operators listed below Subfield 3 A nine character data entry to be read in F9 0 format It will be converted to an integer if the array is a array or if a special array operator such as Q is being used Note that an exponent is permissible but not required Likewise a decimal is permissible but not required If no decimal is supplied it is assumed to be immediately to the left of the exponent if any and otherwise to the right of the last column This entry is referred to as N in the following discussion 3 A list of data operators and their effect on the array being input follows Blank indicates a single entry of data The data entry in the third subfield is entered in the location indicated by the pointer and the pointer is advanced by one However an entirely blank field is ignored 4 or indicates exponentiation The data entry in the third field is entered and multiplied by 10 Ni where Ny is the data numerator in the first subfield with the sign indicated by the data 83 operator itself The pointer advances by one In cases where an exponent is needed this option allows the entering of more significant figures than the blank option a has the same effect as R indicates that the
61. o No of of of of of 2 New WIMS Tape Control 1 2 of of of of of of of of resonance tables on new WIMS tape burnup chain edits nuclide data file edits resonance table edits Pl matrix edits Must be either 0 or 4 Record 8 nuclides energy groups NG N1 N2 N3 69 groups into which there is fission source 27 fast groups 14 resonance groups 13 thermal groups 42 fissile nuclides fission product nuclides N2 N3 are included only for completeness of edit LNIN No NG No NO No N1 No N2 No N3 No NNF No NNFP No NG N1 facilities the program cannot change the group structure These values must be identical to the old WIMS tape values T Terminate Block 1 38 Block 2 8 List of nuclide identification numbers NIN s to be placed on new WIMS tape LNIN 9 List of resonance table identification numbers RIN s to be placed on new WIMS tape LRIN Note 1 Arrays 8 and 9 define the data to be placed on the new WIMS tape Any data on the old tape but not in these lists will be deleted RIN s must be in same order as NIN s which must be in the same order as the old tape New NIN s can be placed in any order on the new tape Note 2 In the following description arrays marked with an underscore e g 10 are used for both data input and data edit If the data is input from an old WIMS tape or a W
62. ocks 5 6 and 7 are repeated NES times Block 8 NER gt 0 5 Resonance Table Edit Control 5 1 RINS RIN of nuclide to be edited Order must match RIN list of 9 array 2 ID5 Data source lt 0 old WIMS tape 0 card input gt O WIMSTAR tape ID5 is idendification no of data 3 ISF Cross section type 1 Absorption 2 Fission yield 1 4 NTE No of temperatures Required only if ID5 0 5 NSP No of p values Required only if ID5 0 T Terminate Block 8 Note Each resonance file contains one absorption or two absorption and fission yield records for each resonance tabulation depending on NF in the nuclide data file of the nuclide to which the tabulation refers If a nuclide contains a fission yield tabulation and is to be edited two sets of Blocks 8 9 and 10 are required the first edits absorption ISF 1 and the second edits fission yield ISF 2 both with the same value of RIN5 Each set counts towards the total NER A Block 9 NER gt 0 10 11 12 Group edit flags N2 entries 0 no edit 1 edit corresponding to each of the N2 resonance energy groups Block 10 is input for each 1 in 10 array N2 List of temperatures increasing order Required only if ID5 0 NTE List of LE values increasing order Required only if ID5 0 NSP Terminate Block 9 Block 10 Repeat for each 1 in 10 array and NER gt O 13 T Resonance cross sections
63. of thermal temperatures 20 exceeded 3 array read error NTE exceeded NSP exceeded WIMS tape problem fission flag NF not consistent with resonance tables WIMS tape problem resonance tables not complete Available core exhausted 8 9 10 11 array read error Group boundaries on WIMSTAR tape don t match old WIMS tape See footnote on page 79 79 Routine Error Description WIMLIB 1 0 1 2 array read error 2 3 4 5 Input energy structure does not match old WIMS tape 2 3 4 5 Total Fast Resonance Thermal 6 Program error WIMSTAR tape structure invalid 7 8 9 10 Input energy structure does not match WIMSTAR tape 7 8 9 10 Total Fast Resonance Thermal WIMSXS 1 2 4 5 6 Available core exhausted 3 Insufficient data on AMPX Master Library tape 7 Energy boundaries don t match 8 Can t find requested nuclide on AMPX Master Library tape 9 9 array read error WIMFMT 1 Available core exhausted WIAPE2 1 4 Energy boundaries don t match on all input WIMSTAR tapes 2 Available core exhausted Program error WIMSTAR tape structure invalid WXCAL 1 2 Available core exhausted 3 4 Energy group structure input does not match AMPX Master Library tape 5 7 8 array read error Increase the REGION size on the job step EXEC card and rerun the step Not all uses of dynamic core can be checked beforehand Eg READ 1 K A I 1 1 K If there is not enough space in one of these cases an IHO24
64. omic weight atomic r number trigger see below f number of temperatures at which thermal data tabulated number of resonance tabulations associ ated with this nuclide NF is the fissile and resonance trigger and may take the following values non fissile no resonance tabulation non fissile resonace absorption tabulation 240 fissile tabulation of resonance absorption only e g Pu fissile tabulation of resonance absorption and fission wh FO fissile no resonance tabulation Note that the atomic weight should be exact because it ig used in WIMS for calculating number densities from physical densities Record 2 PSCAT J J 1 N2 Le potential scattering cross section for resonance groups XISS J J 1 N2 Eo lT slowing down power divided by lethargy width for resonance groups TR J J 1 N1 N2 Ter transport cross section for fast and resonace groups Record 3 Record 4 Record 5 59 ABS J J 1 N1 N2 o gt absorption cross section for fast and resonance groups CHI J J 1 N2 x for resonance groups not used beyond WIMSB ALAMDA J J 1 N2 Goldstein Cohen parameter for resonance groups XNUFIS J J 1 N1 N2 vo and o f f FIS J J 1 N1 N2 and fission cross sections fission yield for fast and resonance groups Record 3 is present only for fissile nuclides i e NF gt 2 For nuclides with resonance tabulation s
65. oss section weighting 10 000 If NPTXS produced point cross section data are to be used arrays dimensioned by LUNR are also used to accommodate the point area Consequently 6 8 MSN NPE NPEP IDTAP MODE NGMA 19 LUNR must be the maximum of 10 000 or the following 2 x number of points in weighting function 2 x maximum number of points from point library for MT 1 2 18 or 102 data 700 The number of points in the weighting function is described above The maximum number of interpolation regions needed to describe any ENDF B File 3 array or to describe any TABl1 record produced in XLACS2 250 To date the default value has been adequate for all XLACS2 problems run Identification number of arbitrary function This entry required only if IW 6 0 Logical unit for XLACS2 produced master cross section library 23 Number of reactions for which cross sections are to be punched ENDF B File 3 cross sections averaged by XLACS2 multigroup representation of the weighting function x etc can be punched in free form FIDO format Data to be punched are identified in the 4 array 0 Format of the ENDF library on logical unit NGMA 0 0 Binary formatted library 1 BCD formatted library Logical unit number of the device which contains ENDF B neutron cross section data i e the ENDF library 11 Data for the thermal energy range may or may not be on this
66. p where N is the number density of the jth component of the mixture in atoms barn cm N is the number density of the resonance nuclide in the mixture Sa is the approximate scattering cross section in the resonance region for the jth com ponent in barns atom M is the number of components in the mixture The default value of o is 1 0 x 108 and is the P recommended value for the WIMSTAR application 2 AJIN The 0 j state for which unresolved parameters will be passed to NITAWL The NITAWL unresolved treatment is restricted to one unresolved sequence For the infinite dilution case op This can be approximated by op 1 0 x 108 A table of scattering cross sections is given in Table IV 1 of Reference 20 3 4 5 RFACT SFACT MATPT 5 or 7 5 or 7 26 The most inportant sequence should be identified and passed via this input parameter The resonance with the smallest mean level spacing is usually the most important 21 resolved resonance The r factor in the Nordheim treatment This factor is used to determine the integration mesh spacing lethargy through q py tE 0 where is the mesh size D is the Doppler width at 273 K and Eg is the resonance peak energy A value of 5 0 has been found adequate for most cases 5 0 The s factor in the Ncrdheim resolved resonance treatment This is the number of practical widths on eithe
67. parameters The resonance data can be Doppler broadened using either the numerical integration techniques of Cullen SIGMAL method 12 13 or the conventional PSI CHI method Note that the background ENDF File 3 data will be broadened oniy by the SIGMA1 method making it the preferred method The unresolved processing uses the same techniques as XLACS2 originally developed for the MC 2 code 14 The identification number for the output point cross sections ID19 input by the user must be different for each temperature of each nuclide to allow GXWIMS to differentiate between different temperature sets of the same nuclide NPTXS must be run once for each temperature desired for a given nuclide k The o in this report and in the NPTXS description in the AMPX II user s manual refer to the same quantity 2 2 2 XLACS2 XLACS2 is the AMPX module which caiculates weighted multigroup neutron cross sections from ENDF B data Its calculational routines 15 ona FLANGE 11 16 codes and were originally taken from the SUPERTOG produce full energy range neutron cross section libraries Provisions are included for treating fast resonance and thermal ENDF B data in a single calculation Energy group structure and expansion orders used to represent differential cross sections can be averaged over an arbitrary user supplied weighting function for example the wet and dry spectra used by Winfrith Section 3 2 or over any of several built in we
68. r side of a resonance peak to which the Nordheim treatment is applied A value of 5 0 has proven sufficient 5 0 The meaning of this parameter depends on the value of IW as follows Meaning of MATPT Enter the identification number usually the ENDF MAT number of the point data set on logical unit NUNIT used for the Oy in the weight function Point data sets are usually produced by the NPTXS module To bypass all resonance processing for a resonance nuclide enter the identification number of the point data set on logical unit NUNIT The Top Tes and on y data are taken from the point data set and bd the corresponding ENDF File 2 and File 3 data for 5 or 7 5 or 7 NUNIT MATEL 27 MATNO are ignored The ID19 master data set will contain no resonance parameters Consequently NITAWL can do no further resonance processing For a resonance nuclide enter zero to trigger XLACS2 to do the unresolved resonance processing and to place resonance parameters in the master data set for subsequent resolved resonance processing by NITAWL For a non resonance nuclide set MATPT equal to zero Default value for MATPT is zero 0 If MATPT and or MATEL 0 enter the logical number of the device which contains the point cross section library Otherwise enter zero 0 Number of angles for the Lobatto quadrature used in the calculation of the elastic scattering transfer array The default value 0 whi
69. ram prints the Run Completion Code when a module terminates and begins execution of the next module requested regardless of the value of the previous Run Completion Code The user is warned that modules executed subse quent to an abnormally terminated module will probably not run properly if they require results from the errant module If a system error occurs illegal operation illegal core location reference etc everything grinds to a halt and the AMPX II driver cannot print the module termination message 4 PROGRAM MAINTENANCE INFORMATION This chapter is intended to aid the maintenance programmer in making future modifications to WIMSTAR 4 1 OVERLAY STRUCTURE The recommended overlay structure is illustrated in Figure 2 and the subroutines and common blocks contained in each overlay module are listed in Table 3 The function of each module is described below Module 1 Contains the mainline program input output routines error routines and other service routines required by the rest of the program This module controls segment execution Module 2 segment WIMLIB Module 3 control for segment GXWIMS Module 4 multigroup phase of GXWIMS Module 5 resonance phase of GXWIMS 53 Module 6 segment SCAN Module 7 segment TRANS Module 8 segment TAPMAN Module 9 segment CONVRT 4 2 PROGRAMMING DETAILS The mainline program reads the segment request cards and calls the appropriate segment driver subroutine
70. rary Brookhaven National Laboratory Report BNL 50274 1970 C J Taubman The WIMS 69 Group Library Tape 166259 U K Atomic Energy Authority Winfrith Report AEEW M1324 1975 N M Greene J L Lucius L M Petrie W E Ford III J E White R Q Wright AMPX II Modular Code System for Generating Coupled Multigroup Neutron Gamma Ray Cross Section Libraries from Data in ENDF Format Oak Ridge National Laboratory Report ORNL TM 3706 1978 Revised for AMPX II D E Cullen Program LINEAR Version 77 1 Linearize Data in the Evaluated Nuclear Data File Version B ENDF B Format Lawrence Livermore Laboratory Report UCRL 50400 Vol 17 Part A 1977 R Ruehle RSYST An Integrated Modular System with a Data Basis for Automated Calculation of Nuclear Reactors Oak Ridge National Laboratory Report ORNL TR 2796 1973 J D Macdougali unpublished data Winfrith User Note WIMS 56 Appendix 7 1967 F J Fayers P B Kemshell unpublished data Winfrith User Note WIMS 40 P B Kemshell M Hardcastle unpublished data Winfrith User Note WIMS 73 1969 P B Kemshell Pl Data in the WIMS Library private commu nication 1976 R Kinsey C Dunford What You Should Know About ENDF B Version V Brookhaven National Laboratory Report BNL NCS 22865 R 1977 D E Cullen Program SIGMA1 Version 77 1 Doppler Broaden Evaluated Cross Sections in the Evaluated Nuclear Data File Version B ENDF B Forma
71. reases the GRGN specification and reruns the job WIMSTAR error messages are listed in Appendix A Each segment of WIMSTAR informs the user of the amount of dynamic core available and at termination of the amount of that core that was not used It is suggested that GRGN initially be set to 400 K and subsequently adjusted as necessary NPTXS and XLACS2 use similar core allocation but the size must be set when the load module is created For the present implemen tation the container arrays for NPTXS and XLACS2 were set at 90 000 and 70 000 words respectively and the modules require 640 K and 540 K regions respectively The user is directed to the AMPX II user s manual 4 if these values require adjustment AJAX and RADE use the same system as WIMSTAR 3 6 ERROR PROCESSING The user is directed to the AMPX II user s manual to inter pret error messages from the AMPX modules When WIMSTAR detects an error the following message is printed kKERROR ERROR CODE n IN SUBROUTINE name where name is the subprogram detecting the error and n is the error code number These codes are explained in Appendix A A FORTRAN IH0220I error is then intentionally committed to obtain a subroutine traceback to aid in solving the problem Finally WIMSTAR is terminated with a Run Completion Code of 16 Some routines print additional infor mation and others print self explanatory messages 52 The current version of the AMPX II driver prog
72. s The setup of each segment is essentially the same Variables to be read through the FIDO input System are collected in common blocks and read by calling FIDAS The input output devices must be initialized to assign core space to the buffers Then all remaining available core is assigned to the segment via the call to the assembler routine ALOCAT If a routine to allocate core dynamically is not available at a particular installation the following FORTRAN subroutine may be substituted SUBROUTINE ALOCAT SUB COMMON SPACE D 50000 LIMIT 50000 CALL SUB D LIMIT RETURN END where the length of D is set at compile time The called subroutine then partitions the D vector into smaller parts as required by the particular task Upon entry each segment prints the length of the D vector at termination it prints the amount that was not used during execution of the segment Subroutine GMR with its entry points START END CLOCK ENTER and ERROR prints entry exit timing and error messages 54 The resonance table calculation phase of GXWIMS is performed by RESCAL Four large work spaces W1 W2 W3 W4 two smaller work spaces W5 W6 and three scratch tapes NT1 NT2 NT3 are used in manipulating point data The lethargy mesh is generated using Tr and Sop maintaining AU ax as defined in Section 3 3 1 and eliminating all zero value cross sections The slowing down equation is solved by RECH2 of module RESPU from the RSYST Code
73. s all of the instructions that follow reflect this environment 13 3 1 INPUT DECK SETUP The input data cards required to run each module are grouped in the SYSIN data file as GO SYSIN DD module data cards Each module is initiated by a module request card with an sign punched in column 1 and the module name following starting in column 2 For example to initiate XLACS2 use XLACS2 Data cards specific to the module follow the module request card Any number of sets of request card and data cards can be placed in the SYSIN file to execute a series of modules All AMPX modules including WIMSTAR 4 employ the FIDO input 4 system Appendix B contains a user s guide to this system The input instructions use the following conventions indicates the number of input items expected for an array indicates conditions under which a block or array is expected to be input indicates the default value for an item supplied by the program if the item is not input The default value is zero if not explicitly indicated The If a block is not expected by the program omit the terminating T as well If the block is expected but no data arrays are input the ter minating T must still be input If the do not appear for a block or array that block or array must be input under all conditions unless explicitly indicated otherwise See the sample case in Appendix C 3 2 AMPX MODULE INPUT
74. system Historically WIMS users have been obliged to accept the library tapes provided by Wiafrith O with littie real control over the contents The task of checking the data or preparing new data was arduous time consuming and error prone WIMSTAR 4 combined with the ENDF processing system amPx 11 4 now provides a much simplified method of producing the WIMS data files and updating the WIMS library tape WIMSTAR 4 is designed to be incorporated as a module of the AMPX II system employing the same input conventions scratch tapes and dynamic core utilization techniques Thus in a single computer run it is possible to read the ENDF B tape generate the required WIMS data for any nuclide and place this data into the WIMS library Inspection and editing facilities are provided for 1 the AMPX generated point and multigroup cross sections 2 the WIMSTAR generated multigroup cross sections and resonance tables and 3 all data files contained on an existing WIMS tape In keeping with the design of the AMPX system each step in the process reads data from one or several input tapes and writes the results to another tape thus allowing the user to proceed one step at a time verifying and if necessary correcting the results before continuing WIMSTAR 4 incorporates the linear point cross section tech niques of the TCP system A module from the RSYST system is used to solve the slowing down equation k WIMSTAR is an acronym
75. t Lawrence Livermore Laboratory Report UCRL 50400 Vol 17 Part B 1977 13 14 15 16 17 18 19 20 21 69 D E Cullen C R Weisbin Exact Doppler Broadening of Tabu lated Cross Sections Nucl Sci Eng 60 199 1976 B J Toppel A L Rago D M O Shea MC A Code to Calculate Multigroup Cross Sections Argonne National Laboratory Report ANL 7318 1967 R Q Wright J L Lucius N M Greene C W Craven Jr SUPERTOG A Program to Generate Fine Group Constants and P Scattering Matrices from ENDF B Oak Ridge National Laboratory Report ORNL TM 2679 1969 H C Honeck D R Finch FLANGE II Version 77 1 A Code to Process Thermal Neutron Data from an ENDF B Tape Savannah River Laboratory Report DP 1278 1971 W Gulden RESPU Ein Schnelles Programm zur Berechnung von Gruppenkonstanten im Resonanzbereich Reaktortagung 1973 des Deutschen Atomforums KTG Karlsruhe R Goldstein E R Cohen Theory of Resonance Absorption of Neutrons Nucl Sci Eng 13 132 1962 H C Honeck THERMOS A Thermalization Transport Theory Code for Reactor Lattice Calculations Brookhaven National Labora tory Report BNL 5826 1961 L M Petrie N F Cross KENO IV An Improved Monte Carlo Criticality Program Oak Ridge National Laboratory Report ORNL 4938 1975 L W Nordheim The Theory of Resonance Adsorption Proc Symp Appl Math Vol XI p 58 Am Math Soc
76. t be modified via the 14 array the 10 array is used for this Changes via the 10 array are completed before those of the 14 array Only if LENB lt 0 ID3 0 LENB x 3 Terminate Block 4 Sets of Blocks 3 and 4 are repeated NEB times Block 5 NES gt 0 4k Nuclide Data File Edit Control 9 1 NIN4 oF NIN of nuclide to be edited Order must match NIN list of 8 array 40 2 ID4 Data source lt 0 old WIMS tape 0 card input gt 0 WIMSTAR tape ID4 is identification no of data 3 TFR Temperature of fast and resonance data 4 AW Atomic weight 5 IAN Atomic number 6 NF Fissile resonance table trigger 0 non fissile no res tables 1 non fissile res absorption tables 2 fissile res absorption tables 3 fissile res absorption and fission yield tables 4 fissile no res tables 7 NT No of temperatures at which thermal data is tabu lated 8 NZZ No of resonance tabulations 9 NP4 Length of condensed PO scattering matrix Required only if ID4 0 Note Items 4 8 are required only if ID4 gt 0 If ID4 lt 0 they are taken from the old WIMS tape unless overridden via 4 array e g addition of a new resonance table would require increasing NZZ of old tape T Terminate Block 5 Block 6 NES gt 0 10 Potential scattering cross section oi N2 Eo 11 Slowing down power divided by lethargy width a N2 12 13 14 24 25 SW
77. the user follows the array identifier in the array origi nator field with the character U or V the input format is to be specified by the user If U is specified the FORTRAN FORMAT to be used must be supplied in columns 1 72 of the next card The format must be enclosed by the usual parentheses The data for the entire array must follow on successive cards The rules of ordinary FORTRAN input as to exponents blanks etc apply If the array data do not fill the last card the remainder must be left blank V has the same effect as U except that the format read in the last preceding U array is used APPENDIX C SAMPLE CASE This appendix contains the input for a sample case to illus trate some of the features of the AMPX WIMSTAR system The output from WIMSTAR is available from the author on request A sample of the output from NPTXS and XLACS2 can be found in the AMPX II user s manual from Reference 4 of main text NPTXS and XLACS2 were run to generate point and multigroup cross sections for fication no 129601 was used Next GXWIMS calculated the WIMS data and placed it on a new WIMSTAR tape 800226 under identification no 23202 TRANS transferred current 2320 data from NIN 1232 of an old WIMS tape to identification no 23201 of the same WIMSTAR tape Then WIMLIB added the new version of 23205 to the WIMS tape using NIN 2232 and RIN 2232 1 232 at 300 K and 9 values of op Nuclide identi using the old burnup c
78. ubroutines of a package not listed here NPTXS o XLACS2 GXWIMS OLD WIMS TRANS TAPMAN List IB NEW WIMLIB WIMS SCAN WIMS CONVRT LIB BCD LI FIGURE 1 AMPX II WIMSTAR 4 Organization AMPX II WIMSTAR 4 i Numbers refer to modules of Table 3 FIGURE 2 Overlay Structure of WIMSTAR 4 of EL APPENDIX A WIMSTAR ERROR MESSAGES The error message format is ERROR ERROR CODE n IN SUBROUTINE name where name is the routine detecting the error and n is the error code described below A FORTRAN IHO2201I error is then intentionally committed to obtain a subroutine traceback And finally the program is terminated with a Run Completion code of 16 Some routines print additional infor mation and others print self explanatory messages Routine WIMSTR ADDXS ANUSET BURNUP Error Description 1 Illegal segment name 1 Available core exhausted 1 Coefficient limit of 10 exceeded for v E on ENDF 2 Interpolation limit of 5 exceeded for v E on ENDF 3 Point limit of 5C exceeded for v E on ENDF 4 Request for v E outside valid E range 5 Invalid interpolation table for v E on ENDF 6 ENDF read error MF 1 MT 452 LNU out of range 1 3 array read error 2 ID3 lt O not allowed for new nuclide 3 4 10 14 array read error 5 Problem with burnup chain edit data check order and content of NIN lists See footnote on page 79 Routine BURN1 CONVRT EMESH FIND FLAGS
79. ulation 14 IFR Resonance table calculation flag no resonance tables Les tables only no Op tape data P tape data 0 1 2 oa and VO amp tables only no o 3 di tables including Op tape data 4 LE and VO es including Op tape data 15 IFP Aux input from Block 5 0 no 1 yes 16 NWK Size of the auxiliary work spaces 500 T Terminate Block 2 Block 3 3 Calculation Parameters 5 1 DUMAX Maximum lethargy difference 2 l 2 AU ax 7 3 2 _ acl o GT A atomic no of nuclide This is calculated internally if DUMAX 0 0 is entered SIGPA AWA APA 3 amp v lt Average potential scattering cross section of absorber Required only if IFR gt 0 Atomic weight of absorber Effective scattering radius in units of 1071 om from ENDF B file 2 data 2 If APA is entered as 0 0 the program takes the value from the END B tape Relative interpolation error used in processing point cross sections 0 001 4 List of temperatures increasing order NTE 5 List of Required only if IFR gt 0 NSP P values for resonance calculation increasing order T Terminate Block 3 Block 4 Repeat NTE times 6 Temperature Dependent Parameters 8 l ID19 NAM NSD ITFP AMPX nuclide identification no must be different for each temperature AMPX master library tape unit no from XLACS2 module of AMPX Tape unit no of point cross sectio
80. unique weighting function determined from the point values of the total cross section of the material This option should not be used for resonance nuclides Instead the IW 5 option is provided for 1 07E weighting of resonance nuclides 4 Fission arbitrary Maxvellian The arbitrary weighting function is input in the 5 and 6 arrays 5 1 0E where Cr is read from a point cross section library which was produced by the AMPX module NPTXS The logical unit for the point cross section library and the identification number for the point data set to be used are specified in the 71 array 6 Fission arbitrary Maxwellian The arbitrary function is read from an existing library of weighting functions re siding on logical unit 46 The arbitrary function is selected with the NPE parameter in the 2 array 7 1 E 0 5 7 where Le is the SIGP parameter of the 71 array and Sr is read from a point cross section library produced by NPTXS The logical unit for the point cross The weighting spectrum is a composite spectrum composed of a fission spectrum in the high energy range a constant 1 E or arbitrary spectrum IW 1 2 4 or 6 in the intermediate energy range and a Maxwellian spectrum in the low energy range The energies at which the three spectra are joined the temperature of the Maxwellian spectrum and the fission spectrum temperature are specified in the 8 array Note that a constant 1 E or arbitrary spe
81. xample of this technique is shown in the sample case Appendix C for 2320n The alternative is simply to input the burnup chain via cards when running WIMLIB 12 2 3 5 TAPMAN Segment TAPMAN provides editing and listing facilities for WIMSTAR data tapes Records on a WIMSTAR tape contain values for a specific data type burnup chain thermal smooth Pl etc for a given nuclide These records are individually manipulated with TAPMAN allow ing the user increased flexibility in handling the data If segment GXWIMS terminates abnormally and is rerun some records on the WIMSTAR tape may be repeated TAPMAN cun be used to remove the duplicate copies Full listings or simply control record summaries can be printed See Section 4 3 1 to interpret the output 2 3 6 CONVRT Segments WIMLIB TRANS and SCAN handle WIMS tapes in binary format only If the new WIMS tape is to be sent to another installation CONVRT is used to convert it to BCD i e card image format When receiving a tape from elsewhere CONVRT converts it back to binary Two small FORTRAN programs are maintained on cards to perform similar oper ations at the other computer site The attributes of the card image tape should be specified as unlabelled blocked with fixed length 80 character records Seven digits of accuracy are maintained for floating point data 3 INSTRUCTIONS FOR USING THE PROGRAMS AMPX II and WIMSTAR 4 have been implemented on an IBM system Thu

Download Pdf Manuals

image

Related Search

Related Contents

Severin KA4566 coffee maker  Pfister F-049-DECC Installation Guide  X3 User Guide  1 User manual for the 3-D macro lens model 2002 General  Accédez à la fiche technique du produit  MOTOROKR U9 Guia do usuário  Manuale di amministrazione e configurazione  IMPORTANT SAFETY INSTRUCTIONS  2204e/2208e  Kicker KMT67 Owner's Manual  

Copyright © All rights reserved.
Failed to retrieve file