database_impurity.f90

Go to the documentation of this file.

subroutine database_impurity (coreprof_in, coresource_db, coreimpur_db, coresource_new, coreimpur_new)

    USE allocate_deallocate
    USE deallocate_structures

    USE itm_constants
    USE euitm_routines
    USE euitm_schemas
    USE interpolate_cpo
    USE copy_structures

    IMPLICIT NONE


! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! +++ CPOs:
   

    TYPE (type_coreprof),    POINTER :: coreprof_in(:)  
    TYPE (type_coresource),    POINTER :: coresource_db(:),coresource_new(:)    
    TYPE (type_coreimpur),    POINTER :: coreimpur_db(:),  coreimpur_new(:)   

    INTEGER,              PARAMETER  :: nslice = 1             !number of CPO ocurancies in the work flow
    INTEGER                          :: nrho, nrho1
    INTEGER                          :: nnucl
    INTEGER                          :: nion, iion
    INTEGER                          :: nimp       
    INTEGER,             ALLOCATABLE :: nzimp(:)
    INTEGER                          :: nneut
    INTEGER,             ALLOCATABLE :: ncomp(:)
    INTEGER,             ALLOCATABLE :: ntype(:)
    REAL(R8),allocatable,dimension(:) :: rho,rho1,fun
    integer    :: nvalues, ivalues,impvalue,iimp

!get dimensions, composition from coreprof
write(*,*) 'time used',coreprof_in(1)%time
nrho                  =  SIZE(coreprof_in(1)%rho_tor,      dim=1)
write(*,*) 'nrho=', nrho
    CALL get_comp_dimensions(coreprof_in(1)%COMPOSITIONS, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp)
!allocate coreimpur and coresource
CALL allocate_coresource_cpo(nslice, nrho, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp, coresource_new)
CALL allocate_coreimpur_cpo(nslice, nrho, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp, coreimpur_new)

!allocate temperature as it is not allocated by default
  do iimp=1,nimp
     if (.not.associated(coreimpur_new(1)%impurity(iimp)%tz)) &
         allocate(coreimpur_new(1)%impurity(iimp)%tz(nrho,nzimp(iimp)))
  enddo


!allocate source_id as it is not allocated by default
if (.not.associated(coresource_new(1)%VALUES(1)%sourceid%id)) &
    ALLOCATE            (coresource_new(1)%VALUES(1)%sourceid%id(1))
if (.not.associated(coresource_new(1)%VALUES(1)%sourceid%description)) &
    ALLOCATE            (coresource_new(1)%VALUES(1)%sourceid%description(1))

call copy_cpo(coreprof_in(1)%compositions, coresource_new(1)%compositions)
call copy_cpo(coreprof_in(1)%compositions, coreimpur_new(1)%compositions)

coreimpur_new(1)%rho_tor=coreprof_in(1)%rho_tor
coreimpur_new(1)%rho_tor_norm=coreprof_in(1)%rho_tor_norm
coresource_new(1)%values(1)%rho_tor=coreprof_in(1)%rho_tor
coresource_new(1)%values(1)%rho_tor_norm=coreprof_in(1)%rho_tor_norm
coresource_new(1)%VALUES(1)%sourceid%id          = 'impurity'
coresource_new(1)%VALUES(1)%sourceid%flag        = 29
coresource_new(1)%VALUES(1)%sourceid%description = 'Impurity source'

allocate (fun(nrho),rho(nrho))
rho=coreprof_in(1)%rho_tor

!get data from coresource (qe,se)
IF(ASSOCIATED(coresource_db)) then
   nvalues=size(coresource_db(1)%values)
   !write(*,*) 'nvalues',nvalues
   write(*,*) 'time from database',coresource_db(1)%time
   if (nvalues.ge.0) then
      impvalue=0
      do ivalues=1,nvalues
         if (coresource_db(1)%values(ivalues)%sourceid%flag.eq.29) impvalue=ivalues
      enddo
   else
      write(*,*) 'no sources are given on input'
   end if
   write(*,*) 'impvalue',impvalue
   if (associated(coresource_db(1)%values(impvalue)%rho_tor)) then
     nrho1=size(coresource_db(1)%values(impvalue)%rho_tor) 
   else 
     write(*,*) 'rho_tor from database (coresource) is not found'
     nrho1=0
   end if
   write(*,*) 'nrho1',nrho1
   allocate(rho1(nrho1))
   rho1(:)=coresource_db(1)%values(impvalue)%rho_tor(:)
!Se
   if (associated(coresource_db(1)%values(impvalue)%se%exp)) then
       fun(:)                                            =  0.0_r8
       IF(ASSOCIATED(coresource_db(1)%VALUES(impvalue)%Se%exp)) THEN
          !write(*,*) 'se', CORESOURCE_db(1)%VALUES(impvalue)%Se%exp
          if (all(coresource_db(1)%VALUES(impvalue)%Se%exp==0)) then
              coresource_new(1)%VALUES(1)%Se%exp(:)    = 0.0
          else
              CALL l3interp(coresource_db(1)%VALUES(impvalue)%Se%exp,     rho1,   nrho1,  &
                        fun,                                   rho,   nrho)
              coresource_new(1)%VALUES(1)%Se%exp(:)    = fun(:)
          end if
       END IF
   end if
!Qe
   if (associated(coresource_db(1)%values(impvalue)%qe%exp)) then
       fun(:)                                            =  0.0_r8
       IF(ASSOCIATED(coresource_db(1)%VALUES(impvalue)%qe%exp)) THEN  
          !write(*,*) 'qe',CORESOURCE_db(1)%VALUES(impvalue)%qe%exp
          CALL l3interp(coresource_db(1)%VALUES(impvalue)%qe%exp,     rho1,   nrho1,  &
                        fun,                                   rho,   nrho)
          coresource_new(1)%VALUES(1)%Qe%exp(:)    = fun(:)
       END IF
   end if
end if
write(*,*) 'after coresource'
!get data from coreimpur (nz)
call interpolate_impur(coreimpur_db(1),coreimpur_new(1))
write(*,*) 'after impurity'
CALL deallocate_cpo(coresource_db)
call deallocate_cpo(coreimpur_db)
deallocate(fun,rho,rho1)


end subroutine database_impurity