qmomk_read_restart.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: QMOMK_READ_RESTART                                     C
!  Purpose: Reading QMOMK data for restart                             C
!                                                                      C
!                                                                      C
!  Author: Alberto Passalacqua                        Date:            C
!  Reviewer:                                          Date:            C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
SUBROUTINE qmomk_read_restart

  USE param
  USE param1
  USE constant
  USE fldvar
  USE cont
  USE geometry
  USE indices
  USE run
  USE compar
  USE physprop
  USE qmom_kinetic_equation
  USE qmomk_quadrature
  USE functions

  IMPLICIT NONE

  INTEGER :: M, IJK

  OPEN (unit=901, file=trim(run_name)//'_QMOMK.RES', form='Unformatted', status='unknown')

  rewind(901)

  print *,'QMOMK: Reading restart info...'

  READ (901) qmomk_n1
  READ (901) qmomk_u1
  READ (901) qmomk_v1
  READ (901) qmomk_w1

  print *,'QMOMK: Updating moments after restart...'

  DO m = 1, mmax
    DO ijk = ijkstart3, ijkend3
     IF (fluid_at(ijk))  THEN
       CALL moments_twenty_eight_nodes (qmomk_n1(:,ijk,m), &
            qmomk_u1(:,ijk,m), qmomk_v1(:,ijk,m), qmomk_w1(:,ijk,m), qmomk_m1(:,ijk,m))
       CALL eight_node_3d (qmomk_m1(:,ijk,m), qmomk_n1(:,ijk,m), &
            qmomk_u1(:,ijk,m), qmomk_v1(:,ijk,m), qmomk_w1(:,ijk,m))
       CALL moments_twenty_eight_nodes (qmomk_n1(:,ijk,m), &
            qmomk_u1(:,ijk,m), qmomk_v1(:,ijk,m), qmomk_w1(:,ijk,m), qmomk_m1(:,ijk,m))
     END IF
    END DO
  END DO

  qmomk_n0 = qmomk_n1
  qmomk_u0 = qmomk_u1
  qmomk_v0 = qmomk_v1
  qmomk_w0 = qmomk_w1
  qmomk_m0 = qmomk_m1

  print *,'QMOMK: Restart successful!'

END SUBROUTINE qmomk_read_restart