! NBI SHOT !---- JETtransp Update Mon Jun 7 09:33:24 BST 2004 ---------------------- !HOST = pppl !UFILE = 0 FTIME = 0.343 KMDSPLUS = 1 MDS_PATH = '\top.INPUTS' MDS_SERVER = 'rca.fusion.org.uk' MDS_SHOT = 07085C01 MDS_TREE = 'trmast' NSHOT = 07085 TINIT = 0.031 ! End JETtransp Update Mon Jun 7 09:33:24 BST 2004 ---------------------- !---- JETtransp Update Thu Feb 12 16:24:43 GMT 2004 ---------------------- ! MDS SHOT FLAG !------------------------------------------------------------------------------------ ! Run coordination, time steps etc.: !------------------------------------------------------------------------------------ ! PARAMETERS HERE CHANGE FOR EACH SHOT TQMODA(1) = 0.1 ! Time to switch from intrvl 1 (EFIT q) to 2 (PFDE) [ s ] !------------------------------------------------------------------------------------ DTMAXG = 0.001 ! Time step for equilibrium and geometry [ s ] NQMODA = 4,1 ! q-profile method for intervals 1 and 2 [ - ] NQMODB = 1,1 ! -"- [ - ] ! nlncsmoo may be important if you have trouble NLNCSMOO =.false. ! Provide in code smoothing of Te and ne? [ - ] NLSQUEEZ =.false. ! Apply NCLASS "squeezing factor"? [ - ] NLEBAL =.true. ! .true. to perform e- energy balance calculation [ - ] NLTIBL =.true. ! .true. to perform ion energy balance calculation [ - ] NLPBAL =.true. ! .true. to perform particle balance calculation [ - ] LEVTRK = 2 ! Straight line tracking (2=full generality) [ - ] LEVGEO = 7 ! Eqm. code (eg. 7=ESC, 10=RZSolver) [ - ] RZSO_CONTROL(6)= 600 ! if eq solver defaults to RZsolver, adds points round the axis to help contouring round the axis NMOM = 8 ! # moments for flux surface representation (max=8) [ - ] NLBCCW =.FALSE. ! BT direction co-false, ctr-true [ - ] NLJCCW =.TRUE. ! Ip direction co-true, ctr-false [ - ] LEVGB = 3 ! Level of tracking for beam code (3=full gyro) [ - ] MHDOUT = 1 ! MHD diagnostic O/P setting [ - ] ! Miscellaneous data to be read in: NLDFLX =.false. ! Read in diamagnetic flux data [ - ] NLNTX =.true. ! Read in neutron data [ - ] !------------------------------------------------------------------------------------ ! Diagnostic output: !------------------------------------------------------------------------------------ SEDIT = 0.001 ! Time interval for SCALAR output [ s ] STEDIT = 0.001 ! Time interval for PROFILE output [ s ] !------------------------------------------------------------------------------------ ! Plasma composition + miscellaneous: !------------------------------------------------------------------------------------ AIMP = 12.0 ! Impurity species atomic mass (eg; carbon=12) [amu] XZIMP = 6.0 ! Impurity species atomic charge (eg. carbon=6 ) [ e ] NG = 2 ! Initial # of hydrogenic species [ - ] NGMAX = 2 ! Maximum number of hydrogenic species [ - ] APLASM = 2.0, 1.0 ! Bulk plasma atomic weight(s) (species 1 & 2) [amu] BACKZ = 1.0, 1.0 ! Bulk plasma atomic charge(s) (species 1 & 2) [ e ] FTEMIN = 10.0 ! Minimum allowed Te [eV ] !------------------------------------------------------------------------------------ ! Thermal plasma input data symmetrization parameters: !------------------------------------------------------------------------------------ NRITER =-4 ! -4 = vs. r/a [ - ] NRINER =-4 ! -"- [ - ] NRITI2 =-4 ! -"- [ - ] NRIOMG =-4 ! -"- [ - ] NRIZF2 =-4 ! cab 4/11/02 rotation vs r and t NRIQPR =-6 ! -6 = vs. normalized poloidal flux [ - ] NSYTER = 0 ! Do not symmetrize [ - ] NSYNER = 0 ! -"- [ - ] NSYQPR = 0 ! -"- [ - ] NSYTI2 = 0 ! -"- [ - ] NSYOMG = 0 ! -"- [ - ] NSYZF2 = 0 !------------------------------------------------------------------------------------ ! Zeff model: !------------------------------------------------------------------------------------ NLZEFM =.false. ! .true. to use resisitivity Zeff for composition [ - ] NLZFI2 =.true. ! .true. to read in 2D Zeff data [ - ] !NZEFMOD = 0 ! Set=0 for flat Zeff if UFILE not read in [ - ] !XZEFFI = 1.0 ! Manually set "composition" Zeff to this value [ - ] NLZFIM =.false. ! .false. to use composition Zeff for resistivity [ - ] !XPZEFRAT = 1.0 ! If NLZFIM =.true., resistivity Zeff shape control [ - ] !XPZEFF = 0.0 ! If NLZFIM =.true., resistivity Zeff shape control [ - ] !XZEFIM = 1.0 ! Manually set "resistivity" Zeff to this value [ - ] !------------------------------------------------------------------------------------ ! Example, options: !------------------------------------------------------------------------------------ ! ! eg. Set composition Zeff to 1.0 and resistivity Zeff to 1.1 (both flat): ! ! NLZEFM =.false. ! NLZFI2 =.false. ! NZEFMOD = 0 ! XZEFFI = 1.0 ! NLZFIM =.true. ! XPZEFRAT = 1.0 ! XPZEFF = 0.0 ! XZEFIM = 1.1 ! ! Other options: ! ! NLZFI2 = .true., don't set NZEFMOD or XZEFFI : ! Read in 2D Zeff data for composition ! ! NLZFIM = .false., don't set XPZEFRAT, XPZEFF or XZEFIM : ! Use composition Zeff for resistivity !------------------------------------------------------------------------------------ ! Particle balance model and fuelling: !------------------------------------------------------------------------------------ NMODEL = 1 ! Use measured Zeff and quasi-neutrality [ - ] NLTAUP =.false. ! Use .RCY recycling source for edge flux [ - ] NLRCYC =.true. ! -"- [ - ] FRAC = 0.95, 0.05 ! D/H species fractions (if NG=2) [ - ] GFRAC = 0.95, 0.05 ! D/H species fractions for gas puffing [ - ] RFRAC = 0.90, 0.10 ! Ratio of hydrogenic species for recycling [ - ] !------------------------------------------------------------------------------------ ! Ti data + transport calculations: !------------------------------------------------------------------------------------ NKIMOD = 4 ! Ti forward soltn. (unless NLTIPRO=.true.) [ - ] NLXKIE =.false. ! Do NOT set chi_i = chi_e (unless NLTIPRO=.true.) [ - ] NLTIPRO =.true. ! Read experimental Ti data? [ - ] NLTI2 =.true. ! Skip predictive model with NLTIPRO=.true. [ - ] NLTI2TX =.true. ! Interpret T_i(meas) as the impurity temperature [ - ] NLTINC =.true. ! Read in Ti and TZ individually into NCLASS [ - ] TIXLIM = 0.9 ! Max XI of valid input Ti data [ - ] TIFACX = 1.0 ! Edge Ti/Te factor (use with above) [ - ] GIEFAC = 0.0 ! Ti = (1-FIEFAC)*(GIEFAC*Te) + FIEFAC*Ti [ - ] FIEFAC = 1.0 ! Ti = (1-FIEFAC)*(GIEFAC*Te) + FIEFAC*Ti [ - ] ALPH0I = 1.0 ! Convective energy loss multiplier (i) [ - ] ALPH0E = 1.0 ! Convective energy loss multiplier (e) [ - ] !------------------------------------------------------------------------------------ ! Poloidal Field Diffusion model: !------------------------------------------------------------------------------------ ! a) Resistivity model: ! model switches: NLSPIZ NLETAW NLRESTSC NLRES_SAU ! --------------------------------------------------------------- ! Spitzer (default) T F F F ! Spitzer (Sauter vsn) T F F T ! TRANSP neoclassical F F F F ! NCLASS neoclassical F T F F ! TSC neoclassical F F T F ! Sauter neoclassical F F F T NLSPIZ =.false. NLETAW =.false. NLRESTSC =.false. NLRES_SAU =.true. ! b) Bootstrap model: ! model switches: NLBOOT NLBOOTW NLBOOT_SAU ! --------------------------------------------------------------- ! (default: NONE) F T/F T/F ! TRANSP aspect ratio approx. T F F ! NCLASS neoclassical T T F ! Sauter neoclassical T F T NLBOOT =.true. NLBOOTW =.false. NLBOOT_SAU =.true. !------------------------------------------------------------------------------------ !Use only for predictive Ti !XTILIM = 0.9 ! Max r/a for matching conductivity [ - ] !XTIKIM = 0.1 ! Min r/a for matching conductivity [ - ] !------------------------------------------------------------------------------------ ! Plasma rotation: !------------------------------------------------------------------------------------ NLVPHI =.true. ! Turn on plasma rotation [ - ] NLVWNC =.true. ! Use NCLASS to analyse rotation data [ - ] NLOMGVTR =.true. ! Set .true. if impurity rotation supplied [ - ] !------------------------------------------------------------------------------------ ! Neutrals model: !------------------------------------------------------------------------------------ NSOMOD = 1 ! 1=FRANTIC [ - ] MOD0ED = 1 ! 1=Use Temperature values given by E0IN [ - ] ! Fuelling species temperatures: E0IN = 20.0, 3.0, 0.55, 20.0, 3.0, 0.55 ! [eV ] FH0ESC = 0.0 ! Fraction of CX neutrals returning as "warm" [ - ] RECYCB = 0.0 ! Fraction of .RCY neutrals returning as "warm" [ - ] NLRECO = .true. ! Turn on "recombination" [ - ] !------------------------------------------------------------------------------------ ! NBI: !------------------------------------------------------------------------------------ !NSIGEXC = 1 ! Enable excited states deposition model [ - ] !------------------------------------------------------------------------------------ ! TRDAT section of the NAMELIST: !------------------------------------------------------------------------------------ $TRDATA LFIXUP = 0 ! Set=0 to ignore labels [ - ] RMJCHK =-1.0 ! Range checking disabled (set < 1) [ - ] PRECUR ='F' ! Prefix for .CUR file [ - ] EXTCUR ='CUR' PRERBZ ='F' ! Prefix for .RBZ file [ - ] EXTRBZ ='RBZ' PREMRY ='M' ! Prefix for .MRY file [ - ] EXTMRY ='MRY' PREVSF ='F' ! Prefix for .VSF file [ - ] EXTVSF ='VSF' PRETER ='F' ! Prefix for .TER file [ - ] EXTTER ='TER' PRENER ='F' ! Prefix for .NER file [ - ] EXTNER ='NER' PRETI2 ='F' ! Prefix for .TIR file [ - ] EXTTI2 ='TI2' PREQPR ='F' ! Prefix for .QPR file [ - ] EXTQPR ='QPR' PREBOL ='F' ! Prefix for .BOL file [ - ] EXTBOL ='BOL' PRENTX ='F' ! Prefix for .NTR file [ - ] EXTNTX ='NTX' PRERCY ='F' ! Prefix for .RCY file [ - ] EXTRCY ='RCY' PREOMG ='F' ! Prefix for .OMG file [ - ] EXTOMG ='OMG' PRENB2 ='P' ! Prefix for .NB2 file [ - ] EXTNB2 ='NB2' PREZF2 ='F' ! Prefix for .ZFT file [ - ] EXTZF2 ='ZF2' $END ! %NBLIST.FOR-- BEAM NAMELIST DATA WRITTEN TO END OF FILE !NBLIST.FOR ! beam constants, configuration = MAST2002 NLBEAM = .T DTBEAM = 0.001 ! beam code timestep NPTCLS = 2000 ! no. of MC ptcls NDEP0 = 500 ! no. of deposition tracks GOOMAX = 5000.0 ! max goosing factor DTN = 0.05e-6 ! orbit timestep control REDGE = 17.0 ! aperture dimensions XZEDGE = 0.0 NBSHAP = 2 ! 1=rectangular source grid FOCLR = 4.E2 ! focal lengths FOCLZ = 4.E2 DIVR = 1.49E-2 ! divergences DIVZ = 1.49E-2 BMWIDR = 15.0 ! source grid dimensions BMWIDZ = 0.0 DN0OUT = 5.E11 ! "external" neutral density NCIRLM = 0 ! circle limiters CRLMR1 = 0.0 CRLMY1 = 0.0 CRLMRD = 0.0 NLINLM = 6 ! line limiters ALNLMR = 20.952,20.952,105.8,199.0,199.0,105.8 ALNLMY = -145.789,144.709,144.709,42.85,-42.49,-145.789 ALNLMT = 90.0,180.0,225.0,-90.0,-45.0,0.0 NLBFLR = .TRUE. ! see p.9 of flr.tex, enables FLR NLBGFLR = .TRUE. ! enables upgrade to the FLR GFLR_MIN = 0.0 ! guiding centre limit GFLR_LL = 30.0 ! max search size of algorithm in a,b, s GFLR_RL = 1.0E6 ! no of particles within this larmor rad ! BEAMLINE 1M %NBLIST.FOR GENERATED OUTPUT NLCO(1) = .T EINJA(1) = 0.000000E+00 PINJA(1) = 0.000000E+00 TBONA(1) = 8.899999E-02 TBOFFA(1) = 1.946000E-01 ABEAMA(1) = 1.000000E+00 XZBEAMA(1) = 1.000000E+00 PDELTA(1) = 1.000000E-03 RTCENA(1) = 7.000000E+01 XLBAPA(1) = 4.690000E+02 XLBTNA(1) = 6.643000E+02 FFULLA(1) = 0.000000E+00 FHALFA(1) = 0.000000E+00 XYBAPA(1) = 0.000000E+00 XYBSCA(1) = 0.000000E+00 NBAPSHA(1) = 2 RAPEDGA(1) = 1.700000E+01 ! BEAMLINE 2M %NBLIST.FOR GENERATED OUTPUT NLCO(2) = .T EINJA(2) = 0.000000E+00 PINJA(2) = 0.000000E+00 TBONA(2) = 8.899999E-02 TBOFFA(2) = 2.172000E-01 ABEAMA(2) = 2.000000E+00 XZBEAMA(2) = 1.000000E+00 PDELTA(2) = 1.000000E-03 RTCENA(2) = 7.000000E+01 XLBAPA(2) = 4.690000E+02 XLBTNA(2) = 6.643000E+02 FFULLA(2) = 0.000000E+00 FHALFA(2) = 0.000000E+00 XYBAPA(2) = 0.000000E+00 XYBSCA(2) = 0.000000E+00 NBAPSHA(2) = 2 RAPEDGA(2) = 1.700000E+01 NBEAM= 2 ! NBLIST.FOR: NO. OF BEAMS ! NBLIST.FOR: 20% .GT. MAX(INJECTION ENERGY): EBDMAX= 4.911562E+04 NLBDAT=.T ! ** BEAM DATA UFILE FLAG **