37
38
39
42 USE multi_fvm_mod
45 USE ebcs_mod
47
48
49
50#include "implicit_f.inc"
51
52
53
54#include "units_c.inc"
55#include "com01_c.inc"
56#include "com04_c.inc"
57
58
59
60 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
61 INTEGER ID,UID
62 TYPE (MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM
63 TYPE (SURF_) ,TARGET, DIMENSION(NSURF) :: IGRSURF
64 CHARACTER(LEN=NCHARTITLE), INTENT(IN) :: TITR
65 TYPE(SUBMODEL_DATA) LSUBMODEL(NSUBMOD)
66 TYPE(t_ebcs_monvol), INTENT(INOUT) :: EBCS
67
68
69
70 INTEGER NGR2USR,ISU,NSEG
71 INTEGER SURF,SENS,MONVOL
73 LOGICAL IS_AVAILABLE,IS_ENCRYPTED
74 INTEGER, DIMENSION(:), POINTER :: INGR2USR
75 INTEGER IFLAGUNIT,J
76
77
78
79
80 ebcs%title = trim(titr)
81
82 iflagunit=0
83 DO j=1,unitab%NUNITS
84 IF (unitab%UNIT_ID(j) == uid) THEN
85 iflagunit = 1
86 EXIT
87 ENDIF
88 ENDDO
89 IF (uid/=0.AND.iflagunit==0) THEN
90 CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,
91 . i2=uid,
93 . c1='EBCS',
94 . c2='EBCS',
95 . c3=titr)
96 ENDIF
97
99 CALL hm_get_intv(
'entityid', surf ,is_available,lsubmodel)
100 CALL hm_get_intv(
'sensorid', sens ,is_available,lsubmodel)
101 CALL hm_get_intv(
'monvol_ID', monvol ,is_available,lsubmodel)
102
103 isu=0
104 ingr2usr => igrsurf(1:nsurf)%ID
105 IF (surf/=0) isu=
ngr2usr(surf,ingr2usr,nsurf)
106 nseg=0
107 IF (isu/=0) nseg=igrsurf(isu)%NSEG
108 IF (surf==0)THEN
109 ierr=ierr+1
110 WRITE(istdo,'(6X,A)')' ** A SURFACE SHOULD BE INPUT'
111 WRITE(iout, '(6X,A)')' ** A SURFACE SHOULD BE INPUT'
112 ELSEIF(isu==0)THEN
113 ierr=ierr+1
114 WRITE(istdo,*)' ** ERROR SURFACE NOT FOUND, ID=',surf
115 WRITE(iout,*) ' ** ERROR SURFACE NOT FOUND, ID=',surf
116 ELSEIF(nseg==0)THEN
117 ierr=ierr+1
118 WRITE(istdo,*)' ** ERROR EMPTY SURFACE, ID=',surf
119 WRITE(iout,*) ' ** ERROR EMPTY SURFACE, ID=',surf
120 ELSE
121 WRITE(iout,1008)
id,trim(titr),monvol,surf,sens,nseg
122 ENDIF
123
124 ebcs%title = titr
125 ebcs%monvol_id = monvol
126 ebcs%sensor_id = sens
127
128 IF (multi_fvm%IS_USED) THEN
129 CALL ancmsg(msgid = 1602, msgtype = msgerror, anmode = aninfo,
130 . i1 =
id, c1 = trim(titr), c2 =
"NOT COMPATIBLE WITH LAW 151")
131 ENDIF
132
133 IF(nspmd > 1 ) THEN
134 CALL ancmsg(msgid=755,msgtype=msgerror,anmode=aninfo,c1=
'ELEMENTARY BOUNDARY CONDITIONS')
135 ENDIF
136
137
138 RETURN
139
140
141 1008 FORMAT(//
142 . 'EBCS NUMBER :',i8,1x,a,/,
143 . ' CONNECTION TO MONITORED VOLUME. . . . . . ',i8,/,
144 . ' THROUGH SURFACE . . . . . . . . . . . . . ',i8,/,
145 . ' ACTIVATED BY SENSOR . . . . . . . . . . . ',i8,/,
146 . ' NUMBER OF SEGMENTS FOUND. . . . . . . . . ',i8,/)
147
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
subroutine hm_option_is_encrypted(is_encrypted)
integer, parameter nchartitle
integer function ngr2usr(iu, igr, ngr)
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)