hm_read_analy.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_analy (nanaly, iparith, ipari0, lsubmodel)

Function/Subroutine Documentation

hm_read_analy()

subroutine hm_read_analy	(	integer, intent(out)	nanaly,
		integer, intent(out)	iparith,
		integer, intent(out)	ipari0,
		type(submodel_data), dimension(*), intent(in)	lsubmodel )

Definition at line 36 of file hm_read_analy.F.

C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   PREREAD /ANALY USING HM_READER
C-----------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C     NAME            DESCRIPTION    
C     NANALY          N2D3D INTEGER
C     IPARITH         FLAG FOR PARAL. ARITH.
C     LSUBMODEL       SUBMODEL STRUCTURE     
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE hm_option_read_mod
      USE submodel_mod , ONLY : submodel_data
      USE eos_param_mod , ONLY : analy_temp
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C INPUT ARGUMENTS
      INTEGER, INTENT(OUT)          ::NANALY,IPARITH,IPARI0
      TYPE(SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
C OUTPUT ARGUMENTS
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: NB_ANALY
      INTEGER :: ITEMP
      LOGICAL :: IS_AVAILABLE
C============================================================================
      is_available = .false.
C----------------------------------------------
      CALL hm_option_count('/ANALY',nb_analy)
      IF (nb_analy>0) THEN
        CALL hm_option_start('/ANALY')
        CALL hm_option_read_key(lsubmodel)
        CALL hm_get_intv('N2D3D',nanaly,is_available,lsubmodel)      
        CALL hm_get_intv('IPARITH',iparith,is_available,lsubmodel)
        CALL hm_get_intv('ANALY_TEMP',itemp,is_available,lsubmodel)
        IF(iparith == 0)THEN
          ipari0  = 1
          iparith = 1
        ELSEIF(iparith == 1)THEN
          ipari0 = 1
        ELSEIF(iparith == 2)THEN
          ipari0 = 0
        ELSE
          ipari0=iparith-1
        ENDIF
        ! ANALY_TEMP : 0 No temperature evolution in expansion for solid EoS
        ! ANALY_TEMP : 1 Temperature evolution also in expansion for solid EoS (cutoff removed)
        analy_temp = 0
        IF(itemp == 1)analy_temp = 1
      ELSE
        nanaly  = 0
        iparith = 1
        analy_temp = 0
      ENDIF
 
C--------------------------------------------------------------
      RETURN