joiner.f90

Go to the documentation of this file.

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

MODULE joiner 

  IMPLICIT NONE 

CONTAINS



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

!-------------------------------------------------------!
     subroutine join_neoclassic(NEOCLASSICin1,          &
                                neoclassicin2,          &
                                neoclassicout)
!-------------------------------------------------------!
!     This routine 
!     
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   
!     Kontacts:     
!                                                       !
!     Comments:     
!                   
!                                                       !
!-------------------------------------------------------!


     USE euitm_schemas
     USE itm_constants
     USE copy_structures
     USE deallocate_structures
     USE allocate_deallocate

     IMPLICIT NONE

     TYPE (type_neoclassic),  pointer :: neoclassicin1(:)
     TYPE (type_neoclassic),  pointer :: neoclassicin2(:)
     TYPE (type_neoclassic),  POINTER :: neoclassicout(:)
     TYPE (type_compositions_type)    :: compositions 

     INTEGER,              PARAMETER  :: nslice = 1             !number of CPO ocurancies in the work flow
     INTEGER                          :: nrho1, nrho2 
     INTEGER                          :: nnucl
     INTEGER                          :: nion, iion
     INTEGER                          :: nimp, iimp, izimp      
     INTEGER,             ALLOCATABLE :: nzimp(:)
     INTEGER                          :: nneut
     INTEGER,             ALLOCATABLE :: ncomp(:)
     INTEGER,             ALLOCATABLE :: ntype(:)



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


     nrho1 = size(neoclassicin1(1)%rho_tor)
     nrho2 = size(neoclassicin2(1)%rho_tor)
!     CALL GET_COMP_DIMENSIONS  (NEOCLASSICin(1)%COMPOSITIONS, NNUCL, NION,  NIMP,  NZIMP, NNEUT, NTYPE, NCOMP)
     

     allocate(neoclassicout(1))        
     call copy_cpo(neoclassicin1(1), neoclassicout(1))   


! ++ sigma 
     if (associated(neoclassicin2(1)%sigma)) then
        if (associated(neoclassicout(1)%sigma)) deallocate(neoclassicout(1)%sigma)
        allocate(neoclassicout(1)%sigma(nrho1))
        call l3interp(neoclassicin2(1)%sigma, neoclassicin2(1)%rho_tor, nrho2,  &
                      neoclassicout(1)%sigma, neoclassicout(1)%rho_tor, nrho1)
     end if 

! ++ bootstrap current 
     if (associated(neoclassicin2(1)%jboot)) then
        if (associated(neoclassicout(1)%jboot)) deallocate(neoclassicout(1)%jboot)
        allocate(neoclassicout(1)%jboot(nrho1))
        call l3interp(neoclassicin2(1)%jboot, neoclassicin2(1)%rho_tor, nrho2,  &
                      neoclassicout(1)%jboot, neoclassicout(1)%rho_tor, nrho1)
     end if 


!+++ 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          = 'neoclassical'
    !CORETRANSP(1)%VALUES(1)%transportid%flag        = 2
    !CORETRANSP(1)%VALUES(1)%transportid%description = 'Neoclassical'
     




    RETURN     


    END SUBROUTINE join_neoclassic
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  


END MODULE joiner