Wed Feb 7 14:52:43 CET 2007 lrwxrwxrwx 1 mab users 56 2007-02-07 14:42 /home/mab/toric5_svn/base/Linux/Ser/Storic.e -> /home/mab/toric5_svn/src/../base/Linux/Ser/Storic_f95f.e macmab 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 /home/mab/runs_toric5/newInput/equigs_19314 Interpolation is made with 7 Chebichev polynomials n and T profiles read from file /home/mab/runs_toric5/newInput/profnt_19314 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 = 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 Global power balance Power using the Parseval formula Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 5.7661E-02 MW/kA^2 ( 18.89 %) Power to ions 2 fundamental 1.9781E-01 MW/kA^2 ( 64.82 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 4.9700E-02 MW/kA^2 ( 16.29 %) - from the fast wave 4.6547E-02 MW/kA^2 ( 15.25 %) - from the ib wave 3.1533E-03 MW/kA^2 ( 1.03 %) Total absorbed power 3.0517E-01 MW/kA^2 Total radiated power (J*E) 3.1173E-01 MW/kA^2 Total radiated power (Poynting) 3.1247E-01 MW/kA^2 - towards the wall 5.0078E-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. 5.7661E-02 MW/kA^2 ( 18.90 %) Power to ions 2 fundamental 1.9781E-01 MW/kA^2 ( 64.82 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 4.9694E-02 MW/kA^2 ( 16.28 %) - from the fast wave 4.6540E-02 MW/kA^2 ( 15.25 %) - from the ib wave 3.1533E-03 MW/kA^2 ( 1.03 %) Total absorbed power 3.0517E-01 MW/kA^2 Accuracy of the 2dim abs. power -2.2335E-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.3937E-01 (A/W) Current per incident power(NT) 3.8009E-01 (A/W) Current drive figure of merit 1.9000E-01 (A/W/M^2) Current drive figure of merit(NT) 3.0170E-01 (A/W/M^2) Current drive estimate, fast Wave: RF current per incident power 2.3603E-01 (A/W) Current per incident power(NT) 3.7421E-01 (A/W) Current drive figure of merit 1.8735E-01 (A/W/M^2) Current drive figure of merit(NT) 2.9703E-01 (A/W/M^2) Current drive estimate, IB Wave: RF current per incident power 3.3328E-03 (A/W) Current per incident power(NT) 5.8815E-03 (A/W) Current drive figure of merit 2.6455E-03 (A/W/M^2) Current drive figure of merit(NT) 4.6685E-03 (A/W/M^2) The output has been stored in netcdf format toric5 was called 1 times, total 318.74 secs gloini was called 1 times, total 4.75 secs fftl_ini was called 1 times, total 4.56 secs equil was called 1 times, total 0.03 secs mhdcoe was called 9150 times, total 0.16 secs fftl_inv was called 2870621 times, total 18.94 secs profnt was called 1 times, total 0.00 secs coords was called 1919 times, total 0.47 secs fftl_dir was called 27558 times, total 0.27 secs rescof was called 1 times, total 0.15 secs solver was called 1 times, total 261.31 secs zetaRe was called 6585377 times, total 37.56 secs eldibw was called 39401 times, total 49.16 secs direl_flr was called 39401 times, total 20.95 secs cibess was called 1876958 times, total 15.17 secs elepwr was called 1 times, total 52.57 secs elefld was called 453 times, total 0.56 secs powtwd was called 453 times, total 50.10 secs Wed Feb 7 14:58:02 CET 2007