User Input for NJOY 97 Search down for the desired module in caps: NJOY, RECONR, BROADR, UNRESR, HEATR, THERMR, GROUPR, GAMINR, MODER, DTFR, MATXSR, ACER, WIMSR, PLOTR, VIEWR, MIXR, PURR, LEAPR, GASPR NJOY c ****************************************************************** c * * c * njoy nuclear data processing system * c * version 97.0 * c * 31 Oct 97 * c * * c ****************************************************************** c * * c * njoy is a system of processing modules intended to convert * c * evaluated nuclear data in the endf format into forms useful * c * for practical applications. * c * * c * reconr...reconstruct pointwise cross sections from endf/b * c * resonance parameters and interpolation schemes. * c * * c * broadr...doppler broaden and thin pointwise cross sections. * c * * c * unresr...compute effective pointwise self-shielded cross * c * sections in the unresolved energy range. * c * * c * heatr...compute heat production cross sections (kerma) * c * and damage energy production. * c * * c * thermr...generate neutron scattering cross sections and * c * point-to-point scattering kernels in the thermal range * c * for free or bound atoms. * c * * c * groupr...generate self-shielded multigroup cross sections and * c * group-to-group scattering and photon production matrices. * c * * c * gaminr...compute multigroup photon interaction cross sections, * c * scattering matrices, and heat production. * c * * c * errorr...construct multigroup covariance matrices. * c * * c * covr...process covariance data from errorr * c * * c * moder...convert between endf/b standard coded mode and the * c * njoy blocked binary mode. * c * * c * dtfr...output and plot multigroup data for discrete ordinates * c * transport codes. * c * * c * ccccr...format multigroup data into the cccc standard * c * interface files isotxs, brkoxs, and dlayxs. * c * * c * matxsr...convert multigroup data into the comprehensive matxs * c * cross section interface format. * c * * c * resxsr...prepare a cccc-like file of pointwise resonance * c * cross sections for thermal flux calculations * c * * c * acer...prepare library for the los alamos continuous energy * c * monte-carlo code mcnp. * c * * c * powr...convert multigroup data into libraries for the thermal * c * powr reactor codes epri-cell and epri-cpm. * c * * c * wimsr...convert multigroup data into libraries for the * c * reactor assembly codes wims-d or wims-e. * c * * c * plotr...plot endf, pendf, gendf, or exp. cross sections, * c * distributions, or matrices. * c * * c * viewr...view plots from plotr, dtfr, covr, etc. in postscript * c * * c * mixr...mix file 3 cross sections (for example, to make * c * elemental cross sections) for plotting, etc. * c * * c * purr...generate unresolved-resonance probability tables * c * for the mcnp monte carlo code * c * * c * leapr...generate s(alpha,beta) for thermal moderators * c * * c * gaspr...add gas production (mt203-207) to pendf * c * * c * each module is an independently compiled code segment. the * c * main program controls the order in which modules are used * c * and contains the utility subroutines used by all modules. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 module option * c * * c * module six character module name, e.g., reconr. * c * it is not necessary to use quotes. * c * * c * repeat card 1 for each module desired, and * c * use the name "stop" to terminate program. * c * * c * see the comments at the start of each module for * c * its specific input instructions. * c * * c *--- code conversion --------------------------------------------* c * * c * code conversion is handled using the "*if" constructs of * c * the upd code (which is a code for maintaining fortran codes * c * similar to update). the logical variables that can be * c * set are the following: * c * * c * sw -- set for short-word machines (32 bit). leave * c * unset for 64 bit machines. * c * * c * wordio -- set to activate word-addressable io * c * (lanl cray unicos only. see acer) * c * * c ****************************************************************** RECONR c ****************************************************************** c * * c * reconstruct pointwise cross sections * c * * c * resonance cross sections are calculated using the methods of * c * resend, with modifications to the method of generating the * c * energy grid and the method of combining resonance and * c * background cross sections. * c * * c * this program generates an energy grid which is the union of * c * an input grid (if any), the resonance energies (if any), and * c * the energies of cross sections in mf3 and mf13 (or mf23). * c * the pointwise cross sections are then computed on this grid * c * and points are added so that the resonance cross sections and * c * any cross sections represented by non-linear interpolation * c * are reproduced within a specified tolerance by linear inter- * c * polation. psi-chi reconstruction can be used if desired. * c * sections which are not cross sections (mu,nu) and photon * c * multiplicities (mf12) are not processed. redundant reactions * c * are reconstructed to be the sum of their parts. the pendf * c * tape contains point cross sections in mf3 and mf13 (or mf23) * c * and a description of the processing in mf1. the mf1 diction- * c * ary is updated. the c1 and c2 fields of the second card in * c * mf1 contain the temperature and reconstruction tolerance * c * respectively. an mf2 appropriate to no resonance parameters * c * is constructed with the potential scattering length added. * c * * c * if unresolved parameters are present, the infinitely dilute * c * cross sections are computed on a special energy grid chosen * c * to preserve the required interpolation properties. this * c * table is added to the pendf tape using a special format in * c * mf2/mt152, and the table is also used to compute the * c * uresolved contributions in mf3. this allows resolved * c * resonance cross sections which overlap the resolved range * c * to be recovered by subtraction. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * nendf unit for endf/b tape * c * npend unit for pendf tape * c * card 2 * c * tlabel 66 character label for new pendf tape * c * delimited with quotes, ended with /. * c * card 3 * c * mat material to be reconstructed * c * ncards number of cards of descriptive data for new mf1 * c * (default=0.) * c * ngrid number of user energy grid points to be added. * c * (default=0.) * c * card 4 * c * err fractional reconstruction tolerance used when * c * resonance-integral error criterion (see errint) * c * is not satisfied. * c * tempr reconstruction temperature (deg kelvin) * c * (default=0.) * c * errmax fractional reconstruction tolerance used when * c * resonance-integral error criterion is satisfied * c * (errmax.ge.err, default=20.*err) * c * errint maximum resonance-integral error (in barns) * c * per grid point (default=err/10000) * c * (note: the max cross section difference for * c * linearization, errlim, and for reconstruction, * c * errmin, are also tied to errint. to get maximum * c * accuracy, set errint to a very small number. * c * for economical production, use the defaults.) * c * card 5 * c * cards ncards of descriptive comments for mt451 * c * each card delimited with quotes, ended with /. * c * card 6 * c * enode users energy grid points * c * * c * cards 3, 4, 5, 6 must be input for each material desired * c * mat=0/ terminates execution of reconr. * c * * c ****************************************************************** BROADR c ****************************************************************** c * * c * doppler broaden and thin neutron point cross sections * c * * c * a modified version of the kernel broadening method developed * c * for sigma1 (d.e.cullen, llnl) is used. cross sections * c * for low threshold reactions are unionized on the grid of the * c * total cross section, then broadened and thinned in parallel. * c * high threshold reactions are not broadened. the total and * c * and nonelastic are reconstructed to equal the sum of parts. * c * * c * the output energy grid for broadened cross sections is * c * constructed adaptively (as in reconr) so that the results * c * represent the true function within the given tolerance. * c * energy values are either removed from the original grid, or * c * new values are added between the original points. thus, a * c * given energy range can have more or fewer points than the * c * original energy grid. * c * * c * for high temperatures and low energies where the original * c * sigma1 breaks down, a new direct expansion of the doppler * c * integral is used. * c * * c * while doing the first temperature, broadr computes a number * c * of thermal quantities and resonance integrals. values for * c * both 0.0253 ev and the first temperature in ev (tev) are * c * printed out. integrals are calculated at tev only. * c * * c * the results are written out in pendf format with each * c * temperature represented as a different mat. dictionaries * c * are corrected to reflect the effects of thinning. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * nendf input endf tape (for thermal nubar only) * c * nin input pendf tape * c * nout output pendf tape * c * card 2 * c * mat1 material to be processed * c * ntemp2 number of final temperatures (maximum=10) * c * istart restart (0 no, 1 yes, default 0) * c * istrap bootstrap (0 no, 1 yes, default 0) * c * temp1 starting temperature from nin (default=0.k) * c * card 3 * c * errthn fractional tolerance for thinning * c * thnmax max. energy for broadening and thinning * c * (default=1 mev) * c * errmax fractional tolerance used when integral criterion * c * is satisfied (same usage as in reconr) * c * (errmax.ge.errthn, default=20*errthn) * c * errint parameter to control integral thinning * c * (usage as in reconr) (default=errthn/10000) * c * set very small to turn off integral thinning. * c * (a good choice for the convergence parameters * c * errthn, errmax, and errint is the same set of * c * values used in reconr) * c * card 4 * c * temp2 final temperatures (deg kelvin) * c * card 5 * c * mat1 next mat number to be processed with these * c * parameters. terminate with mat1=0. * c * * c *---input options------------------------------------------------* c * * c * the output tape will contain the ntemp2 final temperatures * c * specified. it is necessary to have temp1.le.temp2(1). * c * if temp2.eq.temp1, the data will be thinned only. * c * * c * restart continue broadening an existing pendf tape. all * c * temperatures are copied through temp1. additional * c * final temperatures are added by starting with the * c * data at temp1. * c * * c * bootstrap if bootstrap is not requested, each final tempera- * c * ture is generated by broadening directly from * c * temp1 to temp2. if bootstrap is requested, each * c * final temperature is broadened from the preceding * c * temperature. this option is faster due to the * c * thinning in the previous step. however, errors * c * accumulate. * c * * c * thnmax the upper limit for broadening and thinning is the * c * lowest of the input value of thnmax, the lowest * c * reaction threshold, or the start of the unresolved * c * range. if there is resolved-unresolved overlap, * c * the overlap region is included in the broadening. * c * a negative value of thnmax will override the * c * resolved and threshold limits. This allows the * c * first few threshold reactions to be broadened if * c * desired. the magnitude of thnmax must be chosen * c * to keep the number of broadenable reactions less * c * than or equal to the maximum of ntt (10). * c * * c ****************************************************************** UNRESR c ****************************************************************** c * * c * compute unresolved resonance cross-sections * c * * c * the method of etox is used to compute self-shielded * c * unresolved resonance cross-sections on the energy grid of * c * the unresolved parameters. subsequent interpolation is * c * to be on the cross-sections and not on the parameters. * c * additional energy grid points are added at quarter lethargy * c * intervals if only three or fewer grid points are found. * c * the accurate hwang quadrature set is used for the integrals. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * nendf unit for endf/b tape * c * nin unit for input pendf tape * c * nout unit for output pendf tape * c * card 2 * c * matd material to be processed * c * ntemp no. of temperatures (10 max) * c * nsigz no. of sigma zeroes (10 max) * c * iprint print option (0=min, 1=max) (default=0) * c * card 3 * c * temp temperatures in kelvin (including zero) * c * card 4 * c * sigz sigma zero values (including infinity) * c * cards 2, 3, 4 must be input for each material desired * c * matd=0/ terminates execution of unresr. * c * * c ****************************************************************** HEATR c ****************************************************************** c * * c * compute heating kerma (kinetic energy release in material) * c * and radiation damage energy production * c * * c * the prompt kerma is computed pointwise on the grid of the * c * total cross section from the input pendf tape and written * c * onto the output pendf tape at infinite dilution using the * c * 300 series of mt numbers. all temperatures on the input pendf * c * tape for the desired material are processed. the dictionary * c * is revised. reaction q values are obtained from the endf/b * c * tape unless the user enters his own value. partial kermas * c * can be requested for self-shielding calculations or other * c * purposes. the code uses the energy balance method where * c * photon files are available and deposits all photon energy * c * locally when files are not available. this assures * c * consistency between neutron heating and energy deposition by * c * subsequent photon interactions. an exception is made for * c * capture where recoil is computed by momentum conservation. * c * photon files are used to estimate the average photon momentum * c * when available. a diagnostic message is printed if the * c * momentum calculation leads to a significant error in * c * energy conservation. * c * * c * if desired, the energy-balance kerma factors can be compared * c * with conservative kinematic limits (set iprint=2). * c * a plot file for viewr can be automatically prepared. * c * * c * damage energy is computed using the lindhard electronic * c * screening damage function with a 25 ev displacement * c * threshold. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * nendf unit for endf/b tape * c * nin unit for input pendf tape * c * nout unit for output pendf tape * c * nplot unit for graphical check output * c * card 2 * c * matd material to be processed * c * npk number of partial kermas desired (default=0) * c * nqa number of user q values (default=0) * c * ntemp number of temperatures to process * c * (default=0, meaning all on pendf) * c * local 0/1=gamma rays transported/deposited locally * c * (default=0) * c * iprint print (0 min, 1 max, 2 check) (default=0) * c * card 3 for npk gt 0 only * c * mtk mt numbers for partial kermas desired * c * total (mt301) will be provided automatically. * c * partial kerma for reaction mt is mt+300 * c * and may not be properly defined unless * c * a gamma file for mt is on endf tape. * c * special values allowed-- * c * 303 non-elastic (all but mt2) * c * 304 inelastic (mt51 thru 91) * c * 318 fission (mt18 or mt19, 20, 21, 38) * c * 401 disappearance (mt102 thru 120) * c * 443 total kinematic kerma (high limit) * c * damage energy production values-- * c * 444 total * c * 445 elastic (mt2) * c * 446 inelastic (mt51 thru 91) * c * 447 disappearance (mt102 thru 120) * c * cards 4 and 5 for nqa gt 0 only * c * card 4 * c * mta mt numbers for users q values * c * card 5 * c * qa user specified q values (ev) * c * (if qa.ge.99.e6, read in variable qbar * c * for this reaction) * c * card 5a variable qbar (for reactions with qa flag only) * c * qbar tab1 record giving qbar versus e (1000 words max) * c * * c ****************************************************************** THERMR c ****************************************************************** c * * c * generate neutron scattering cross sections and point-to-point * c * scattering kernels in the thermal range * c * * c * add pointwise scattering cross sections and scattering * c * matrices to an existing pendf tape. cross sections are * c * added to mf3 and matrices are written in mf6 (using a * c * modified format). both using mtref for inelastic and * c * mtref+1 for elastic (if any). * c * * c * multiple scattering types (ie, h free and h in h2o) can be * c * written on one pendf tape by using different values of mtref * c * for each thermr run. if data for one mtref is already on * c * the tape, it will be replaced with the new cross sections. * c * the energy grid for coherent scattering is determined * c * adaptively so as to represent the sharp bragg edges to * c * a specified tolerance using linear interpolation. the * c * secondary energy grid for inelastic scattering is also deter- * c * mined adaptively. angular dependence is represented as * c * equally probable cosines. the initial energy grid is wired in* c * (see egrid in calcem). a special projection interpolation * c * scheme is used in groupr to integrate this relatively * c * coarse grid. * c * * c * current capabilities... * c * 1.) compute free-gas scattering matrices and normalize * c * to the elastic cross section on the old pendf tape. * c * 2.) compute incoherent matrices from read-in * c * s(alpha,beta) data. * c * 3.) compute coherent scattering from hexagonal lattices. * c * 4.) compute incoherent elastic scattering. * c * * c * for endf 3 to 5 formats, the constants used for coherent * c * elastic, incoherent elastic, and short-collision-time calcu- * c * lations are obtained from internal data statements based on * c * the original general atomic report on the evaluations * c * (ga-8774 revised, endf-269, july 1978). * c * * c * for endf6 format libraries, these constants are included * c * in the format. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * nendf endf/b tape for mf7 data * c * nin old pendf tape * c * nout new pendf tape * c * card 2 * c * matde material desired on endf tape * c * matdp material desired on pendf tape * c * nbin number of equi-probable angles * c * ntemp number of temperatures * c * iinc inelastic options * c * 0 none * c * 1 compute as free gas * c * 2 reserved * c * 3 reserved * c * 4 read s(a,b) and compute matrix * c * icoh elastic options * c * 0 none * c * 1 graphite * c * 2 beryllium * c * 3 beryllium oxide * c * 11 polyethylene * c * 12 h(zrh) * c * 13 zr(zrh) * c * natom number of principal atoms * c * mtref mt for inelastic reaction (201-250 only) * c * iprint print option (0=minimum, 1=maximum, * c * 2=max. normal + intermediate results) * c * (default=0) * c * card 3 * c * tempr temperatures (kelvin) * c * card 4 * c * tol tolerance * c * emax maximum energy for thermal treatment * c * (for temperatures greater than 3000, * c * emax and the energy grid are scaled by * c * temp/300. free gas only.) * c * * c * nendf can be endf-6 format (e.g., from leapr) while * c * nin and nout are endf-4 or 5 format, if desired. * c * * c ****************************************************************** GROUPR c ****************************************************************** c * * c * compute self-shielded group-averaged cross sections * c * * c * produces self-shielded cross sections, neutron scattering * c * matrices, and photon production matrices. scattering and * c * photon matrices may be self-shielded if desired (see init). * c * bondarenko weighting is normally used. optionally, the flux * c * can be computed for an infinite mixture of heavy absorber * c * and light moderator. delayed neutron data and thermal * c * scattering matrices are handled specially. * c * * c * the integration over initial energy is handled in the same * c * way for all reaction types by using the integrand * c * feed*xsec*flux . * c * feed is the source into final energy group gprime and * c * legendre order l from initial energy e (see getff). for * c * vectors, the feed is 1. or a yield (nubar, mubar). for two * c * body scattering, a center-of-mass gaussian integration is used * c * to obtain accurate results even for small legendre components * c * of the group-to-group scattering. additional initial energy * c * quadrature points are added to integrate the known polynomial * c * order of this feed function. feed for tabulated continuum * c * reactions is computed exactly on the endf/b grid points and * c * then interpolated at e. a special projection interpolation * c * scheme is used for thermal matrices (see getaed). the feed * c * for analytic continuum reactions is exact. * c * * c *---input specifications (free format)---------------------------* c * * c * card1 * c * nendf unit for endf/b tape * c * npend unit for pendf tape * c * ngout1 unit for input gout tape (default=0) * c * ngout2 unit for output gout tape (default=0) * c * card2 * c * matb material to be processed * c * matb lt 0 is a flag to add mts to and/or replace * c * individual mts on ngout1. * c * ign neutron group structure option * c * igg gamma group structure option * c * iwt weight function option * c * lord legendre order * c * ntemp number of temperatures * c * nsigz number of sigma zeroes * c * iprint long print option (0/1=minimum/maximum) * c * (default=1) * c * card3 * c * title run label (up to 80 characters delimited by *, * c * ended with /) (default=blank) * c * card4 * c * temp temperatures in kelvin * c * card5 * c * sigz sigma zero values (including infinity) * c * * c * if ign=1, read neutron group structure (6a and 6b) * c * card6a * c * ngn number of groups * c * card6b * c * egn ngn+1 group breaks (ev) * c * * c * if igg=1, read gamma group structure (7a and 7b) * c * card7a * c * ngg number of groups * c * card7b * c * egg ngg+1 group breaks (ev) * c * * c * weight function options (8a,8b,8c,8d) * c * card8a flux calculator parameters (iwt.lt.0 only) * c * ehi break between computed flux and bondarenko flux * c * (must be in resolved range) * c * sigpot estimate of potential scattering cross section * c * nflmax maximum number of computed flux points * c * ninwt tape unit for new flux parameters (default=0) * c * jsigz index of reference sigma zero in sigz array * c * (default=0) * c * alpha2 alpha for admixed moderator (def=o=none) * c * sam admixed moderator xsec in barns per absorber * c * atom (def=0=none) * c * beta heterogeniety parameter (def=0=none) * c * alpha3 alpha for external moderator (def=0=none) * c * gamma fraction of admixed moderator cross section in * c * external moderator cross section (def=0) * c * card8b tabulated (iwt=1 or -1 only) * c * wght read weight function as tab1 record. * c * end with a /. * c * card8c analytic flux parameters (iwt=4 or -4 only) * c * eb thermal break (ev) * c * tb thermal temperature (ev) * c * ec fission break (ev) * c * tc fission temperature (ev) * c * card8d input resonance flux (iwt=0 only) * c * ninwt tape unit for flux parameters * c * * c * card9 * c * mfd file to be processed * c * mtd section to be processed * c * mtname description of section to be processed * c * repeat for all reactions desired * c * mfd=0/ terminates this temperature/material. * c * card10 * c * matd next mat number to be processed * c * matd=0/ terminates groupr run. * c * * c *---options for input variables----------------------------------* c * * c * ign meaning * c * --- ------- * c * 1 arbitrary structure (read in) * c * 2 csewg 239-group structure * c * 3 lanl 30-group structure * c * 4 anl 27-group structure * c * 5 rrd 50-group structure * c * 6 gam-i 68-group structure * c * 7 gam-ii 100-group structure * c * 8 laser-thermos 35-group structure * c * 9 epri-cpm 69-group structure * c * 10 lanl 187-group structure * c * 11 lanl 70-group structure * c * 12 sand-ii 620-group structure * c * 13 lanl 80-group structure * c * 14 eurlib 100-group structure * c * 15 sand-iia 640-group structure * c * 16 vitamin-e 174-group structure * c * 17 vitamin-j 175-group structure * c * * c * igg meaning * c * --- ------- * c * 0 none * c * 1 arbitrary structure (read in) * c * 2 csewg 94-group structure * c * 3 lanl 12-group structure * c * 4 steiner 21-group gamma-ray structure * c * 5 straker 22-group structure * c * 6 lanl 48-group structure * c * 7 lanl 24-group structure * c * 8 vitamin-c 36-group structure * c * 9 vitamin-e 38-group structure * c * 10 vitamin-j 42-group structure * c * * c * iwt meaning * c * --- ------- * c * 1 read in smooth weight function * c * 2 constant * c * 3 1/e * c * 4 1/e + fission spectrum + thermal maxwellian * c * 5 epri-cell lwr * c * 6 (thermal) -- (1/e) -- (fission + fusion) * c * 7 same with t-dep thermal part * c * 8 thermal--1/e--fast reactor--fission + fusion * c * 9 claw weight function * c * 10 claw with t-dependent thermal part * c * 11 vitamin-e weight function (ornl-5505) * c * 12 vit-e with t-dep thermal part * c * -n compute flux with weight n * c * 0 read in resonance flux from ninwt * c * * c * mfd meaning * c * --- ------- * c * 3 cross section or yield vector * c * 5 fission chi by short-cut method * c * 6 neutron-neutron matrix (mf4/5) * c * 8 neutron-neutron matrix (mf6) * c * 12 photon prod. xsec (photon yields given, mf12) * c * 13 photon prod. xsec (photon xsecs given, mf13) * c * 16 neutron-gamma matrix (photon yields given) * c * 17 neutron-gamma matrix (photon xsecs given) * c * 18 neutron-gamma matrix (mf6) * c * note: if necessary, mfd=13 will automatically change * c * to 12 and mfd=16 will automatically change to 17 or 18.* c * 21 proton production matrix (mf6) * c * 22 deuteron production (mf6) * c * 23 triton production (mf6) * c * 24 he-3 production (mf6) * c * 25 alpha production (mf6) * c * 26 residual nucleus (a>4) production (mf6) * c * 31 proton production matrix (mf4) * c * 32 deuteron production (mf4) * c * 33 triton production (mf4) * c * 34 he-3 production (mf4) * c * 35 alpha production (mf4) * c * 36 residual nucleus (a>4) production (mf4) * c * note: if necessary, mfd=21-26 will * c * automatically change to 31-36. * c * 90+n radioactive nuclide production from nth * c * subsection for mtd in file 9. * c * 100+n radioactive nuclide production from nth * c * subsection for mtd in file 10. * c * * c * mtd meaning * c * --- ------- * c * -n process all mt numbers from the previous * c * entry to n inclusive * c * 221-250 reserved for thermal scattering * c * 258 average lethargy * c * 259 average inverse velocity (m/sec) * c * * c * automatic reaction processing options * c * ------------------------------------- * c * 3/ do all reactions in file3 on input pendf * c * 6/ do all matrix reactions in endf dictionary * c * 10/ do all radioactive isotope productions * c * 13/ do all photon production cross sections * c * 16/ do all photon production matrices * c * 21/ do all proton production matrices * c * 22/ do all deuteron production matrices * c * 23/ do all triton production matrices * c * 24/ do all he-3 production matrices * c * 25/ do all alpha production matrices * c * 26/ do all a>4 production matrices * c * * c ****************************************************************** GAMINR c ****************************************************************** c * * c * compute multigroup photon cross sections * c * * c * produce multigroup photon interaction cross sections * c * and heating kerma factors using endf/b cross sections * c * and coherent and incoherent form factors. initial energy * c * quadrature techiques are identical to those used in groupr. * c * secondary energy-angle quadrature is performed using gaussian * c * integration. * c * * c *---input specifications (free format)---------------------------* c * * c * card1 * c * nendf unit for endf/b tape * c * npend unit for pendf tape * c * ngam1 unit for input ngam tape (default=0) * c * ngam2 unit for output ngam tape (default=0) * c * card2 * c * matb material to be processed * c * input materials in ascending order * c * igg gamma group structure option * c * iwt weight function option * c * lord legendre order * c * iprint print option (0/1=minimum/maximum) (default=1) * c * card3 * c * title run label up to 80 characters (delimited by *, * c * ended with /) * c * card4 (igg=1 only) * c * ngg number of groups * c * egg ngg+1 group bounds (ev) * c * card5 (iwt=1 only) * c * wght weight function as tab1 record * c * card6 * c * mfd file to be processed * c * mtd section to be processed * c * mtname description of section to be processed * c * repeat for all reactions desired * c * mfd=0/ terminates this material * c * mfd=-1/ is a flag to process all sections present * c * for this material (termination is automatic) * c * card7 * c * matd next mat number to be processed * c * terminate gaminr run with matd=0. * c * * c *---options for input variables----------------------------------* c * * c * igg meaning * c * --- ------- * c * 0 none * c * 1 arbitrary structure (read in) * c * 2 csewg 94-group structure * c * 3 lanl 12-group structure * c * 4 steiner 21-group gamma-ray structure * c * 5 straker 22-group structure * c * 6 lanl 48-group structure * c * 7 lanl 24-group structure * c * 8 vitamin-c 36-group structure * c * 9 vitamin-e 38-group structure * c * 10 vitamin-j 42-group structure * c * * c * iwt meaning * c * --- ------- * c * 1 read in * c * 2 constant * c * 3 1/e + rolloffs * c * * c ****************************************************************** MODER c ****************************************************************** c * * c * change the mode of an endf/b tape. * c * * c * also works for pendf, gendf and covariance tapes. * c * moder can also be used to select materials from an endf-type * c * tape, or to merge several materials into a new tape. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 unit numbers * c * nin input unit * c * nout output unit * c * * c * a positive unit is coded (mode 3). * c * a negative unit is blocked binary (njoy mode). * c * * c * note: abs(nin) ge 1 and le 19 is a flag to select various * c * materials from one or more input tapes, with or * c * without mode conversion. the kind of data to be * c * processed is keyed to nin as follows: * c * nin=1, for endf or pendf input and output, * c * 2, for gendf input and output, * c * 3, for errorr-format input and output. * c * * c * cards 2 and 3 for abs (nin) ge 1 and le 19 only. * c * * c * card 2 * c * tpid tapeid for nout. 66 characters allowed * c * (delimited with *, ended with /) * c * card 3 * c * nin input unit * c * terminate moder by setting nin=0 * c * matd material on this tape to add to nout * c ****************************************************************** DTFR c ****************************************************************** c * * c * convert output of groupr to dtf format. * c * * c * processes neutron and gamma production cross sections and * c * matrices. the neutron tables can have reduced table length. * c * up-scatter is allowed. the absorption reaction is computed * c * from the total cross section and total scattering. any edits * c * can be produced which are either given in the endf/b file * c * or are linear combinations of endf/b cross sections. the * c * fission nu*sigf and chi are computed from the fission matrices * c * for all partial fission reactions. chi includes source * c * weighting. the pl tables for l.gt.0 contain the pl weighted * c * total in the total position and the pl transport cross section * c * in the absorption position. the gamma tables have gamma group * c * 1 in position 1, 2 in position 2, etc, with a table length * c * equal to the number of gamma groups. * c * * c * plots can be prepared in viewr format. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 units * c * nin input unit with data from groupr (binary). * c * nout output unit containing dtf tables (coded). * c * (default=0=none) * c * npend input unit with pendf tape for point plots. * c * (default=0=none) * c * nplot output plot info for plotr module * c * (default=0=none) * c * card 2 options * c * iprint print control (0 minimum, 1 maximum) * c * ifilm film control (0/1/2=no/yes with 1 plot per frame/ * c * yes with 4 plots per frame (default=0) * c * iedit edit control (0/1=in table/separate) (default=0) * c * * c * cards 3 through 5 only for iedit=0 * c * * c * card 3 neutron tables * c * nlmax number of neutron tables desired. * c * ng number of neutron groups * c * iptotl position of total cross section * c * ipingp position of in-group scattering cross section. * c * itabl neutron table length desired. * c * ned number of entries in edit table (default=0). * c * ntherm number of thermal groups (default=0). * c * card 3a only for ntherm ne 0 * c * card 3a thermal incoherent and coherent mts * c * mti mt for thermal incoherent data * c * mtc mt for thermal coherent data (default=0) * c * nlc no. coherent legendre orders (default=0) * c * card 4 edit names * c * six character hollerith names for edits for as many * c * cards as needed. there will be iptotl-3 names read. * c * each name is delimited with *. * c * card 5 edit specifications * c * ned triplets of numbers on as many cards as needed. * c * positions can appear more than once. * c * reaction types can appear more than once. * c * jpos position of edit quantity. * c * mt endf/b reaction number. * c * mult multiplicity to be used when adding this mt. * c * * c * card 6 for iedit=1 * c * * c * card 6 claw-format tables * c * nlmax number of neutron tables (def=5) * c * ng number of neutron groups (def=30) * c * (number of thermal groups is zero) * c * * c * card 7 gamma ray tables * c * nptabl number of gamma tables desired (default=0) * c * ngp number of gamma groups (default=0) * c * card 8 material description * c * one card for each table set desired. * c * empty card (/) terminates execution of dtfr. * c * hisnam 6-character isotope name * c * mat material number as in endf/b (default=0) * c * jsigz index number of sigma-zero desired (default=1) * c * dtemp temperature desired (default=300) * c * * c ****************************************************************** MATXSR c ****************************************************************** c * * c * produce matxs interface file from njoy intermediate cross * c * section data from group or gaminr * c * * c * the matxs file uses a generalized, flexible format based on * c * the cccc interface conventions. working from groupr and/or * c * gaminr output tapes, this module can process neutron, * c * thermal, photon, and charged-particle data into a * c * single output file. this file can then be used by the * c * transx code to prepare data libraries for transport codes. * c * c * a matxs file specification may be found following the * c * input instructions. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 units * c * ngen1 input unit for data from groupr * c * ngen2 input unit for data from gaminr * c * nmatx output unit for matxs * c * ngen3 incident proton data from groupr (default=0) * c * ngen4 incident deuteron data from groupr (default=0) * c * ngen5 incident triton data from groupr (default=0) * c * ngen6 incident he3 data from groupr (default=0) * c * ngen7 incident alpha data from groupr (default=0) * c * card 2 user identification * c * ivers file version number (default=0) * c * huse user id (up to 16 characters, delimited by *, * c * ended by /) (default=blank) * c * card 3 file control * c * npart number of particles for which group * c * structures are given * c * ntype number of data types in set * c * nholl number of cards to be read for hollerith * c * identification record. * c * nmat number of materials desired * c * card 4 set hollerith identification * c * hsetid hollerith identification of set * c * (each line can be up to 72 characters, * c * delimited with *, ended by /) * c * card 5 particle identifiers * c * hpart hollerith identifiers for particles * c * (up to 8 characters each) * c * card 6 energy groups * c * ngrp number of groups for each particle * c * card 7 data type identifiers * c * htype hollerith identifiers for data types * c * (up to 8 characters each) * c * card 8 input particle ids * c * jinp input particle id for each data type * c * card 9 output particle ids * c * joutp output particle id for each data type * c * card 10 material data (one card per material) * c * hmat hollerith material identifier * c * (up to 8 characters each) * c * matno integer material identifier * c * (endf/b mat number) * c * matgg mat number for photoatomic data * c * (default=100*(matno/100) as in endf-6) * c * * c ****************************************************************** ACER c ****************************************************************** c * * c * prepare a data library for mcnp, * c * the los alamos continuous energy monte carlo code * c * * c * --- continuous (fast) data --- * c * * c * reaction cross sections are reconstructed on the grid of the * c * total cross section from the input pendf tape (assumed to be * c * linearized and unionized). redundant reactions (except for * c * mt1, mt452, and reactions needed for photon yields) are * c * removed. mt18 is considered redundant if mt19 is present. * c * angular distributions are converted into equal probability * c * bins. tabulated energy distributions are converted into * c * "law 4" probability distributions. analytic secondary energy * c * distributions are converted into their ace formats. * c * energy-angle distributions are converted to the "law 44" * c * cummulative probability distributions form with parameters * c * for kalbach systematics. all photon production cross * c * sections are combined on the cross section energy grid. * c * if provided, multigroup photon production data is summed * c * and converted into a set of equally probable emission * c * energies for each input group. detailed photon production * c * data can be generated directly from files 12, 13, 14, 15, * c * and 16 from the input endf tapes and written out using the * c * "law 4" cummulative energy distribution format. * c * * c * --- thermal data --- * c * * c * the data from the pendf tape as prepared by the thermr * c * module is read in. inelastic and incoherent elastic cross * c * sections are stored directly. coherent elastic cross * c * sections are converted to a cummulative "stair step" form * c * and stored. the angular representation for incoherent * c * elastic is stored directly. none is needed for coherent * c * elastic. the incoherent inelastic energy distributions * c * are converted into probability bins with the equally * c * probable angles left unchanged. the bins can have equal * c * probabilities or variable probabilities. in the latter * c * case, outlying bins with smaller probabilities are provided * c * to extend the sampling to rare events. * c * * c * --- dosimetry data --- * c * * c * endf cross sections for dosimetry reactions are simply * c * stored in ace format without changing the energy grid. * c * the endf interpolation law is also provided. * c * * c * --- photoatomic data --- * c * * c * photon interaction cross sections are stored on the grid of * c * the total cross section. the coherent form factor is * c * stored together with an integral over the form factor that * c * is used inb sampling for coherent scattering. the * c * incoherent scattering function is simply stored. photon * c * heating is calculated from the incoherent scattering data, * c * the pair production data, and the photoelectric absorption * c * data. flourescence data. * c * * c * --- output --- * c * * c * the ace output file can be type 1 (formatted), type 2 * c * (f77 binary), or type 3 (lanl binary). a line of file * c * directory information is also written. it must normally * c * be edited to tell the system the path to the file. acer * c * can also be used to print, edit, or convert the mode of * c * existing ace-format files. * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * nendf unit for input endf/b tape * c * npend unit for input pendf tape * c * ngend unit for input multigroup photon data * c * nace unit for output ace tape * c * ndir unit for output mcnp directory * c * card 2 * c * iopt type of acer run option * c * 1 fast data * c * 2 thermal data * c * 3 dosimetry data * c * 4 photo-atomic data * c * 7 read type 1 ace files to print or edit * c * 8 read type 2 ace files to print or edit * c * 9 read type 3 ace files to print or edit * c * iprint print control (0 min, 1 max, default=1) * c * ntype ace output type (1, 2, or 3, default=1) * c * suff id suffix for zaid (default=.00) * c * nxtra number of iz,aw pairs to read in (default=0) * c * card 3 * c * hk descriptive character string (70 char max) * c * delimited by quotes * c * card 4 (nxtra.gt.0 only) * c * iz,aw nxtra pairs of iz and aw * c * * c * --- fast data (iopt=1 only) --- * c * * c * card 5 * c * matd material to be processed * c * tempd temperature desired (kelvin) (default=300) * c * card 6 * c * err tolerance for thinning distributions * c * iopp detailed photons (0=no, 1=yes, default=1) * c * card 7 * c * type of thinning is determined by sign of thin(1) * c * (pos. or zero/neg.=energy skip/integral fraction) * c * (all entries defaulted=no thinning) * c * thin(1) e1 energy below which to use all energies (ev) * c * or iwtt weighting option (1=flat,2=1/e) * c * (1/e actually has weight=10 when e lt .1) * c * thin(2) e2 energy above which to use all energies * c * or target number of points * c * thin(3) iskf skip factor--use every iskf-th energy * c * between e1 and e2, or rsigz reference sigma zero * c * * c * --- thermal data (iopt=2 only) --- * c * * c * card 8 * c * matd material to be processed * c * tempd temperature desired (kelvin) (default=300) * c * tname thermal zaid name ( 6 char max, def=za) * c * card 8a * c * iza01 moderator component za value * c * iza02 moderator component za value (def=0) * c * iza03 moderator component za value (def=0) * c * card 9 * c * mti mt for thermal incoherent data * c * nbint number of bins for incoherent scattering * c * mte mt for thermal elastic data * c * ielas 0/1=coherent/incoherent elastic * c * nmix number of atom types in mixed moderator * c * (default=1, not mixed) * c * (example, 2 for beo or c6h6) * c * emax maximum energy for thermal treatment (ev) * c * (default=1000.=determined from mf3, mti) * c * iwt weighting option * c * 0/1=variable/constant (default=variable) * c * * c * --- dosimetry data (iopt=3 only) --- * c * * c * card 10 * c * matd material to be processed * c * tempd temperature desired (kelvin) (default=300) * c * * c * --- photo-atomic data (iopt=4 only) --- * c * * c * card 11 * c * matd material to be processed * c * * c * --- print or edit existing files (iopt=7-9) --- * c * * c * no additional input cards are required. mount the old * c * ace tape on "npend". the code can modify zaid, hk, * c * the (iz,aw) list, and the type of the file. use suff<0 * c * to leave the old zaid unchanged. use just "/" on * c * card 3 to leave the comment field hk unchanged. use * c * nxtra=0 to leave the old iz,aw list unchanged. * c * the code can modify zaid, hk, and type of file. * c * * c * exhaustive consistency checks are automatically made on * c * the input file. if ngend.ne.0, a set of standard ace plots * c * are prepared on unit ngend as plotr input instructions. * c * * c ****************************************************************** WIMSR c ****************************************************************** c * * c * format multigroup cross sections from groupr for wims * c * * c *---input specifications (free format)---------------------------* c * * c * card 1 * c * ngendf unit for input gendf tape * c * nout unit for output wims tape * c * * c * card 2 * c * iprint print option * c * 0=minimum (default) * c * 1=regular * c * 2=1+intermediate results * c * iverw wims version * c * 4=wims-d (default) * c * 5=wims-e * c * igroup group option * c * 0=69 groups (default) * c * 9=user's choice * c * * c * card 2a (igroup.eq.9 only) * c * ngnd number of groups * c * nfg number of fast groups * c * nrg number of resonance groups * c * igref reference group (default is last fast group) * c * * c * card 3 * c * mat endf mat number of the material to be processed * c * nfid not used * c * rdfid identification of material for the wims library * c * iburn burnup data option * c * -1=suppress printout of burnup data * c * 0=no burnup data is provided (default) * c * 1=burnup data is provided in cards 5 and 6 * c * * c * card 4 * c * ntemp no. of temperatures to process for this material * c * in the thermal energy range c * (0=all found on input tape) * c * nsigz no. of sigma zeroes to process for this material * c * (0=all found on input tape) * c * sgref reference sigma zero * c * (.ge. 1.e10 to select all cross sect. at inf.dil.* c * but fully shielded elastic x-sect, * c * .lt. 1.e10 to select all x-sect at inf.dil. * c * =sig0 from the list on groupr input to * c * select all x-sect. at that sig0) * c * ires resonance absorber indicator * c * 0=no resonance tables * c * >0=ires temperatures processed * c * sigp potential cross section from endf/b. * c * (if zero, replace by the elastic cross section) * c * mti thermal inelastic mt (default=0=none) * c * mtc thermal elastic mt (default=0=none) * c * ip1opt include p1 matrices * c * 0=yes * c * 1=no, correct p0 ingroups (default) * c * inorf resonance fission (if found) * c * 0=include resonance fission (default) * c * 1=do not include * c * isof fission spectrum * c * 0=do not include fission spectrum (default) * c * 1=include fission spectrum * c * ifprod fission product flag * c * 0=not a fission product (default) * c * 1=fission product, no resonance tables * c * 2=fission product, resonance tables * c * jp1 transport correction neutron current spectrum flag * c * 0=use p1-flux for transport correction (default) * c * >0=read in jp1 values of the neutron current * c spectrum from input * c * * c * following two cards are for iburn gt 0 only * c * card 5 * c * ntis no. of time-dependent isotopes * c * efiss energy released per fission * c * * c * card 6 (repeat this card ntis times) * c * identa ident of fission product isotope * c * yield fission yield of identa from burnup of mat * c * * c * card 7 * c * lambda resonance-group goldstein lambdas (13 for * c * default 69-group structure, nrg otherwise). * c * * c * card 8 (only when jp1>0) c * p1flx current spectrum (jp1 entries read, the rest are * c * set with the default p1-flux calculated by njoy). * c ****************************************************************** PLOTR c ****************************************************************** c * * c * plot cross sections * c * * c * handles endf/b data, pendf or gendf data at specified temper- * c * atures, or experimental input data. several plots can be * c * given on each set of axes, with both left and right scales. * c * also, several graphs can be given on each page or display. * c * error bars may be included for input data. flexible titles * c * and legend blocks are allowed. all standard combinations of * c * log and linear axes are supported, either grids or tick marks * c * can be requested, and scales can be chosen automatically * c * or set by the user. in some cases, the x axis is thinned. * c * in other cases, extra points are added so that, for example, * c * linear-linear data plots correctly on a log-log graph. * c * a limited capability for 3-d plots of angle and energy * c * is included, and the endf-vi file 6 format is supported. * c * * c * plotr writes plot commands on an output file for later use * c * by the viewr module or an external graphics program. * c * * c *---input--------------------------------------------------------* c * * c * card 0 * c * nplt unit for output plot commands * c * nplt0 unit for input plot commands * c * default=0=none * c * output plot commands are appended * c * to the input plot commands, if any. * c * card 1 * c * lori page orientation (def=1) * c * 0 portrait (7.5x10in) * c * 1 landscape (10x7.5in) * c * istyle character style (def=2) * c * 1 = roman * c * 2 = swiss * c * size character size option * c * pos = height in page units * c * neg = height as fraction of subplot size * c * (default=0.30) * c * ipcol page color (def=white) * c * 0=white * c * 1=navajo white * c * 2=blanched almond * c * 3=antique white * c * 4=very pale yellow * c * 5=very pale rose * c * 6=very pale green * c * 7=very pale blue * c * * c * -----repeat cards 2 through 13 for each curve----- * c * * c * card 2 * c * iplot plot index * c * 99 = terminate plotting job * c * 1 = new axes, new page * c * -1 = new axes, existing page * c * n = nth additional plot on existing axes * c * -n = start a new set of curves using * c * the alternate y axis * c * default = 1 * c * iwcol window color (def=white) * c * color list same as for ipcol above * c * factx factor for energies (default=1.) * c * facty factor for cross-sections (default=1.) * c * xll,yll lower-left corner of plot area * c * ww,wh,wr window width, height, and rotation angle * c * (plot area defaults to one plot per page) * c * * c * -----cards 3 thru 7 for iplot = 1 or -1 only----- * c * * c * card 3 * c * t1 first line of title * c * 60 characters allowed. * c * default=none * c * * c * card 3a * c * t2 second line of title * c * 60 characters allowed. * c * default=none * c * * c * card 4 * c * itype type for primary axes * c * 1 = linear x - linear y * c * 2 = linear x - log y * c * 3 = log x - linear y * c * 4 = log x - log y * c * set negative for 3d axes * c * default=4 * c * jtype type for alternate y axis or z axis * c * 0 = none * c * 1 = linear * c * 2 = log * c * default=0 * c * igrid grid and tic mark control * c * 0 = no grid lines or tic marks * c * 1 = grid lines * c * 2 = tic marks on outside * c * 3 = tic marks on inside * c * default=2 * c * ileg option to write a legend. * c * 0 = none * c * 1 = write a legend block with upper left * c * corner at xtag,ytag (see below) * c * 2 = use tag labels on each curve with * c * a vector from the tag to the curve * c * default=0 * c * xtag x coordinate of upper left corner * c * of legend block * c * ytag y coord of upper left corner * c * default=upper left corner of plot * c * * c * card 5 * c * el lowest energy to be plotted * c * eh highest energy to be plotted * c * xstep x axis step * c * default = automatic scales * c * (default all 3, or none) * c * (the actual value of xstep is * c * ignored for log scales) * c * * c * card 5a * c * xlabl label for x axis * c * 60 characters allowed. * c * default="energy (ev)" * c * * c * card 6 * c * yl lowest value of y axis. * c * yh highest value of y axis. * c * ystep step for y ayis (linear scales only) * c * default = automatic scales * c * (default all 3, or none) * c * (the actual value of ystep is * c * ignored for log scales) * c * * c * card 6a * c * ylabl label for y axis * c * 60 characters allowed. * c * default="cross section (barns)" * c * * c * card 7 (jtype.gt.0 only) * c * rbot lowest value of secondary y axis or z axis * c * rtop highest value of secondary y axis or z axis * c * rstep step for secondary y axis or z axis * c * default for last three = automatic * c * * c * card 7a (jtype.gt.0 only) * c * rl label for alternate y axis or z axis * c * 60 characters allowed. * c * default=blank * c * * c * -----cards 8 thru 9 are always given----- * c * * c * card 8 * c * iverf version of endf tape * c * set to zero for data on input file * c * and ignore rest of parameters on card * c * set to 1 for gendf data * c * nin input tape * c * can change for every curve if desired. * c * matd desired material * c * mfd desired file * c * mtd desired section * c * mtd=0 means loop over all reactions in mfd * c * (usually one page per mt, but for mf=3, * c * resonance reactions may have several pages) * c * temper temperature for endf/b data (k) * c * default=0. * c * nth,ntp,nkh see below (defaults=1) * c * * c * special meanings for nth,ntp,nkh for file 3 or 5 data * c * nth number of subsection to plot * c * (works for isomer prod, delayed n, etc.) * c * ntp not used * c * nkh not used * c * * c * special meanings for nth,ntp,nkh for file 6 data * c * nth index for incident energy * c * ntp number of dep. variable in cyle to plot * c * (or angle number for law 7) * c * nkh number of outgoing particle to plot * c * * c * special meanings for nth,ntp,nkh for gendf mf=3 data * c * nth=0 for flux per unit lethargy * c * nth=1 for cross section (default) * c * ntp=1 for infinite dilution (default) * c * ntp=2 for next lowest sigma-zero values, etc. * c * nkh=1 for p0 weighting (default) * c * nkh=2 for p1 weighting (total only) * c c * special meaning for nth for gendf mf=6 data * c * nth=1 plot 2-d spectrum for group 1 * c * nth=2 plot 2-d spectrum for group 2 * c * etc. * c * no special flags are needed for mf=6 3d plots * c * * c * special meanings for nth and ntp for mf7 plots * c * nth is index for indep. variable (alpha or beta) * c * ntp=1 selects alpha as indep. variable (default) * c * ntp=2 selects beta as indep. variable * c * nkh=1 selects normal s(alpha,beta) * c * nkh=2 selects script s(alpha,-beta) * c * nkh=3 selects script s(alpha,beta) * c * * c * -----cards 9 and 10 for 2d plots only----- * c * * c * card 9 * c * icon symbol and connection option * c * 0 = points connected, no symbols * c * -i = points not connected, symbol at every * c * ith point * c * i = points connected, symbol at every ith * c * points * c * default=0 * c * isym no. of symbol to be used * c * 0 = square * c * 1 = octagon * c * 2 = triangle * c * 3 = cross * c * 4 = ex * c * 5 = diamond * c * 6 = inverted triangle * c * 7 = exed square * c * 8 = crossed ex * c * 9 = crossed diamond * c * 10 = crossed octagon * c * 11 = double triangle * c * 12 = crossed square * c * 13 = exed octagon * c * 14 = triangle and square * c * 15 = filled circle * c * 16 = open circle * c * 17 = open square * c * 18 = filled square * c * default=0 * c * idash type of line to plot * c * 0 = solid * c * 1 = dashed * c * 2 = chain dash * c * 3 = chain dot * c * 4 = dot * c * default=0 * c * iccol curve color (def=black) * c * 0=black * c * 1=red * c * 2=green * c * 3=blue * c * 4=magenta * c * 5=cyan * c * 6=brown * c * 7=purple * c * ithick thickness of curve (def=1) * c * 0 = invisible (for shaded areas) * c * ishade shade pattern * c * 0 = none * c * 1 to 10 = 10% to 100% gray * c * 11 to 20 = 45 deg right hatching * c * 21 to 30 = 45 deg left hatching * c * 31 to 40 = 45 deg cross hatching * c * 41 to 50 = shades of green * c * 51 to 60 = shades of red * c * 61 to 70 = shades of brown * c * 71 to 80 = shades of blue * c * default=0 * c * * c * card 10 ---ileg.ne.0 only--- * c * aleg title for curve tag or legend block * c * 60 characters allowed. * c * default=blank * c * * c * card 10a ---ileg.eq.2 only--- * c * xtag x position of tag title * c * ytag y position of tag title * c * xpoint x coordinate of vector point * c * (.le.0 to omit vector) * c * * c * -----card 11 for 3d plots only----- * c * * c * card 11 * c * xv,yv,zv abs. coords of view point * c * defaults=15.,-15.,15. * c * x3,y3,z3 abs. sides of work box volume * c * defaults=2.5,6.5,2.5 * c * * c * set x3 or y3 negative to flip the order of the * c * axis on that side of the work box. * c * * c * -----cards 12 thru 13 for iverf = 0 only----- * c * * c * card 12 * c * nform format code for input data * c * 0 = free format input with * c * optional x and y error bars * c * * c * card 13 ---nform = 0 only--- * c * xdata dependent value * c * terminate with empty card (/) * c * ydata independent value * c * yerr1 lower y error limit * c * no y error bar if zero * c * yerr2 upper y error limit * c * if zero, equals yerr1 * c * xerr1 x left error limit * c * no x error bar if zero * c * xerr2 x right error limit * c * if zero, equals xerr1 * c * * c * * c * all curves contain at least 10 points per decade (see delta). * c * code can plot curves containing fewer than 2000 points (see * c * max) without thinning. curves with more points are thinned * c * based on a minimum spacing determined from max and the * c * length of the x axis. * c * * c ****************************************************************** VIEWR c ****************************************************************** c * * c * view plots generated by plotr, covr, or dtfr. * c * postscript version * c * * c * character specifications are similar to disspla, except that * c * the default case is lower instead of upper. this allows * c * mixed-case strings to be used in postscript mode. the * c * following shift characters are used: * c * < = upper-case standard * c * > = lower-case or mixed-case standard * c * [ = upper-case greek * c * ] = lower-case or mixed-case greek * c * # = instructions * c * give one of the shift characters twice to get it instead of * c * its action. the following instructions are supported: * c * Ev = elevate by v as a fraction of the height * c * if v is missing or D is given, use .5 * c * Lv = lower by v as a fraction of the height * c * if v is missing or D is given, use .5 * c * Hv = change height by v as a fraction of the height * c * if v is missing or D is given, use .5 * c * Fi = change to font number i * c * Mi = change mode number, where mode 0 is the lower 128 * c * postscript characters and mode 1 is the upper 128 * c * X = reset E, L, or H. Font and Mode must be * c * reset explicitly. * c * c is a real number, i is an integer. * c * * c *---input--------------------------------------------------------* c * * c * card 1 * c * infile input file * c * nps postscript output file * c * * c *---data file format---------------------------------------------* c * * c * card 1 * c * lori page orientation (def=1) * c * 0 portrait (7.5x10in) * c * 1 landscape (10x7.5in) * c * istyle character style (def=2) * c * 1 roman * c * 2 swiss * c * size character size option * c * pos = height in page units * c * neg = height as fraction of subplot size * c * (default=.30) * c * ipcol page color (def=white) * c * 0=white * c * 1=navajo white * c * 2=blanched almond * c * 3=antique white * c * 4=very pale yellow * c * 5=very pale rose * c * 6=very pale green * c * 7=very pale blue * c * * c * -----repeat cards 2 through 13 for each curve----- * c * * c * card 2 * c * iplot plot index * c * 99 = terminate plotting job * c * 1 = new axes, new page * c * -1 = new axes, existing page * c * n = nth additional plot on existing axes * c * -n = start a new set of curves using * c * the alternate y axis * c * default = 1 * c * iwcol window color (def=white) * c * color list same as for ipcol above * c * factx factor for energies (default=1.) * c * facty factor for cross-sections (default=1.) * c * xll,yll lower-left corner of plot area * c * ww,wh,wr window width, height, and rotation angle * c * (plot area defaults to one plot per page) * c * * c * -----cards 3 thru 7 for iplot = 1 or -1 only----- * c * * c * card 3 * c * t1 first line of title * c * 60 characters allowed. * c * default=none * c * * c * card 3a * c * t2 second line of title * c * 60 characters allowed. * c * default=none * c * * c * card 4 * c * itype type for primary axes * c * 1 = linear x - linear y * c * 2 = linear x - log y * c * 3 = log x - linear y * c * 4 = log x - log y * c * set negative for 3d axes * c * 0 = no plot, titles only * c * default=4 * c * jtype type for alternate y axis or z axis * c * 0 = none * c * 1 = linear * c * 2 = log * c * default=0 * c * igrid grid and tic mark control * c * 0 = no grid lines or tic marks * c * 1 = grid lines * c * 2 = tic marks on outside * c * 3 = tic marks on inside * c * default=2 * c * ileg option to write a legend. * c * 0 = none * c * 1 = write a legend block with upper left * c * corner at xtag,ytag (see below) * c * 2 = use tag labels on each curve with * c * a vector from the tag to the curve * c * default=0 * c * xtag x coordinate of upper left corner * c * of legend block * c * ytag y coord of upper left corner * c * default=upper left corner of plot * c * * c * card 5 * c * xmin lowest energy to be plotted * c * xmax highest energy to be plotted * c * xstep x axis step * c * default = automatic scales * c * (for linear, give all 3, or none) * c * (for log, give first 2, or none) * c * * c * card 5a * c * xlabl label for x axis * c * 60 characters allowed. * c * (default = no label, no numbering) * c * * c * card 6 * c * ymin lowest value of y axis. * c * ymax highest value of y axis. * c * ystep step for y ayis (linear scales only) * c * default = automatic scales * c * (for linear, give all 3, or none) * c * (for log, give first 2, or none) * c * * c * card 6a * c * ylabl label for y axis * c * 60 characters allowed. * c * (default = no label, no numbering) * c * * c * card 7 (jtype.gt.0 only) * c * rmin lowest value of secondary y axis or z axis * c * rmax highest value of secondary y axis or z axis * c * rstep step for secondary y axis or z axis * c * (default = automatic scale) * c * (for linear, give all 3, or none) * c * (for log, give first 2, or none) * c * * c * card 7a (jtype.gt.0 only) * c * rl label for alternate y axis or z axis * c * 60 characters allowed. * c * (default = no label, no numbering) * c * * c * card 8 -- dummy input card for consistency with plotr * c * it always should be 0/ * c * * c * -----cards 9 and 10 for 2d plots only----- * c * * c * card 9 * c * icon symbol and connection option * c * 0 = points connected, no symbols * c * -i = points not connected, symbol at every * c * ith point * c * i = points connected, symbol at every ith * c * points * c * default=0 * c * isym no. of symbol to be used * c * 0 = square * c * 1 = octagon * c * 2 = triangle * c * 3 = cross * c * 4 = ex * c * 5 = diamond * c * 6 = inverted triangle * c * 7 = exed square * c * 8 = crossed ex * c * 9 = crossed diamond * c * 10 = crossed octagon * c * 11 = double triangle * c * 12 = crossed square * c * 13 = exed octagon * c * 14 = triangle and square * c * 15 = filled circle * c * 16 = open circle * c * 17 = open square * c * 18 = filled square * c * 19 = filled diamond * c * 20 = filled triangle * c * 21 = filled inverted triangle * c * 22 = crossed circle * c * 23 = exed circle * c * 24 = exed diamond * c * default=0 * c * idash type of line to plot * c * 0 = solid * c * 1 = dashed * c * 2 = chain dash * c * 3 = chain dot * c * 4 = dot * c * 5 = invisible * c * default=0 * c * iccol curve color (def=black) * c * 0=black * c * 1=red * c * 2=green * c * 3=blue * c * 4=magenta * c * 5=cyan * c * 6=brown * c * 7=purple c * ithick controls thickness of curve * c * 0 = invisible (for shaded areas) * c * (default=1) * c * ishade shade pattern * c * 0 = none * c * 1 to 10 = 10% to 100% gray * c * 11 to 20 = 45 deg right hatching * c * 21 to 30 = 45 deg left hatching * c * 31 to 40 = 45 deg cross hatching * c * 41 to 50 = shades of green * c * 51 to 60 = shades of red * c * 61 to 70 = shades of brown * c * 71 to 80 = shades of blue * c * default=0 * c * * c * card 10 ---ileg.ne.0 only--- * c * aleg title for curve tag or legend block * c * 60 characters allowed. * c * default=blank * c * * c * card 10a ---ileg.eq.2 only--- * c * xtag x position of tag title * c * ytag y position of tag title * c * xpoint x coordinate of vector point * c * (.le.0 to omit vector) * c * * c * -----card 11 for 3d plots only----- * c * * c * card 11 * c * xv,yv,zv abs. coords of view point * c * defaults= 15.,-15.,15. * c * x3,y3,z3 abs. sides of work box volume * c * defaults=2.5,6.5,2.5 * c * * c * set x3 negative to flip the order of the axis on * c * that side of the box (secondary energy, cosine). * c * * c * card 12 * c * nform format code for input data * c * 0 = free format input with * c * optional x and y error bars * c * 1 = free format input for a * c * 3d family of curves z(x) vs y * c * * c * card 13 ---nform = 0 only--- 2-d data * c * xdata independent value * c * terminate with empty card (/) * c * ydata dependent value * c * yerr1 lower y error limit * c * no y error bar if zero * c * yerr2 upper y error limit * c * if zero, equals yerr1 * c * xerr1 x left error limit * c * no x error bar if zero * c * xerr2 x right error limit * c * if zero, equals xerr1 * c * * c * card 14 ---nform = 1 only--- 3-d data * c * y y value for curve * c * repeat cards 13 and 13a for each curve * c * terminate with empty card (/) * c * * c * card14a ---nform = 1 only--- * c * x x value * c * z z value * c * repeat card 13a for each point in curve * c * terminate with empty card (/) * c * disspla version requires same x grid * c * for each value of y. * c * * c ****************************************************************** MIXR c ****************************************************************** c * * c * mixr * c * * c * construct a