hm_read_upwind.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_upwind (jupw, eta1, eta2, eta3, lsubmodel, unitab)

Function/Subroutine Documentation

hm_read_upwind()

subroutine hm_read_upwind	(	integer, intent(inout)	jupw,
		intent(inout)	eta1,
		intent(inout)	eta2,
		intent(inout)	eta3,
		type(submodel_data), dimension(nsubmod), intent(in)	lsubmodel,
		type(unit_type_), intent(in)	unitab )

Definition at line 38 of file hm_read_upwind.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE hm_option_read_mod
      USE submodel_mod
      USE message_mod
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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: JUPW
      my_real, INTENT(INOUT) :: eta1, eta2, eta3
      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 :: NUPWIND
      LOGICAL :: IS_AVAILABLE
 
      jupw = 0
      eta1 = one
      eta2 = one
      eta3 = one
      
      CALL hm_option_count('/UPWIND', nupwind)
      
      IF (nupwind > 0) THEN
         CALL hm_option_start('/UPWIND')
         CALL hm_option_next()
         jupw = 1
         CALL hm_get_floatv('ENG_UPWIND_eta1', eta1, is_available, lsubmodel, unitab)
         CALL hm_get_floatv('ENG_UPWIND_eta2', eta2, is_available, lsubmodel, unitab)
         CALL hm_get_floatv('ENG_UPWIND_eta3', eta3, is_available, lsubmodel, unitab)
         IF (eta1 == zero) eta1 = one
         IF (eta2 == zero) eta2 = one
         IF (eta3 == zero) eta3 = one
         
         IF (eta3 /= one) THEN
            CALL ancmsg(msgid = 1628, msgtype = msgwarning, anmode = aninfo)
         ENDIF
      ENDIF