mod_profiles.f90

Go to the documentation of this file.

module mod_profiles

  use itm_types
  use helena_spline

  implicit none

  integer(itm_i4), parameter :: max_points_profile = 2001

  integer(itm_i4) :: npts

  real(r8), allocatable, dimension(:) :: dp, fdf, j_tor, q_in
  target :: dp, fdf, j_tor, q_in

  real(r8), allocatable :: dpres(:), p_int(:), dgam(:), gam_int(:)
  target :: dpres, p_int, dgam, gam_int

  real(r8), allocatable :: psi_in(:)

  real(r8), allocatable :: psivec(:)

  type(spline_coefficients) :: gam_spline, dp_spline, j_spline


  contains

  subroutine allocate_spline_coefficients(spline, n)
    
    implicit none

    type(spline_coefficients) :: spline
    integer(itm_i4) :: n, alloc_error
    
    allocate(spline%sp1(n), spline%sp2(n), spline%sp3(n), spline%sp4(n), &
     stat = alloc_error)
    if (alloc_error /= 0) then
      print *, "Allocation Error"
      stop
    end if

  end subroutine allocate_spline_coefficients

  subroutine deallocate_spline_coefficients(spline)

    implicit none

    type(spline_coefficients) :: spline
    integer(itm_i4) :: dealloc_error
    
    deallocate(spline%sp1, spline%sp2, spline%sp3, spline%sp4, &
     stat = dealloc_error)
    if (dealloc_error /= 0) then
      print *, "Couldn't deallocate space for spline"
      stop
    end if

  end subroutine deallocate_spline_coefficients

end module mod_profiles