| item0 | item1 |
|---|---|
| avgtim(i) | i=1,2..n Avgtim is a list of pivot times used to linearize the transport equations (similar to the Crank Nicholson central pivot time). It is required because some quantities are determined only at the pivot times and not at the full time (given in time(i)). It is defined as avgtim(i) = time - (1.0-theta)*dtt (The association of the dependent variables with the appropriate time arrays is done automatically so users should need not be concerned about this) |
| beamon | inone namelist 2 input Time (sec) at which beam heating is turned on |
| btime | inone namelist 2 input Time interval (sec) during which beam heating is on |
| chieinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The electron thermal diffusivity,cm**2/sec. In analysis mode this quantity is obtained from the conductive flux divided by the negative electron temperature gradient and the electron density. The conductive flux is obtained from the volume integral of qconde (see qconde). |
| chietrinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max Chietrinv is defined only in analysis mode. The transport electron chi obtained by neglecting losses due to convection. (The energy associated with convection is counted as being part of the conductive channel in analysis mode,cm**2/sec. |
| chiiinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ion thermal diffusivity,cm**2/sec. In analysis mode this quantity is obtained from the conductive flux divided by the negative ion temperature gradient and the total primary ion density. The conductive flux is obtained from the volume integral of qcondi (see qcondi). |
| chiineo(j,i) | j=1,2..kj,i=1,2..nprtlst_max The neoclassical ion thermal Diffusivity,cm**2/sec. The factor wneo(3,3) is NOT included here! |
| chiitrinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max Chiitrinv is defined only in analysis mode. The transport ion chi obtained by neglecting losses due to convection. (The energy associated with convection is counted as being part of the conductive channel in analysis mode,cm**2/sec. |
| curboot_ohml(j,i) | j=1,2..kj,i=1,2..ki3d The experimental bootstrap current profile, amps/cm**2 Uses eta parallel and models of beam and rf driven current (as stored in curdri) to separate out the bootstrap current. Hence this value of the bootstrap current DOES NOT depend on a theoretical model of the bootstrap current. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| curdrive_ohml(j,i) | j=1,2..kj,i=1,2..ki3d The experimental total driven current profile, amps/cm**2 Uses eta parallel and models of bootstrap and rf driven current (as stored in curdri) to separate out the bootstrap current. Hence this value of the bootstrap current DOES NOT depend on a theoretical model of the bootstrap current. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| curni_ohml(j,i) | j=1,2..kj,i=1,2..ki3d The experimental total noninductive current profile, amps/cm**2. Uses whatever model of resistivity is currently stored in etap to separate the product of eta*curden_ohmic . Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| deetot(i) | i=1,...n The total rate of change of electron stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| deitot(i) | i=1,...n The total rate of change of ion stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| detot(i) | i=1,...n The total rate of change of plasma stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| dpedt(j,i) | j=1,2..kj,i=1,2..nprtlst_max The local rate of change in the stored electron energy density (ie 1.5*pressure),watts/cm**3 . The derivative is obtained from a simple forward difference approximation! |
| dpidt(j,i) | j=1,2..kj,i=1,2..nprtlst_max The local rate of change in the stored ion energy density (ie 1.5*pressure),watts/cm**3 . The derivative is obtained from a simple forward difference approximation! |
| eatot(i) | i=1,...n The (fusion progenated) alpha stored energy,joules. |
| ebtot(i) | i=1,...n The beam ion stored energy,joules. |
| eetot(i) | i=1,...n The electron energy,joules. |
| eitot(i) | i=1,...n The thermal ion energy,joules. |
| eneiv(j,i) | j=1,2..kj, i=1,2..nprtlst_max The electron density at timeiv(i),#/cm**3 as a function of rho(j) |
| eniiv(j,i) | j=1,2..kj, i=1,2..nprtlst_max The ion densities at timeiv(i),#/cm**3, as a function of rho(j) |
| entaue(i) | i=1,...n The n*tau confinement parameter. Total number of electrons in the plasma times taue energy confinement time divided by plasma volume,sec/cm**3. |
| entot(i) | i=1,...n The total number of electrons in the plasma. |
| eta_par_ohml(j,i) | j=1,2..kj,i=1,2..ki3d. The experimental resistivity, ohm-cm. Given the beam,rf,and bootstrap current profiles (as determined by the theoretical input models selected by the user) get the resistivity that satisfies Faraday's Law. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| etot(i) | i=1,...n The total (kinetic) plasma energy,joules. |
| exptl_neutron_rate(i) | i=1,...n This is just a simple linear interpolation of the input experimental neutron rate onto the time(i) grid.,#/sec. |
| fast_ion_target | inone namelist 2 input integer,if set to 1 will let the incoming neutral beam see the stored fast ion density no attempt to modify the stopping rates due to the fact that the fast ion distribution is not Maxwellian has been made! |
| gconde(j,i) | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to conduction, watts/cm**2. |
| gcondi(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to conduction, watts/cm**2. |
| gconvde(j,i) | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to conduction plus convection,watts/cm**2. |
| gconvdi(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to conduction plus convection,watts/cm**2. |
| gconve(j,i) | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to convection, watts/cm**2. |
| gconvi(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to convection, watts/cm**2. |
| hdepsmth | inone namelist 2 input set this param to a positive value (gt.0.0) to turn off the smoothing of hibrz and hdepz in subroutine POSTNUB. if this option is used then enough zones must be specified for adequate resolution (zones = number of radial grid points) and enough injected neutrals must be followed to minimize the statistical noise enough so that no greatly uneven profiles result! this option was added because the smoothing of the profiles by subroutine SMOOTH can lead to unphysical peaking of the birth and deposition profiles. |
| imaxbeamconvg | inone namelist 2 input iterations are done to get consistency in the thermal and fast ion densities (ONLY). default imaxbeamconvg = 5 |
| ishot | inone namelist 3 input ishot should be set to shot number in netCDF file The tdem mode NEVER solves the MHD (i.e., Grad-Shafranov) equation. Instead it reads a precomputed file of MHD data (in netCDF format) as a function of SPACE AND TIME and interpolates this data for the values required at each time step. |
| iterate_beam | inone namelist 2 input logical, if .true., allow up to imaxbeamconvg iterations default iterate_beam = .false. |
| mhdmethd | inone namelist 3 input set mhdmethd = 'green' to get solution of Grad-Shafranov equation using Green's function for entire MHD grid (this is the slowww method.It should not be used routinely. It was included here primarily as an aid to verify solutions obtained with the new cyclic reduction solver used in subroutine CMPLTCYR.) set mhdmethd = 'cycred' to use cyclic reduction method of inverting del-star(default). This method is much faster and should be used routinely. mhdmode = 'no coils' for fixed boundary cases mhdmethd = 'cycred' yields the cyclic reduction solution of the fixed boundary problem. mhdmethd = 'sorpicrd' yields the successive overrelaxation method. mhdmethd ='tdem' for time-dependent eqdsk mode if mhdmethd ='tdem' (i.e., the multiple eqdsk case) then this list of times is not used (the netCDF file shot_12.cdf supplies the data instead) The tdem mode also requires that you set the following switches (these are set automatically for you even if you set them differently or didn't set them at all) if mhdmethd = 'tdem' then the following settings are FORCED ieqdsk = 0 mhdmode = 'no coils' ifixshap = 0 also eqdskin must be set to an appropriate netCDF file (shot_12.cdf) which we CANNOT set for the user |
| ONETWOPAGES | An area on the hydra where online postscript information regarding Onetwo is kept. At this time it is very incomplete but some useful stuff is there. |
| pbeame(j,i) | j=1,2..kj,i=1,2..nprtlst_max The integrated power to the electrons due to beam ion friction,watts. |
| pbeami(j,i) | j=1,2..kj,i=1,2..nprtlst_max The integrated power to the ions due to beam ion friction,watts. |
| pconden(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the electron energy flux due to conduction,watts, as a function of rho(j). (see qconde for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pconden(j,i)/(4*pi**2*R0*H*rho). |
| pcondin(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the ion energy flux due to conduction,watts, as a function of rho(j). (see qcondi for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pcondin(j,i)/(4*pi**2*R0*H*rho). |
| pconvden(j,i) | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for electrons,watts pconvden(j,i) =pconven(j,i)+pconden(j,i) |
| pconvdin(j,i) | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for ions,watts pconvdin(j,i) = pconvin(j,i)+pcondin(j,i) |
| pconven(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the electron energy flux due to convection,watts, as a function of rho(j). (see qconve for more info). The (negative) convective flux,watts/cm**2,can be (should be) obtained as pconven(j,i)/(4*pi**2*R0*H*rho). |
| pconvin(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the ion energy flux due to convection,watts, as a function of rho(j). (see qconvi for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pconvin(j,i)/(4*pi**2*R0*H*rho). |
| pdelten(j,i) | j=1,2..kj,i=1,2..nprtlst_max The volume integral of qdelten out to the value rho(j). |
| pe2d(j,i) | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qe2d out to rho(j),watts. |
| pfusi(j,i) | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qfusi out to rho(j),watts. |
| pi2d(j,i) | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qi2d out to rho(j),watts. |
| pneuti(j,i) | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qneuti out to rho(j),watts. |
| pohm(j,i) | j=1,2..kj,i=1,2..nprtlst_max The volume integrated ohmic power delivered to the electrons,watts. |
| prfe(j,i) | j=1,2..kj,i=1,2..nprtlst_max The volume integrated rf power delivered to the electrons,watts. |
| prfi(j,i) | j=1,2..kj,i=1,2..nprtlst_max The volume integrated rf power delivered to the ions,watts. |
| proti(j,i) | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qroti out to rho(j),watts. |
| qbeame(j,i) | j=1,2..kj,i=1,2..nprtlst_max The power to the electrons due to beam ion friction,watts/cm*3. |
| qbeami(j,i) | j=1,2..kj,i=1,2..nprtlst_max The power to the ions due to beam ion friction,watts/cm*3. |
| qconden(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the electron energy flux due to conduction,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! qconden = -qconde and qconde is chosen (in analysis mode ) to guarantee energy conservation |
| qcondin(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the ion energy flux due to conduction,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. qcondin = - qcondi and qcondi is determined (in analysis mode) so as to guarantee energy conservation |
| qconvden(j,i) | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for electrons,watts/cm**3 qconvden(j,i) = qconven(j,i)+qconden(j,i) |
| qconvdin(j,i) | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for ions,watts/cm**3 qconvdin(j,i) = qconvin(j,i)+qcondin(j,i) |
| qconven(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the electron energy flux due to convection,watts/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. |
| qconvin(j,i) | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the ion energy flux due to convection,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. |
| qdelten(j,i) | j=1,2..kj,i=1,2..nprtlst_max The negative of the neoclassical electron ion energy exchange term,watts/cm**3 qdelten = - qdelt with qdelt = c*(te-ti) where c is a POSITIVE parameter. |
| qe2d(j,i) | j=1,2..kj,i=1,2..nprtlst_max Heating of electron distribution function due to two dimensional MHD effects,watts/cm**3. This term is identically zero in the usual analysis mode runs. It is NOT zero when the time depedent eqdsk mode (TDEM) is used. See the ONETWOPAGES for detailed description. |
| qfusi(j,i) | j=1,2..kj,i=1,2..nprtlst_max Represents heating of ion distribution due to alpha particle slowing down,watts/cm**3. |
| qi2d(j,i) | j=1,2..kj,i=1,2..nprtlst_max Heating of thermal ion distribution function due to two dimensional MHD effects,watts/cm**3. This term is identically zero in the usual analysis mode runs. It is NOT zero when the time depedent eqdsk mode (TDEM) is used. See the ONETWOPAGES for detailed description. |
| qneut(j,i) | j=1,2..kj,i=1,2..nprtlst_max This parameter is defined as qneut = qioni - qcx It represents the net power input due to electron impact ionization of neutrals and recombination of thermal ions (qioni) and charge exchange of thermal ions with beam neutrals ,thermal neutrals and secondary charge exchange of beam ions with thermal neutrals to form a thermal ion (if ibcx is on). See the ONETWOPAGES for more details. |
| qohm(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ohmic power delivered to the electrons, watts/cm**3. |
| qrfe(j,i) | j=1,2..kj,i=1,2..nprtlst_max The (total) rf power delivered to the electrons, watts/cm**3. |
| qrfi(j,i) | j=1,2..kj,i=1,2..nprtlst_max The (total) rf power delivered to the ions, watts/cm**3. |
| qroti(j,i) | j=1,2..kj,i=1,2..nprtlst_max Heating of ion distribution associated with input of toroidal momentum watts/cm**3 |
| relative error | inone namelist 2 input for beam convergence (The relative error is measured in the change in electron density if icenez=0 and as a relative change in the thermal ion species corresponding to the beam if icenez=1) |
| relaxden | inone namelist 2 input relaxation parameter for beam iteration,must be in (0,1] with 1 meaning no relaxation. |
| relaxden_err | inone namelist 2 input |
| taue(i) | i=1,...n The total energy confinement time in sec (includes electron,primary and impurity ions,beamions and alphas) etot+eatot+ebtot taue= ------------------ qtot-detot |
| tauec(i) | i=1,...n The central(ie rho=0) energy confinement time in sec (includes electron,primary and impurity ions,beamions and alphas) |
| taup(i) | i=1,...n The total particle (ie electron) confinement time in sec |
| teiv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The electron temperature at timeiv(i),kev, as a function of rho(j) |
| tiiv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ion temperature at timeiv(i),kev, as a function of rho(j) . |
| timbplt | inone namelist 2 input Times (up to 5) to produce data for Freya-like plots of beam deposition. output is processed by the NUBPLT code. Defaulted to OFF. timbplt(1) .le. time0 .and. beamon .lt. time0 gives o/p for initial time. |
| time(i) | i=1,2....n This is the basic time list. An entry to this list is made at the end of each time step that the code makes sucessfully. Thus n is the total number of time steps taken. Each time step is dynamically chosen by the code so time is not an equispaced vector. |
| timeiv(i) | i=1,2..nprtlst_max timeiv(i) is the time,in sec, for plots of interest in analysis mode (power balance,diffusivities,as well as densities and temperatures,etc). timeiv(i) is APPROXIMATELY controlled by the user through the input time list prtlst,see cray102.f for explanation. The actual times are timeiv(i)=prtlst(i)-(1-theta)*(current time step) where theta is the Crank Nicholson pivot parameter and the current time step is dynamically adjusted by the code ( you control the maximum time step with dtmax and Crank-Nicholson parameter theta in inone,see cray102.f). |
| tniv(j,k,i) | j=1,2..kj,i=1,2..nprtlst_max The neutral temperature at timeiv(i),kev as a function of rho(j) . |
| what is it | A description of the Onetwo output variables used in PLOT12, TRPLOT and other codes. Reference to ONETWOPAGES means the online documents in /u/stjohn/ONETWOPAGES/. |
| xchietot(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode electron thermal diffusivity, cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| xchiitot(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion thermal diffusivity, cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| xdchitot(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion particle diffusivity, of primary ion #1,cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| xkangtot(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode momentum diffusivity, g*cm/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| xkeinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The electron thermal conductivity,1./(cm*sec). In analysis mode this quantity is obtained from the conductive flux divided by the negative electron temperature gradient. The conductive flux is obtained from the volume integral of qconei (see qconde). |
| xketot(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode electron thermal conductivity, 1./(cm*sec). (total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| xkiinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The ion thermal conductivity,1./(cm*sec). In analysis mode this quantity is obtained from the conductive flux divided by the negative ion temperature gradient. The conductive flux is obtained from the volume integral of qcondi (see qcondi). |
| xkineo(j,i) | j=1,2..kj,i=1,2..nprtlst_max The neoclassical ion thermal coductivity,1/(cm*sec)cray102.f. The factor wneo(3,3) is NOT included here! |
| xkitot(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion thermal conductivity, 1./(cm*sec). (total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| xndinv(j,i) | j=1,2..kj,i=1,2..nprtlst_max The total simulation OR ANALYSIS mode ion particle diffusivity of that primary ion species which is also present as a neutral species,cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc. in simulation mode. In analysis mode it is just the primary ion particle flux divided by the density and density gradient). |
| item0: | i=1,2..n Avgtim is a list of pivot times used to linearize the transport equations (similar to the Crank Nicholson central pivot time). It is required because some quantities are determined only at the pivot times and not at the full time (given in time(i)). It is defined as avgtim(i) = time - (1.0-theta)*dtt (The association of the dependent variables with the appropriate time arrays is done automatically so users should need not be concerned about this) |
| item1: | i=1,2..n Avgtim is a list of pivot times used to linearize the transport equations (similar to the Crank Nicholson central pivot time). It is required because some quantities are determined only at the pivot times and not at the full time (given in time(i)). It is defined as avgtim(i) = time - (1.0-theta)*dtt (The association of the dependent variables with the appropriate time arrays is done automatically so users should need not be concerned about this) |
| item0: | inone namelist 2 input Time (sec) at which beam heating is turned on |
| item1: | inone namelist 2 input Time (sec) at which beam heating is turned on |
| item0: | inone namelist 2 input Time interval (sec) during which beam heating is on |
| item1: | inone namelist 2 input Time interval (sec) during which beam heating is on |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The electron thermal diffusivity,cm**2/sec. In analysis mode this quantity is obtained from the conductive flux divided by the negative electron temperature gradient and the electron density. The conductive flux is obtained from the volume integral of qconde (see qconde). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The electron thermal diffusivity,cm**2/sec. In analysis mode this quantity is obtained from the conductive flux divided by the negative electron temperature gradient and the electron density. The conductive flux is obtained from the volume integral of qconde (see qconde). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Chietrinv is defined only in analysis mode. The transport electron chi obtained by neglecting losses due to convection. (The energy associated with convection is counted as being part of the conductive channel in analysis mode,cm**2/sec. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Chietrinv is defined only in analysis mode. The transport electron chi obtained by neglecting losses due to convection. (The energy associated with convection is counted as being part of the conductive channel in analysis mode,cm**2/sec. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ion thermal diffusivity,cm**2/sec. In analysis mode this quantity is obtained from the conductive flux divided by the negative ion temperature gradient and the total primary ion density. The conductive flux is obtained from the volume integral of qcondi (see qcondi). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ion thermal diffusivity,cm**2/sec. In analysis mode this quantity is obtained from the conductive flux divided by the negative ion temperature gradient and the total primary ion density. The conductive flux is obtained from the volume integral of qcondi (see qcondi). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The neoclassical ion thermal Diffusivity,cm**2/sec. The factor wneo(3,3) is NOT included here! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The neoclassical ion thermal Diffusivity,cm**2/sec. The factor wneo(3,3) is NOT included here! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Chiitrinv is defined only in analysis mode. The transport ion chi obtained by neglecting losses due to convection. (The energy associated with convection is counted as being part of the conductive channel in analysis mode,cm**2/sec. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Chiitrinv is defined only in analysis mode. The transport ion chi obtained by neglecting losses due to convection. (The energy associated with convection is counted as being part of the conductive channel in analysis mode,cm**2/sec. |
| item0: | j=1,2..kj,i=1,2..ki3d The experimental bootstrap current profile, amps/cm**2 Uses eta parallel and models of beam and rf driven current (as stored in curdri) to separate out the bootstrap current. Hence this value of the bootstrap current DOES NOT depend on a theoretical model of the bootstrap current. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item1: | j=1,2..kj,i=1,2..ki3d The experimental bootstrap current profile, amps/cm**2 Uses eta parallel and models of beam and rf driven current (as stored in curdri) to separate out the bootstrap current. Hence this value of the bootstrap current DOES NOT depend on a theoretical model of the bootstrap current. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item0: | j=1,2..kj,i=1,2..ki3d The experimental total driven current profile, amps/cm**2 Uses eta parallel and models of bootstrap and rf driven current (as stored in curdri) to separate out the bootstrap current. Hence this value of the bootstrap current DOES NOT depend on a theoretical model of the bootstrap current. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item1: | j=1,2..kj,i=1,2..ki3d The experimental total driven current profile, amps/cm**2 Uses eta parallel and models of bootstrap and rf driven current (as stored in curdri) to separate out the bootstrap current. Hence this value of the bootstrap current DOES NOT depend on a theoretical model of the bootstrap current. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item0: | j=1,2..kj,i=1,2..ki3d The experimental total noninductive current profile, amps/cm**2. Uses whatever model of resistivity is currently stored in etap to separate the product of eta*curden_ohmic . Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item1: | j=1,2..kj,i=1,2..ki3d The experimental total noninductive current profile, amps/cm**2. Uses whatever model of resistivity is currently stored in etap to separate the product of eta*curden_ohmic . Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item0: | i=1,...n The total rate of change of electron stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| item1: | i=1,...n The total rate of change of electron stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| item0: | i=1,...n The total rate of change of ion stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| item1: | i=1,...n The total rate of change of ion stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| item0: | i=1,...n The total rate of change of plasma stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| item1: | i=1,...n The total rate of change of plasma stored energy, watts. Simple forward differencing is used to get this quantity. You may want to control the time step with dtmax to get a smooth value. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The local rate of change in the stored electron energy density (ie 1.5*pressure),watts/cm**3 . The derivative is obtained from a simple forward difference approximation! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The local rate of change in the stored electron energy density (ie 1.5*pressure),watts/cm**3 . The derivative is obtained from a simple forward difference approximation! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The local rate of change in the stored ion energy density (ie 1.5*pressure),watts/cm**3 . The derivative is obtained from a simple forward difference approximation! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The local rate of change in the stored ion energy density (ie 1.5*pressure),watts/cm**3 . The derivative is obtained from a simple forward difference approximation! |
| item0: | i=1,...n The (fusion progenated) alpha stored energy,joules. |
| item1: | i=1,...n The (fusion progenated) alpha stored energy,joules. |
| item0: | i=1,...n The beam ion stored energy,joules. |
| item1: | i=1,...n The beam ion stored energy,joules. |
| item0: | i=1,...n The electron energy,joules. |
| item1: | i=1,...n The electron energy,joules. |
| item0: | i=1,...n The thermal ion energy,joules. |
| item1: | i=1,...n The thermal ion energy,joules. |
| item0: | j=1,2..kj, i=1,2..nprtlst_max The electron density at timeiv(i),#/cm**3 as a function of rho(j) |
| item1: | j=1,2..kj, i=1,2..nprtlst_max The electron density at timeiv(i),#/cm**3 as a function of rho(j) |
| item0: | j=1,2..kj, i=1,2..nprtlst_max The ion densities at timeiv(i),#/cm**3, as a function of rho(j) |
| item1: | j=1,2..kj, i=1,2..nprtlst_max The ion densities at timeiv(i),#/cm**3, as a function of rho(j) |
| item0: | i=1,...n The n*tau confinement parameter. Total number of electrons in the plasma times taue energy confinement time divided by plasma volume,sec/cm**3. |
| item1: | i=1,...n The n*tau confinement parameter. Total number of electrons in the plasma times taue energy confinement time divided by plasma volume,sec/cm**3. |
| item0: | i=1,...n The total number of electrons in the plasma. |
| item1: | i=1,...n The total number of electrons in the plasma. |
| item0: | j=1,2..kj,i=1,2..ki3d. The experimental resistivity, ohm-cm. Given the beam,rf,and bootstrap current profiles (as determined by the theoretical input models selected by the user) get the resistivity that satisfies Faraday's Law. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item1: | j=1,2..kj,i=1,2..ki3d. The experimental resistivity, ohm-cm. Given the beam,rf,and bootstrap current profiles (as determined by the theoretical input models selected by the user) get the resistivity that satisfies Faraday's Law. Calculated only in tdem mode and then only if the poloidal magnetic field is run in analysis mode. See the ONETWOPAGES for detailed description. |
| item0: | i=1,...n The total (kinetic) plasma energy,joules. |
| item1: | i=1,...n The total (kinetic) plasma energy,joules. |
| item0: | i=1,...n This is just a simple linear interpolation of the input experimental neutron rate onto the time(i) grid.,#/sec. |
| item1: | i=1,...n This is just a simple linear interpolation of the input experimental neutron rate onto the time(i) grid.,#/sec. |
| item0: | inone namelist 2 input integer,if set to 1 will let the incoming neutral beam see the stored fast ion density no attempt to modify the stopping rates due to the fact that the fast ion distribution is not Maxwellian has been made! |
| item1: | inone namelist 2 input integer,if set to 1 will let the incoming neutral beam see the stored fast ion density no attempt to modify the stopping rates due to the fact that the fast ion distribution is not Maxwellian has been made! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to conduction, watts/cm**2. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to conduction, watts/cm**2. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to conduction, watts/cm**2. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to conduction, watts/cm**2. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to conduction plus convection,watts/cm**2. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to conduction plus convection,watts/cm**2. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to conduction plus convection,watts/cm**2. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to conduction plus convection,watts/cm**2. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to convection, watts/cm**2. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The electron energy flux due to convection, watts/cm**2. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to convection, watts/cm**2. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ion energy flux due to convection, watts/cm**2. |
| item0: | inone namelist 2 input set this param to a positive value (gt.0.0) to turn off the smoothing of hibrz and hdepz in subroutine POSTNUB. if this option is used then enough zones must be specified for adequate resolution (zones = number of radial grid points) and enough injected neutrals must be followed to minimize the statistical noise enough so that no greatly uneven profiles result! this option was added because the smoothing of the profiles by subroutine SMOOTH can lead to unphysical peaking of the birth and deposition profiles. |
| item1: | inone namelist 2 input set this param to a positive value (gt.0.0) to turn off the smoothing of hibrz and hdepz in subroutine POSTNUB. if this option is used then enough zones must be specified for adequate resolution (zones = number of radial grid points) and enough injected neutrals must be followed to minimize the statistical noise enough so that no greatly uneven profiles result! this option was added because the smoothing of the profiles by subroutine SMOOTH can lead to unphysical peaking of the birth and deposition profiles. |
| item0: | inone namelist 2 input iterations are done to get consistency in the thermal and fast ion densities (ONLY). default imaxbeamconvg = 5 |
| item1: | inone namelist 2 input iterations are done to get consistency in the thermal and fast ion densities (ONLY). default imaxbeamconvg = 5 |
| item0: | inone namelist 3 input ishot should be set to shot number in netCDF file The tdem mode NEVER solves the MHD (i.e., Grad-Shafranov) equation. Instead it reads a precomputed file of MHD data (in netCDF format) as a function of SPACE AND TIME and interpolates this data for the values required at each time step. |
| item1: | inone namelist 3 input ishot should be set to shot number in netCDF file The tdem mode NEVER solves the MHD (i.e., Grad-Shafranov) equation. Instead it reads a precomputed file of MHD data (in netCDF format) as a function of SPACE AND TIME and interpolates this data for the values required at each time step. |
| item0: | inone namelist 2 input logical, if .true., allow up to imaxbeamconvg iterations default iterate_beam = .false. |
| item1: | inone namelist 2 input logical, if .true., allow up to imaxbeamconvg iterations default iterate_beam = .false. |
| item0: | inone namelist 3 input set mhdmethd = 'green' to get solution of Grad-Shafranov equation using Green's function for entire MHD grid (this is the slowww method.It should not be used routinely. It was included here primarily as an aid to verify solutions obtained with the new cyclic reduction solver used in subroutine CMPLTCYR.) set mhdmethd = 'cycred' to use cyclic reduction method of inverting del-star(default). This method is much faster and should be used routinely. mhdmode = 'no coils' for fixed boundary cases mhdmethd = 'cycred' yields the cyclic reduction solution of the fixed boundary problem. mhdmethd = 'sorpicrd' yields the successive overrelaxation method. mhdmethd ='tdem' for time-dependent eqdsk mode if mhdmethd ='tdem' (i.e., the multiple eqdsk case) then this list of times is not used (the netCDF file shot_12.cdf supplies the data instead) The tdem mode also requires that you set the following switches (these are set automatically for you even if you set them differently or didn't set them at all) if mhdmethd = 'tdem' then the following settings are FORCED ieqdsk = 0 mhdmode = 'no coils' ifixshap = 0 also eqdskin must be set to an appropriate netCDF file (shot_12.cdf) which we CANNOT set for the user |
| item1: | inone namelist 3 input set mhdmethd = 'green' to get solution of Grad-Shafranov equation using Green's function for entire MHD grid (this is the slowww method.It should not be used routinely. It was included here primarily as an aid to verify solutions obtained with the new cyclic reduction solver used in subroutine CMPLTCYR.) set mhdmethd = 'cycred' to use cyclic reduction method of inverting del-star(default). This method is much faster and should be used routinely. mhdmode = 'no coils' for fixed boundary cases mhdmethd = 'cycred' yields the cyclic reduction solution of the fixed boundary problem. mhdmethd = 'sorpicrd' yields the successive overrelaxation method. mhdmethd ='tdem' for time-dependent eqdsk mode if mhdmethd ='tdem' (i.e., the multiple eqdsk case) then this list of times is not used (the netCDF file shot_12.cdf supplies the data instead) The tdem mode also requires that you set the following switches (these are set automatically for you even if you set them differently or didn't set them at all) if mhdmethd = 'tdem' then the following settings are FORCED ieqdsk = 0 mhdmode = 'no coils' ifixshap = 0 also eqdskin must be set to an appropriate netCDF file (shot_12.cdf) which we CANNOT set for the user |
| item0: | An area on the hydra where online postscript information regarding Onetwo is kept. At this time it is very incomplete but some useful stuff is there. |
| item1: | An area on the hydra where online postscript information regarding Onetwo is kept. At this time it is very incomplete but some useful stuff is there. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The integrated power to the electrons due to beam ion friction,watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The integrated power to the electrons due to beam ion friction,watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The integrated power to the ions due to beam ion friction,watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The integrated power to the ions due to beam ion friction,watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the electron energy flux due to conduction,watts, as a function of rho(j). (see qconde for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pconden(j,i)/(4*pi**2*R0*H*rho). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the electron energy flux due to conduction,watts, as a function of rho(j). (see qconde for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pconden(j,i)/(4*pi**2*R0*H*rho). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the ion energy flux due to conduction,watts, as a function of rho(j). (see qcondi for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pcondin(j,i)/(4*pi**2*R0*H*rho). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the ion energy flux due to conduction,watts, as a function of rho(j). (see qcondi for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pcondin(j,i)/(4*pi**2*R0*H*rho). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for electrons,watts pconvden(j,i) =pconven(j,i)+pconden(j,i) |
| item1: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for electrons,watts pconvden(j,i) =pconven(j,i)+pconden(j,i) |
| item0: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for ions,watts pconvdin(j,i) = pconvin(j,i)+pcondin(j,i) |
| item1: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for ions,watts pconvdin(j,i) = pconvin(j,i)+pcondin(j,i) |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the electron energy flux due to convection,watts, as a function of rho(j). (see qconve for more info). The (negative) convective flux,watts/cm**2,can be (should be) obtained as pconven(j,i)/(4*pi**2*R0*H*rho). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the electron energy flux due to convection,watts, as a function of rho(j). (see qconve for more info). The (negative) convective flux,watts/cm**2,can be (should be) obtained as pconven(j,i)/(4*pi**2*R0*H*rho). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the ion energy flux due to convection,watts, as a function of rho(j). (see qconvi for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pconvin(j,i)/(4*pi**2*R0*H*rho). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the integrated divergence of the ion energy flux due to convection,watts, as a function of rho(j). (see qconvi for more info). The negative convective flux,watts/cm**2,can be (should be) obtained as pconvin(j,i)/(4*pi**2*R0*H*rho). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The volume integral of qdelten out to the value rho(j). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The volume integral of qdelten out to the value rho(j). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qe2d out to rho(j),watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qe2d out to rho(j),watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qfusi out to rho(j),watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qfusi out to rho(j),watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qi2d out to rho(j),watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qi2d out to rho(j),watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qneuti out to rho(j),watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qneuti out to rho(j),watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The volume integrated ohmic power delivered to the electrons,watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The volume integrated ohmic power delivered to the electrons,watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The volume integrated rf power delivered to the electrons,watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The volume integrated rf power delivered to the electrons,watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The volume integrated rf power delivered to the ions,watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The volume integrated rf power delivered to the ions,watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qroti out to rho(j),watts. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Volume integrated qroti out to rho(j),watts. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The power to the electrons due to beam ion friction,watts/cm*3. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The power to the electrons due to beam ion friction,watts/cm*3. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The power to the ions due to beam ion friction,watts/cm*3. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The power to the ions due to beam ion friction,watts/cm*3. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the electron energy flux due to conduction,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! qconden = -qconde and qconde is chosen (in analysis mode ) to guarantee energy conservation |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the electron energy flux due to conduction,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! qconden = -qconde and qconde is chosen (in analysis mode ) to guarantee energy conservation |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the ion energy flux due to conduction,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. qcondin = - qcondi and qcondi is determined (in analysis mode) so as to guarantee energy conservation |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the ion energy flux due to conduction,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. qcondin = - qcondi and qcondi is determined (in analysis mode) so as to guarantee energy conservation |
| item0: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for electrons,watts/cm**3 qconvden(j,i) = qconven(j,i)+qconden(j,i) |
| item1: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for electrons,watts/cm**3 qconvden(j,i) = qconven(j,i)+qconden(j,i) |
| item0: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for ions,watts/cm**3 qconvdin(j,i) = qconvin(j,i)+qcondin(j,i) |
| item1: | j=1,2..kj,i=1,2..nprtlst_max This is the sum of conduction plus convection terms for ions,watts/cm**3 qconvdin(j,i) = qconvin(j,i)+qcondin(j,i) |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the electron energy flux due to convection,watts/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the electron energy flux due to convection,watts/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the ion energy flux due to convection,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The NEGATIVE of the divergence of the ion energy flux due to convection,watt/cm**3, as a function of rho(j). In the power balance tables sources of power are positive and sinks are negative. Since conduction or convection out of a volume element is a sink of power these two quantities appear negative in these tables whenever the flow of power is out of the local volume element! That is why the NEGATIVE of the divergence of the flux is given. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The negative of the neoclassical electron ion energy exchange term,watts/cm**3 qdelten = - qdelt with qdelt = c*(te-ti) where c is a POSITIVE parameter. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The negative of the neoclassical electron ion energy exchange term,watts/cm**3 qdelten = - qdelt with qdelt = c*(te-ti) where c is a POSITIVE parameter. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Heating of electron distribution function due to two dimensional MHD effects,watts/cm**3. This term is identically zero in the usual analysis mode runs. It is NOT zero when the time depedent eqdsk mode (TDEM) is used. See the ONETWOPAGES for detailed description. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Heating of electron distribution function due to two dimensional MHD effects,watts/cm**3. This term is identically zero in the usual analysis mode runs. It is NOT zero when the time depedent eqdsk mode (TDEM) is used. See the ONETWOPAGES for detailed description. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Represents heating of ion distribution due to alpha particle slowing down,watts/cm**3. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Represents heating of ion distribution due to alpha particle slowing down,watts/cm**3. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Heating of thermal ion distribution function due to two dimensional MHD effects,watts/cm**3. This term is identically zero in the usual analysis mode runs. It is NOT zero when the time depedent eqdsk mode (TDEM) is used. See the ONETWOPAGES for detailed description. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Heating of thermal ion distribution function due to two dimensional MHD effects,watts/cm**3. This term is identically zero in the usual analysis mode runs. It is NOT zero when the time depedent eqdsk mode (TDEM) is used. See the ONETWOPAGES for detailed description. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max This parameter is defined as qneut = qioni - qcx It represents the net power input due to electron impact ionization of neutrals and recombination of thermal ions (qioni) and charge exchange of thermal ions with beam neutrals ,thermal neutrals and secondary charge exchange of beam ions with thermal neutrals to form a thermal ion (if ibcx is on). See the ONETWOPAGES for more details. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max This parameter is defined as qneut = qioni - qcx It represents the net power input due to electron impact ionization of neutrals and recombination of thermal ions (qioni) and charge exchange of thermal ions with beam neutrals ,thermal neutrals and secondary charge exchange of beam ions with thermal neutrals to form a thermal ion (if ibcx is on). See the ONETWOPAGES for more details. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ohmic power delivered to the electrons, watts/cm**3. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ohmic power delivered to the electrons, watts/cm**3. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The (total) rf power delivered to the electrons, watts/cm**3. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The (total) rf power delivered to the electrons, watts/cm**3. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The (total) rf power delivered to the ions, watts/cm**3. |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The (total) rf power delivered to the ions, watts/cm**3. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max Heating of ion distribution associated with input of toroidal momentum watts/cm**3 |
| item1: | j=1,2..kj,i=1,2..nprtlst_max Heating of ion distribution associated with input of toroidal momentum watts/cm**3 |
| item0: | inone namelist 2 input for beam convergence (The relative error is measured in the change in electron density if icenez=0 and as a relative change in the thermal ion species corresponding to the beam if icenez=1) |
| item1: | inone namelist 2 input for beam convergence (The relative error is measured in the change in electron density if icenez=0 and as a relative change in the thermal ion species corresponding to the beam if icenez=1) |
| item0: | inone namelist 2 input relaxation parameter for beam iteration,must be in (0,1] with 1 meaning no relaxation. |
| item1: | inone namelist 2 input relaxation parameter for beam iteration,must be in (0,1] with 1 meaning no relaxation. |
| item0: | inone namelist 2 input |
| item1: | inone namelist 2 input |
| item0: | i=1,...n The total energy confinement time in sec (includes electron,primary and impurity ions,beamions and alphas) etot+eatot+ebtot taue= ------------------ qtot-detot |
| item1: | i=1,...n The total energy confinement time in sec (includes electron,primary and impurity ions,beamions and alphas) etot+eatot+ebtot taue= ------------------ qtot-detot |
| item0: | i=1,...n The central(ie rho=0) energy confinement time in sec (includes electron,primary and impurity ions,beamions and alphas) |
| item1: | i=1,...n The central(ie rho=0) energy confinement time in sec (includes electron,primary and impurity ions,beamions and alphas) |
| item0: | i=1,...n The total particle (ie electron) confinement time in sec |
| item1: | i=1,...n The total particle (ie electron) confinement time in sec |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The electron temperature at timeiv(i),kev, as a function of rho(j) |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The electron temperature at timeiv(i),kev, as a function of rho(j) |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ion temperature at timeiv(i),kev, as a function of rho(j) . |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ion temperature at timeiv(i),kev, as a function of rho(j) . |
| item0: | inone namelist 2 input Times (up to 5) to produce data for Freya-like plots of beam deposition. output is processed by the NUBPLT code. Defaulted to OFF. timbplt(1) .le. time0 .and. beamon .lt. time0 gives o/p for initial time. |
| item1: | inone namelist 2 input Times (up to 5) to produce data for Freya-like plots of beam deposition. output is processed by the NUBPLT code. Defaulted to OFF. timbplt(1) .le. time0 .and. beamon .lt. time0 gives o/p for initial time. |
| item0: | i=1,2....n This is the basic time list. An entry to this list is made at the end of each time step that the code makes sucessfully. Thus n is the total number of time steps taken. Each time step is dynamically chosen by the code so time is not an equispaced vector. |
| item1: | i=1,2....n This is the basic time list. An entry to this list is made at the end of each time step that the code makes sucessfully. Thus n is the total number of time steps taken. Each time step is dynamically chosen by the code so time is not an equispaced vector. |
| item0: | i=1,2..nprtlst_max timeiv(i) is the time,in sec, for plots of interest in analysis mode (power balance,diffusivities,as well as densities and temperatures,etc). timeiv(i) is APPROXIMATELY controlled by the user through the input time list prtlst,see cray102.f for explanation. The actual times are timeiv(i)=prtlst(i)-(1-theta)*(current time step) where theta is the Crank Nicholson pivot parameter and the current time step is dynamically adjusted by the code ( you control the maximum time step with dtmax and Crank-Nicholson parameter theta in inone,see cray102.f). |
| item1: | i=1,2..nprtlst_max timeiv(i) is the time,in sec, for plots of interest in analysis mode (power balance,diffusivities,as well as densities and temperatures,etc). timeiv(i) is APPROXIMATELY controlled by the user through the input time list prtlst,see cray102.f for explanation. The actual times are timeiv(i)=prtlst(i)-(1-theta)*(current time step) where theta is the Crank Nicholson pivot parameter and the current time step is dynamically adjusted by the code ( you control the maximum time step with dtmax and Crank-Nicholson parameter theta in inone,see cray102.f). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The neutral temperature at timeiv(i),kev as a function of rho(j) . |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The neutral temperature at timeiv(i),kev as a function of rho(j) . |
| item0: | A description of the Onetwo output variables used in PLOT12, TRPLOT and other codes. Reference to ONETWOPAGES means the online documents in /u/stjohn/ONETWOPAGES/. |
| item1: | A description of the Onetwo output variables used in PLOT12, TRPLOT and other codes. Reference to ONETWOPAGES means the online documents in /u/stjohn/ONETWOPAGES/. |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode electron thermal diffusivity, cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode electron thermal diffusivity, cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion thermal diffusivity, cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion thermal diffusivity, cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion particle diffusivity, of primary ion #1,cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion particle diffusivity, of primary ion #1,cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode momentum diffusivity, g*cm/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode momentum diffusivity, g*cm/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The electron thermal conductivity,1./(cm*sec). In analysis mode this quantity is obtained from the conductive flux divided by the negative electron temperature gradient. The conductive flux is obtained from the volume integral of qconei (see qconde). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The electron thermal conductivity,1./(cm*sec). In analysis mode this quantity is obtained from the conductive flux divided by the negative electron temperature gradient. The conductive flux is obtained from the volume integral of qconei (see qconde). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode electron thermal conductivity, 1./(cm*sec). (total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode electron thermal conductivity, 1./(cm*sec). (total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The ion thermal conductivity,1./(cm*sec). In analysis mode this quantity is obtained from the conductive flux divided by the negative ion temperature gradient. The conductive flux is obtained from the volume integral of qcondi (see qcondi). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The ion thermal conductivity,1./(cm*sec). In analysis mode this quantity is obtained from the conductive flux divided by the negative ion temperature gradient. The conductive flux is obtained from the volume integral of qcondi (see qcondi). |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The neoclassical ion thermal coductivity,1/(cm*sec)cray102.f. The factor wneo(3,3) is NOT included here! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The neoclassical ion thermal coductivity,1/(cm*sec)cray102.f. The factor wneo(3,3) is NOT included here! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion thermal conductivity, 1./(cm*sec). (total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation mode ion thermal conductivity, 1./(cm*sec). (total meaning sum of all models that were selected such as neoclassical+RLW+..etc.) In analysis mode this quantity is never calculated ! |
| item0: | j=1,2..kj,i=1,2..nprtlst_max The total simulation OR ANALYSIS mode ion particle diffusivity of that primary ion species which is also present as a neutral species,cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc. in simulation mode. In analysis mode it is just the primary ion particle flux divided by the density and density gradient). |
| item1: | j=1,2..kj,i=1,2..nprtlst_max The total simulation OR ANALYSIS mode ion particle diffusivity of that primary ion species which is also present as a neutral species,cm**2/sec.(total meaning sum of all models that were selected such as neoclassical+RLW+..etc. in simulation mode. In analysis mode it is just the primary ion particle flux divided by the density and density gradient). |