OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fxbodfp.F File Reference
#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "param_c.inc"
#include "parit_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine fxbodfp1 (fxbipm, fxbrpm, fxbnod, fxbmod, fxbdep, fxbvit, fxbacc, a, ar, nme, nmod, itn, fxbelm, fxbsig, elbuf, partsav, x, d, iparg, nfx, nsn, mfext, iad_elem, fr_elem, nsnt, fskyfxb, iadn, iadsky, elbuf_tab)
subroutine fxbodfp2 (fxbipm, fxbrpm, fxbglm, fxbcpm, fxbcps, fxblm, fxbfls, fxbdls, fxbdep, fxbvit, nme, nmod, mvn, mcd, se, sv, fsav, fxbfp, wfext, fxbfc, fxbgrvi, fxbgrvr, nlgrav, igrv, npc, tf, fxbgrp, tfgrav, sensor_tab, nsensor, mfext, agrv, python)
subroutine fxbmajp1 (fxbdep, fxbvit, fxbacc, fxbrpm, dt1, nme, nmod, tmrot)
subroutine fxbmajp2 (mc, rc, mvn, nme, nmod, dt1, fxbcpm, fxbcps, fxblm, fxbrpm)

Function/Subroutine Documentation

◆ fxbmajp1()

subroutine fxbmajp1 ( fxbdep,
fxbvit,
fxbacc,
fxbrpm,
dt1,
integer nme,
integer nmod,
tmrot )

Definition at line 453 of file fxbodfp.F.

455C-----------------------------------------------
456C I m p l i c i t T y p e s
457C-----------------------------------------------
458#include "implicit_f.inc"
459C-----------------------------------------------
460C D u m m y A r g u m e n t s
461C-----------------------------------------------
462 INTEGER NME,NMOD
463 my_real
464 . fxbdep(*),fxbvit(*),fxbacc(*),fxbrpm(*),dt1,tmrot(*)
465C-----------------------------------------------
466C L o c a l V a r i a b l e s
467C-----------------------------------------------
468 INTEGER I,J,II
469 my_real
470 . ddep(nme+nmod)
471C
472C-----------------------------------------------
473C Update displacements and velocities
474C-----------------------------------------------
475 DO i=1,nme+nmod
476 fxbdep(i)=fxbdep(i)+dt1*fxbvit(i)
477 ddep(i)=dt1*fxbvit(i)
478 ENDDO
479C-----------------------------------------------
480C Transfer matrix
481C-----------------------------------------------
482 DO i=1,3
483 DO j=1,3
484 fxbrpm(1+3*(i-1)+j)=fxbrpm(1+3*(i-1)+j)+
485 . ddep(3*(j-1)+i)-ddep(9+i)
486 tmrot(3*(j-1)+i)=fxbrpm(1+3*(i-1)+j)
487 END DO
488 END DO
489C
490 RETURN
my_real
#define my_real
Definition cppsort.cpp:32

◆ fxbmajp2()

subroutine fxbmajp2 ( mc,
rc,
mvn,
integer nme,
integer nmod,
dt1,
fxbcpm,
fxbcps,
fxblm,
fxbrpm )

Definition at line 497 of file fxbodfp.F.

499C-----------------------------------------------
500C I m p l i c i t T y p e s
501C-----------------------------------------------
502#include "implicit_f.inc"
503C-----------------------------------------------
504C D u m m y A r g u m e n t s
505C-----------------------------------------------
506 INTEGER NME, NMOD
507 my_real
508 . mc(*), rc(*), mvn(*), dt1, fxbcpm(*), fxbcps(*), fxblm(*),
509 . fxbrpm(*)
510C-----------------------------------------------
511C L o c a l V a r i a b l e s
512C-----------------------------------------------
513 INTEGER I, II, J
514 my_real
515 . alpha, fac, tmrot(9)
516C
517 DO i=1,3
518 DO j=1,3
519 tmrot(3*(j-1)+i)=fxbrpm(1+3*(i-1)+j)
520 ENDDO
521 ENDDO
522C-----------------------------------------------
523C Mass and stiffness coupling
524C-----------------------------------------------
525 DO i=1,nme*nmod
526 rc(i)=zero
527 mc(i)=zero
528 END DO
529C
530 DO i=1,9
531 ii=nmod*nme*(i-1)
532 DO j=1,nme*nmod
533 mc(j)=mc(j)+fxbrpm(1+i)*fxbcpm(ii+j)
534 rc(j)=rc(j)+tmrot(i)*fxbcps(ii+j)
535 END DO
536 END DO
537C
538 alpha=fxbrpm(13)
539 fac=one+half*dt1*alpha
540 DO i=1,nme
541 DO j=1,nmod
542 ii=nmod*(i-1)+j
543 mvn(ii)=mc(ii)/fxblm(j)/fac
544 ENDDO
545 ENDDO
546C
547 RETURN
alpha
#define alpha
Definition eval.h:35

