hm_read_ale_close.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "com01_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_ale_close (unitab, lsubmodel, geo)

Function/Subroutine Documentation

hm_read_ale_close()

subroutine hm_read_ale_close	(	type(unit_type_), intent(in)	unitab,
		type(submodel_data), dimension(nsubmod), intent(in)	lsubmodel,
		intent(inout)	geo )

Definition at line 34 of file hm_read_ale_close.F.

C-----------------------------------------------
C   Description
C-----------------------------------------------
C This subroutine is reading option /ALE/CLOSE
C  in Starter input file
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE submodel_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-----------------------------------------------
#include      "param_c.inc"
#include      "com01_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
      TYPE(SUBMODEL_DATA), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
      my_real, DIMENSION(NPROPG), INTENT(INOUT) :: geo
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real :: htest, hclose
      LOGICAL :: IS_AVAILABLE1, IS_AVAILABLE2
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      iclose = 0
      CALL hm_get_floatv('htest', htest, is_available1, lsubmodel, unitab)
      CALL hm_get_floatv('hclose', hclose, is_available2, lsubmodel, unitab)
      IF (is_available1 .OR. is_available2) THEN
         iclose = 1
         IF (htest == zero) htest = em20
         IF (hclose == zero) hclose = htest * em01
         hclose = min(htest, hclose)
         geo(129) = hclose
         geo(130) = htest
         WRITE(iout, 2010) htest, hclose
      ENDIF
      RETURN
C-----------------------------------------------
 2010 FORMAT(
     & 5x,'    ALE ELEMENT CLOSURE  ',/,
     & 5x,'    -------------------  ',/,
     & 5x,'MAXIMUM LENGTH FOR ACTIVATION . . . . .=',1pg20.13/,
     & 5x,'MAXIMUM LENGTH FOR RESISTANCE . . . . .=',1pg20.13/)
C-----------------------------------------------