ignore_transport.f90

Go to the documentation of this file.

!-------------------------------------------------------!
!-------------------------------------------------------!

  SUBROUTINE ignore_transport    (COREPROF, CORETRANSP) 

!-------------------------------------------------------!
!     This routine is used by the ETS workflow          !
!     in case when transport should be ignored          !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     Denis.Kalupin@efda.org              !
!                                                       !
!     Comments:     output CORETRANSP CPO is            !
!                   allocated inside the module         !
!                                                       !
!-------------------------------------------------------!


      USE  euitm_schemas
      USE  euitm_routines
      USE  allocate_deallocate
      USE  itm_types
      USE  copy_structures
      USE  deallocate_structures


      IMPLICIT NONE


      INTEGER                           :: ifail
 

! +++ CPO derived types:
      TYPE (type_coreprof),    POINTER  :: coreprof(:)         !input CPO with internal ETS parameters profiles from previous time
      TYPE (type_coretransp),  POINTER  :: coretransp(:)       !output CPO with transport


! +++ Dimensions:
      INTEGER,               PARAMETER  :: nslice = 1          !number of CPO ocurancies in the work flow
      INTEGER                           :: nrho                !number of radial points     (input, determined from COREPROF CPO)
      INTEGER                           :: nion                !number of ion species       (input, determined from COREPROF CPO)
      INTEGER                           :: nimp                !number of impurities                (input)
      INTEGER                           :: nnucl               !number of nuclei 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



! +++ Set dimensions:
      nrho                     = SIZE (coreprof(1)%rho_tor, dim=1)

      CALL get_comp_dimensions(coreprof(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)
      call deallocate_cpo(coretransp(1)%COMPOSITIONS)
      CALL copy_cpo(coreprof(1)%COMPOSITIONS, coretransp(1)%COMPOSITIONS)


! +++ Save output in CPO:
      coretransp(1)%time                     =  coreprof(1)%time     !time     [s]
      coretransp(1)%values(1)%rho_tor        =  coreprof(1)%rho_tor  !rho      [m]
      coretransp(1)%values(1)%rho_tor_norm   =  coreprof(1)%rho_tor/coreprof(1)%rho_tor(nrho)  !rho_norm [m]



!+++ ADD IDENTIFIER TO OUTPUT CPO VALUES(1):
     ALLOCATE            (coretransp(1)%VALUES(1)%transportid%id(1))
     ALLOCATE            (coretransp(1)%VALUES(1)%transportid%description(1))
     coretransp(1)%VALUES(1)%transportid%id          = 'not_provided'
     coretransp(1)%VALUES(1)%transportid%flag        = 13
     coretransp(1)%VALUES(1)%transportid%description = 'No data provided'

      RETURN


      END SUBROUTINE ignore_transport  

!-------------------------------------------------------!
!-------------------------------------------------------!