convert_neutrals.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
      MODULE convert_neutrals

      CONTAINS


! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  

      SUBROUTINE convert_cpo_to_neutrals_types                                  &
                 (equilibrium, coreprof, control_integer, control_double, &
                  geometry,    profiles, control)

!-------------------------------------------------------!
!     This routine converts CPOs into NEUTRALS          !
!     derived types.                                    !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     D.Kalupin@fz-juelich.de             !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!


      USE euitm_schemas

      USE ets_plasma

      IMPLICIT NONE

! +++ CPO derived types:
      TYPE (type_equilibrium),    POINTER       :: equilibrium(:)
      TYPE (type_coreprof),       POINTER       :: coreprof(:)
      INTEGER,                    INTENT(IN)    :: control_integer(2)  !integer control parameters
      REAL (R8),                  INTENT(IN)    :: control_double(5)   !real control parameters

! +++ NEUTRALS derived types:
      TYPE (magnetic_geometry),   INTENT(INOUT) :: geometry          !contains all geometry quantities
      TYPE (plasma_profiles),     INTENT(INOUT) :: profiles          !contains profiles of plasma parameters
      TYPE (run_control),         INTENT(INOUT) :: control           !contains all parameters required by run

      INTEGER                                   :: nrho, neq

      nrho                     = SIZE(coreprof(1)%rho_tor)
      neq                      = SIZE(equilibrium(1)%profiles_1d%rho_tor)


! +++ Convert geometry:
      geometry%RHO             = coreprof(1)%rho_tor

      CALL l3deriv(equilibrium(1)%profiles_1d%volume, equilibrium(1)%profiles_1d%rho_tor, neq,  &
               geometry%VPR, geometry%RHO, nrho)


! +++ Convert profiles:
      profiles%ZION            = coreprof(1)%composition%zion 
      profiles%MION            = coreprof(1)%composition%amn

      profiles%NI              = coreprof(1)%ni%value
      profiles%NE              = coreprof(1)%ne%value 
      profiles%TI              = coreprof(1)%ti%value
      profiles%TE              = coreprof(1)%te%value

! +++ Convert control parameters:
      control%SOLVER_TYPE      = control_integer(1)
      control%SIGMA_SOURCE     = control_integer(2)
      control%TAU              = control_double(1)
      control%AMIX             = control_double(2) 
      control%AMIXTR           = control_double(3)
      control%CONV             = control_double(4) 
      control%CONVREC          = control_double(5)


      RETURN


      END SUBROUTINE convert_cpo_to_neutrals_types


! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  




! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  

      SUBROUTINE convert_neutrals_to_cpo_types                  &
                 (sources, coresource)

!-------------------------------------------------------!
!     This routine converts NEUTRALS into the CPO       !
!     derived types.                                    !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     D.Kalupin@fz-juelich.de             !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!


      USE euitm_schemas

      USE ets_plasma

      IMPLICIT NONE

! +++ CPO derived types:
      TYPE (type_coresource), POINTER  :: coresource(:)    !input CPO with sources

! +++ NEUTRALS derived types:
      TYPE (sources_and_sinks)         :: sources          !contains profiles of sources


! +++ Convert profiles of sources:
      coresource(1)%VALUES(1)%qe%exp    = sources%QE_EXP
      coresource(1)%VALUES(1)%si%exp    = sources%SI_EXP
      coresource(1)%VALUES(1)%qi%exp    = sources%QI_EXP



      RETURN


      END SUBROUTINE convert_neutrals_to_cpo_types


! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  



      END MODULE convert_neutrals