synchrotron.f90

Go to the documentation of this file.

MODULE synchrotron

CONTAINS


! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + SYNCHROTRON + + + + + + + + + +   

  SUBROUTINE synchrotron_radiation(COREPROF, CORESOURCE)

!-------------------------------------------------------!
!     This routine calculate sybchrotron radiation      !
!                                                       !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   I.M.Ivanova-Stanik                  !
!     Kontacts:     irena.ivanova-stanik@ifpilm.pl      !
!                                                       !
!     Comments:     might change after the ITM          !
!                   data stucture is finalized          !
!                                                       !
!-------------------------------------------------------!


    USE itm_types
    USE itm_constants
    USE euitm_schemas
    USE copy_structures
    USE allocate_deallocate
      
      
    IMPLICIT NONE
 
! +++ CPO derived types:
    TYPE (type_coreprof),    POINTER :: coreprof(:)          !input CPO main plasma
    TYPE (type_coresource),  POINTER :: coresource(:)        !output CPO with syncrhotron radiation

      

! +++ Dimensions:
    INTEGER                          :: nrho                 !number of radial points     (input, determined from COREPROF CPO)
    INTEGER                          :: nnucl                !number of ion species       (input, determined from COREPROF CPO)
    INTEGER                          :: nion                 !number of ion species       (input, determined from COREPROF CPO)
    INTEGER                          :: nimp                 !number of impurity species  (input, determined from COREIMPUR CPO)
    INTEGER                          :: nneut                !number of neutrals species                (input)
    INTEGER,             ALLOCATABLE :: ntype(:)             !number of impurity ionization states (input)     
    INTEGER,             ALLOCATABLE :: ncomp(:)             !max_number of distinct atoms enter the composition-"1" wich is neutral
    INTEGER,             ALLOCATABLE :: nzimp(:)
  
    INTEGER,               PARAMETER :: nocur = 1            !number of CPO ocurancies in the work flow
  
     

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

! +++ Allocate output CPOs:
    CALL allocate_coresource_cpo(nocur, nrho, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp, coresource)
    CALL copy_cpo(coreprof(1)%COMPOSITIONS, coresource(1)%COMPOSITIONS)
     

    coresource(1)%datainfo%cocos           = 13
    coresource(1)%time                     = coreprof(1)%time
    coresource(1)%VALUES(1)%rho_tor        = coreprof(1)%rho_tor
    coresource(1)%VALUES(1)%rho_tor_norm   = coreprof(1)%rho_tor/coreprof(1)%rho_tor(nrho)


! +++ Syncrhotron radiation    
    coresource(1)%VALUES(1)%qe%exp = 0.01_r8*6.2e-20_r8*coreprof(1)%toroid_field%B0**2.0*coreprof(1)%ne%value*coreprof(1)%te%value


! +++ ADD IDENTIFIER TO OUTPUT CPO VALUES(1):
    ALLOCATE            (coresource%VALUES(1)%sourceid%id(1))
    ALLOCATE            (coresource(1)%VALUES(1)%sourceid%description(1))
    coresource(1)%VALUES(1)%sourceid%id          = 'syncrotronradiation'
    coresource(1)%VALUES(1)%sourceid%flag        = 17
    coresource(1)%VALUES(1)%sourceid%description = 'Source from syncrotron radiation'


    IF(ALLOCATED(ntype)) DEALLOCATE (ntype)
    IF(ALLOCATED(ncomp)) DEALLOCATE (ncomp)
    IF(ALLOCATED(nzimp)) DEALLOCATE (nzimp)



    RETURN


    END SUBROUTINE synchrotron_radiation
! + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 
! + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +




END MODULE synchrotron
! + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 
! + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +