33
34
35
36 USE sensor_mod
37
38
39
40#include "implicit_f.inc"
41
42
43
44#include "units_c.inc"
45#include "comlock.inc"
46#include "task_c.inc"
47
48
49
50 my_real ,
INTENT(IN) :: time,timestep
51 TYPE (SENSOR_STR_) ,INTENT(INOUT) :: SENSOR
52
53
54
55 INTEGER :: NSENS
56 my_real :: tdelay,tstop,infinity
57 DOUBLE PRECISION :: TDELAY_DP,TSTOP_DP,NEXT_CYCLE
58 parameter(infinity = 1.0e20)
59
60 nsens = sensor%SENS_ID
61 tdelay = sensor%TDELAY
62 tstop = sensor%VALUE
63 next_cycle = time + timestep
64 tdelay_dp = tdelay
65 tstop_dp = tstop
66
67 IF (next_cycle >= tdelay_dp .AND. sensor%STATUS == 0) THEN
68 sensor%STATUS = 1
69 sensor%TSTART = tdelay
70 IF (ispmd == 0) THEN
71#include "lockon.inc"
72 WRITE (iout ,1000) sensor%SENS_ID,tdelay
73 WRITE (istdo,1000) sensor%SENS_ID,tdelay
74#include "lockoff.inc"
75 ENDIF
76 ELSE IF (next_cycle >= tstop_dp .AND. sensor%STATUS == 1) THEN
77 sensor%STATUS = 0
78 sensor%TSTART = infinity
79 sensor%TDELAY = infinity
80 IF (ispmd == 0) THEN
81#include "lockon.inc"
82 WRITE (iout ,2000) sensor%SENS_ID,tstop
83 WRITE (istdo,2000) sensor%SENS_ID,tstop
84#include "lockoff.inc"
85 ENDIF
86 ENDIF
87
881000 FORMAT(' SENSOR (TIME) NUMBER ',i10,' ACTIVATED AT TIME ',1pe12.5)
892000 FORMAT(' SENSOR (TIME) NUMBER ',i10,' DESACTIVATED AT TIME ',1pe12.5)
90
91 RETURN