29
30
31
32#include "implicit_f.inc"
33
34
35
36#include "com04_c.inc"
37#include "param_c.inc"
38#include "units_c.inc"
39#include "task_c.inc"
40
41
42
43 INTEGER NFVMODI, MONVOL(*)
45 . volmon(*)
46
47
48
49 INTEGER I, ID, K1, KK1, N, ITYP, IGMERG, ILVOUT, IFOUND
51 . cgmerg, cnmerg
52
53 IF(ispmd==0) WRITE (iout,1000)
54 DO i=1,nfvmodi
56 READ(iin,'(I10,2F20.0,I10)') igmerg, cgmerg, cnmerg,ilvout
57
58 ifound=0
59 k1=1
60 kk1=1
61 DO n=1,nvolu
62 IF (monvol(k1) ==
id)
THEN
63 ityp=monvol(k1-1+2)
64 ifound=1
65 IF (ityp/=6.AND.ityp/=8) THEN
66 IF(ispmd==0)THEN
67 WRITE(iout,'(A)') ' '
68 WRITE(iout,'(A)') ' /FVMBAG/MODIF : '
69 WRITE(iout,
'(A,I10,A)')
' MONITORED VOLUME ID ',
id,
70 . ' IS NOT OF TYPE FVMBAG OR'
71 WRITE(iout,'(A)') ' FVMBAG1 INPUT IS IGNORED'
72 WRITE(iout,'(A)') ' '
73 ENDIF
74 ELSE
75 IF(igmerg==1.OR.igmerg==2) monvol(k1-1+60)=igmerg
76 IF(cgmerg > zero ) volmon(kk1-1+31)=cgmerg
77 IF(cnmerg > zero ) volmon(kk1-1+34)=cnmerg
78 monvol(k1-1+44)=ilvout
79 IF(ispmd==0)THEN
80 WRITE(iout,1100)
id,monvol(k1-1+60),volmon(kk1-1+31),
81 . volmon(kk1-1+34),monvol(k1-1+44)
82 ENDIF
83 ENDIF
84 GO TO 100
85 ENDIF
86 k1=k1+nimv
87 kk1=kk1+nrvolu
88 ENDDO
89 IF(ifound == 0) THEN
90 IF(ispmd==0)THEN
91 WRITE(iout,'(A)') ' '
92 WRITE(iout,'(A)') ' /FVMBAG/MODIF : '
93 WRITE(iout,
'(A,I10,A)')
' MONITORED VOLUME ID ',
id,
94 . ' DOES NOT EXIST '
95 WRITE(iout,'(A)') ' INPUT IS IGNORED '
96 WRITE(iout,'(A)') ' '
97 ENDIF
98 ENDIF
99 100 ENDDO
100
101 1000 FORMAT(/' NEW FVMBAG PARAMETERS ')
102 1100 FORMAT(/' FVMBAG ID. . . . . . . . . . . . . . .=',i10,
103 . /' GLOBAL MERGING TYPE. . . . . . . . .=',i10,
104 . /' GLOBAL MERGING COEFFICIENT . . . . .=',1pg20.13,
105 . /' NEIGHBORHOOD MERGING COEFFICIENT . .=',1pg20.13,
106 . /' OUTPUT LEVEL . . . . . . . . . . . .=',i10/)
107
108 RETURN