fc2k_changeocc.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +   
  SUBROUTINE changeocc                                      &
                      (coreprof_in,     coreprof_out,       &
                       coretransp_in,   coretransp_out,     &
                       coresource_in,   coresource_out,     &
                       coreimpur_in,    coreimpur_out,      &
                       coreneutrals_in, coreneutrals_out,   &
                       neoclassic_in,   neoclassic_out,     &
                       equilibrium_in,  equilibrium_out,    &
                       toroidfield_in,  toroidfield_out) 

! +++ Declaration of variables: 
    USE itm_types
    USE euitm_schemas
    USE copy_structures


! +++ CPO derived types:
    TYPE (type_coreprof),     POINTER  :: coreprof_in(:)        !input CPO  
    TYPE (type_coreprof),     POINTER  :: coreprof_out(:)       !output CPO  
    TYPE (type_coretransp),   POINTER  :: coretransp_in(:)      !input CPO  
    TYPE (type_coretransp),   POINTER  :: coretransp_out(:)     !output CPO  
    TYPE (type_coresource),   POINTER  :: coresource_in(:)      !input CPO  
    TYPE (type_coresource),   POINTER  :: coresource_out(:)     !output CPO  
    TYPE (type_coreimpur),    POINTER  :: coreimpur_in(:)       !input CPO  
    TYPE (type_coreimpur),    POINTER  :: coreimpur_out(:)      !output CPO  
    TYPE (type_coreneutrals), POINTER  :: coreneutrals_in(:)    !input CPO  
    TYPE (type_coreneutrals), POINTER  :: coreneutrals_out(:)   !output CPO  
    TYPE (type_neoclassic),   POINTER  :: neoclassic_in(:)      !input CPO
    TYPE (type_neoclassic),   POINTER  :: neoclassic_out(:)     !output CPO
    TYPE (type_equilibrium),  POINTER  :: equilibrium_in(:)     !input CPO
    TYPE (type_equilibrium),  POINTER  :: equilibrium_out(:)    !output CPO
    TYPE (type_toroidfield),  POINTER  :: toroidfield_in(:)     !input CPO
    TYPE (type_toroidfield),  POINTER  :: toroidfield_out(:)    !output CPO

! +++ Allocate output CPOs:
!    IF(ASSOCIATED(COREPROF_OUT))     DEALLOCATE(COREPROF_OUT)
!    IF(ASSOCIATED(CORETRANSP_OUT))   DEALLOCATE(CORETRANSP_OUT)
!    IF(ASSOCIATED(CORESOURCE_OUT))   DEALLOCATE(CORESOURCE_OUT)
!    IF(ASSOCIATED(COREIMPUR_OUT))    DEALLOCATE(COREIMPUR_OUT)
!    IF(ASSOCIATED(CORENEUTRALS_OUT)) DEALLOCATE(CORENEUTRALS_OUT)
!    IF(ASSOCIATED(NEOCLASSIC_OUT))   DEALLOCATE(NEOCLASSIC_OUT)
!    IF(ASSOCIATED(EQUILIBRIUM_OUT))  DEALLOCATE(EQUILIBRIUM_OUT)
!    IF(ASSOCIATED(TOROIDFIELD_OUT))  DEALLOCATE(TOROIDFIELD_OUT)

    ALLOCATE(coreprof_out(1))
    ALLOCATE(coretransp_out(1))
    ALLOCATE(coresource_out(1))
    ALLOCATE(coreimpur_out(1))
    ALLOCATE(coreneutrals_out(1))
    ALLOCATE(neoclassic_out(1))
    ALLOCATE(equilibrium_out(1))
    ALLOCATE(toroidfield_out(1))

! +++ Copy input CPOs to output CPOs:

    CALL copy_cpo(coreprof_in(1),     coreprof_out(1))
    CALL copy_cpo(coretransp_in(1),   coretransp_out(1))
    CALL copy_cpo(coresource_in(1),   coresource_out(1))
    CALL copy_cpo(coreimpur_in(1),    coreimpur_out(1))
    CALL copy_cpo(coreneutrals_in(1), coreneutrals_out(1))
    CALL copy_cpo(neoclassic_in(1),   neoclassic_out(1))
    CALL copy_cpo(equilibrium_in(1),  equilibrium_out(1)) 
    CALL copy_cpo(toroidfield_in(1),  toroidfield_out(1)) 

    RETURN 

    END SUBROUTINE changeocc 
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +