qmomk_make_arrays.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: QMOMK_MAKE_ARRAYS                                      C
!  Purpose: DES - Initialize/Restart QMOMK arrays                      C
!                                                                      C
!                                                                      C
!  Author: Alberto Passalacqua                        Date:            C
!  Reviewer:                                          Date:            C
!  Comments:                                                           C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

SUBROUTINE qmomk_make_arrays

  USE compar
  USE exit, only: mfix_exit
  USE funits
  USE geometry
  USE param1
  USE qmom_kinetic_equation
  USE run

  IMPLICIT NONE

  INTEGER CHECK_MPI

  IF(coordinates == 'CYLINDRICAL') THEN
     WRITE (unit_log, *) ' '
     WRITE (unit_log, *) 'Cylindrical coordinates are being used. STOP'
     WRITE (unit_log, *) 'QMOMK should only be run using cartesian coordinates.'
     WRITE (*, *) ' '
     WRITE (*, *) 'Cylindrical coordinates are being used. STOP'
     WRITE (*, *) 'QMOMK should only be run using cartesian coordinates.'
     CALL mfix_exit(mype)
  END IF

  check_mpi = nodesi * nodesj * nodesk
  IF((check_mpi.NE.1).AND.(qmomk)) THEN
     WRITE (unit_log, *) ' '
     WRITE (unit_log, *) 'QMOMK being run on multiple processors. STOP'
     WRITE (unit_log, *) 'QMOMK should only be run serially on one processor.'
     WRITE (*, *) ' '
     WRITE (*, *) 'QMOMK being run on multiple processors. STOP'
     WRITE (*, *) 'QMOMK should only be run serially on one processor.'
     CALL mfix_exit(mype)
  END IF

  IF(run_type == 'RESTART_1') THEN !  Read Restart
     CALL qmomk_read_restart
     WRITE(*,*) 'QMOMK_RES file read at Time= ', time
     WRITE(unit_log,*) 'QMOMK_RES file read at Time= ', time
  ELSE IF (run_type == 'RESTART_2') THEN
     WRITE(unit_log,*) 'Restart 2 is not implemented with QMOMK'
     WRITE(*,*) 'Restart 2 is not implemented with QMOMK'
     CALL mfix_exit(mype)
  END IF

  RETURN
END SUBROUTINE qmomk_make_arrays