&toricainp ! ! Read from subroutine hfread, "contained" in t4_torica ! isol = 1, ! mode of operation of the TORIC code ! = 0 Explorative(diagnostic) mode: ! Magnetic configuration ! Density & temperature profiles, ! Dielectric properties vs X ! Dielectric properties vs TH (IPLTHT = 1) ! = 1 Run mode, Maxwellian plasmas ! Construction and inversion of the stiffness ! matrix; results are printed, plotted, and ! written on disk if required. ! = 2 Run mode, with non-maxwellian minority ! In this case, the namelist /qlmininp/ must specify ! the minority species and the file where the ! data from the run of SSFPQL are to be read. ! ! Note: always check vacuum extrapolation with ! isol = -1 before using a new configuration! nvrb = 3, ! -------------------------------------------------------------------- ! Definition of the variables and mesh for the run ! -------------------------------------------------------------------- ! Number of components of the electric field: ! NVRB = 2 - parallel field omitted (for comparison ! with old codes omitting E_//.) ! NVRB = 3 - parallel field evaluated (normal option) nelm = 450, ! Number of radial finite elements ntt = 512, ! Number of points in the poloidal angle (theta) mesh: ! must be a power of 2 (FFT) and satisfy NTT >= 2*NMOD+1 nmod = 255, ! Number of modes in the solution (preferably odd): ! NOTE: to optimize use of the information in the poloidal ! mesh the best choice is NMOD = NTT/2-1 (see below). ! A value of NMOD << NTT/2-1 can lead to inaccuracies. ! ! Example: With MASTCH=1 the largest job accepted by the ! IPP IBMaix or a normal laptop has (about): ! NTT = 64, NMOD = 31, NELM = 450. ! ! nelm = 200, ntt = 256, nmod = 127 can perform well ! on the TianHe cluster. ! -------------------------------------------------------------------- nptvac = 50, mxmvac = 15, ! It is recommented not to exceed MXMVAC = 15 ! -------------------------------------------------------------------- ! H.F. and antenna parameters ! -------------------------------------------------------------------- freqcy = 16.0E6, ! Applied frequency, in hz. nphi = 18, ! Toroidal mode number. antlen = 50, ! Length of the antenna in the poloidal direction (cm) antlc = 1.2, ! Effective propagation constant of the antenna. ! ANTLC > 0 - the current ends abruptly at the feeders ! ANTLC < 0 - the current goes to zero at the feeders ! (works only with a large value of ntt). ! ANTLC = -999. - poloidally uniform excitation theant = 0., ! poloidal position of the central point of the antenna ! (degrees from the outer equatorial plane) ! -------------------------------------------------------------------- iflr = 1, ibweld = 2,
ibpol = 1, iqtor = 1, icoll = 0, enhcol = 1., ! Must be very large (10^5) to have any effect imdedg = 2, iezvac = 1, icosig = 0, ! Use ICOSIG <> 0 only in exceptional cases! iregax = 1, ! IREGAX=1 and intchb > 0 recommended when ! IBW are propagative near the magnetic axis mastch = 1, lenwrd = 8, ! ---------------------------------------------------------------------- ! Input of parameters for the control of the output: ! ---------------------------------------------------------------------- iout = 1, ! printed output control: ! = 0 normal output. ! = 1 detailed output. ! = 2 complete output: includes radial behaviour of the ! Fourier coefficients of the electric field. iwdisk = 1, ! output on disk (fort.9, only with ISOL=1) ! = 0 skipped ! = 1 done (specify file name in the JCL job) idlout = 1, ! output for IDL plots ! IDLOUT = 0 skipped ! IDLOUT > 0 executed (normal setting) ipltht = 1, ! plot of dielectric tensor as function of theta ! on selected magnetic surfaces io_ncdf = 1, ! The output for plots is made in netcdf format ! At present only available for ISOL = 1 or 2 ! ---------------------------------------------------------------------- zeff = 1.5, ! If Zeff is >= 1._r8 it will be used throughout ! < 1._r8 it will be reevaluated. iclres = 0 ! iclres = 1 adds collisions around isolated ion-ion ! resonances. Use with care! ! dnures = Width (cm) / tnures = strength (a.u.) dnures = 5.0 tnures = 10.0 ! scratchpath = '/tmp/' timing_on = .true. pcblock = 4 /
! ------------------------------------------------------------------------ ! The following data are used only with a default analytic equilibrium ! (igsmhd = 0). ! rtor = 167.0 ! Major radius (cm) rplasm = 44.0 ! Plasma radius (cm) ! ashift = 5.0 ! Shafranov Shift aellip = 1.00 ! Ellipticity on axis sellip = 1.60 ! Ellipticity at separatrix strian = -0.22 ! Triangularity at separatrix iudsym = 1 azvert = 0.0 strigz = 0.0 sp_xs1 = 0.0 ! bzero = 1.97 ! Magnetic field at the geometric center aicurr = 830. ! Toroidal current (kA) ppjte = 2.0 ppjti = 1.5 ! ! ----------------------------------------------------------------------- ! The following data are used only with analytic profiles (iprodf = 0) ! nspec = 2 ! Number of ion species ! Should not exceed parameter nspmx mainsp = 1 ! Reference species (for charge neutrality)
denec = 5.0E13 ! Central electron density (cm^-3) tempec = 4 ! Central electron temperature (keV)
tempic(1) = 3.5 ! Central ion temperature (keV) tempic(2) = 3.5 ! tempic(3) = 8.0
ppnei = 2.0 ! Inner power in the density profile ppnee = 0.500 ! Outer power in the density profile
pptei = 2.0 ! Inner power in the electron temperature profile pptee = 0.5 ! Outer power in the electron temperature profile
pptii(1) = 2.0 ! Inner power in the ion temperature profile pptie(1) = 0.5 ! Outer power in the ion temperature profile pptii(2) = 2.0 pptie(2) = 0.5 ! pptii(3) = 4.0 ! pptie(3) = 1.5
dnsepr = 0.5E13, ! Density at the last closed surface tesepr = 0.20, ! Electron temperature at the last closed surface tisepr(1) = 0.20, ! Ion temperature at the last closed surface tisepr(2) = 0.20, ! tisepr(3) = 10.,
so_thickness = 0., ! Scrape-off thickness (cm) gldn = 3.0, ! Density gradient length in scrape-off glte = 3., ! Electron temperature gradient length in scrape-off glti(1) = 3., ! Ion temperature gradient length in scrape-off glti(2) = 3., /
gfreadinp
用于设置 gfile 文件的位置, 只在 toricmode = "equil" 时会被读取.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
&gfreadinp ! ! Read from subroutine gfread in module t4_mod_toi2mex ! gfile= 'H2DNv3_2g033332.00000' ! nmhd = 165 ! number of radial psi points default=201 ntheta = 129 ! number of "theta" points default = 129
i2mex_LAST_NORM_SURFACE_IS = 0.95 ! ! Name speaks for itself ic1 = 2 ! index of first (non-singular) surface inc = -1 ! index of last surface ! ! (-1: makes it equal to npsi) imom1 = 10 ! Number of Fourier modes per mag. surface /
var_name (for legibility) nprodt number of points in the radial variable psi nspec number of ion species mainsp species to be used to impose charge neutrality kdiff_itemp Control of individual ion temperature profiles = 0 - All ion species have the same temperature (the filecontainsonly one T_i profile) = 1 - Each ion species has a different temperature (the file must contain NSPEC t_i profiles)
var_name (for legibility) nprodt number of points in the radial variable psi nspec number of ion species mainsp species to be used to impose charge neutrality kdiff_idens Control of individual ion density profiles kdiff_itemp Control of individual ion temperature profiles
! ********************************************************************* ! Namelist entries needed by SSFPQL (quasilinear solver for the ions) ! ********************************************************************* ! ! namelist /ssfpqlinp/ ! & power, iex_sumnph, joutql, nrdpsi, legmod, idlout, iwrdql, ! & path_to_t4data, t4data_file, timing_on ! ! -------------------------------------------------------------------- ! &ssfpqlinp ! ! Read by subroutine t4_ssfpql in module t4_mod_ssfpql ! power = 3.6, ! Total absorbed power (MW) emxfac(1) = 200., emxfac(2) = 100., ! multipliers of T_i in the definition ! of the max energy in the velocity mesh iex_sumnph = 0, ! Control of input from TORIC ! = 0 - only one toroidal mode ! = 1 - weighted superposition of several toroidal modes ! (performed by t4_sumnph, with its own namelist) joutql = 1, ! Control of printed output: ! = 0 - essential results (recommended) ! = 1 - detailed results nrdpsi = 46, ! Number of radial points (uniform in psi) where ! the qlfp equation for ions is to be solved. legmod = 6, ! number of Legendre polynomials in the ion ! distribution functions ! Recommended: 5 <= legmod <= 7 nptu = 201,
idlout = 1, ! Plots of distribution functions, power profiles, etc. ! = 1 - done ! = 0 - skipped. iwdisk = 1, ! 1 - netCDF format ! -1 - binary path_to_t4data = './' ! output of toricmode = "toric" (named "fort.9") t4data_file = 'fort.9' iwrdql = 1, ! Output for iteration by TORIC with nonMaxwellian ! minority distribution. timing_on = .true. ! Control of time trace: ! .true. - time tracing on ! .false. - time tracing off / ! -------------------------------------------------------------------- !