write_parameters.f90

Go to the documentation of this file.

subroutine write_parameters
!-----------------------------------------------------------------------
! subroutine to write the output to unit 20
!-----------------------------------------------------------------------
  
  use itm_types
  use mod_dat
  use mod_map
  use mod_mesh
  use mod_output
  use mod_profiles
  use mod_helena_io, only : out_he
  
  implicit none
  
  write(out_he, *) '*****************************************'
  write(out_he, 10)
  write(out_he, *) '*****************************************'
  write(out_he, *) '*          HELENA                       *'
  write(out_he, *) '*****************************************'
  
  write(out_he, *)
  write(out_he, *)
  write(out_he, 20) ellip, tria, quad
  write(out_he, 21) mharm, isol, ias
  write(out_he, *)
  write(out_he, 30) hbt
  write(out_he, 31) input_type
  write(out_he, 32) radial_coordinate
  if (trim(input_type) == "p' and j_tor" &
   .or. trim(input_type) == "p and j_tor" &
   .or. trim(input_type) == "p' and q") then
    write(out_he,36) current_averaging
  end if
  write(out_he, *)
  write(out_he, 40) eps, alfa, b
  write(out_he, 41) ip, betap
  write(out_he, *)
  write(out_he, 50) nr, np, nrmap, npmap, nchi, niter
  write(out_he, *)
  if (trim(input_type) == "p' and j_tor" &
   .or. trim(input_type) == "p and j_tor" &
   .or. trim(input_type) == "p' and q") then
    write(out_he, 51) nmesh, npts, errcur
  end if
  write(out_he, *)
  write(out_he, *)
  
 10 format(' *  helena : grad-shafranov equilibrium  *')
 20 format(' $shape    ellip =', e16.8, ', tria  =', e16.8, ', quad  =', &
     e16.8)
 21 format('           mharm =', i3, ', isol =', i2, ', ias = ', i2)
 30 format(' $profile  hbt =', l7, ',')
 31 format('           input_type = ', a12, ',')
 32 format('           radial_coordinate = ', a6, ',')
 36 format('           averaging type: current_averaging = ', a6)
 40 format(' $phys     eps  =', e16.8, ', alfa =', e16.8, &
     ', b =', e16.8, ',')
 41 format('           Ip =', e16.8, ', betap = ', e16.8)
 50 format(' $num      nr =', i3, ', np =', i3, ', nrmap =', i3, &
     ', npmap =', i3, ', nchi = ', i3, ', niter = ', i3)
 51 format('           nmesh =', i3, ', npts =', i4, ', errcur =', e16.8)
  
  return
end subroutine write_parameters