fc2k_workaround_impurity.f90

Go to the documentation of this file.

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

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

  SUBROUTINE itmimpurity (EQUILIBRIUM_ITER, COREPROF_ITER,  CORETRANSP_ITER,      &
                          coreimpur_old,    coreimpur_iter, coreneutrals_iter,    &
                          coresource_iter,  coresource_new, coreimpur_new,        &
                          control_integer,  control_double)
 

  USE euitm_schemas 
  USE itm_types
  USE impurity
  USE neutrals
  

  
  IMPLICIT NONE
  
  TYPE (type_equilibrium), POINTER  :: equilibrium_iter(:)
  TYPE (type_coreprof),    POINTER  :: coreprof_iter(:)
  TYPE (type_coretransp),  POINTER  :: coretransp_iter(:)
  TYPE (type_coresource),  POINTER  :: coresource_new(:), coresource_iter(:)
  TYPE (type_coreimpur),   POINTER  :: coreimpur_old(:),coreimpur_iter(:),coreimpur_new(:)
  TYPE (type_coreneutrals),POINTER  :: coreneutrals_iter(:)
 
  INTEGER,               INTENT(IN) :: control_integer(2)    
  REAL (R8),             INTENT(IN) :: control_double(5)   
   
  ! PIS 2016-08-18 simple check to see that we indeed have impurities loaded if not -return          
  IF (.NOT. ASSOCIATED (coreimpur_old(1)%impurity) ) RETURN

  CALL impurity_ets(equilibrium_iter, coreprof_iter,  coretransp_iter,      &
                      coreimpur_old,    coreimpur_iter, coreneutrals_iter,    &
                      coresource_iter,  coresource_new, coreimpur_new,        &
                      control_integer,  control_double)
     
  
  RETURN
  
  END SUBROUTINE itmimpurity