convergence_check.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
MODULE convergence_check

CONTAINS



! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +   
  SUBROUTINE check_convergence (COREPROF_ITER, COREPROF_NEW, CONTROL_DOUBLE) 

!-------------------------------------------------------!
!     This routinechecks the convergence of plasma      !
!     profiles.                                         !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     D.Kalupin@fz-juelich.de             !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!


! +++ Declaration of variables: 
    USE euitm_schemas
    USE copy_structures
    USE ets_plasma
    USE copy_cpo_ets


    IMPLICIT NONE

! +++ CPO derived types:
    TYPE (type_coreprof), POINTER  :: coreprof_iter(:) !input/output CPO with internal ETS parameters profiles 
    TYPE (type_coreprof), POINTER  :: coreprof_new(:)  !input/output CPO with internal ETS parameters profiles 


! +++ Internal ETS derived types:
    REAL (R8)                      :: control_double(5)!real control parameters
    REAL (R8)                      :: conv

! +++ local

    REAL (R8)                      :: err_psi, err_ni, err_te, err_ti, err_vtor
    REAL (R8), PARAMETER           :: psi_0 = 1.0e-3_r8
    REAL (R8), PARAMETER           :: ni_0 = 1.0e10_r8
    REAL (R8), PARAMETER           :: te_0 = 1.0e-3_r8
    REAL (R8), PARAMETER           :: ti_0 = 1.0e-3_r8
    REAL (R8), PARAMETER           :: vtor_0 = 1.0e-3_r8
 
! +++ Check convergence:

    err_psi  = maxval(abs(1.0_r8 - (abs(coreprof_iter(1)%psi%value)+psi_0)/(abs(coreprof_new(1)%psi%value)+psi_0)))
    err_ni   = maxval(abs(1.0_r8 - (abs(coreprof_iter(1)%ni%value)+ni_0)/(abs(coreprof_new(1)%ni%value)+ni_0)))
    err_te   = maxval(abs(1.0_r8 - (abs(coreprof_iter(1)%te%value)+te_0)/(abs(coreprof_new(1)%te%value)+te_0)))
    err_ti   = maxval(abs(1.0_r8 - (abs(coreprof_iter(1)%ti%value)+ti_0)/(abs(coreprof_new(1)%ti%value)+ti_0)))
    err_vtor = maxval(abs(1.0_r8 - (abs(coreprof_iter(1)%vtor%value)+vtor_0)/(abs(coreprof_new(1)%vtor%value)+vtor_0)))

    !write(*,*) ' psi err ', err_psi
    !write(*,*) '  ni err ', err_ni
    !write(*,*) '  te err ', err_te
    !write(*,*) '  ti err ', err_ti
    !write(*,*) 'vtor err ', err_vtor

    conv = max(err_psi, err_ni, err_te, err_ti, err_vtor)

    control_double(4)   = conv

    RETURN 



  END SUBROUTINE check_convergence
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  



! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +   
  SUBROUTINE check_convergence_neutrals (CORENEUTRALS_ITER, CORENEUTRALS_NEW, CONV) 

!-------------------------------------------------------!
!     This routinechecks the convergence of plasma      !
!     profiles.                                         !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D. Coster                           !
!     Kontacts:     David.Coster@ipp.mpg.de             !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!


! +++ Declaration of variables: 
    USE euitm_schemas
    USE copy_structures
    USE ets_plasma
    USE copy_cpo_ets


    IMPLICIT NONE

    INTEGER                        :: ineut, itype

! +++ CPO derived types:
    TYPE (type_coreneutrals), POINTER  :: coreneutrals_iter(:) !input/output CPO with internal ETS parameters profiles 
    TYPE (type_coreneutrals), POINTER  :: coreneutrals_new(:)  !input/output CPO with internal ETS parameters profiles 


! +++ Internal ETS derived types:
    REAL (R8)                      :: conv

! +++ local

    REAL (R8)                      :: err_neut
    REAL (R8), PARAMETER           :: ni_0 = 1.0e10_r8
 
! +++ Check convergence:

    err_neut = 0.0

    do ineut = 1, size(coreneutrals_iter(1)%PROFILES)
       do itype = 1, size(coreneutrals_iter(1)%PROFILES(ineut)%NEUTRALTYPE)
          err_neut = max(err_neut, &
               maxval(abs(1.0_r8 - (abs(coreneutrals_iter(1)%PROFILES(ineut)%NEUTRALTYPE(itype)%n0%value)+ni_0)/  &
                                   (abs(coreneutrals_new(1)%PROFILES(ineut)%NEUTRALTYPE(itype)%n0%value)+ni_0))))
       enddo
    enddo

    !write(*,*) '  n0 err ', err_neut

    conv = err_neut

    RETURN 



  END SUBROUTINE check_convergence_neutrals
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  



! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +   
  SUBROUTINE check_convergence_impurities (COREIMPUR_ITER, COREIMPUR_NEW, CONV) 

!-------------------------------------------------------!
!     This routinechecks the convergence of plasma      !
!     profiles.                                         !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D. Coster                           !
!     Kontacts:     David.Coster@ipp.mpg.de             !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!


! +++ Declaration of variables: 
    USE euitm_schemas
    USE copy_structures
    USE ets_plasma
    USE copy_cpo_ets


    IMPLICIT NONE

    INTEGER                        :: iimp

! +++ CPO derived types:
    TYPE (type_coreimpur), POINTER  :: coreimpur_iter(:) !input/output CPO with internal ETS parameters profiles 
    TYPE (type_coreimpur), POINTER  :: coreimpur_new(:)  !input/output CPO with internal ETS parameters profiles 


! +++ Internal ETS derived types:
    REAL (R8)                      :: control_double(5)!real control parameters
    REAL (R8)                      :: conv

! +++ local

    REAL (R8)                      :: err_imp
    REAL (R8), PARAMETER           :: nz_0 = 1.0e10_r8

! +++ Check convergence:

    err_imp = 0.0_r8

    do iimp = 1, size(coreimpur_iter(1)%IMPURITY)
       err_imp = max(err_imp, &
            maxval(abs(1.0_r8 - (abs(coreimpur_iter(1)%IMPURITY(iimp)%NZ(:,:))+nz_0)/  &
                                (abs(coreimpur_new(1)%IMPURITY(iimp)%NZ(:,:))+nz_0))))
       
       
    enddo

    !write(*,*) 'nimp err ', err_imp

    conv = err_imp

    RETURN 



  END SUBROUTINE check_convergence_impurities
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  



END MODULE convergence_check