◆ fxbodfp1()

subroutine fxbodfp1 ( integer, dimension(*) fxbipm,
fxbrpm,
integer, dimension(*) fxbnod,
fxbmod,
fxbdep,
fxbvit,
fxbacc,
a,
ar,
integer nme,
integer nmod,
integer itn,
integer, dimension(*) fxbelm,
fxbsig,
elbuf,
partsav,
x,
d,
integer, dimension(nparg,*) iparg,
integer nfx,
integer nsn,
mfext,
integer, dimension(2,*) iad_elem,
integer, dimension(*) fr_elem,
integer nsnt,
fskyfxb,
integer iadn,
integer, dimension(*) iadsky,
type (elbuf_struct_), dimension (ngroup) elbuf_tab )

Definition at line 34 of file fxbodfp.F.

40C-----------------------------------------------
41C M o d u l e s
42C-----------------------------------------------
43 USE elbufdef_mod
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48C-----------------------------------------------
49C C o m m o n B l o c k s
50C-----------------------------------------------
51#include "com01_c.inc"
52#include "com04_c.inc"
53#include "com08_c.inc"
54#include "param_c.inc"
55#include "parit_c.inc"
56C-----------------------------------------------
57C D u m m y A r g u m e n t s
58C-----------------------------------------------
59 INTEGER FXBIPM(*), FXBNOD(*), NME, NMOD, ITN, FXBELM(*),
60 . IPARG(NPARG,*), NFX, NSN, IAD_ELEM(2,*), FR_ELEM(*),
61 . NSNT, IADN, IADSKY(*)
62 my_real
63 . fxbrpm(*), fxbmod(*), fxbdep(*), fxbvit(*), fxbacc(*),
64 . a(3,*), ar(3,*), fxbsig(*), elbuf(*), partsav(npsav,*),
65 . x(3,*), d(3,*), mfext(*), fskyfxb(nsnt,*)
66 TYPE (ELBUF_STRUCT_), DIMENSION (NGROUP) :: ELBUF_TAB
67C-----------------------------------------------
68C L o c a l V a r i a b l e s
69C-----------------------------------------------
70 INTEGER IFILE, LMOD, ISH, NELS, NELC, NELT, NELP, NELTG, LVSIG,
71 . IRCS, DN, IRCM, I, IAD, II, N, IFAC(NUMNOD), J, JJ
72 my_real
73 . rt(9)
74C=======================================================================
75C
76 CALL fxbmajp1(fxbdep, fxbvit, fxbacc, fxbrpm, dt1 ,
77 . nme , nmod , rt )
78C
79 ifile=fxbipm(29)
80 IF (ifile==0) THEN
81 lmod=fxbipm(3)*6
82 ELSEIF (ifile==1) THEN
83 lmod=nsn*6
84 ENDIF
85 ish=fxbipm(16)
86C
87 nels=fxbipm(21)
88 nelc=fxbipm(22)
89 nelt=fxbipm(34)
90 nelp=fxbipm(35)
91 neltg=fxbipm(23)
92 lvsig=nels*7+nelc*10+nelt*2+nelp*8+neltg*10
93 ircs=fxbipm(31)
94 CALL fxbsgmaj(elbuf, fxbelm, fxbsig, fxbdep, fxbipm,
95 . fxbrpm(15), partsav, rt , itn , iparg ,
96 . nfx , lvsig , ircs ,elbuf_tab)
97 IF (itn==0) THEN
98 dn=fxbipm(3)-fxbipm(18)
99 ircm=fxbipm(30)
100 CALL fxbdepla(fxbdep, fxbrpm, x, d, dn,
101 . nsn, fxbnod, nme, nmod, fxbmod,
102 . ish , ifile , nfx, ircm)
103 ENDIF
104C-----------------------------------------------
105C External modal forces
106C-----------------------------------------------
107 DO i=1,numnod
108 ifac(i)=1
109 ENDDO
110 IF (nspmd>1) THEN
111 DO i=1,nspmd
112 DO j=iad_elem(1,i),iad_elem(1,i+1)-1
113 jj=fr_elem(j)
114 ifac(jj)=ifac(jj)+1
115 ENDDO
116 ENDDO
117 ENDIF
118C
119 IF (iparit==0) THEN
120 DO i=1,12
121 mfext(i)=zero
122 iad=(i-1)*lmod
123 DO ii=1,nsn
124 n=fxbnod(ii)
125 mfext(i)=mfext(i)+(a(1,n)*fxbmod(iad+1)
126 . +a(2,n)*fxbmod(iad+2)
127 . +a(3,n)*fxbmod(iad+3))/ifac(n)
128 iad=iad+6
129 ENDDO
130 ENDDO
131 IF (ish>0) THEN
132 DO i=13,nme
133 mfext(i)=zero
134 iad=(i-1)*lmod
135 DO ii=1,nsn
136 n=fxbnod(ii)
137 mfext(i)=mfext(i)+(ar(1,n)*fxbmod(iad+4)
138 . +ar(2,n)*fxbmod(iad+5)
139 . +ar(3,n)*fxbmod(iad+6))/ifac(n)
140 iad=iad+6
141 ENDDO
142 ENDDO
143 ENDIF
144C
145 IF (nmod>0) THEN
146 DO i=1,nmod
147 mfext(nme+i)=zero
148 iad=(nme+i-1)*lmod
149 DO ii=1,nsn
150 n=fxbnod(ii)
151 mfext(nme+i)=mfext(nme+i)+
152 . (a(1,n)*(fxbrpm(2)*fxbmod(iad+1)+fxbrpm(3)*fxbmod(iad+2)+
153 . fxbrpm(4)*fxbmod(iad+3))
154 . +a(2,n)*(fxbrpm(5)*fxbmod(iad+1)+fxbrpm(6)*fxbmod(iad+2)+
155 . fxbrpm(7)*fxbmod(iad+3))
156 . +a(3,n)*(fxbrpm(8)*fxbmod(iad+1)+fxbrpm(9)*fxbmod(iad+2)+
157 . fxbrpm(10)*fxbmod(iad+3))
158 . +ar(1,n)*(fxbrpm(2)*fxbmod(iad+4)+fxbrpm(3)*fxbmod(iad+5)+
159 . fxbrpm(4)*fxbmod(iad+6))
160 . +ar(2,n)*(fxbrpm(5)*fxbmod(iad+4)+fxbrpm(6)*fxbmod(iad+5)+
161 . fxbrpm(7)*fxbmod(iad+6))
162 . +ar(3,n)*(fxbrpm(8)*fxbmod(iad+4)+fxbrpm(9)*fxbmod(iad+5)+
163 . fxbrpm(10)*fxbmod(iad+6)))/ifac(n)
164 iad=iad+6
165 ENDDO
166 ENDDO
167 ENDIF
168 ELSE
169 DO i=1,nsn
170 fskyfxb(iadn+i,1)=iadsky(i)
171 ENDDO
172C
173 DO i=1,12
174 iad=(i-1)*lmod
175 DO ii=1,nsn
176 n=fxbnod(ii)
177 fskyfxb(iadn+ii,1+i)=a(1,n)*fxbmod(iad+1)
178 . +a(2,n)*fxbmod(iad+2)
179 . +a(3,n)*fxbmod(iad+3)
180 iad=iad+6
181 ENDDO
182 ENDDO
183 IF (ish>0) THEN
184 DO i=13,nme
185 iad=(i-1)*lmod
186 DO ii=1,nsn
187 n=fxbnod(ii)
188 fskyfxb(iadn+ii,1+i)=ar(1,n)*fxbmod(iad+4)
189 . +ar(2,n)*fxbmod(iad+5)
190 . +ar(3,n)*fxbmod(iad+6)
191 iad=iad+6
192 ENDDO
193 ENDDO
194 ENDIF
195C
196 IF (nmod>0) THEN
197 DO i=1,nmod
198 iad=(nme+i-1)*lmod
199 DO ii=1,nsn
200 n=fxbnod(ii)
201 fskyfxb(iadn+ii,1+nme+i)=
202 . a(1,n)*(fxbrpm(2)*fxbmod(iad+1)+fxbrpm(3)*fxbmod(iad+2)+
203 . fxbrpm(4)*fxbmod(iad+3))
204 . +a(2,n)*(fxbrpm(5)*fxbmod(iad+1)+fxbrpm(6)*fxbmod(iad+2)+
205 . fxbrpm(7)*fxbmod(iad+3))
206 . +a(3,n)*(fxbrpm(8)*fxbmod(iad+1)+fxbrpm(9)*fxbmod(iad+2)+
207 . fxbrpm(10)*fxbmod(iad+3))
208 . +ar(1,n)*(fxbrpm(2)*fxbmod(iad+4)+fxbrpm(3)*fxbmod(iad+5)+
209 . fxbrpm(4)*fxbmod(iad+6))
210 . +ar(2,n)*(fxbrpm(5)*fxbmod(iad+4)+fxbrpm(6)*fxbmod(iad+5)+
211 . fxbrpm(7)*fxbmod(iad+6))
212 . +ar(3,n)*(fxbrpm(8)*fxbmod(iad+4)+fxbrpm(9)*fxbmod(iad+5)+
213 . fxbrpm(10)*fxbmod(iad+6))
214 iad=iad+6
215 ENDDO
216 ENDDO
217 ENDIF
218 ENDIF
219C
220 RETURN
fxbdepla
subroutine fxbdepla(fxbdep, fxbrpm, x, d, dn, nsn, fxbnod, nme, nmod, fxbmod, ish, ifile, nfx, ircm)
Definition fxbdepla.F:31
fxbmajp1
subroutine fxbmajp1(fxbdep, fxbvit, fxbacc, fxbrpm, dt1, nme, nmod, tmrot)
Definition fxbodfp.F:455
fxbsgmaj
subroutine fxbsgmaj(elbuf, fxbelm, fxbsig, fxbdep, fxbipm, eiel, partsav, rt, itn, iparg, nfx, lvsig, ircs, elbuf_tab)
Definition fxbsgmaj.F:35

◆ fxbodfp2()

subroutine fxbodfp2 ( integer, dimension(*) fxbipm,
fxbrpm,
fxbglm,
fxbcpm,
fxbcps,
fxblm,
fxbfls,
fxbdls,
fxbdep,
fxbvit,
integer nme,
integer nmod,
mvn,
mcd,
se,
sv,
fsav,
fxbfp,
wfext,
fxbfc,
integer, dimension(*) fxbgrvi,
fxbgrvr,
integer nlgrav,
integer, dimension(nigrv,*) igrv,
integer, dimension(*) npc,
tf,
fxbgrp,
tfgrav,
type (sensor_str_), dimension(nsensor) sensor_tab,
integer, intent(in) nsensor,
mfext,
agrv,
type(python_), intent(inout) python )

Definition at line 233 of file fxbodfp.F.

240C-----------------------------------------------
241C M o d u l e s
242C-----------------------------------------------
243 USE python_funct_mod, ONLY : python_
244 use finter_mixed_mod, ONLY : finter_mixed
245 USE sensor_mod
246C-----------------------------------------------
247C I m p l i c i t T y p e s
248C-----------------------------------------------
249#include "implicit_f.inc"
250C-----------------------------------------------
251C C o m m o n B l o c k s
252C-----------------------------------------------
253#include "com04_c.inc"
254#include "com08_c.inc"
255#include "param_c.inc"
256C-----------------------------------------------
257C D u m m y A r g u m e n t s
258C-----------------------------------------------
259 INTEGER ,INTENT(IN) :: NSENSOR
260 INTEGER FXBIPM(*), NME, NMOD, FXBGRVI(*), NLGRAV, IGRV(NIGRV,*),
261 . NPC(*)
262 my_real
263 . fxbrpm(*), fxbglm(*), fxbcpm(*), fxbcps(*), fxblm(*),
264 . fxbfls(*), fxbdls(*), fxbdep(*), fxbvit(*), mvn(*),
265 . mcd(nme,*), se(*), sv(*), fsav(*), fxbfp(*), wfext,
266 . fxbfc(*), fxbgrvr(*), tf(*), fxbgrp(*), tfgrav,
267 . mfext(*), agrv(lfacgrv,*)
268 TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) :: SENSOR_TAB
269 TYPE(python_), intent(inout) :: python
270C-----------------------------------------------
271C L o c a l V a r i a b l e s
272C-----------------------------------------------
273 INTEGER IBLO, I, IADG, IADGR, IG, NL, IFUNC, ISENS, K, IIM, NMST,
274 . IAD, II, III, IIAD
275 my_real
276 . mc(nme*nmod), rc(nme*nmod), ts, fi, alpha, beta,
277 . mfint(nme+nmod), fdalpha(nmod), enint,
278 . mfgrav(nmod+nme)
279C
280 iblo=fxbipm(28)
281 IF (nmod>0) CALL fxbmajp2(mc, rc, mvn, nme, nmod,
282 . dt1, fxbcpm, fxbcps, fxblm, fxbrpm)
283C-----------------------------------------------
284C Projection of gravity forces
285C-----------------------------------------------
286 DO i=1,nme+nmod
287 mfgrav(i)=zero
288 ENDDO
289C
290 IF (nlgrav>0) THEN
291 iadg=0
292 iadgr=0
293 DO ig=1,nlgrav
294 nl=fxbgrvi(iadg+1)
295 ifunc=igrv(3,nl)
296 isens=0
297 DO k=1,nsensor
298 IF (igrv(6,nl)==sensor_tab(k)%SENS_ID) isens=k ! do it in starter !!!
299 ENDDO
300 IF (isens==0) THEN
301 ts=tt
302 ELSE
303 ts=tt-sensor_tab(isens)%TSTART
304 IF (ts<zero) cycle
305 ENDIF
306 IF (ifunc > 0) THEN
307 fi=agrv(1,nl)*finter_mixed(python,nfunct,ifunc,ts,npc,tf)
308 ELSE
309 fi=agrv(1,nl)
310 ENDIF
311 DO i=1,nme
312 mfgrav(i)=mfgrav(i)+fi*fxbgrvr(iadgr+i)
313 ENDDO
314 iim=0
315 DO i=1,nmod
316 mfgrav(nme+i)=mfgrav(nme+i)
317 . +fxbrpm(2)*fi*fxbgrvr(iadgr+nme+iim+1)
318 . +fxbrpm(3)*fi*fxbgrvr(iadgr+nme+iim+2)
319 . +fxbrpm(4)*fi*fxbgrvr(iadgr+nme+iim+3)
320 . +fxbrpm(5)*fi*fxbgrvr(iadgr+nme+iim+4)
321 . +fxbrpm(6)*fi*fxbgrvr(iadgr+nme+iim+5)
322 . +fxbrpm(7)*fi*fxbgrvr(iadgr+nme+iim+6)
323 . +fxbrpm(8)*fi*fxbgrvr(iadgr+nme+iim+7)
324 . +fxbrpm(9)*fi*fxbgrvr(iadgr+nme+iim+8)
325 . +fxbrpm(10)*fi*fxbgrvr(iadgr+nme+iim+9)
326 iim=iim+9
327 ENDDO
328 iadg=iadg+2+fxbgrvi(iadg+2)
329 iadgr=iadgr+nme+9*nmod
330 ENDDO
331 ENDIF
332C-----------------------------------------------
333C External forces work
334C-----------------------------------------------
335 tfgrav=zero
336 DO i=1,nme+nmod
337 tfgrav=tfgrav+half*(fxbgrp(i)+mfgrav(i))*fxbvit(i)*dt1
338 mfext(i)=mfext(i)+mfgrav(i)
339 wfext=wfext+half*(fxbfp(i)+mfext(i))*fxbvit(i)*dt1
340 fxbgrp(i)=mfgrav(i)
341 fxbfp(i)=mfext(i)
342 ENDDO
343C-----------------------------------------------
344C Internal modal forces
345C-----------------------------------------------
346 nmst=fxbipm(5)
347 alpha=fxbrpm(13)
348 beta=fxbrpm(14)
349 IF (nmod>0.AND.iblo==0) THEN
350 DO i=1,nme
351 mfint(i)=zero
352 iad=nmod*(i-1)
353 DO ii=1,nmod
354 mfint(i)=mfint(i)+rc(iad+ii)*fxbdep(nme+ii)
355 ENDDO
356 ENDDO
357 ELSE
358 DO i=1,nme
359 mfint(i)=zero
360 ENDDO
361 ENDIF
362 IF (nmst>0) THEN
363 DO i=1,nmst
364 mfint(nme+i)=zero
365 fxbfc(i)=zero
366 DO ii=1,i-1
367 iad=(ii-1)*(2*nmod-ii+2)/2+(i-ii+1)
368 mfint(nme+i)=mfint(nme+i)+fxbfls(iad)*fxbdep(nme+ii)
369 fxbfc(i)=fxbfc(i)+fxbfls(iad)*beta*fxbvit(nme+ii)
370 ENDDO
371 DO ii=i,nmod
372 iad=(i-1)*(2*nmod-i+2)/2+(ii-i+1)
373 mfint(nme+i)=mfint(nme+i)+fxbfls(iad)*fxbdep(nme+ii)
374 fxbfc(i)=fxbfc(i)+fxbfls(iad)*beta*fxbvit(nme+ii)
375 ENDDO
376 fdalpha(i)=alpha*fxblm(i)*fxbvit(nme+i)
377 ENDDO
378 ENDIF
379 IF ((nmod-nmst)>0) THEN
380 DO i=1,nmod-nmst
381 iii=nmst+i
382 mfint(nme+iii)=zero
383 fxbfc(iii)=zero
384 IF (nmst>0) THEN
385 DO ii=1,nmst
386 iad=(ii-1)*(2*nmod-ii+2)/2+(iii-ii+1)
387 mfint(nme+iii)=mfint(nme+iii)+fxbfls(iad)*
388 . fxbdep(nme+ii)
389 fxbfc(iii)=fxbfc(iii)+fxbfls(iad)*beta*fxbvit(nme+ii)
390 ENDDO
391 ENDIF
392 mfint(nme+iii)=mfint(nme+iii)+fxbdls(i)*fxbdep(nme+iii)
393 fxbfc(iii)=fxbfc(iii)+fxbdls(i)*beta*fxbvit(nme+iii)
394 fdalpha(iii)=alpha*fxblm(iii)*fxbvit(nme+iii)
395 ENDDO
396 ENDIF
397 enint=zero
398 IF (nmod>0) THEN
399 DO i=1,nmod
400 enint=enint+half*fxbdep(nme+i)*mfint(nme+i)
401 ENDDO
402 ENDIF
403 fxbrpm(11)=enint
404 fsav(1)=enint
405 fsav(3)=wfext
406C-----------------------------------------------
407C Condensation of global unknowns
408C-----------------------------------------------
409 DO i=1,nme
410 se(i)=mfext(i)-mfint(i)
411 ENDDO
412 IF (nmod>0) THEN
413 DO i=1,nmod
414 sv(i)=mfext(nme+i)-mfint(nme+i)-fxbfc(i)-fdalpha(i)
415 ENDDO
416 IF (iblo==0) THEN
417 DO i=1,nme
418 iad=nmod*(i-1)
419 DO ii=1,nmod
420 se(i)=se(i)-mvn(iad+ii)*sv(ii)
421 ENDDO
422 ENDDO
423 ENDIF
424 ENDIF
425 IF (iblo==1) RETURN
426 iii=0
427 DO i=1,nme
428 DO ii=i,nme
429 iii=iii+1
430 mcd(i,ii)=fxbglm(iii)
431 IF (i/=ii) mcd(ii,i)=mcd(i,ii)
432 ENDDO
433 ENDDO
434 IF (nmod>0) THEN
435 DO i=1,nme
436 iad=nmod*(i-1)
437 DO ii=1,nme
438 iiad=nmod*(ii-1)
439 DO iii=1,nmod
440 mcd(i,ii)=mcd(i,ii)-mvn(iad+iii)*mc(iiad+iii)
441 ENDDO
442 ENDDO
443 ENDDO
444 ENDIF
445C
446 RETURN
fxbmajp2
subroutine fxbmajp2(mc, rc, mvn, nme, nmod, dt1, fxbcpm, fxbcps, fxblm, fxbrpm)
Definition fxbodfp.F:499
nl
character *2 function nl()
Definition message.F:2354