Advanced Input

An example full input called “Input_TestModBCFull” is located in the folder “input_files”, to understand what most of the “non-standard” variables do, a good understanding of tomcat is required. the variables are:

  • namecase: The name of the machine

  • shortname: short name of the machine

  • N0: core density

  • N1: edge density

  • nspeci: number of species present in the plasma

  • z: species charge numbers

  • a: species mass numbers

  • concentr: species concentrations

  • ichargeneut: speci where you enforce charge neutrality (most often majority ion)

  • iforcechargeneut: boolean forcing charge neutrality computation

  • IAUTO: boolean dispersion equation automatic root ordering

  • IL & IR: if IAUTO zero, specify order of root ordering.

  • ap: minor radius

  • R0: major radius

  • IDSHAPE,xkap,del0,delt: machine definitions, a boolean to enable D shape, kappa (elongation), Shafranov shift of the magnetic axis and triangularity

  • T0: species core temperature

  • T1: species edge temperature

  • exponn= exponential in density profile

  • expont= exponential in temperature profile

  • frek: frequency in hertz

  • ntor: toroidal mode number

  • xleft/right: start and stop point simulation

  • xleftin/xrigthin: inner domain, between xleft and xleftin is a detailed region with high resolution solving, similar for xright and xrightin.

  • numLFS: number of elements on the LFS space inbetween xright and xrightin

  • numHFS: number of elements on the LFS space inbetween xleft and xleftin

  • zleft/right: Integration interval is from (xleft,zleft) to (xright,zright)

  • iautobreak: boolean to enable automatic determination of density, temperature and their gradients to be continuous across the matching point near the edge.

  • xDensL/R, DensL/R ,rhoLN/RN, xTempR/L ,TempR/L, RhoLT/RT: If iautobreak=0, the x position where density continuity is enforced, the corresponding density, the magnetic surface labeling parameter, x position where temperature continuity is enforced, the corresponding temperature and magnetic surface labeling on the right side.

  • LambdaT: exponential decay length of the temperature in the scrape off layer

  • LambdaN: exponential decay length of the density in the scrape off layer

  • icollis= boolean to enable collisions

  • colfreq= collision frequency, per species

  • iasymp= asymptotic expansion of plasma dispersion function used or not

  • ky= wavevector in the y (Z) direction

  • IDISP: Enable root finding everywhere in the plasma.

  • IHFS: Antenna at the HFS or not

  • ILAUNCH: wave excitation based on incoming root number (standard:1)

  • inewBC: 1 (standard) when boundary conditions are directly imposed as equations, 0 when they are imposed via the edge system equations and 2 as 1 but allowing to launch serveral modes at the same time. Two vectors SOURCEL and SOURCER are included as input which give the amplitudes of the six ‘launched’ modes. These modes are pair nrs for HFS and impair numbers for LFA ‘launch’. The remaining 6 modes’ amplitudes follow from the integration itself. By manipulating the vectors IL and IR, the amplitude of up to 6 modes can be prescribed at one side.

  • xANT: x position of antenna

  • JyANT: current density on antenna in the poloidal direction

  • JfiANT: current density on antenna in the toroidal direction.

  • SOURCEL: eigenvector complex amplitude of corresponding excited wave, used if IHFS is 1

  • SOURCEL: eigenvector complex amplitude of corresponding excited wave, used if IHFS is 0

  • ialpha: consider self-produced alpha particle profile (assumes DT plasma is being studied)

  • ibeam: simulate a neutral beam profile (should be left on when you don’t want charge neutrality violated, can keep NBEAMMAX and TBEAMMAX low for no effect)

  • W0N: half width of gaussian profile for beam density

  • W0T: half width of gaussian profile for beam temperature

  • RHOBEAMN: amplitude of gaussian profile for beam density

  • RHOBEAMT: amplitude of gaussian profile for beam temperature

  • NBEAMMAX: maximum gaussian beam density per species

  • TBEAMMAX: maximum gaussian beam temperature per species

  • irhophi: switch between cartesian and polar coordinates

  • ixors: switch between real cartesian x or normalized variable s along integration path

  • iforceunityS: force flux at excitation point to one.