reinitialize.f

Go to the documentation of this file.

MODULE reinit
   CONTAINS
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  SUBROUTINE: REINITIALIZE                                            !
!  Purpose: read and verify input data, open files                     !
!                                                                      !
!  Author: P. Nicoletti                               Date: 04-DEC-91  !
!  Reviewer: M.SYAMLAL, W.ROGERS, P.NICOLETTI         Date: 24-JAN-92  !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE reinitialize(filename)

      use run, only: reinitializing

      use error_manager

      IMPLICIT NONE

      CHARACTER(LEN=*), intent(in) :: filename

      INTEGER :: IER

      ier = 0
      reinitializing = .true.

! Clean up reaction data if needed
      CALL reinit_rxn_data

! Read in the namelist variables from the ascii input file.
      CALL read_namelist(2, filename)

      CALL reinitialize0(ier)

      reinitializing = .false.

      IF(ier /=0) THEN
         WRITE(err_msg, 2000)
      ELSE
         WRITE(err_msg, 2100)
      ENDIF

 2000 FORMAT(2/70('*'),/'Reinitialization failed!',/'Correct all ',    &
         'reported errors and reinitialize again.',/70('*'))

 2100 FORMAT(2/,70('*'),/'Successfully reinitialized!'/70('*'))
      CALL flush_err_msg(header=.false., footer=.false.)

      RETURN
      END SUBROUTINE reinitialize

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  SUBROUTINE: REINITIALIZE0                                           !
!  Purpose: read and verify input data, open files                     !
!                                                                      !
!  Author: P. Nicoletti                               Date: 04-DEC-91  !
!  Reviewer: M.SYAMLAL, W.ROGERS, P.NICOLETTI         Date: 24-JAN-92  !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE reinitialize0(pIER)

      USE cutcell, only: cartesian_grid
      use coeff, only: init_coeff

      use error_manager

      IMPLICIT NONE

      INTEGER, INTENT(OUT) :: pIER
      INTEGER :: IER

! Set the default error flag to ERROR state.
      pier = 1

      CALL check_run_control()
      IF(reinit_error()) RETURN

      CALL check_numerics()
      IF(reinit_error()) RETURN

      CALL check_output_control()
      IF(reinit_error()) RETURN

      CALL check_gas_phase
      IF(reinit_error()) RETURN

      CALL check_solids_phases
      IF(reinit_error()) RETURN

      CALL check_initial_conditions
      IF(reinit_error()) RETURN
      CALL check_boundary_conditions
      IF(reinit_error()) RETURN
      CALL check_internal_surfaces
      IF(reinit_error()) RETURN
      CALL check_point_sources

      CALL check_chemical_rxns
      IF(reinit_error()) RETURN
      CALL check_odepack_stiff_chem
      IF(reinit_error()) RETURN

! Convert (mass, volume) flows to velocities.
      CALL set_bc_flow
      IF(reinit_error()) RETURN

! This is all that happens in SET_L_SCALE so it needs moved, maybe
! this should go in int_fluid_var.?
!     CALL SET_L_SCALE
!      L_SCALE(:) = L_SCALE0

! Set constant physical properties
      CALL set_constprop
      IF(reinit_error()) RETURN

! Set initial conditions
      CALL set_ic
      IF(reinit_error()) RETURN

! Set point sources.
      CALL set_ps
      IF(reinit_error()) RETURN

! Set boundary conditions
      CALL zero_norm_vel
      IF(reinit_error()) RETURN
      CALL set_bc0
      IF(reinit_error()) RETURN
      IF(cartesian_grid) CALL cg_set_bc0
      IF(reinit_error()) RETURN

! Set gas mixture molecular weight
      CALL set_mw_mix_g
      IF(reinit_error()) RETURN

! Initialize densities.
      CALL set_ro_g
      IF(reinit_error()) RETURN
      CALL set_ro_s
      IF(reinit_error()) RETURN

! Initialize time dependent boundary conditions
      CALL set_bc1
      IF(reinit_error()) RETURN

! Check the field variable data and report errors.
      IF(.NOT.cartesian_grid) CALL check_data_20
      IF(reinit_error()) RETURN

! Parse residual strings
      CALL parse_resid_string ()
      IF(reinit_error()) RETURN

      CALL rrates_init(ier)
      IF(reinit_error()) RETURN

! Calculate all the coefficients once before entering the time loop
      CALL init_coeff(ier)
      IF(reinit_error()) RETURN

! After reinitialization, the field vars should pass these checks too
      CALL check_data_30()
      IF(reinit_error()) RETURN

      CALL write_res0()
      CALL write_out0()

! Made it here without error.
      pier = 0

      RETURN
      END SUBROUTINE reinitialize0


!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  SUBROUTINE: REINIT_RXN_DATA                                         !
!  Purpose: read and verify input data, open files                     !
!                                                                      !
!  Author: P. Nicoletti                               Date: 04-DEC-91  !
!  Reviewer: M.SYAMLAL, W.ROGERS, P.NICOLETTI         Date: 24-JAN-92  !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE reinit_rxn_data

      use parse, only: rxn_name, des_rxn_name
      use parse, only: rxn_chem_eq, des_rxn_chem_eq
      use parse, only: usrdh, des_usrdh
      use parse, only: usrfdh, des_usrfdh

      use rxns, only: no_of_rxns
      use rxns, only: reaction
      use des_rxns, only: no_of_des_rxns
      use des_rxns, only: des_reaction

      use error_manager

      IMPLICIT NONE

      INTEGER :: LC

! Reaction Names: Allocate/Initialize
      IF(allocated( rxn_name )) deallocate(rxn_name)
! Chemical Equations: Allocate/Initialize
      IF(allocated( rxn_chem_eq )) deallocate(rxn_chem_eq)
! User defined heat of reaction: Allocate/Initialize
      IF(allocated( usrdh ))deallocate(usrdh)
! User defined heat of reaction partitions: Allocate/Initialize
      IF(allocated( usrfdh )) deallocate(usrfdh)


      IF(allocated(reaction)) THEN
         DO lc=1,no_of_rxns
            IF(allocated(reaction(lc)%HoR)) &
               deallocate(reaction(lc)%HoR)
            IF(allocated(reaction(lc)%rPhase)) &
               deallocate(reaction(lc)%rPhase)
            IF(allocated(reaction(lc)%Species)) &
               deallocate(reaction(lc)%Species)
         ENDDO
         deallocate(reaction)
      ENDIF
      no_of_rxns = 0

! Reaction Names: Allocate/Initialize
      IF(allocated( des_rxn_name )) deallocate(des_rxn_name)
! Chemical Equations: Allocate/Initialize
      IF(allocated( des_rxn_chem_eq )) deallocate(des_rxn_chem_eq)
! User defined heat of reaction: Allocate/Initialize
      IF(allocated( des_usrdh )) deallocate( des_usrdh)
! User defined heat of reaction partitions: Allocate/Initialize
      IF(Allocated( des_usrfdh )) deallocate( des_usrfdh)

      IF(allocated(des_reaction)) THEN
         DO lc=1,no_of_des_rxns
            IF(allocated(des_reaction(lc)%HoR)) &
               deallocate(des_reaction(lc)%HoR)
            IF(allocated(des_reaction(lc)%rPhase)) &
               deallocate(des_reaction(lc)%rPhase)
            IF(allocated(des_reaction(lc)%Species)) &
               deallocate(des_reaction(lc)%Species)
         ENDDO
         deallocate(des_reaction)
      ENDIF
      no_of_des_rxns = 0


      RETURN
      END SUBROUTINE reinit_rxn_data
END MODULE reinit