Homec4science
Diffusion COSOlver (new_coulomb) Browse

cosolver/src/control.f90
83b8f38f6c51new_coulomb

Last Change
Raw File

control.f90
View Options

SUBROUTINE control
USE basic
use basic_mpi
USE self_colls
USE ei_colls
USE ie_colls
USE prec_const
USE operator_model
USE basis_transformation
USE output
USE timer
IMPLICIT NONE
integer :: ik
! 0. Init
!
CALL ppinit
CALL daytim('Started at ')
!
IF (me .eq. 0 ) WRITE(*,*) 'MPI initialized ...'
CALL init_prec_const
! 1. Prologue
! 1.1 Read inputs/MPI
CALL prologue
! 2. Main loop - Collisional Moment
IF(me .eq. 0) WRITE(*,*) '============= Collision Operator Model ============='
! 2.01 Initialize collision operator models
IF( me .eq. 0) WRITE(*,*) 'Initialize collision operator models ...'
CALL set_operator_model
IF(me .eq. 0) WRITE(*,*) '.... operator models initialized'
! 2.02 Allocate Collision Matrices
IF( me .eq. 0) WRITE(*,*) 'Allocate collisional matrices ...'
CALL memory
IF( me .eq. 0) WRITE(*,*) '...Allocation done'
! ..... MPI WAIT ....
CALL mpi_barrier(MPI_COMM_WORLD,ierr)
! 2.03 Load Basis transformation Coefficients
IF( me .eq. 0) WRITE(*,*) 'Load basis Transformation coefficients ....'
CALL load_basis_transformation
IF( me .eq.0 ) WRITE(*,*) '... done '
! ..... MPI WAIT ....
CALL mpi_barrier(MPI_COMM_WORLD,ierr)
! start clock
call clock(1)
do ik= ikp_start, ike
if( kpscan) call update_kperp( ik)
! 1.2 Initialize output file
CALL diagnose(0)
!----------------------------------------------------------------
! 2.1 Electron-Ion Collisions
IF(eicolls .and. is_e) THEN
call ei_colls_control
ENDIF
! 2.2 Electron-electron Collisions
IF(eecolls .and. is_self) THEN
call ee_colls_control
ENDIF
! 2.3 Ion-Electron Collisions
IF (iecolls .and. is_i) THEN
call ie_colls_control
ENDIF
! 2.2 Ion-Ion Collisions
IF(iicolls .and. is_self) THEN
call ii_colls_control
ENDIF
!----------------------------------------------------------------
! 3. Epilogue
! 3.2 MPI sync
CALL ppsync
! End clock
call clock(-1)
! 3.2 Write output
CALL mpi_barrier(MPI_COMM_WORLD,ierr)
CALL diagnose(-1)
CALL mpi_barrier(MPI_COMM_WORLD,ierr)
!
enddo
IF(me .eq. 0) CALL daytim('Done at ')
IF(me .eq. 0) WRITE(*,*) 'Ellapsed Time [s]: ',system_ellapsed_time
CALL ppexit
!
END SUBROUTINE control
c4science · Help