mod_output.f90

Go to the documentation of this file.

module mod_output

  use itm_types

  implicit none

  character(len = 100) :: txtout(40)
  logical :: diagnostics_on
  logical :: standard_output, elite_output
  logical :: profiles_output, additional_output
  logical :: xmgrace_output
  logical :: vec_file, eqdsk_file

  contains

  subroutine write_equidistant_profiles(equilibrium, pskn1, nr)
!-----------------------------------------------------------------------------
! subroutine to write out high resolution equidistant profiles for a re-run
!  uses equilibrium CPO and transfer_1d
!-----------------------------------------------------------------------------

    use mod_dat, only : output, cpsurfin
    use plot_data
    use mod_profiles
    use mod_helena_io
    use mod_interpolation
    use euitm_schemas

    implicit none

!-- number of grid points for equidistant radial grid
    integer(itm_i4), parameter :: npts_out = 2000
    integer(itm_i4), parameter :: profile_out = 22

    type (type_equilibrium), intent(in) :: equilibrium
    real(r8), intent(in) :: pskn1(nr)
    integer(itm_i4), intent(in) :: nr

    real(r8), dimension(npts_out) :: ps_equi, dp_equi, fdf_equi, p_equi, j_equi
    integer(itm_i4) :: i

    do i = 1, npts_out
      ps_equi(i) = pskn1(1) + dble(i - 1) / dble(npts_out - 1) &
       * (pskn1(size(pskn1)) - pskn1(1))
    end do

!-- transfer profiles to equidistant psi grid
    call transfer_1d_profile(pskn1, equilibrium%profiles_1d%jphi, &
     ps_equi, j_equi, spline_type = 1)
    call transfer_1d_profile(pskn1, equilibrium%profiles_1d%pprime, &
     ps_equi, dp_equi, spline_type = 1)
    call transfer_1d_profile(pskn1, equilibrium%profiles_1d%ffprime, &
     ps_equi, fdf_equi, spline_type = 1)
    call transfer_1d_profile(pskn1, equilibrium%profiles_1d%pressure, &
     ps_equi, p_equi, spline_type = 1)

!-- write out equidistant profiles for re-run
    if (trim(output) == 'profiles' .or. trim(output) == 'full') then
      open (unit = profile_out, file = trim(adjustl(path)) // 'psi.out', &
       status = 'replace', form = 'formatted', action = 'write', &
       iostat = i_error)
      if (i_error /= 0) then
        write(iu6, *) ' ERROR writing psi.out '
      end if
      write(profile_out, 26) (ps_equi(i) * cpsurfin, i = 1, npts_out)
      close(profile_out)

      open (unit = profile_out, file = trim(adjustl(path)) // 'dp.out', &
       status = 'replace', form = 'formatted', action = 'write', &
       iostat = i_error)
      if (i_error /= 0) then
        write(iu6, *) ' ERROR writing dp.out '
      end if
      write(profile_out, 26) (dp_equi(i), i = 1, npts_out)
      close(profile_out)

      open (unit = profile_out, file = trim(adjustl(path)) // 'p.out', &
       status = 'replace', form = 'formatted', action = 'write', &
       iostat = i_error)
      if (i_error /= 0) then
        write(iu6, *) ' ERROR writing p.out '
      end if
      write(profile_out, 26) (p_equi(i), i = 1, npts_out)
      close(profile_out)

      open (unit = profile_out, file = trim(adjustl(path)) // 'fdf.out', &
       status = 'replace', form = 'formatted', action = 'write', &
       iostat = i_error)
      if (i_error /= 0) then
        write(iu6, *) ' ERROR writing fdf.out '
      end if
      write(profile_out, 26) (fdf_equi(i), i = 1, npts_out)
      close(profile_out)

      open (unit = profile_out, file = trim(adjustl(path)) // 'j_tor.out', &
       status = 'replace', form = 'formatted', action = 'write', &
       iostat = i_error)
      if (i_error /= 0) then
        write(iu6, *) ' ERROR writing j_tor.out '
      end if
      write(profile_out, 26) (j_equi(i), i = 1, npts_out)
      close(profile_out)

    end if

 26 format(e24.16)

  end subroutine write_equidistant_profiles

end module mod_output