107
108
109
110
111
112
113
114#include "implicit_f.inc"
115
116
117
118#include "com04_c.inc"
119#include "param_c.inc"
120#include "scr16_c.inc"
121#include "scr17_c.inc"
122#include "task_c.inc"
123#include "units_c.inc"
124#include "com01_c.inc"
125
126
127
128 INTEGER NPBY(NNPBY,*),NOM_SECT(LNOPT1,*),NOM_WALL(LNOPT1,*)
130
131
132
133 INTEGER I,N
134
135 IF(nspmd > 1)THEN
137 IF (ispmd/=0) THEN
138 DO n=1,nsect+nrbody+nrwall
139 DO i=1,6
140 fopt(i,n) = zero
141 ENDDO
142 ENDDO
143 RETURN
144 ENDIF
145 END IF
146
147
148 WRITE(iugeo,'(A)')'/RBODY /VECTOR /FORCE'
149 WRITE(iugeo,'(A)')
150 . 'Rigid bodies : forces and moments'
151 IF (outyy_fmt==2) THEN
152 WRITE(iugeo,'(A)')'#FORMAT: (I8,1P3E16.9,/,8X,,1P3E16.9) '
153 WRITE(iugeo,'(2A,/,2A)')'#RBODYID',
154 . ' FX FY FZ',
155 . '# ',
156 . ' MX MY MZ'
157 ELSE
158 WRITE(iugeo,'(A)')'#FORMAT: (I10,1P3E20.13,/,10X,1P3E20.13) '
159 WRITE(iugeo,'(2A,/,2A)')'# RBODY_ID',
160 . ' FX FY FZ',
161 . '# ',
162 . ' MX MY MZ'
163
164 END IF
165
166
167 IF (outyy_fmt==2) THEN
168 WRITE(iugeo,'(I8,1P3E16.9,/,8X,1P3E16.9)')
169 . (npby(6,n),(fopt(i,nsect+n),i=1,3),
170 . (fopt(i,nsect+n),i=4,6),n=1,nrbody)
171 ELSE
172 WRITE(iugeo,'(I10,1P3E20.13,/,10X,1P3E20.13)')
173 . (npby(6,n),(fopt(i,nsect+n),i=1,3),
174 . (fopt(i,nsect+n),i=4,6),n=1,nrbody)
175 END IF
176
177
178
179
180 WRITE(iugeo,'(A)')'/SECTION /VECTOR /FORCE'
181 WRITE(iugeo,'(a)')
182 . 'sections forces and moments'
183 IF (OUTYY_FMT==2) THEN
184 WRITE(IUGEO,'(a)')'#FORMAT: (I8,1P3E16.9,/,16X,,1P3E16.9) '
185 WRITE(iugeo,'(2A,/,2A)')'#SECT_ID',
186 . ' FX FY FZ',
187 . '# ',
188 . ' MX MY MZ'
189 ELSE
190 WRITE(iugeo,'(A)')'#FORMAT: (I10,1P3E20.13,/,10X,1P3E20.13) '
191 WRITE(iugeo,'(2A,/,2A)')'# SECT_ID',
192 . ' FX FY FZ',
193 . '# ',
194 . ' MX MY MZ'
195
196 ENDIF
197
198
199 IF (outyy_fmt==2) THEN
200 WRITE(iugeo,'(I8,1P3E16.9,/,8X,1P3E16.9)')
201 . (nom_sect(1,n),(fopt(i,n),i=1,3),
202 . (fopt(i,n),i=4,6),n=1,nsect)
203 ELSE
204 WRITE(iugeo,'(I10,1P3E20.13,/,10X,1P3E20.13)')
205 . (nom_sect(1,n),(fopt(i,n),i=1,3),
206 . (fopt(i,n),i=4,6),n=1,nsect)
207 ENDIF
208
209 WRITE(iugeo,'(A)')'/RWALL /VECTOR /FORCE'
210 WRITE(iugeo,'(A)')
211 . 'Rigid wall forces'
212 IF (outyy_fmt==2) THEN
213 WRITE(iugeo,'(A)')'#FORMAT: (2I8,1P3E16.9) '
214 WRITE(iugeo,'(2A)')'#RWAL_ID',
215 . ' FX FY FZ'
216 ELSE
217 WRITE(iugeo,'(A)')'#FORMAT: (2I10,1P3E20.13) '
218 WRITE(iugeo,'(2A)')'# RWAL_ID',
219 . ' FX FY FZ'
220
221 ENDIF
222
223 IF (outyy_fmt==2) THEN
224 WRITE(iugeo,'(I8,1P3E16.9)')
225 . (nom_wall(1,n),(fopt(i,nsect+nrbody+n),i=1,3),n=1,nrwall)
226 ELSE
227 WRITE(iugeo,'(I10,1P3E20.13)')
228 . (nom_wall(1,n),(fopt(i,nsect+nrbody+n),i=1,3),n=1,nrwall)
229 ENDIF
230
231
232 RETURN
subroutine spmd_glob_dsum9(v, len)