read_sensor_gauge.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	read_sensor_gauge (sensor_ptr, sens_id, titr, lgauge, unitab, lsubmodel)

Function/Subroutine Documentation

read_sensor_gauge()

subroutine read_sensor_gauge	(	type (sensor_str_)	sensor_ptr,
		integer, intent(in)	sens_id,
		character(len=nchartitle)	titr,
		integer, dimension(3,*), intent(in)	lgauge,
		type (unit_type_), intent(in)	unitab,
		type (submodel_data), dimension(nsubmod)	lsubmodel )

Definition at line 38 of file read_sensor_gauge.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------  
      USE unitab_mod
      USE message_mod
      USE submodel_mod
      USE sensor_mod
      USE hm_option_read_mod
      USE names_and_titles_mod , ONLY : nchartitle
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      "com04_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ,INTENT(IN) :: SENS_ID
      INTEGER ,DIMENSION(3,*)    ,INTENT(IN)  :: LGAUGE
      CHARACTER(LEN=NCHARTITLE)::TITR
      TYPE (SENSOR_STR_) :: SENSOR_PTR
      TYPE (SUBMODEL_DATA) ,DIMENSION(NSUBMOD) :: 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 :: I,J,NGAU,NPARIS,NPARRS,NVAR,SENS_TYPE
      my_real :: tdel
      INTEGER ,DIMENSION(6) :: IGAU,G_ID
      my_real ,DIMENSION(6) :: pres,tmin
      LOGICAL :: IS_AVAILABLE
C--------------------------------
C     SENSOR : PRESSURE GAUGE LIST
C=======================================================================
      is_available = .false.
      sens_type = 10
c--------------------------------------------------
card1
      CALL hm_get_floatv('Tdelay' ,tdel  ,is_available,lsubmodel,unitab)
card2
      CALL hm_get_intv  ('NIP'    ,ngau  ,is_available,lsubmodel)
card3...
      DO i=1,min(ngau,6)
        CALL hm_get_int_array_index  ('GaugeArrayIds',g_id(i),i,is_available,lsubmodel)
        CALL hm_get_float_array_index('Fporp' ,pres(i),i,is_available,lsubmodel,unitab)
        CALL hm_get_float_array_index('Fport' ,tmin(i),i,is_available,lsubmodel,unitab)
      ENDDO
c--------------------------------------------------
      IF (ngau > 6) THEN
         CALL ancmsg(msgid=833, msgtype=msgerror, anmode=aninfo, i1=sens_id,c1=titr)
      ENDIF
c
      DO i=1,ngau
        igau(i) = 0
        DO j = 1,nbgauge
          IF (g_id(i) == lgauge(2,j)) THEN
            igau(i) = j
            EXIT
          ENDIF
        ENDDO
        IF (igau(i) == 0) THEN
          CALL ancmsg(msgid=834, msgtype=msgerror, anmode=aninfo_blind_2,
     .                i1=sens_id, c1=titr)
          EXIT
        END IF
      ENDDO
c--------------------------------------------------
c
      sensor_ptr%TYPE    = sens_type
      sensor_ptr%SENS_ID = sens_id
      sensor_ptr%STATUS  = 0            ! status = deactivated
      sensor_ptr%TSTART  = infinity
      sensor_ptr%TCRIT   = infinity
      sensor_ptr%TMIN    = zero         ! TMIN global
      sensor_ptr%TDELAY  = tdel         ! time delay before activation
      sensor_ptr%VALUE   = zero
 
      nparis = ngau + 1
      nparrs = ngau * 3
      nvar   = 0 
c
      sensor_ptr%NPARI  = nparis
      sensor_ptr%NPARR  = nparrs
      sensor_ptr%NVAR   = nvar
c
      ALLOCATE (sensor_ptr%IPARAM(nparis))
      ALLOCATE (sensor_ptr%RPARAM(nparrs))
      ALLOCATE (sensor_ptr%VAR(nvar))
      sensor_ptr%VAR(:) = zero
c
      sensor_ptr%IPARAM(1) = ngau
      DO i = 1,ngau
        sensor_ptr%IPARAM(i+1) = igau(i)
      END DO
      j = 0
      DO i = 1,ngau
        sensor_ptr%RPARAM(j+1) = pres(i)
        sensor_ptr%RPARAM(j+2) = tmin(i)
        sensor_ptr%RPARAM(j+3) = infinity       
        j = j+3
      END DO    
c------------------------------------------------------------
      WRITE(iout, 1000) sens_id,tdel
      WRITE(iout, 2000) ngau
      DO i = 1,ngau
        WRITE(iout, 3000) g_id(i),pres(i),tmin(i)
      END DO
c------------------------------------------------------------
 1000 FORMAT(
     . 5x,'    SENSOR TYPE 10: PRESSURE GAUGE        '/,
     . 5x,'    ------------------------------        '/,
     . 5x,'SENSOR ID. . . . . . . . . . . . . . . . .=',i10/
     . 5x,'TIME DELAY BEFORE ACTIVATION . . . . . . .=',e12.4)
 2000 FORMAT(
     . 5x,'NUMBER OF GAUGES . . . . . . . . . . . . .=',i10)
 3000 FORMAT(
     . 5x,'   GAUGE ID .  . . . . . . . . . . . . . .=',i10/
     . 5x,'   MIN PRESSURE FOR ACTIVATION . . . . . .=',e12.4/
     . 5x,'   MINIMUM PRES. DURATION FOR ACTIVATION .=',e12.4/)
c-----------
      RETURN