calculate_globals.f90

Go to the documentation of this file.

SUBROUTINE calculate_globals (GEOMETRY, PROFILES, SOURCES, GLOBAL, DIAG)
!-------------------------------------------------------!
!     This routine calculates global plasma             !
!     parameters based on profiles of plasma            !
!     parameters and equilibrium quantities.            !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     Denis.Kalupin@euro-fusion.org       !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!


    USE itm_types
    USE itm_constants
    USE ets_plasma
    USE ets_math

    IMPLICIT NONE


! +++ Internal Transport Solver derived types:
    TYPE (magnetic_geometry)           :: geometry            !contains all geometry quantities
    TYPE (plasma_profiles)             :: profiles            !contains profiles of plasma parameters
    TYPE (sources_and_sinks)           :: sources             !contains profiles of sources
    TYPE (global_param)                :: global              !contains global plasma parameters
    TYPE (diagnostic)                  :: diag                !contains error messages and warnings

! +++ Temporary functions:
    REAL(R8), DIMENSION(PROFILES%NRHO) :: fun, intfun
    INTEGER                            :: nrho


    nrho               =  geometry%NRHO

    fun                =  geometry%VPR * profiles%CURR_PAR / 2.0e0_r8 / itm_pi * geometry%BGEO / geometry%FDIA**2
    CALL integr_fun(nrho, geometry%RHO, fun, intfun)
    global%CURRENT     =  intfun(nrho) * geometry%FDIA(nrho)


    fun                =  geometry%VPR * (sources%CURR_EXP + sources%CURR_IMP * profiles%PSI) / 2.0e0_r8 / itm_pi * geometry%BGEO / geometry%FDIA**2
    CALL integr_fun(nrho, geometry%RHO, fun, intfun)
    global%CURRENT_NI  =  intfun(nrho) * geometry%FDIA(nrho)
    
    fun                =  profiles%QOH * geometry%VPR
    CALL integr_fun(nrho, geometry%RHO, fun, intfun)
    global%QOH         =  intfun(nrho)     

    global%VLOOP       = profiles%PSI(nrho)-profiles%PSIM(nrho)  +  2.0_r8*itm_pi*geometry%BGEO*geometry%RHO_BND*geometry%RHO_BND_PRIME/profiles%QSF(nrho)


    RETURN

END SUBROUTINE calculate_globals
!-------------------------------------------------------!
!-------------------------------------------------------!