The number of poloidal Fourier modes should not exceed half the number of points in the poloidal mesh minus one: NMOD -> NTT/2-1 Reading mhd equilibrium from file ./equigs_19314 Interpolation is made with 7 Chebichev polynomials n and T profiles read from file ./profnt_19314a t4_mod_profnt : Profiles nprodt, nspec, mainsp, kdiff_itemp 201 2 1 0 t4_mod_profnt Ne: n_e t4_mod_profnt Te: T_e t4_mod_profnt : Ni 2 n_i 2 t4_mod_profnt : Ti 1 T_i 1 TOROIDAL IC FULL-WAVE CODE TORIC $Id: t4_output.for,v 1.7 2006/06/20 14:44:22 Doug_McCune 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 Jaing-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 = 90.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 : pwanor, elenor, pwnorm, qtnorm 2.039358713369091E-02 376.9911184307751 8.472222222222223E-04 2.136196186865549E-05 Global power balance Power using the Parseval formula Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 4.2715E-01 MW/kA^2 ( 18.39 %) Power to ions 2 fundamental 1.5162E+00 MW/kA^2 ( 65.28 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 3.7939E-01 MW/kA^2 ( 16.33 %) - from the fast wave 3.5419E-01 MW/kA^2 ( 15.25 %) - from the ib wave 2.5194E-02 MW/kA^2 ( 1.08 %) Total absorbed power 2.3228E+00 MW/kA^2 Total radiated power (J*E) 2.3777E+00 MW/kA^2 Total radiated power (Poynting) 2.3799E+00 MW/kA^2 - towards the wall 1.9030E-07 MW/kA^2 Powers from the explicit theta integration Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 4.2715E-01 MW/kA^2 ( 18.39 %) Power to ions 2 fundamental 1.5162E+00 MW/kA^2 ( 65.28 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 3.7934E-01 MW/kA^2 ( 16.33 %) - from the fast wave 3.5414E-01 MW/kA^2 ( 15.25 %) - from the ib wave 2.5194E-02 MW/kA^2 ( 1.08 %) Total absorbed power 2.3227E+00 MW/kA^2 Accuracy of the 2dim abs. power -2.2920E-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 Gauss, 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 9.7473E-03 (A/W) Current per incident power(NT) 1.5901E-02 (A/W) Current drive figure of merit 7.7370E-03 (A/W/M^2) Current drive figure of merit(NT) 1.2622E-02 (A/W/M^2) Current drive estimate, fast Wave: RF current per incident power 9.5997E-03 (A/W) Current per incident power(NT) 1.5643E-02 (A/W) Current drive figure of merit 7.6199E-03 (A/W/M^2) Current drive figure of merit(NT) 1.2417E-02 (A/W/M^2) Current drive estimate, IB Wave: RF current per incident power 1.4755E-04 (A/W) Current per incident power(NT) 2.5841E-04 (A/W) Current drive figure of merit 1.1712E-04 (A/W/M^2) Current drive figure of merit(NT) 2.0512E-04 (A/W/M^2) The output has been stored in netcdf format %subroutine t4_tofpp: call t4_tofpp_drive... Toric data will be read from file ./toric.data.1001 Toric data read from file ./toric.data.1001 The data have been read succefully Reading mhd equilibrium from file ./equigs_19314 Interpolation is made with 7 Chebichev polynomials n and T profiles read from file ./profnt_19314a t4_mod_profnt : Profiles nprodt, nspec, mainsp, kdiff_itemp 201 2 1 0 t4_mod_profnt Ne: n_e t4_mod_profnt Te: T_e t4_mod_profnt : Ni 2 n_i 2 t4_mod_profnt : Ti 1 T_i 1 Inconsistency detected in the power from toric Total power = 2.4053E+01 to electrons = 3.9176E+00 to ions 1 f = 0.0000E+00 h = 4.5435E+00 to ions 2 f = 1.5592E+01 h = 0.0000E+00 TOROIDAL IC FULL-WAVE CODE TORIC Version 4.10 - Sept. 2004 Preliminary interface to ffpmod TOKAMAK PARAMETERS: Major radius = 166.896 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 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 PLASMA PARAMETERS Central electron density = 6.531E+13 cm-3 Central electron temperature = 4.353 keV Electron density at the separat. = 0.152E+14 cm-3 Electron temperature at the sep. = 0.197 keV PLASMA MODEL: Ions: finite Larmor radius approximation ELD of IB waves taken into account Electrons: Landau damping and TTMP with mixed term Collisions omitted Poloidal magnetic field taken into account, sign = 1 Toroidal broadening of IC resonances taken into account Jaing-Wu-Povinelli formulation in vacuum Ez suppressed in vacuum WAVE AND ANTENNA PARAMETERS Frequency = 30.500 Mhz Half-length of the antenna = 90.000 cm (LC)-constant of the antenna = 1.200 MESH PARAMETERS Number of poloidal modes 31 Number of poloidal mesh points 64 Number of radial mesh points 453 Fields and power deposition will be evaluated for a launched power of 4.000 MW Summary: total power to electrons (W) 1.5597E+07 to ions 1 fund. 0.0000E+00 harm. 1.8040E+07 to ions 2 fund. 6.1916E+07 harm. 0.0000E+00 Total accounted power (MW) 9.5553E+01 Nominal power 4.0000E+00 Total accounted rf curr. (A) 1.4651E+07 Nominal current 9.5747E+05 toric was called 1 times, total 191.56 secs gloini was called 1 times, total 5.53 secs fftl_ini was called 1 times, total 5.41 secs equil was called 1 times, total 0.03 secs mhdcoe was called 9236 times, total 0.41 secs fftl_inv was called 2870621 times, total 5.64 secs profnt was called 1 times, total 0.00 secs coords was called 1919 times, total 0.24 secs fftl_dir was called 27558 times, total 0.24 secs rescof was called 1 times, total 0.08 secs solver was called 1 times, total 159.62 secs zetaRe was called 6585377 times, total 12.97 secs eldibw was called 39401 times, total 18.21 secs direl_flr was called 39401 times, total 9.42 secs cibess was called 1876932 times, total -8.34 secs elepwr was called 1 times, total 20.36 secs elefld was called 453 times, total 0.41 secs powtwd was called 453 times, total 18.74 secs