write_equilibrium.f90

Go to the documentation of this file.

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

  use itm_types
  use mod_map
  use mod_mesh, only : nr

  use euitm_schemas
  use mod_helena_io

  implicit none

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

  integer(itm_i4) :: i, j

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

!-- dimensions
  write(out_map, 1) nr, nchi

!-- global parameters
  write(out_map, 2) equilibrium_out%eqgeometry%geom_axis%r
  write(out_map, 2) equilibrium_out%eqgeometry%a_minor
  write(out_map, 2) equilibrium_out%global_param%mag_axis%position%r
  write(out_map, 2) equilibrium_out%global_param%mag_axis%bphi

!-- 1d profiles
  write(out_map, 3) (equilibrium_out%profiles_1d%psi(i), i = 1, nr)
  write(out_map, 3) (equilibrium_out%profiles_1d%pressure(i), i = 1, nr)
  write(out_map, 3) (equilibrium_out%profiles_1d%F_dia(i), i = 1, nr)
  write(out_map, 3) (equilibrium_out%profiles_1d%q(i), i = 1, nr)
  write(out_map, 3) (equilibrium_out%profiles_1d%pprime(i), i = 1, nr)

!-- boundary
  write(out_map, 3) (equilibrium_out%eqgeometry%boundary(1)%r(i), i = 1, &
   nchi)
  write(out_map, 3) (equilibrium_out%eqgeometry%boundary(1)%z(i), i = 1, &
   nchi)

!-- grid
  write(out_map, 3) ((equilibrium_out%coord_sys%position%r(i, j), &
   i = 1, nr), j = 1, nchi)
  write(out_map, 3) ((equilibrium_out%coord_sys%position%z(i, j), &
   i = 1, nr), j = 1, nchi)
  write(out_map, 3) ((equilibrium_out%coord_sys%g_11(i, j), &
   i = 1, nr), j = 1, nchi)
  write(out_map, 3) ((equilibrium_out%coord_sys%g_12(i, j), &
   i = 1, nr), j = 1, nchi)
  write(out_map, 3) ((equilibrium_out%coord_sys%g_33(i, j), &
   i = 1, nr), j = 1, nchi)

  close(out_map)

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

end subroutine write_equilibrium