!List of files which contain default values for FPP namelist variables: $fpdflt deffiles = 'fpbasic.dft' $end $FPPIN COMMENTS = 'Ignitor D(He-3), VV moments, R2=0', 'Changes from ignpao22.in:', ' Carbon, with Zeff=1.3, lowered n_he3/ne=.10->.07', ' increased kpar from 0.05 to 0.08', ' reduced nshell from 20 to 10 to speed up...' ! Multiple ion species information: ! The first ion in this list is assumed to be the "main" ion species, ! for which FRANTIC will calculate a neutral density. nions=4 !actual number of ions species used nAions=2,3,3,12 !list of Atomic mass A_i of each ion species nZions=1,2,1,6 !list of Atomic charge Z_i of each ion species namions='D','He3','T','C' !label for each ion species nfasti=2 !index i of the fast ion to be Fokker-Planck'd ! nfasti is the ion species for which we will solve the bounce-averaged ! Fokker-Planck/Quasilinear equation. This could be the fast ions ! produced by neutral beam injection or ICRF minority heating, or could ! be the Maxwellian plus energetic tail for second harmonic ICRF ! heating of a majority ion species. ! ! If Te0>0 then FPP will use parabolic-to-a-power formulas for the ! temperature and density profiles. The parabolic formulas are of ! the form: ! Te(r) = (Te0-Tea)*(1-r**2/a**2)**Teexp + Tea ! Te0 = 20.0e3 !central Te(0) for parabolic formula (eV) Tea = 0.5e3 !edge Te(a) for parabolic formula (eV) Teexp = 1.0 !exponent for parabolic formula for Te Delec0 = 3.0e14 !central ne(0) for parabolic formula (/cm**3) Deleca = 1.e13 !edge ne(a) for parabolic formula (/cm**3) Delecexp= 1.0 !exponent for parabolic formula for ne Ti0=4*15.0e3 !list of central ion temperatures (eV) Tia=4*0.5e3 !list of edge ion temperatures (eV) Tiexp=4*1.0 !list of ion temperature parabolic exponents rations=0.75,0.07,0.05 ! fraction n_i(r)/n_e(r) for each ion species rZeff=1.30 ! >0 used to set the impurity and main ion fractions nshell=10 !# of radial shells in FPP calculation (10=default) !impurity info: AN00 = 1.e09 TEDG = 0.2500000 , !machine parameters: RMIN = 47. , RMAJ = 132.0 , RINWAL = 85.0, FPELLIP = 1.85, FPSHIFT0= 0.0, VLOOPFP = 0.2 , BTOR = 1.3e5 PLACUR = 1.2e4 , !Plasma current (KiloAmps) QAX = 0.8, !beam parameters: NUMNBI = 0.0 RBEAMS = 34.80000 , 11*0.0000000E+00, VBEAMS = 48000.00 , 11*0.0000000E+00, PBEAMS = 12*0.0000000E+00, BMIXFULL = 0.2650000 , 11*0.0000000E+00, BMIXHALF = 0.3290000 , 11*0.0000000E+00, TRISE = 1.0000000E-03, TOFF = -1.0000000E-03, FALLR = 1000.000 , FALLZ = 1000.000 , BMEDGHW = 5.000000 , BMEDGHH = 20.00000 , ! description of CX analyzers: LDIR = -1 NLINE = 0 ! # of different Rtan's for CX analyzers TANS = 13.00000 , 70.00000 , 102.0000 , 47*0.0000000E+00, RCRSFP = 185.0000 , PHICRSFP = 0.0000000E+00, Xmax=250.e3 !maximum energy on CX f(E) plot Ymin=1.e-4 !minimum Y (relative to Maximum Y) on CX f(E) plot ELIM2 = 250000.0 , BKGMAX = 0.0000000E+00, NENCX = 0, NSHCX = 0, LLNFPL = 1, FMIN = 28.00000 , FMAX = 38.00000 , EMAX = 250000.0 , LCXRTP = 0, !----------------------------------------------------------------------- ! ICRF parameters: ! prftot = 18.0e6 ! total input rf power (watts) FREQRF = 132.e6 ! frequency of RF generator in Hz RKPAR = 0.08 ! abs(kpar) (+/- used in code) (cm^-1) NHARM = 1 !1=fundamental 2=2cd harmonic ICRF heating CONCMI = 0.10 ! nmin/ne for the fundamental minority species LSHOOT = 1 !.1. = Smithe gen'l geometry code, 0=hand guess ! LSHOOT = 1 requires the following parameters ! (also, SPRUCE currently assumes two out-of-phase antenna straps, ! on the low field side). ! antRcen=184.0 ! major radius R of the antenna center (cm) antHH=38.0 ! antenna half height (cm) ! The antenna is assumed to be curved to lie parallel to the ! vacuum vessel. Its length (along its curvature) may be shorter ! than its height (along a straight line from tip to tip). antwidth=12.5 !toroidal width of each antenna (cm), single strap antsep=110. !toroidal separation between the two antennas (cm), next port ! FPP/SPRUCE currently assumes that the antenna current has a ! poloidal dependence of cos(theta/theta0), so the current ! goes to zero at the tips like a simple half-wave antenna. !************************************************************* !VV shape description in Fourier moments (in cm) for i=0 to mxvvmoms: ! (mxvvmoms=5 at the moment). ! ! A Lao-Hirschman type fit to the CMOD VV: VVRmoms(0)=132.0,47.0,0.0,0.,0. ! R(th) = Sum VVRmoms(i) cos(i th) VVZmoms(0)= 0.,87.0,-9.8,0.,0. ! Z(th) = Sum VVZmoms(i) sin(i th) ! VVZmoms(0) must be 0.0 ! ! For comparison, the Lao-Hirschman 2-moment representation is of ! the form: ! ! R(th) = R0 + a cos(th) + R2 cos(2 th) ! Z(th) = E ( a sin(th) - R2 sin(2 th) ) ! ! Describing the VV shape with the 4 parameters Rmin, Rmax, ! Rtop, and Ztop (these last two give the position of the top of ! the VV) we can calculate the Lao-Hirschman coefficients via: ! ! a = (Rmax-Rmin)/2 ! Rx = (Rmax+Rmin)/2 ! ! The next two equations can be solved by combining them together ! to make a cubic equation. Or for small d they can be solved ! iteratively by using R2=0.0 as an initial guess: ! ! d = (Rx-R2-Rtop)/a ! a measure of "D-ness" or triangularity ! R2 = 3 a d / (9 - 8 d**2) ! ! Once the above two equations are solved, we can then find: ! ! R0 = Rx-R2 ! E = 3 (Ztop/a) (9 - 8 d**2) / (9 - 4 d**2)**(3/2) ! !**************************************************************** ! LSHOOT = .f. requires the following guesses at ICRF fields: !cp i LRFmod ERFpro is actually 1=E-plus 2=Power profile shape !cp r ERFpro(maxshe) E_plus(r) or Prf(r) depending of LRFmod !c ERFpro is just the profile shape. Its value gets normalized to !prftot. !cp i LKmodel k_perp model (0= cold plasma k_perp 2=input rkperp(r)) !cp r rkperp RF k-perp(r) profile (cm^-1) !c LRFMOD = 1, Erfpro=20*1.0, Lkmodel=0 ! ! Time (all in seconds) step information: TSTART = 0.0 ! start time of calculation TCALC = 0.5 ! finish time of calculation DTFPP = 0.01 ! time step of main FPP output NMULTI = 2 ! (half) # of FP time steps to take each DTFPP ! Energy grid specification: LGRID = 0, EGRMAX = 200000.0 , EMINEV = 0.0000000E+00, DEOVE = 0.1000000 , DEMIN = 1000.0000 , DEMAX = 1000000. , NSKIPCX = 2, NPLTE = 1, NPLTT = 2, NIANGLS = 45, 35, 25, 16, 6, EFFOA = 100000.0 , 60000.00 , 30000.00 , 10000.00 , 500.0000 , LPlotE = 1, LPlotA = 1, LPlotF = 1, LFSAVE = F, LFGET = F, CXONLY = F, lcxsrc=.t. !turn on the source part of the CX operator lnsrc=.t. !turn on a source to maintain the specified density !c !c------------------------------------------------------------------------ !c Fast Ion Radial transport due to asymmetric drag and ripple: !c !c FPP will now calculate radial transport of fast ions. The two !c important components of this transport are a convective velocity v_r due !c due to the asymmetric edge drag mechanism !c discussed in my thesis and my 1986 APS invited talk, and D_r due to the !c Goldston and Towner banana drift ripple diffusion process. !c !cp r convmult convection multiplier (1=standard, 0=off) !cp r XRipMult ripple diffusion multiplier (1=standard,0=off) !c !c parameters for exponential ripple model: !c delta(R,y) = Ripdmin * exp( sqrt((R-RipRmin)**2+y**2)/Ripexpl ) !cp r RipRmin major radius of minimum ripple (cm) !cp r Ripdmin ripple (peak to average) at that point !cp r RipexpL exponential scale length increase (cm) !cp i NumTF # of TF coils !c !c Rob's fit to PLT ripple: !c data RipRmin /112/ !major radius of minimum ripple (cm) !c data Ripdmin /1.1e-5/ !ripple (peak to average) at that point !c data RipexpL /9.97/ !exponential scale length increase (cm) !c data NumTF /18/ ! # of TF coils ! radial transport model: MRLOOP = 10 !Mrloop .ne. 0 to do RF feedback of Power. But turn off all transport ! models because: 1. don't have CIT ripple parameters, ! 2. apparently need to include GWB stochasticity thresholds on ripple ! transport, and 3. neoclassical transport not debugged yet! CONVMULT = 0.0, XRIPMULT = 0.000000 , XneoMult = 0 RIPRMIN = 228.0000 , RIPDMIN = 1.8950001E-05, RIPEXPL = 19.17000 , NUMTF = 20 $END