/
parameters
executable file
·55 lines (55 loc) · 2.97 KB
/
parameters
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
#****************************************************************************
#Parameter file for SPINNAKER
#****************************************************************************
#Setup of the 2-dimensional grid
#----------------------------------------------------------------------------
grid_size = 200#<=200; Number of grip points in z
nu_channel = 400#<=400; Number of grid points in frequency
grid_delta = 10#Every x_th grid point is written to output file
z_halo = 8.#Halo size [kpc]
first_data_point_at_0kpc = 1#First data point at z=0 kpc? 1:yes -1: no
normalize_intensities = -1#At z=0, I(nu_i)=1, 1:yes -1: no
#----------------------------------------------------------------------------
#Output
#----------------------------------------------------------------------------
nu_1 = 1.37e9# [Hz]
nu_2 = 4.86e9# [Hz]
nu_3 = 6.20e9# [Hz]
nu_4 = 8.40e9# [Hz]
mode = 1#Advection: 1; Diffusion: 2; Cylindrical diffusion: 3; Radial diffusion: 4
epsilon = -1#1:write emissivities instead of intensities (only if model=1)
FWHM_effective_beam = 1.2#Effective beam size as FWHM [kpc]
#----------------------------------------------------------------------------
#Setup of the advection and diffusion model
#----------------------------------------------------------------------------
gamma_in = 2.7#Injection CRe index
rad_field = 0.11#Radiation energy density = U_IRF/U_B
V0 = 300.0e5#Advection speed [cm s^-1]
velocity_field = 0#0: constant, -1: radius power-law, 1: exponential V = V0*exp(z/h_V)
#velocity_field = 2: power-law with z, V = V0 *(1 + (z/h_V)^beta)
h_V = 1.#[kpc]
adiabatic_losses = -1#Yes: 1, No: -1 (only for advection)
D0 = 3.0e28#Diffusion coefficient [cm^2 s^-1]
mu_diff = 0.5#Energy dependence D=D0*E^mu_diff
V_rot = 187.e5#rotational velocity [cm s^-1]
#----------------------------------------------------------------------------
#Magnetic field setup
#----------------------------------------------------------------------------
galaxy_mode = 1#1: thin and thick disc, -1: piece-wise exponential
#2: constant magnetic field (isotropic turbulent)
z1 = 10.0#if galaxy_mode=-1, transition from zone 1 to zone 2 [kpc]
B0 = 13.5e-6#Magnetic field strength [G], 1mikroGauss=1.0e-10Tesla
B1 = 8.6e-6#Magnetic field strength [G], 1mikroGauss=1.0e-10Tesla
h_B1 = 0.6#Scaleheight of the B-field in zone 1 [kpc]
h_B2 = 5.0#Scaleheight of the B-field in zone 2 [kpc]
beta = 1.0#Power-law index for the velocity, V = V0 * (R/R0)^-beta
R0 = 1.04e22#Jet radius at the base [cm]
#----------------------------------------------------------------------------
#Use a magnetic field model (needs edit of the source files)
#----------------------------------------------------------------------------
model = -1#-1: no model, 1: model (which needs edit of the source files)
initialize_model = -1#Initialze magnetic field model
model_north = 1#Use northern part of the model
update_model = -1#Rescale magnetic field model
xi = 3.65#Model normalization
#EOF*************************************************************************