coredelta2coreprof.f90

Go to the documentation of this file.

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

  !-------------------------------------------------------!
  !     This routine adds values from coredelta CPO       !
  !     to a coreprof CPO                                 !
  !-------------------------------------------------------!
  !     Delta:       ---                                  !
  !     Developers:   D.Kalupin                           !
  !     Kontacts:     Denis.Kalupin@efda.org              !
  !                                                       !
  !     Comments:     created for  ETS workflow           !
  !                                                       !
  !-------------------------------------------------------!

    SUBROUTINE coredelta2coreprof                         &
              (coreprof_pivot, coredelta, coreprof_new)

      USE itm_constants
      USE euitm_schemas
      USE copy_structures
      USE interpolate_cpo
      USE deallocate_structures
      USE allocate_deallocate

      
      IMPLICIT NONE



! +++ CPOs  
      TYPE (type_coreprof),   POINTER :: coreprof_pivot(:)  
      TYPE (type_coreprof),   POINTER :: coreprof_new(:)  
      TYPE (type_coredelta),  POINTER :: coredelta_interp(:)  
      TYPE (type_coredelta),  POINTER :: coredelta(:)    
      TYPE (type_coredelta),  POINTER :: coredelta_v1(:)    

      INTEGER                         :: nrho
      INTEGER,              PARAMETER :: nslice = 1             
      INTEGER                         :: nrho1, nrho2
      INTEGER                         :: nnucl
      INTEGER                         :: nion, iion
      INTEGER                         :: nimp       
      INTEGER,            ALLOCATABLE :: nzimp(:)
      INTEGER                         :: nneut
      INTEGER,            ALLOCATABLE :: ncomp(:)
      INTEGER,            ALLOCATABLE :: ntype(:)
 

! +++ Interpolate coredelta CPO to the grid of pivot coreprof CPO  

      nrho                       = SIZE(coreprof_pivot(1)%rho_tor)
      CALL get_comp_dimensions(coreprof_pivot(1)%COMPOSITIONS, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp)
      CALL allocate_coredelta_cpo(nslice, nrho, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp,  coredelta_interp)
      CALL copy_cpo(coreprof_pivot(1)%COMPOSITIONS,   coredelta_interp(1)%COMPOSITIONS)
      coredelta_interp(1)%VALUES(1)%rho_tor = coreprof_pivot(1)%rho_tor



      ALLOCATE                    (coredelta_v1(1))
      ALLOCATE                    (coredelta_v1(1)%VALUES(1))
      CALL copy_cpo(coredelta(1)%COMPOSITIONS,   coredelta_v1(1)%COMPOSITIONS)
      CALL copy_cpo(coredelta(1)%VALUES(1),      coredelta_v1(1)%VALUES(1))


      CALL interpolate_delta(coredelta_v1(1),     coredelta_interp(1))



! +++ Allocate output CPO and internal derived types:
      ALLOCATE                    (coreprof_new(1))
      CALL copy_cpo(coreprof_pivot(1),   coreprof_new(1))



! +++ Fill output CPO with data:
      coreprof_new(1)%psi%value  = coreprof_pivot(1)%psi%value  + coredelta_interp(1)%VALUES(1)%delta_psi
      coreprof_new(1)%te%value   = coreprof_pivot(1)%te%value   + coredelta_interp(1)%VALUES(1)%delta_te
      coreprof_new(1)%ne%value   = coreprof_pivot(1)%ne%value   + coredelta_interp(1)%VALUES(1)%delta_ne
      coreprof_new(1)%ti%value   = coreprof_pivot(1)%ti%value   + coredelta_interp(1)%VALUES(1)%delta_ti
      coreprof_new(1)%ni%value   = coreprof_pivot(1)%ni%value   + coredelta_interp(1)%VALUES(1)%delta_ni
      coreprof_new(1)%vtor%value = coreprof_pivot(1)%vtor%value + coredelta_interp(1)%VALUES(1)%delta_vtor


      CALL deallocate_cpo(coredelta_interp)
      CALL deallocate_cpo(coredelta_v1)
 

      RETURN


    END SUBROUTINE coredelta2coreprof
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +