2896.head The master node of this job is node008 The working directory is /home/jwright/ttt2 The node file is /usr/spool/PBS/aux/2896.head =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- This job runs on the following nodes: node008 node008 node017 node017 node002 node002 node016 node016 =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- This job has allocated 8 nodes =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- command to EXE: /usr/local/mpich-gm.lhf/bin/mpirun -np 8 -machinefile /usr/spool/PBS/aux/2896.head /usr/local/toric/toric5/toric5/base/Linux/Par/Ptoric.e using TMPDIR pwd /home/jwright pwd /home/jwright/ttt2 AAMAIN: I am 0 node008. 8 MPI context: 91 &EQUIDATA IGSMHD=1,INTCHB=7,IDPROF=-1,RTOR=167.0000000000000,RPLASM=45.00000000000000,ASHIFT=5.000000000000000,AELLIP=1.250000000000000,SELLIP=1.600000000000000,STRIAN=-0.2200000000000000,IUDSYM=1,AZVERT=-5.000000000000000,STRIGZ=0.000000000000000E+00,SP_XS1=0.000000000000000E+00,BZERO=2.080000000000000,AICURR=750.0000000000000,PPJTE=2.000000000000000,PPJTI=1.500000000000000,NSPEC=2,MAINSP=1,ATM=2.000000000000000,1.000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,AZI=1.000000000000000,1.000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,ACONC=0.000000000000000E+00,6.000000000000000E-02,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,DENEC=65000000000000.00,TEMPEC=2.500000000000000,TEMPIC=2.000000000000000,2.000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,DNSEPR=18000000000000.00,TESEPR=0.2000000000000000,TISEPR=0.2000000000000000,0.2000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,DENIC=0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,DNISEPR=0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,PPNEI=3.500000000000000,PPNEE=0.5000000000000000,PPTEI=1.500000000000000,PPTEE=1.000000000000000,PPTII=2.000000000000000,2.000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,PPTIE=1.000000000000000,1.000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,PPNII=0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,PPNIE=0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,GLDN=3.000000000000000,GLTE=3.000000000000000,GLTI=3.000000000000000,3.000000000000000,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,SO_THICKNESS=0.000000000000000E+00,GLNI=0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,0.000000000000000E+00,DIST_PLAFARS=0.000000000000000E+00,DIST_PLAANT=8.000000000000000,DIST_PLAWALL=18.00000000000000,35.00000000000000,9.000000000000000,30.00000000000000,INPUTPATH=./ ,EQUIL_FILE=equigs_19314 ,PROFNT_FILE=profnt_19314a ,KDIFF_CONC=0/ Reading mhd equilibrium from file ./equigs_19314 n and T profiles read from file ./profnt_19314a TOROIDAL IC FULL-WAVE CODE TORIC t4_output.F,v 1.8 2005/11/17 12:11:12 cvs Exp Solution of the wave equations TOKAMAK PARAMETERS: Major radius = 166.896 cm Magnetic Axis = 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 PLASMA PARAMETERS Central electron density = 6.531E+13 cm-3 Central electron temperature = 4.353 keV Central electron beta = 2.948E+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 Concentration (100*n_i/n_e) = 94.000 % Central temperature = 4.331 keV Partial beta on axis = 2.757E+00% f/f_cy: inner ~= 1.55 axis = 2.09 outer ~= 2.64 Harmonic resonance at X = -7.278 cm tangent to the surface r/a = 0.173 on the high-field side Ion Species 2 : Charge (atomic units) 1.0 Mass (atomic units) 1.0 Concentration (100*n_i/n_e) = 6.000 % Central temperature = 4.331 keV Partial beta on axis = 1.760E-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.173 on the high-field side Ion-ion resonance at X = -14.062 cm Ion-ion cutoff at X = -11.983 cm 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 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 = 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 AAMAIN: I am 1 node008. 8 MPI context: 91 zgebs2d check 1 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 1 MPI FINALIZE status: 0 AAMAIN: I am 5 node002. 8 MPI context: 91 zgebs2d check 5 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 5 MPI FINALIZE status: 0 AAMAIN: I am 7 node016. 8 MPI context: 91 zgebs2d check 7 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 7 MPI FINALIZE status: 0 AAMAIN: I am 3 node017. 8 MPI context: 91 zgebs2d check 3 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 3 MPI FINALIZE status: 0 AAMAIN: I am 6 node016. 8 MPI context: 91 zgebs2d check 6 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 6 MPI FINALIZE status: 0 AAMAIN: I am 2 node017. 8 MPI context: 91 zgebs2d check 2 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 2 MPI FINALIZE status: 0 AAMAIN: I am 4 node002. 8 MPI context: 91 zgebs2d check 4 0 186 453 186 0 0 F Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 4 MPI FINALIZE status: 0 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 = 205.377 inner side: x = -49.955 R = 116.940 Antenna at psi =: 1.182 (mesh point 378) outer side: x = 46.690 R = 213.586 inner side: x = -53.838 R = 113.058 Wall at psi =: 1.415 (mesh point 453) outer side: x = 57.164 R = 224.060 inner side: x = -58.791 R = 108.105 Two-point boundary condition at the antenna Two-point boundary condition at the plasma edge zgebs2d check 0 0 186 453 186 0 0 T Psi partition: 1 57 113 169 226 283 340 397 454 56 56 56 57 57 57 57 57 0 56 453 8 5 91 0 11 24 24 Global power balance Power using the Parseval formula Power to ions 1 fundamental 0.0000E+00 MW/kA^2 ( 0.00 %) first harm. 1.9901E-01 MW/kA^2 ( 16.45 %) Power to ions 2 fundamental 7.5825E-01 MW/kA^2 ( 62.66 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 2.5280E-01 MW/kA^2 ( 20.89 %) - from the fast wave (*) 2.2661E-01 MW/kA^2 ( 18.73 %) - from the ib wave 2.6186E-02 MW/kA^2 ( 2.16 %) Total absorbed power 1.2101E+00 MW/kA^2 Total radiated power (J*E) 1.2472E+00 MW/kA^2 Total radiated power (Poynting) 1.2482E+00 MW/kA^2 - towards the wall 4.6859E-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. 1.9901E-01 MW/kA^2 ( 16.33 %) Power to ions 2 fundamental 7.6697E-01 MW/kA^2 ( 62.93 %) first harm. 0.0000E+00 MW/kA^2 ( 0.00 %) Power to the electrons (tot) 2.5283E-01 MW/kA^2 ( 20.74 %) - from the fast wave (*) 2.2665E-01 MW/kA^2 ( 18.60 %) - from the ib wave 2.6186E-02 MW/kA^2 ( 2.15 %) Total absorbed power 1.2188E+00 MW/kA^2 Accuracy of the 2dim abs. power 7.2392E-03 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 5.6684E-01 (10^20 M^3) Current drive estimate, total: RF current per incident power 1.4592E-02 (A/W) Current per incident power(NT) 2.8470E-02 (A/W) Current drive figure of merit 1.3804E-02 (A/W/M^2) Current drive figure of merit(NT) 2.6934E-02 (A/W/M^2) Current drive estimate, fast Wave: RF current per incident power 1.4193E-02 (A/W) Current per incident power(NT) 2.7730E-02 (A/W) Current drive figure of merit 1.3427E-02 (A/W/M^2) Current drive figure of merit(NT) 2.6234E-02 (A/W/M^2) Current drive estimate, IB Wave: RF current per incident power 3.9926E-04 (A/W) Current per incident power(NT) 7.3961E-04 (A/W) Current drive figure of merit 3.7772E-04 (A/W/M^2) Current drive figure of merit(NT) 6.9970E-04 (A/W/M^2) The output has been stored on disk toric4 was called 1 times, total 111.97 secs gloini was called 1 times, total 5.07 secs fftl_ini was called 1 times, total 5.01 secs equil was called 1 times, total 0.03 secs mhdcoe was called 2499 times, total 0.06 secs fftl_inv was called 1428041 times, total 5.85 secs profnt was called 1 times, total 0.00 secs coords was called 1513 times, total 0.23 secs fftl_dir was called 18716 times, total 0.09 secs mblock was called 453 times, total 50.16 secs mblock:inside plasma was called 954 times, total 46.66 secs operev was called 954 times, total 42.92 secs zetaRe was called 3257537 times, total 7.99 secs eldibw was called 19831 times, total 14.24 secs direl_flr was called 19831 times, total 7.08 secs pzgetrf was called 453 times, total 12.87 secs pzgetrs,Gk=Ek\Bk was called 1 times, total 0.03 secs pzgetrs,yk=Ek\yk was called 453 times, total 0.58 secs yvec<-yk was called 453 times, total 1.38 secs writeGk was called 452 times, total 0.98 secs pzgemm was called 452 times, total 10.89 secs pzgetrs,Gk=Ek\Gk was called 451 times, total 17.42 secs ykm1<-yvec was called 452 times, total 0.71 secs pzgemm,yk=yk-Gk*ykm1 was called 452 times, total 0.16 secs cibess was called 848544 times, total 3.79 secs mblock:inside vacuum was called 396 times, total 2.31 secs opevac was called 396 times, total 0.79 secs solve: xvec <- x_in was called 1 times, total 0.02 secs yk<-yvec was called 453 times, total 0.73 secs xvec<-xk was called 453 times, total 1.05 secs readGk was called 452 times, total 0.48 secs pzgemm, xk=xk-Gk*xkp was called 452 times, total 0.77 secs solve: x_in <- xvec was called 1 times, total 0.04 secs solve:bcast was called 1 times, total 0.04 secs elepwr was called 1 times, total 7.29 secs elefld was called 56 times, total 0.07 secs powtwd was called 56 times, total 4.64 secs MPI FINALIZE status: 0