fc2k_plasma_composition2.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +   
  SUBROUTINE fc2k_plasma_composition2 (COREPROF_DB, COREPROF_OUT,             &
                                 nion,      nimp,      nneut,           &
                                 amn_ion,   zn_ion,    z_ion,           &
                                 amn_imp,   zn_imp,    maxz_imp,        &
                                 ncomp_fix,  ntype_fix,                   &
                                 ncold,     nthermal,  nfast,   nnbi) 


! +++ Declaration of variables: 
    use itm_types
    USE euitm_schemas
    USE allocate_deallocate
    USE copy_structures
    USE deallocate_structures
    USE plasma_composition



    IMPLICIT NONE

! +++ CPO derived types:
    TYPE (type_coreprof),     POINTER   :: coreprof_db(:)         !input  CPO slice 
    TYPE (type_coreprof),     POINTER   :: coreprof_out(:)        !output CPO slice 
    TYPE (type_compositionc), POINTER   :: compositionc(:)

! +++ Input:
    INTEGER,                  PARAMETER  :: nocur = 1             !number of CPO ocurancies in the work flow
    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
 
    INTEGER                              :: ncold, nthermal, nfast, nnbi     

    REAL (R8)                            :: amn_ion(30),  zn_ion(30),  z_ion(30)
    REAL (R8)                            :: amn_imp(100), zn_imp(100), maxz_imp(100)

    integer                              :: ncomp_fix,ntype_fix
    INTEGER                              :: ncomp_in(130)
    INTEGER                              :: ntype_in(130)
    integer                              :: i1
! +++ Local:
    REAL (R8),              ALLOCATABLE  :: amn(:),    zn(:),    zion(:)
    REAL (R8),              ALLOCATABLE  :: amnimp(:), znimp(:), maxzimp(:)


    ALLOCATE (amn(30),     zn(30),     zion(30))
    ALLOCATE (amnimp(130), znimp(130), maxzimp(130))

    amn      = amn_ion
    zn       = zn_ion
    zion     = z_ion
    amnimp   = amn_imp
    znimp    = zn_imp
    maxzimp  = maxz_imp
    do i1=1,130
       ncomp_in(i1)=ncomp_fix
       ntype_in(i1)=ntype_fix
    enddo



    CALL set_plasma_composition(coreprof_out,                          &
                                 nion,      nimp,      nneut,           &
                                 amn,       zn,        zion,            &
                                 amnimp,    znimp,     maxzimp,         &
                                 ncomp_in,  ntype_in,                   &
                                 ncold,     nthermal,  nfast,   nnbi) 



    RETURN 


    END SUBROUTINE fc2k_plasma_composition2
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +