hm_read_ale_muscl.F File Reference

#include "implicit_f.inc"
#include "scr05_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_ale_muscl (lsubmodel, unitab)

Function/Subroutine Documentation

hm_read_ale_muscl()

subroutine hm_read_ale_muscl	(	type(submodel_data), dimension(nsubmod), intent(in)	lsubmodel,
		type(unit_type_), intent(in)	unitab )

Definition at line 39 of file hm_read_ale_muscl.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C This subroutine is reading /ALE/MUSCL option in user input file
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE hm_option_read_mod
      USE submodel_mod
      USE message_mod
      USE alemuscl_mod , only : alemuscl_param
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE(SUBMODEL_DATA), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
      TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------     
      INTEGER :: NALEMUSCL,IALEMUSCL
      LOGICAL :: IS_AVAILABLE
      my_real :: beta
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------       
      CALL hm_option_count('/ALE/MUSCL', nalemuscl)
      
      IF (nalemuscl > 0) THEN
         CALL hm_option_start('/ALE/MUSCL')
         CALL hm_option_next()
         CALL hm_get_floatv('ALE_MUSCL_BETA', beta, is_available, lsubmodel, unitab)
         CALL hm_get_intv('MAT_Iflag', ialemuscl, is_available, lsubmodel)
         
         !---default BETA
         IF (beta == zero) THEN
           beta = two
         ENDIF
 
         !---default IALEMUSCL
         IF (ialemuscl == 0) THEN
           ialemuscl = 1
         ELSEIF (ialemuscl == 1) THEN
           ialemuscl = 2
         ELSEIF (ialemuscl < -1 .OR. ialemuscl >= 2) THEN
           ialemuscl = 1
         ENDIF
                  
         !---MUSCL off / disabled
         IF(ialemuscl == -1)THEN
           ialemuscl = 0
           beta = two
           nalemuscl=0
         ENDIF
         
         alemuscl_param%IALEMUSCL = ialemuscl
         alemuscl_param%Beta = beta
         
         IF (iresp == 1) THEN
            ! Not available in single precision
            CALL ancmsg(msgid = 1563, anmode = aninfo, msgtype = msgerror)
         ENDIF
      ENDIF
C-----------------------------------------------