43#include "implicit_f.inc"
54 INTEGER,
INTENT(IN) :: NEL
55 my_real,
DIMENSION(MVSIZ) ,
INTENT(OUT) :: LM_MAX
56 my_real,
DIMENSION(MVSIZ) ,
INTENT(IN) ::
57 . X1, X2, X3, X4, X5, X6, , X8,
58 . Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,
59 . Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8
65 . RX, RY, RZ, SX, SY, SZ, TX, TY, TZ,NORM
68 rx = x1(i)+x2(i)+x5(i)+x6(i)-x3(i)-x4(i)-x7(i)-x8(i)
69 ry = y1(i)+y2(i)+y5(i)+y6(i)-y3(i)-y4(i)-y7(i)-y8(i)
70 rz = z1(i)+z2(i)+z5(i)+z6(i)-z3(i)-z4(i)-z7(i)-z8(i)
71 sx = x5(i)+x6(i)+x7(i)+x8(i
72 sy = y5(i)+y6(i)+y7(i)+y8(i)-y1(i)-y2(i)-y3(i)-y4(i)
73 sz = z5(i)+z6(i)+z7(i)+z8(i)-z1(i)-z2(i)-z3(i)-z4(i)
74 tx = x3(i)+x2(i)+x7(i)+x6(i)-x1(i)-x4(i)-x5(i)-x8(i)
75 ty = y3(i)+y2(i)+y7(i)+y6(i)-y1(i)-y4(i)-y5(i)-y8(i)
76 tz = z3(i)+z2(i)+z7(i)+z6(i)-z1(i)-z4(i)-z5(i)-z8(i)
77 norm = rx*rx + ry*ry + rz*rz
79 norm = sx*sx + sy*sy + sz*sz
80 lm_max(i) =
max(lm_max(i),norm)
81 norm = tx*tx + ty*ty + tz*tz
82 lm_max(i) =
max(lm_max(i),norm)
85 lm_max(i) = fourth*sqrt(lm_max(i))
subroutine s20forc3(timers, output, elbuf_tab, ng, pm, geo, ixs, x, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, nloc_dmg, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, ixs20, iads20, nel, fx, fy, fz, volnp, rx, ry, rz, sx, sy, sz, tx, ty, tz, stig, stin, ul, xx, yy, zz, vx, vy, vz, vdxx, vdyy, vdzz, dnidr, dnids, dnidt, px, py, pz, ipm, istrain, temp, fthe, fthesky, iexpan, gresav, grth, igrth, table, igeo, voln, condn, condnsky, itask, ioutprt, mat_elem, h3d_strain, dt, idel7nok, snpc, stf, sbufmat, svis, idtmins, nsvois, iresp, idel7ng, idtmin, maxfunc, imon_mat, userl_avail, glob_therm, impl_s, idyna, sensors)
subroutine s8eforc3(timers, output, elbuf_tab, ng, pm, geo, ixs, x, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, nloc_dmg, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, icp, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, nvc, ipm, itask, istrain, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mssa, dmels, table, igeo, xdp, voln, condn, condnsky, d, sensors, ioutprt, mat_elem, h3d_strain, dt, snpc, stf, sbufmat, svis, nsvois, idtmins, iresp, maxfunc, userl_avail, glob_therm, impl_s, idyna)
subroutine sdlenmax(lm_max, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, nel)
subroutine szforc3(timers, elbuf_tab, ng, pm, geo, ixs, x, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, icp, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, fskym, msnf, nvc, ipm, itask, qmv, istrain, imatvis, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mssa, dmels, table, igeo, xdp, condn, condnsky, d, tagprt_sms, sensors, ioutprt, nale, nloc_dmg, mat_elem, h3d_strain, dt, output, nsvois, sz_bufvois, snpc, stf, sbufmat, svis, iresp, maxfunc, userl_avail, glob_therm, impl_s, idyna)