toricmode is toric Reading mhd equilibrium from file ./equigs_19314 Interpolation is made with 7 Chebichev polynomials n and T profiles read from file ./profnt_19314 Profiles nprodt, nspec, mainsp, kdiff_idens, kdiff_itemp 201 2 1 1 0 profnt: # radial pts = 201 profnt: # species = 2 profnt: # ion temp = 1 profnt: species ( 2, 1)_( 1, 1)_ t4_mod_profnt Ne: n_e t4_mod_profnt Te: T_e t4_mod_profnt Ni: 1 n_i t4_mod_profnt Ti: 1 T_i 1 t4_mod_profnt Ni: 2 n_i 2 TOROIDAL IC FULL-WAVE CODE TORIC Id string: t4_output.F,v 1.15 2006/08/23 22:03:54 jwright Exp Solution of the wave equations TOKAMAK PARAMETERS: Major radius = 166.896 cm Magnetic Axis maior radius = 172.645 cm Separatrix radius = 43.474 cm Plasma radius (edge of s.o.) = 43.474 cm Faraday Shield omitted = Antenna radius = 51.474 cm Vacuum vessel radius = 61.474 cm Central magnetic field = 1.971 Tesla Toroidal plasma current = 836.058 kA Safety factor q on axis = 1.008 Safety factor q at plasma edge = 2.924 Equilibrium from ntcc file: highest poloidal mode retained 5 Shift (axis, sep., wall) = -5.750 0.000 -5.136 cm Ellipticity = 1.348 1.506 1.419 Triangularity = 0.000 0.047 0.054 Vertical shift = 4.811 0.000 -2.793 cm Vertical Triangularity = 0.000 -0.034 0.020 Skewdness = 0.000 1.307 Psi_pol at last closed surface = 0.900 PLASMA PARAMETERS Central electron density = 6.531E+13 cm-3 Central electron temperature = 4.353 keV Central electron beta = 3.125E+00 % Electron density at the separat. = 1.519E+13 cm-3 Electron temperature at the sep. = 0.197 keV Ion Species 1 : Charge (atomic units) 1.0 Mass (atomic units) 2.0 Central density = 6.139E+13 Concentration (100*n_i/n_e) = 94.000 % Central temperature = 4.331 keV Partial beta on axis = 2.923E+00 % f/f_cy: inner ~= 1.55 axis = 2.09 outer ~= 2.64 Fundam. resonance at X = -90.061 cm - outside the plasma on the high-field side Harmonic resonance at X = -7.278 cm tangent to the surface r/a = 0.178 on the high-field side Ion Species 2 : Charge (atomic units) 1.0 Mass (atomic units) 1.0 Central density = 3.919E+12 Concentration (100*n_i/n_e) = 6.000 % Central temperature = 4.331 keV Partial beta on axis = 1.866E-01 % f/f_cy: inner ~= 0.77 axis = 1.05 outer ~= 1.32 Fundam. resonance at X = -7.278 cm tangent to the surface r/a = 0.178 on the high-field side Harmonic resonance at X = 157.691 cm - outside the plasma on the low-field side Ion-ion resonance at X = -14.062 cm tangent to the surface r/a = 0.338 on the high-field side Ion-ion cutoff at X = -11.984 cm tangent to the surface r/a = 0.290 on the high-field side PLASMA MODEL: Ions: finite Larmor radius approximation ELD of IB waves taken into account Parallel electric field taken into account Electrons: Landau damping and TTMP with mixed term Current drive efficiency evaluated with Z_eff = 1.50 Collisions omitted Poloidal magnetic field taken into account, sign = 1 Toroidal broadening of IC resonances taken into account Jiang-Wu-Povinelli formulation in vacuum Ez suppressed in vacuum WAVE AND ANTENNA PARAMETERS Frequency = 30.500 Mhz Toroidal wavenumber NPHI = 12 Equiv. parallel index = 10.873 Equiv. k_parallel (m^-1) = 6.951 Length of the antenna = 50.000 cm (LC)-constant of the antenna = 1.200 ko = 6.3924E-01 lo = 1.5644E+02 ko*R = 1.1036E+00 ko*a = 3.9297E-01 MESH PARAMETERS Number of poloidal modes 31 Number of poloidal mesh points 64 Number of radial mesh points 453 Psi step: 3.145E-03 Number of poloidal modes in vacuum 15 Number of mesh points in vacuum -1 Psi step in vacuum 3.145E-03 Plasma edge at psi =: 1.000 (mesh point 319) outer side: x = 38.481 R = 211.126 inner side: x = -49.955 R = 122.690 Antenna at psi =: 1.182 (mesh point 378) outer side: x = -53.838 R = 118.808 inner side: x = 46.690 R = 219.335 Wall at psi =: 1.415 (mesh point 453) outer side: x = 57.164 R = 229.809 inner side: x = -58.791 R = 113.854 Two-point boundary condition at the antenna Two-point boundary condition at the plasma edge t4_mod_powtwd : edjnor, pwnorm, qtnorm 30.87928070884211 0.1281366870390934 1.343151855755696E-04 Global power balance Power using the Parseval formula Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 2.4372E-02 MW/kA^2 ( 17.89 %) Power to ions 2 fundamental 8.8427E-02 MW/kA^2 ( 64.91 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 2.3437E-02 MW/kA^2 ( 17.20 %) - from the fast wave 2.1856E-02 MW/kA^2 ( 16.04 %) - from the ib wave 1.5809E-03 MW/kA^2 ( 1.16 %) Total absorbed power 1.3624E-01 MW/kA^2 Total radiated power (J*E) 1.3899E-01 MW/kA^2 Total radiated power (Poynting) 1.3916E-01 MW/kA^2 - towards the wall 5.5429E-08 MW/kA^2 Powers from the explicit theta integration Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 2.4372E-02 MW/kA^2 ( 17.89 %) Power to ions 2 fundamental 8.8427E-02 MW/kA^2 ( 64.91 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 2.3433E-02 MW/kA^2 ( 17.20 %) - from the fast wave 2.1852E-02 MW/kA^2 ( 16.04 %) - from the ib wave 1.5809E-03 MW/kA^2 ( 1.16 %) Total absorbed power 1.3623E-01 MW/kA^2 Accuracy of the 2dim abs. power -2.5570E-05 Total power in W (MW) for 1 A (1 kA) in the antenna numerically equal to the loading resistance in Ohms Electric field in V/m, magnetic field in Tesla, power density in W/cm3, for 1 MW coupled to the plasma Volume average electron density 4.7560E-01 (10^20 M^3) Current drive estimate, total: RF current per incident power 2.5412E-01 (A/W) Current per incident power(NT) 4.1024E-01 (A/W) Current drive figure of merit 2.0171E-01 (A/W/M^2) Current drive figure of merit(NT) 3.2563E-01 (A/W/M^2) Current drive estimate, fast Wave: RF current per incident power 2.5023E-01 (A/W) Current per incident power(NT) 4.0351E-01 (A/W) Current drive figure of merit 1.9862E-01 (A/W/M^2) Current drive figure of merit(NT) 3.2029E-01 (A/W/M^2) Current drive estimate, IB Wave: RF current per incident power 3.8932E-03 (A/W) Current per incident power(NT) 6.7278E-03 (A/W) Current drive figure of merit 3.0903E-03 (A/W/M^2) Current drive figure of merit(NT) 5.3402E-03 (A/W/M^2) The output has been stored in netcdf format toric5 was called 1 times, total 296.28 secs gloini was called 1 times, total 5.38 secs fftl_ini was called 1 times, total 5.21 secs equil was called 1 times, total 0.04 secs mhdcoe was called 9150 times, total 0.15 secs fftl_inv was called 2870621 times, total 10.56 secs profnt was called 1 times, total 0.00 secs coords was called 1919 times, total 0.25 secs fftl_dir was called 27558 times, total 0.14 secs rescof was called 1 times, total 0.12 secs solver_stri was called 1 times, total 263.38 secs zetaRe was called 6585377 times, total 10.61 secs eldibw was called 39401 times, total 33.13 secs direl_flr was called 39401 times, total 14.48 secs add_block_row was called 453 times, total 137.01 secs writeGk was called 452 times, total 3.01 secs cibess was called 1876906 times, total 13.17 secs stri_solve was called 1 times, total 2.38 secs readGk was called 452 times, total 1.84 secs elepwr was called 1 times, total 27.00 secs elefld was called 453 times, total 1.08 secs powtwd was called 453 times, total 24.68 secs