write_old_equilibrium.f90

Go to the documentation of this file.

subroutine write_old_equilibrium(path, equilibrium_out)
!----------------------------------------------------------------------------
!  This subroutine writes out the resulting HELENA equilibrium for use
!  by ILSA.
!----------------------------------------------------------------------------

  use itm_types
  use mod_dat, only : zgeo
  use mod_map
  use mod_mesh, only : nr
  use mod_parameter, only : twopi, mu0

  use euitm_schemas
  use mod_helena_io

  implicit none

  type (type_equilibrium) :: equilibrium_out
  character(len = 132) :: path 

  character (len=30), parameter :: file_old_map_out = 'old_equilibrium'

  real(r8) :: rmag, bmag
  integer(itm_i4) :: i, j

  rmag = equilibrium_out%global_param%mag_axis%position%r
  bmag = equilibrium_out%global_param%mag_axis%bphi

!-------------------- write output for ILSA ---------------------------
  open (unit = out_map, file = trim(adjustl(path)) // file_old_map_out, &
   status = 'replace', form = 'formatted', action = 'write', iostat = i_error)

  write(out_map, 3)  nr - 1, nchi
  write(out_map, 1)  (cs(i), i = 1, nr)
  write(out_map, 1)  (equilibrium_out%profiles_1d%q(i), i = 1, nr)
  write(out_map, 2)  0._r8, dqs(nr)
  write(out_map, 1)  (dqs(i), i = 2, nr)
  write(out_map, 1)  (0._r8, i = 1, nr)
  write(out_map, 2)  0._r8, 0._r8
  write(out_map, 1)  (chi(i), i = 1, nchi)
  write(out_map, 1)  ((equilibrium_out%coord_sys%g_11(i, j) &
   / (twopi * rmag * bmag)**2, j = 1, nchi), i = 2, nr)
  write(out_map, 1)  ((equilibrium_out%coord_sys%g_12(i, j) &
   / (twopi * bmag), j = 1, nchi), i = 2, nr)
  write(out_map, 2)  equilibrium_out%profiles_1d%psi(nr) / (twopi &
   * bmag * rmag**2),  equilibrium_out%eqgeometry%a_minor / rmag
  write(out_map, 1)  ((1._r8 / equilibrium_out%coord_sys%g_33(i, j) &
   / rmag**2, j = 1, nchi), i = 2, nr)
  write(out_map, 4)  1._r8
  write(out_map, 1)  (equilibrium_out%profiles_1d%pressure(i) &
   / bmag**2 * mu0, i = 1, nr)
  write(out_map, 2)  0._r8, dp0(nr)
  write(out_map, 1)  (equilibrium_out%profiles_1d%F_dia(i) / (rmag &
   * bmag), i = 1, nr)
  write(out_map, 2)  0._r8, 0._r8
!----------------------------------------- additional data for vacuum --
  write(out_map, 1) ((equilibrium_out%eqgeometry%boundary(1)%r(j) &
   - equilibrium_out%eqgeometry%geom_axis%r) &
   / equilibrium_out%eqgeometry%a_minor, j = 1, nchi)
  write(out_map, 1) ((equilibrium_out%eqgeometry%boundary(1)%z(j) &
   - zgeo) / equilibrium_out%eqgeometry%a_minor, j = 1, nchi)
  write(out_map, 4) equilibrium_out%eqgeometry%a_minor &
   / equilibrium_out%eqgeometry%geom_axis%r
  write(out_map, 1) (((equilibrium_out%coord_sys%position%r(i, j) &
   - equilibrium_out%eqgeometry%geom_axis%r) &
   / equilibrium_out%eqgeometry%a_minor, j = 1, nchi), i = 2, nr)
  write(out_map, 1) (((equilibrium_out%coord_sys%position%z(i, j) &
   - zgeo) / equilibrium_out%eqgeometry%a_minor, j = 1, nchi), &
   i = 2, nr)
  close(out_map)

1 format(4e16.8)
2 format(2e16.8)
3 format(2i6)
4 format(e16.8)

end subroutine write_old_equilibrium