fc2k_changepsi.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +   
  SUBROUTINE fc2k_changepsi   (EQUILIBRIUM_IN,                      &
                               coreprof_in,     coreprof_out) 

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


! +++ CPO derived types:
    TYPE (type_equilibrium),  POINTER  :: equilibrium_in(:)      !input CPO

    TYPE (type_coreprof),     POINTER  :: coreprof_in(:)         !input CPO  
    TYPE (type_coreprof),     POINTER  :: coreprof_out(:)        !output CPO  
    REAL(R8)                           :: rho_tor_rescale


! +++ Allocate output CPOs:
    ALLOCATE(coreprof_out(1))


! +++ Copy input CPOs to output CPOs:
    CALL copy_cpo(coreprof_in(1),     coreprof_out(1))


! +++ Copy input CPOs to output CPOs:
    rho_tor_rescale             = equilibrium_in(1)%profiles_1d%rho_tor(SIZE(equilibrium_in(1)%profiles_1d%rho_tor)) / &
                                  coreprof_out(1)%rho_tor(SIZE(coreprof_out(1)%rho_tor))
    WRITE(*,*) 'RHO_TOR COREPROF rescale factor = ', rho_tor_rescale
    coreprof_out(1)%rho_tor     = coreprof_out(1)%rho_tor * rho_tor_rescale

! +++ Modify psi:
    CALL l3interp(equilibrium_in(1)%profiles_1d%psi, equilibrium_in(1)%profiles_1d%rho_tor, SIZE(equilibrium_in(1)%profiles_1d%rho_tor),&
                        coreprof_out(1)%psi%value,         coreprof_out(1)%rho_tor,               SIZE(coreprof_out(1)%rho_tor)) 

    RETURN 

    END SUBROUTINE fc2k_changepsi
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +