fc2k_workaround_ets.f90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + + 
!-------------------------------------------------------!
!     This routine provides a wrapper for the           !
!     FC2K actor generator to call the actual           !
!     transport solver routine.                         !
!-------------------------------------------------------!
!     Source:       ---                                 !
!     Developers:   D.Kalupin                           !
!     Kontacts:     Denis.Kalupin@euro-fusion.org       !
!                                                       !
!     Comments:     ---                                 !
!                                                       !
!-------------------------------------------------------!
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
SUBROUTINE fc2k_transport_solver                              &
       (coreprof_old,      coreprof_iter,      coreprof_new,  &
        equilibrium_old,   equilibrium_iter,                  &
        coretransp_iter,                                      &
        coresource_iter,                                      &
        coreimpur_iter,                                       &
        corefast_iter,                                        &
        control_integer,   control_double,                    & 
!       CONTROL_INTEGER, CONTROL_DOUBLE, HYPER_DIFF) !AF 25.Apr.2016, 22.Aug.2016
        hyper_diff,user_out_outputflag,user_out_diagnosticinfo)  

  USE itm_types
  USE euitm_schemas
  USE ets_plasma
  USE ets
! module needed when moving to TRANSPORT_SOLVER
!  USE TRANSPORT_SOLVER
!  USE ISO_C_BINDING 
 
  IMPLICIT NONE

  TYPE (type_coreprof),    POINTER    :: coreprof_old(:)     
  TYPE (type_coreprof),    POINTER    :: coreprof_iter(:)    
  TYPE (type_coreprof),    POINTER    :: coreprof_new(:)
  TYPE (type_equilibrium), POINTER    :: equilibrium_old(:)  
  TYPE (type_equilibrium), POINTER    :: equilibrium_iter(:)
  TYPE (type_coretransp),  POINTER    :: coretransp_iter(:)
  TYPE (type_coresource),  POINTER    :: coresource_iter(:)
  TYPE (type_coreimpur),   POINTER    :: coreimpur_iter(:)
  TYPE (type_corefast),    POINTER    :: corefast_iter(:)

  TYPE (type_param)                   :: code_parameters

  TYPE (diagnostic)                   :: diag                

  INTEGER,                 INTENT(IN) :: control_integer(4) !integer control parameters  
  REAL (R8),               INTENT(IN) :: control_double(6)  !real control parameters   
 
! +++ Stabilization scheme !AF 25.Apr.2016, 22.Aug.2016
  REAL (R8), DIMENSION(2)           :: hyper_diff          !hyper diffusivity

  integer       :: ifail
  integer       :: nch_tmp

!----  Diagnostic info  ----
integer, intent(out)     ::     user_out_outputflag
character(len=:), pointer, intent(out)    ::    user_out_diagnosticinfo


!initial values of diagnostic info
user_out_outputflag=0
!reset diag type
diag%ierr=0
diag%error_message=' '

! Call the actual code:
  CALL transport_solver_interface                        &
!  OLD:              ITER:                NEW:
!  Previous          Previous             output
!  time step         iteration
  (coreprof_old,     coreprof_iter,       coreprof_new,  &   !1-D profiles of plasma parameters
   equilibrium_old,  equilibrium_iter,                   &   !Equilibrium quantities / geometry
                     coretransp_iter,                    &   !transport coefficients
                     coresource_iter,                    &   !sources
                     coreimpur_iter,                     &   !impurity species
                     corefast_iter,                      &   !non-thermal components
!
   control_integer, control_double,                      &   !control switches
!  CONTROL_INTEGER, CONTROL_DOUBLE, HYPER_DIFF,  ifail,  &   !AF 25.Apr.2016, 22.Aug.2016
   hyper_diff,                                           &   !AF 25.Apr.2016, 22.Aug.2016
!  code_parameters)
   diag)                                                     !run diagnostic (errors)         

   if (diag%ierr .lt.0) then
      user_out_outputflag=diag%ierr
      nch_tmp=len_trim(diag%error_message)
      allocate(character(nch_tmp)  :: user_out_diagnosticinfo)
      user_out_diagnosticinfo=diag%error_message(1:nch_tmp)
   end if
      
  RETURN
  
END SUBROUTINE fc2k_transport_solver
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +