sensor_gauge.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "com08_c.inc"
#include "units_c.inc"
#include "comlock.inc"
#include "task_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sensor_gauge (sensor, gauge)

Function/Subroutine Documentation

sensor_gauge()

subroutine sensor_gauge	(	type (sensor_str_)	sensor,
		dimension(llgauge,*)	gauge )

Definition at line 31 of file sensor_gauge.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE sensor_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      "com08_c.inc"
#include      "units_c.inc"
#include      "comlock.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real, DIMENSION(LLGAUGE,*) :: gauge
      TYPE (SENSOR_STR_) :: SENSOR
C----------------------------------------------------------
C Local Variables
C----------------------------------------------------------
      INTEGER I,J,NGAUGE,IGAU,ICRIT
c
      my_real :: pmin,tmin,tstart,tdelay,infinity,maxp
      parameter(infinity = 1.0e20)
C=======================================================================
      IF (sensor%STATUS == 1) RETURN   ! already activated
c
      tdelay = sensor%TDELAY
      ngauge = sensor%IPARAM(1)
 
      j = 0
      maxp = - infinity
      DO i = 1,ngauge
        icrit  = 0
        igau   = sensor%IPARAM(i+1)
        pmin   = sensor%RPARAM(j+1)
        tmin   = sensor%RPARAM(j+2)
        tstart = sensor%RPARAM(j+3)
        IF (abs(gauge(30,igau)) > pmin) icrit  = 1
        IF(abs(gauge(30,igau)) > maxp) maxp = abs(gauge(30,igau))
c----------------------------------------------------------------
c       check activation status
c----------------------------------------------------------------
        IF (sensor%TCRIT + tmin > tt) THEN
          IF (icrit == 0) THEN
            sensor%TCRIT = infinity
          ELSE IF (sensor%TCRIT == infinity) THEN
            sensor%TCRIT = min(sensor%TCRIT, tt)
          END IF
        ELSE IF (sensor%RPARAM(j+3) == infinity) THEN
          sensor%RPARAM(j+3) = sensor%TCRIT + tmin + tdelay
        END IF
        IF (sensor%RPARAM(j+3) <= tt) THEN   ! sensor activation
          sensor%STATUS = 1
          sensor%TSTART = sensor%RPARAM(j+3)
        END IF
c  
        IF (sensor%STATUS == 1 .and. ispmd == 0) THEN
#include "lockon.inc"
          WRITE (istdo,1100) sensor%SENS_ID,sensor%TSTART
          WRITE (iout ,1100) sensor%SENS_ID,sensor%TSTART
          WRITE (iout ,1200) pmin,abs(gauge(30,igau))
#include "lockoff.inc"
        ENDIF
        j = j + 3
      END DO
      sensor%RESULTS(1) = maxp
c-----------------------------------------------------------------------      
1100  FORMAT(' SENSOR NUMBER ',i10,' ,ACTIVATED AT TIME ',1pe12.5)
1200  FORMAT('      TARGET PRESSURE = ',1pe12.5,/
     .       '      CURRENT PRESSURE AFTER TMIN AND TDELAY = ',1pe12.5)
c-----------------------------------------------------------------------      
      RETURN