toricmode is toric 01Jul2008: Toric 5.3_1 updated 07-Jul-2008 Reading mhd equilibrium from file /home/pshare/Toric/tests/data/equigs_19314 Interpolation is made with 7 Chebichev polynomials n and T profiles read from file /home/pshare/Toric/tests/data/profnt_19314 TOROIDAL IC FULL-WAVE CODE TORIC Id string: t4_output.F 476 2008-06-30 10:10:22Z mab Solution of the wave equations TOKAMAK PARAMETERS: Major radius = 166.896 cm Magnetic Axis major 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 0.000 cm Vertical Triangularity = 0.000 -0.034 0.020 Skewdness = 0.000 1.307 Psi_pol at last closed surface = 0.900 Total plasma volume = 9.257 m^3 Area of the cross section = 5.592 m^2 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.983 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 Global power balance Power using the Parseval formula Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 6.0658E-02 MW/kA^2 ( 17.68 %) Power to ions 2 fundamental 2.2361E-01 MW/kA^2 ( 65.16 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 5.8908E-02 MW/kA^2 ( 17.17 %) - from the fast wave 5.4824E-02 MW/kA^2 ( 15.98 %) - from the ib wave 4.0839E-03 MW/kA^2 ( 1.19 %) Total absorbed power 3.4317E-01 MW/kA^2 Total radiated power (J*E) 3.5033E-01 MW/kA^2 Total radiated power (Poynting) 3.5063E-01 MW/kA^2 - towards the wall -4.8484E-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. 6.0658E-02 MW/kA^2 ( 17.68 %) Power to ions 2 fundamental 2.2361E-01 MW/kA^2 ( 65.16 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 5.8899E-02 MW/kA^2 ( 17.16 %) - from the fast wave 5.4815E-02 MW/kA^2 ( 15.97 %) - from the ib wave 4.0839E-03 MW/kA^2 ( 1.19 %) Total absorbed power 3.4316E-01 MW/kA^2 Accuracy of the 2dim abs. power -2.5678E-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.4960E-01 (A/W) Current per incident power(NT) 4.0917E-01 (A/W) Current drive figure of merit 1.9812E-01 (A/W/M^2) Current drive figure of merit(NT) 3.2478E-01 (A/W/M^2) Current drive estimate, fast Wave: RF current per incident power 2.4559E-01 (A/W) Current per incident power(NT) 4.0226E-01 (A/W) Current drive figure of merit 1.9494E-01 (A/W/M^2) Current drive figure of merit(NT) 3.1930E-01 (A/W/M^2) Current drive estimate, IB Wave: RF current per incident power 4.0150E-03 (A/W) Current per incident power(NT) 6.9094E-03 (A/W) Current drive figure of merit 3.1869E-03 (A/W/M^2) Current drive figure of merit(NT) 5.4844E-03 (A/W/M^2) The output has been stored in netcdf format t4_mod_tofppmod: Toric file not read, presumed ready. t4_mod_tofppmod.t4_tofpp_init: power_MW, tpw, zpwsq 4.000000000000000 0.3431720703337321 3414.082632093354 Volume of the plasma = 9.350 m3 Total power according to toric (MW) = 9.9994E-01 Summary: total power to electrons (W) 7.1765E-01 to ions 1 fund. 0.0000E+00 harm. 6.4670E-01 to ions 2 fund. 3.0601E+00 harm. 0.0000E+00 Total accounted power (MW) 4.4245E+00 Nominal power 4.0000E+00 Total accounted rf curr. (A) 4.4284E-01 Nominal current 9.9841E-01 toric5 was called 1 times, total 283.24 secs gloini was called 1 times, total 5.46 secs fftl_ini was called 1 times, total 5.15 secs equil was called 1 times, total 0.09 secs mhdcoe was called 9535 times, total 0.16 secs fftl_inv was called 2870621 times, total 0.78 secs profnt was called 1 times, total 0.06 secs coords was called 2641 times, total 0.27 secs fftl_dir was called 36222 times, total -0.59 secs rescof was called 1 times, total 0.15 secs solver_stri was called 1 times, total 250.95 secs zetaRe was called 6585377 times, total -15.50 secs eldibw was called 39401 times, total 21.93 secs direl_flr was called 39401 times, total 10.23 secs add_block_row was called 453 times, total 116.80 secs writeGk was called 452 times, total 3.22 secs cibess was called 1876958 times, total -0.31 secs stri_solve was called 1 times, total 2.50 secs readGk was called 452 times, total 1.97 secs elepwr was called 1 times, total 26.21 secs elefld was called 453 times, total 0.51 secs powtwd was called 453 times, total 24.36 secs