coreprof2coredelta.f90

Go to the documentation of this file.

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

  !-------------------------------------------------------!
  !     This routine converts two coreprof CPOs           !
  !     to a single coredelta CPO by taking               !
  !     difference between them                           !
  !-------------------------------------------------------!
  !     Delta:       ---                                  !
  !     Developers:   D.Kalupin                           !
  !     Kontacts:     Denis.Kalupin@efda.org              !
  !                                                       !
  !     Comments:     created for  ETS workflow           !
  !                                                       !
  !-------------------------------------------------------!

    SUBROUTINE coreprof2coredelta                         &
              (coreprof_pivot, coreprof_new, coredelta)

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

      
      IMPLICIT NONE



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

      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 pivot CPO to the grid of new CPO  
      ALLOCATE                    (coreprof_interp(1))
      CALL copy_cpo(coreprof_new(1),   coreprof_interp(1))
      CALL interpolate_prof(coreprof_pivot(1), coreprof_interp(1))
      nrho                  = SIZE(coreprof_new(1)%rho_tor)


! +++ Allocate output CPO and internal derived types:
      CALL get_comp_dimensions(coreprof_new(1)%COMPOSITIONS, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp)
      CALL allocate_coredelta_cpo(nslice, nrho, nnucl, nion,  nimp,  nzimp, nneut, ntype, ncomp,  coredelta)
      CALL copy_cpo(coreprof_new(1)%COMPOSITIONS,   coredelta(1)%COMPOSITIONS)


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


! +++ ADD IDENTIFIER TO OUTPUT CPO VALUES(1):
      ALLOCATE                                (coredelta(1)%VALUES(1)%deltaid%id(1))
      ALLOCATE                                (coredelta(1)%VALUES(1)%deltaid%description(1))     
      coredelta(1)%VALUES(1)%deltaid%id          = 'unspecified'
      coredelta(1)%VALUES(1)%deltaid%flag        = 0
      coredelta(1)%VALUES(1)%deltaid%description = 'Unspecified coredelta'


10    IF(ALLOCATED(nzimp)) DEALLOCATE         (nzimp)
      IF(ALLOCATED(ncomp)) DEALLOCATE         (ncomp)
      IF(ALLOCATED(ntype)) DEALLOCATE         (ntype)

      CALL deallocate_cpo(coreprof_interp)
 

      RETURN


    END SUBROUTINE coreprof2coredelta
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +