ExoCross keywords - Trovemaster/exocross GitHub Wiki

The position of the keywords is not important. The input is not key-sensitive.

  • GNS: Nuclear statistical weight gns 2.0

  • Temperature: Temperature (K):

Temperature 3000.0

Aliases: Temp, Temperature

  • Temperature can be split into rotational temperature and vibrational temperature, which in input is given as

    Temperature  (Rotational) 700 Vib (vibrational) 2000 (K)

In this case it is also important to provide the addresses of the vibrational columns using the QN section.

  • QN: The designation of the quantum numbers (QN) to columns can be done using the QN section. In this section vib is used for the range of the vibrational QNs; K is for the column with the rotational "K"; Sym is for symmetry; Nsym is number of the symmetries; Nmodes stands for the number of vibrational modes. The column numbering starts after the J column (4th). The vibrational addressing is required for when vibrational temperature is used.

    QN
  Vib 2 7 
  Nmodes 9 
  K  7
end

  • PF: input value of the partition function.

PF 1000.0

  • Range: Wavenumber (cm-1) or wavelength (um or micron) range, i.e. from nu1 to nu2:

Range 0.0 2000.000

Here is example for microns (wavelength), which however does nor work will all profiles:

 Range 0.0 2000.000 um

Aliases: Window, Frequencies, Wavenumbers, Range

  • um or micron: units of wavelength used in combination with Range

  • cm-1: units of wavenumber used in combination with Range, can be omitted as default.

  • abundance: is the molecular abundance, which is the factor used in intensities or cross-sections.

    abundance 0.97

  • Cut-off: to use the HITRAN-cut-off scheme (see HITRAN 2012 paper).
    nu_crit is to give the corresponding parameter in the HITRAN cut-off scheme.

    cutoff HITRAN (S_crit) 1e-29 nu_crit 2000

  • Npoints: Number of grid points (usually an odd number)

Npoints 100001

Aliases: Npoints, Number-of-points

  • Absorption: type of spectra

Absorption

  • Emission: type of spectra`

Emission

  • lifetime: type of spectra, to compute life time of Ej as 1/sum_i A_ji; saved into a new .life file.

Only one of these three types is needed, it can appear outside any sections anywhere in the program. 'Absorption' is default.

  • Partfunc or Partition-function: To compute the partition function and its contributions for a set of temperature; it also starts a section defining the corresponding parameters
  • Ntemps: Number of temperature steps

Ntemps 500

  • TempMax, maxtemps: Maximal temperature

TempMax 5000.0

  • Moment: Moment to compute (0=none,1=Q1,2=Q2,3=CP)

moment 1

  • CP: the same as Moment 3

CP

  • output: File name for crosssections, stick-spectra etc.

output CH4_stick_T300K

  • Cooling is a section name to compute the cooling function on a grid of temperatures. The parameters used in the Cooliing section are TempMax (maxtemps) and Ntemps, see Partfunc.

  • States, StatesFile, StateFile, States_file, States-File: The name of the .states file States CH4_linelist.states

  • Transitions, TransitionFiles: The name of the .trans file or the section name for the list of the .trans files.

Transitions CH4_linelist.trans

or

Transitions
 CH4_100.trans
 CH4_200.trans
 CH4_300.trans
 CH4_400.trans
end

In combination with the keyword histogram-J the transition filename can be followed by the J-value this file is associated with.

  • Threshold: Intensity threshold. A line is skipped from line profile evaluation, or simply from the output if the corresponding absorption coefficient/emissivity is smaller than Threshold.

Threshold 1e-28

  • Enermax: Energy Threshold used to select states below some energy value

enermax 20000.0

  • Gaussian, Gauss, Gauss0, Doppler, Doppl, Doppl0, Box, Bin, Rect, Sticks, Stick, Voigt, pseudo, pse-Rocco, pse-Liu, Voig-Quad: The type of the line profile.

gauss

Doppler

stick

bin

  • HWHM, HalfWidth: Half width at the half maximum, used for Gauss, Lorenz. HWHM 0.1

  • Mass, Masses: Effective molecular mass (u0), used for Doppl. mass 16.0

  • Verbose: verbose level, defines the amount of data to be printed out in the output Verbose 3

  • histogram computing crosssections using intensities saved on (usually) an equidistant grid.

  • histogram-J computing cross-sections using intensities saved on a grid with J-dependent histogram files.

  • NRAM or Ncash: number of transitions to put into RAM; Alias is LINES-TO-CASH

  • Nprocs: Number of OMP processes. Ideally should be the same as the number of omp-processors allocated, but would work with any number. Aliases are OMP_NUM_PROCS and OMP_PROCS.

  • mem: maximal memory allocated for the job. It is used to estimate how many transitions can be put into RAM. Should be less than or equal too the memory of the system.

  • ioffset is a in integer factor to compute the wavenumber offset in units of HWHM: offset = ioffset*hwhm

  • offset is a factor to compute the wavenumber offset in cm-1

    offset 25 (cm-1)

  • pressure value in bar must be specified (otherwise P=1 bar is assumed) used for Voigt.

  • Species or Broadener starts a section to define the Voigt-type broadening parameters

    gamma(Voigt) = \sum_i gamma_i (T0_i/T)^n P/P0_i ratio_i

gamma or gamma0 is the reference HWHM (cm-1) for Voigt used in Species

n is the exponent n_i for Voigt used in Species

delta is the pressure shift.

T0 is the reference T (K), usually 298 for Voigt used in Species

P0 is the reference pressure in bar, usually 1 for Voigt used in Species

ratio is the mixing ratio of the species (unitless) for example the solar mixing ratio of H2 and He is 0.9 and 0.1 (species).

file is the name of the file with broadening parameters

model is the broadening model, e.g. J (alias a0) or JJ (alias a1).

nquad is the number of quadrature points used for Voigt-Quad

The name of the species should be the first thing on the line.

Example:

 mass 16.0
 pressure 0.5 
 Temperature 1300.0 
 Species
   H2  gamma 0.05 n 0.4 t0 298.0 ratio 0.9
   He  gamma 0.04 n 1.0 t0 298.0 ratio 0.1
 end

or

 Species
   H2  gamma 0.0207 n 0.44 t0 298.0 file 1H2-16O__H2__a1.broad model JJ ratio 0.84
   He  gamma 0.043 n 0.02 t0 298.0 file  1H2-16O__He__a1.broad model JJ ratio 0.16
 end

  • Phoenix is the keyword for converting ExoMol line list to the Phoenix format. The species are expected to specify the Voigt parameters of the broadeners. Phoenix should appear anywhere in the main body of the input in the way as a line profile keyword.

  • HITRAN is to use the HITRAN-format of the transition file or output. Reading from hitran (.par) requires also the definition of the partition function pf and the isotopologue number iso. No .states is needed. To read from HITRAN use HITRAN READ

iso is the isotopologue number of a HITRAN molecule, e.g. 261 for 12C2H2. Example

hitran
iso 261
pf 1000.0
output C2H2_ab_g0.5
Transitions  26_hit12.par

HITRAN` can also form a section where there HITRAN-error codes are specified, e.g.

 hitran
  write
  error-E  qn 4 ierr 4  vmax  10 ierr 3  vmax  20 ierr 2  vmax  30  ierr 1  vmax  40 ierr 0  vmax  100 
  error-S  qn 4 ierr 5  vmax  10 ierr 4  vmax  20 ierr 3  vmax  100 
  error-Air   ierr 4
 end

error-E and error-S are used to specify the ranges for the quantum numbers for different Energy and intensity error-codes, respectively. qn is the quantum number (the number of the QN-column after J) used for the error-specification. ierr is followed by the error code, followed by vmax (keyword) followed by the maximum value of qn this error code applies, which is followed by another error code. error-Air ierr to give the error code for the HITRAN air-broadening (one constant value).
error-self ierr to give the error code for the HITRAN self-broadening (one constant value). error-n ierr to give the error code for the HITRAN n-broadening (one constant value). error-delta ierr to give the error code for the HITRAN line-shift (one constant value).