pedestal.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
MODULE pedestal

CONTAINS

   
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
  SUBROUTINE main_plasma_pedestal (FUN, RHO, RHO_PED, PSI, NRHO) 

!-------------------------------------------------------!
!     This routine provides artificial pedestal         !
!     in profiles of plasma parameters which are not    !
!     solved up to the separatrix.                      !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     Denis.Kalupin@efda.org              !
!                                                       !
!     Ref.:         `Snyder, PPCF 46 (2004), A131`      !
!                                                       !
!-------------------------------------------------------!


    USE itm_types
    USE ets_plasma

    IMPLICIT NONE

! +++ Internal ETS derived types:
    TYPE (magnetic_geometry)          :: geometry       !contains all geometry quantities
    TYPE (plasma_profiles)            :: profiles       !contains profiles of plasma parameters


! +++ Dimensions:
    INTEGER                           :: nrho,irho      !number of radial points     
    INTEGER                           :: nped


! +++ Parameters:
    REAL (R8),          INTENT(INOUT) :: fun(nrho)
    REAL (R8),          INTENT(IN)    :: rho(nrho)
    REAL (R8),          INTENT(IN)    :: psi(nrho)
    REAL (R8),          INTENT(IN)    :: rho_ped
    REAL (R8)                         :: psi_mid
    REAL (R8)                         :: delta
!-------------------------------------------------------!

! +++ Determine the index of pedestal top
    DO irho = 1, nrho
       IF (rho(irho).GT.rho_ped) EXIT
       nped = irho
    END DO

    psi_mid =    (psi(nped) + psi(nrho)) / 2.0_r8
    delta   = abs(psi(nped) - psi(nrho))

! +++ Determine the index of pedestal top
    DO irho = nped+1, nrho
       fun(irho) = fun(nped) + (fun(nped)-fun(nrho)) *  &
                   (tanh(2.0_r8*(1.0_r8-psi_mid)/delta)-tanh(2.0_r8*(psi(irho)-psi_mid)/delta))
    END DO

    RETURN



  END SUBROUTINE main_plasma_pedestal
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  





END MODULE pedestal