|
|
|
## Code parameters
|
|
|
|
|
|
|
|
A plain ASCII text file is used to pass as input:
|
|
|
|
- code-specific parameters
|
|
|
|
- names of the other input files
|
|
|
|
- beam definition at launch
|
|
|
|
- magnetic equilibrium
|
|
|
|
- kinetic profiles
|
|
|
|
- the EC beam parameters intended as real-time adjustable:
|
|
|
|
- launch angles
|
|
|
|
- injected power
|
|
|
|
- polarization
|
|
|
|
|
|
|
|
The name of such file defaults to `gray_params.data`.
|
|
|
|
|
|
|
|
The values must be blank separated, and strings must be enclosed in quotes if they contain spaces or slashes. The following Fortran code snippet is used to read the file:
|
|
|
|
|
|
|
|
```fortran
|
|
|
|
use iso_fortran_env, only: real64 => wp
|
|
|
|
integer :: u
|
|
|
|
integer :: nrayr, nrayth, nstep, idst, igrad, ipass, ipol
|
|
|
|
integer :: iwarm, ilarm, imx, ieccd, iox
|
|
|
|
integer :: ibeam, iequil, iprof, ipsinorm, idesc, ifreefmt
|
|
|
|
integer :: icocos, sgnb, sgni, irho, ixp, iscal
|
|
|
|
integer :: ipec, nrho, istpr, istpl
|
|
|
|
real(wp) :: rwmax, dst
|
|
|
|
real(wp) :: alpha, beta, power, psipol, chipol
|
|
|
|
real(wp) :: rwall, psnbnd, ssplps, sspld, factb, factte, factne
|
|
|
|
character(len=256) :: filenmbm, filenmeq, filenmprf
|
|
|
|
|
|
|
|
read(unit=u, *) nrayr, nrayth, rwmax
|
|
|
|
read(unit=u, *) igrad, ipass, ipol
|
|
|
|
read(unit=u, *) dst, nstep, idst
|
|
|
|
read(unit=u, *) iwarm, ilarm, imx
|
|
|
|
read(unit=u, *) ieccd
|
|
|
|
read(unit=u, *) alpha, beta
|
|
|
|
read(unit=u, *) power
|
|
|
|
read(unit=u, *) iox, psipol, chipol
|
|
|
|
read(unit=u, *) ibeam
|
|
|
|
read(unit=u, *) filenmbm
|
|
|
|
read(unit=u, *) iequil
|
|
|
|
read(unit=u, *) filenmeq
|
|
|
|
read(unit=u, *) icocos, ipsinorm, idesc, ifreefmt
|
|
|
|
read(unit=u, *) ixp, ssplps
|
|
|
|
read(unit=u, *) sgnb, sgni, factb
|
|
|
|
read(unit=u, *) rwall
|
|
|
|
read(unit=u, *) iprof, irho
|
|
|
|
read(unit=u, *) filenmprf
|
|
|
|
read(unit=u, *) psnbnd, sspld
|
|
|
|
read(unit=u, *) factte, factne, iscal
|
|
|
|
read(unit=u, *) ipec, nrho
|
|
|
|
read(unit=u, *) istpr, istpl
|
|
|
|
```
|
|
|
|
|
|
|
|
Any text appearing after the last value expected on any line is ignored.
|
|
|
|
|
|
|
|
The following table describes the meaning of each parameter, and the type of value expected.
|
|
|
|
|
|
|
|
| Line | Field | Parameter | Type | Units | Description |
|
|
|
|
| ---: | ----: | :---------- | :------ | :---: | :--------------------------- |
|
|
|
|
| 1 | 1 | `nrayr` | Integer | - | Radial number of rays in the beam |
|
|
|
|
| 1 | 2 | `nrayth` | Integer | - | Angular number of rays in the beam |
|
|
|
|
| 1 | 3 | `rwmax` | Real | - | Beam radius cut-off: outermost rays at `rwmax`$`\cdot w`$, with beam radius $`w`$ |
|
|
|
|
| 2 | 1 | `igrad` | Integer | - | 0: raytracing; 1: beamtracing |
|
|
|
|
| 2 | 2 | `ipass` | Integer | - | |
|
|
|
|
| 2 | 3 | `ipol` | Integer | - | |
|
|
|
|
| 3 | 1 | `dst` | Real | cm | Beam trajectory integration step |
|
|
|
|
| 3 | 2 | `nstep` | Integer | - | Maximum number of steps |
|
|
|
|
| 3 | 3 | `idst` | Integer | - | Integration variable. 0: $`s`$ ; 1: $`ct`$ ; 2: $`S_R`$ |
|
|
|
|
| 4 | 1 | `iwarm` | Integer | - | Absorption model. 0: no absorption; 1: weakly relativistic; 2: fully relativistic |
|
|
|
|
| 4 | 2 | `ilarm` | Integer | - | Absorption: expansion order in Larmor radius |
|
|
|
|
| 4 | 3 | `imx` | Integer | - | Absorption: maximum iteration number for convergence |
|
|
|
|
| 5 | 1 | `ieccd` | Integer | - | CD model. 0: no CD; 1: Cohen; 11: Marushchenko (momentum conservation) |
|
|
|
|
| 6 | 1 | `alpha` | Real | deg | Poloidal injection angle |
|
|
|
|
| 6 | 2 | `beta` | Real | deg | Toroidal injection angle |
|
|
|
|
| 7 | 1 | `power` | Real | MW | Beam power at launch |
|
|
|
|
| 8 | 1 | `iox` | Integer | - | Polarization. 1: O-mode; 2: X-mode |
|
|
|
|
| 8 | 2 | `psipol` | Real | deg | |
|
|
|
|
| 8 | 3 | `chipol` | Real | deg | |
|
|
|
|
| 9 | 1 | `ibeam` | Integer | - | Beam shape variation in `filenmbm`<br>0: fixed;<br>1: $`(\alpha,\beta,x_0,w,R_c)=f`$(1D steering);<br>2: $`(x_0,w,R_c)=f(\alpha, \beta)`$ |
|
|
|
|
| 10 | 1 | `filenmbm` | String | | Input file for beam description |
|
|
|
|
| 11 | 1 | `iequil` | Integer | - | Magnetic equilibrium in `filenmeq`<br> 0: vacuum; 1: analytical; 2: numerical |
|
|
|
|
| 12 | 1 | `filenmeq` | String | | Input file for magnetic equilibrium |
|
|
|
|
| 13 | 1 | `icocos` | Integer | - | [COCOS][cocos] used in `filenmeq`<br>0 or 10 to force $`B_\phi`$ and $`I_\phi`$ direction |
|
|
|
|
| 13 | 2 | `ipsinorm` | Integer | - | |
|
|
|
|
| 13 | 3 | `idesc` | Integer | - | |
|
|
|
|
| 13 | 4 | `ifreefmt` | Integer | - | |
|
|
|
|
| 14 | 1 | `ixp` | Integer | - | |
|
|
|
|
| 14 | 2 | `ssplps` | Real | - | |
|
|
|
|
| 15 | 1 | `sgnb` | Integer | - | $`B_\phi`$, sign in COCOS=3<br>Used if `icocos` == 0 or 10 |
|
|
|
|
| 15 | 2 | `sgni` | Integer | - | $`I_\phi`$, sign in COCOS=3<br>Used if `icocos` == 0 or 10 |
|
|
|
|
| 15 | 3 | `factb` | Real | - | Magnetic field scaling factor<br>$`B'=B \cdot`$`factb` |
|
|
|
|
| 16 | 1 | `rwall` | Real | m | Major radius of inner first wall<br>Must be `rwall` < min($`R`$) along the LCFS |
|
|
|
|
| 17 | 1 | `iprof` | Integer | - | Profiles in `filenmprf`<br> 0: analytical; 1: numerical |
|
|
|
|
| 17 | 2 | `irho` | Integer | - | Minor radius in `filenmprf`<br> 0: $`\rho_{\rm tor}`$; 1: $`\rho_{\rm pol}`$; 2: $`\psi_{\rm norm}`$ |
|
|
|
|
| 18 | 1 | `filenmprf` | String | | Input file for kinetic profiles |
|
|
|
|
| 19 | 1 | `psnbnd` | Real | - | |
|
|
|
|
| 19 | 2 | `sspld` | Real | - | |
|
|
|
|
| 20 | 1 | `factte` | Real | - | Temperature scaling factor<br>$`T_e'=T_e \cdot`$`factte`$`\cdot`$`factb`$`^{\sigma_T}`$ |
|
|
|
|
| 20 | 2 | `factne` | Real | - | Density scaling factor<br>$`n_e'=n_e \cdot`$`factne`$`\cdot`$`factb`$`^{\sigma_n}`$ |
|
|
|
|
| 20 | 3 | `iscal` | Integer | - | $`T_e`$, $`n_e`$ scaling law with $`B`$<br>0: $`\nu^*`$=const ($`\sigma_T`$=2/3, $`\sigma_n`$=4/3);<br>1: $`n/n_{\rm GW}`$=const ($`\sigma_T`$=$`\sigma_n`$=1);<br>2: none ($`\sigma_T`$=$`\sigma_n`$=0) |
|
|
|
|
| 21 | 1 | `ipec` | Integer | - | |
|
|
|
|
| 21 | 2 | `nrho` | Integer | - | |
|
|
|
|
| 22 | 1 | `istpr` | Integer | - | |
|
|
|
|
| 22 | 2 | `istpl` | Integer | - | |
|
|
|
|
|
|
|
|
[cocos]: <http://dx.doi.org/10.1016/j.cpc.2012.09.010> |
|
|
|
\ No newline at end of file |
