helena_initialization.f90

Go to the documentation of this file.

subroutine helena_initialization
!-----------------------------------------------------------------------
! subroutine to initialize the namelist input variables
!-----------------------------------------------------------------------
  
  use itm_types
  use mod_parameter
  use mod_corners
  use mod_dat
  use mod_dete
  use mod_gauss
  use mod_map
  use mod_mesh
  use mod_meshacc
  use mod_output
  use mod_profiles
  
  implicit none

  hbt = .false.

  input_type = "p' and FF'"
  radial_coordinate = "psi"
  current_averaging = "line"
  match = "psi_bound"
  output = "equilibrium"

  rs(1, 1) = -1._r8
  rs(1, 2) = -1._r8
  rs(2, 1) = -1._r8
  rs(2, 2) = +1._r8
  rs(3, 1) = +1._r8
  rs(3, 2) = +1._r8
  rs(4, 1) = +1._r8
  rs(4, 2) = -1._r8
  ij(1, 1) = 0
  ij(1, 2) = 0
  ij(2, 1) = 1
  ij(2, 2) = 0
  ij(3, 1) = 0
  ij(3, 2) = 1
  ij(4, 1) = 1
  ij(4, 2) = 1
  xgauss(1) = -0.861136311594053_r8
  xgauss(2) = -0.339981043584856_r8
  xgauss(3) =  0.339981043584856_r8
  xgauss(4) =  0.861136311594053_r8
  wgauss(1) =  0.347854845137454_r8
  wgauss(2) =  0.652145154862546_r8
  wgauss(3) =  0.652145154862546_r8
  wgauss(4) =  0.347854845137454_r8
  
  ias = 0
  ellip = 1.0_r8
  tria  = 0.0_r8
  quad  = 0.0_r8
  par1 = 0.0_r8
  par2 = 0.0_r8
  par3 = 1.0_r8
  par4 = 0.0_r8
  mfm = 256
  imesh = 0
  amesh = 0._r8
  bmesh = 0._r8
  isol = 0
  iarc = 0
  n_acc_points = 2
  equidistant = 1._r8

  allocate(s_acc(2))
  allocate(sig(2))
  allocate(weights(2))
  s_acc(1) = 0.5_r8
  sig(1) = 2._r8
  weights(1) = 1._r8
  s_acc(2) = 1._r8
  sig(2) = 0.2_r8
  weights(2) = 1._r8

  eps = 0.1_r8
  alfa = 0._r8
  ip = 0._r8
  q95 = 0._r8
  betap = 0._r8
  wmhd = 0._r8
  b = 0._r8
  amix = 0._r8
  bmix = 0._r8

  nr = 51
  np = 65
  nrmap = 101
  npmap = 65
  nchi  = 128
  niter = 100
  nmesh = 20
  nouter = 1

  verbosity = 1
  nbb = 4
  abb = 1._r8
  bbb = 1._r8
  ampl = 1._r8
  diagnostics_on = .false.
  eqdsk_file = .false.
  vec_file = .false.
  standard_output = .false.
  elite_output = .false.
  profiles_output = .true.
  additional_output = .false.
  xmgrace_output = .true.
  errcur = 1.e-5_r8
  errit  = 1.e-8_r8

  ne%shape = 0
  ne%h = 1._r8
  ne%h_0 = 0.1_r8
  ne%slope = 0.1_r8
  ne%psi_m = 1.0_r8
  ne%delta = 0.1_r8
  ne%psi_0 = 0.0_r8
  ne%alpha = (/ 0.0_r8, 1.0_r8, 0.0_r8/)
  ne%fac = 1._r8

  te%shape = 0        
  te%h = 1._r8
  te%h_0 = 0.1_r8
  te%slope = 0.1_r8
  te%psi_m = 1.0_r8
  te%delta = 0.1_r8
  te%psi_0 = 0.0_r8
  te%alpha = (/ 0.0_r8, 1.0_r8, 0.0_r8/)
  te%fac = 1000._r8

  edgelabe = 99.0_r8
  
  rvac = 1.0_r8
  bvac = 1.0_r8
  zeff = 1.0_r8
  zn0  = 1.0_r8
  rpe  = 0.5_r8
  etaei = 1.0_r8

  return
  
end subroutine helena_initialization