spitzer.f90

Go to the documentation of this file.

MODULE spitzer

CONTAINS
!-------------------------------------------------------!
!-------------------------------------------------------!

  SUBROUTINE spitzer_resistivity    (COREPROF_IN, CORETRANSP_OUT) 

!-------------------------------------------------------!
!     This routine is used by the ETS workflow          !
!     to calculate Spitzer resistivity                  !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     Denis.Kalupin@efda.org              !
!                                                       !
!     Comments:     output CORETRANSP CPO is            !
!                   allocated inside the module         !
!                                                       !
!-------------------------------------------------------!

      USE itm_types
      USE itm_constants
      USE euitm_schemas
      USE allocate_deallocate
      USE plasma_collisionality
      

      IMPLICIT NONE

      INTEGER                           :: ifail

! +++ CPO types:
      TYPE (type_coreprof),    POINTER  :: coreprof_in(:)        !input CPO with profiles
      TYPE (type_coretransp),  POINTER  :: coretransp_out(:)     !output CPO with transport (SIGMA)

! +++ Dimensions:
      INTEGER,                PARAMETER :: nslice = 1            !number of CPO ocurancies in the work flow
      INTEGER                           :: nrho, irho            !number of radial points     (input, determined from COREPROF CPO)
      INTEGER                           :: nnucl                 !number of nuclei species
      INTEGER                           :: nion                  !number of ion species
      INTEGER                           :: nimp                  !number of impurity species
      INTEGER,              ALLOCATABLE :: nzimp(:)              !number of ionization states for each impurity
      INTEGER                           :: nneut                 !number of neutrals species
      INTEGER,              ALLOCATABLE :: ncomp(:)              !number of components for each neutral
      INTEGER,              ALLOCATABLE :: ntype(:)              !number of types for each neutral

! +++ Internal variables:
      REAL (R8)                         :: clog                  !Coulomb logarithm      
      REAL (R8),            ALLOCATABLE :: te(:)
      REAL (R8),            ALLOCATABLE :: ne(:)
      REAL (R8),            ALLOCATABLE :: tau_e(:)
      REAL (R8),            ALLOCATABLE :: sigma(:)

! +++ Constants:
      REAL (R8),              PARAMETER :: me  = itm_me*1e3_r8   !electron mass              [g]
      REAL (R8),              PARAMETER :: e   = itm_qe*3e9_r8   !elementary charge,         [esu]
      REAL (R8),              PARAMETER :: cn  = 1.e-6_r8        !density convergence        from [m^-3] to [cm^-3]
      REAL (R8),              PARAMETER :: cs  = 9.e9_r8         !conductivity convergence   from [(Ohm*m)^-1] to [s^-1]


! +++ Set dimensions:
      nrho                   = SIZE (coreprof_in(1)%rho_tor)
      CALL get_comp_dimensions(coreprof_in(1)%COMPOSITIONS, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp)

! +++ Allocate output CPO:
      CALL allocate_coretransp_cpo(nslice,  nrho,  nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp, coretransp_out)        
      coretransp_out(1)%VALUES(1)%rho_tor = coreprof_in(1)%rho_tor

! +++ Set up profiles:
      ALLOCATE (te(nrho))
      ALLOCATE (ne(nrho))
      ALLOCATE (tau_e(nrho))
      ALLOCATE (sigma(nrho))
      rho_loop1: DO irho =1,nrho
         te(irho)           = coreprof_in(1)%TE%value(irho)
         ne(irho)           = coreprof_in(1)%NE%value(irho)


! +++ Electron collisions:
!       determination of Coulomb logarithm:
         IF(te(irho).GE.10) clog = 24.e0_r8 - 1.15e0_r8*log10(ne(irho)*cn) + 2.30e0_r8*log10(te(irho))
         IF(te(irho).LT.10) clog = 23.e0_r8 - 1.15e0_r8*log10(ne(irho)*cn) + 3.45e0_r8*log10(te(irho))


!       electron collision time:
         tau_e(irho)          = (sqrt(2.d0*me)*(te(irho))**1.5) / 1.8d-19 / (ne(irho)*cn) / clog


!       Plasma electrical conductivity:
         sigma(irho)          = 1.96e0_r8 * e**2 *ne(irho)*cn * tau_e(irho) /me /cs

      END DO rho_loop1

      coretransp_out(1)%VALUES(1)%sigma = sigma

! +++ ADD IDENTIFIER TO OUTPUT CPO VALUES(1):
      ALLOCATE            (coretransp_out(1)%VALUES(1)%transportid%id(1))
      ALLOCATE            (coretransp_out(1)%VALUES(1)%transportid%description(1))
      coretransp_out(1)%VALUES(1)%transportid%id          = 'spitzer'
      coretransp_out(1)%VALUES(1)%transportid%flag        = 14
      coretransp_out(1)%VALUES(1)%transportid%description = 'Spitzer Resistivity'


      IF (ALLOCATED(te))    DEALLOCATE (te)
      IF (ALLOCATED(ne))    DEALLOCATE (ne)
      IF (ALLOCATED(tau_e)) DEALLOCATE (tau_e)
      IF (ALLOCATED(sigma)) DEALLOCATE (sigma)

      RETURN

      END SUBROUTINE spitzer_resistivity    
!-------------------------------------------------------!
!-------------------------------------------------------!


END MODULE spitzer