38
39
40
41#include "implicit_f.inc"
42#include "comlock.inc"
43
44
45
46#include "mvsiz_p.inc"
47
48
49
50 INTEGER I_MEM, NSN,IGAP
51 INTEGER J_STOK,MULNSN,NOINT,INACTI,ESHIFT,II_STOK
52 INTEGER IRECT(4,*),CAND_N(*),CAND_E(*)
53 INTEGER PROV_N(MVSIZ),PROV_E(MVSIZ)
54
56 . xloc(3,*), gap_s(*),
57 . marge, gap, gapmin, gapmax, curv_max(*),
58 . depth, xm0(3,*), nod_normal(3,*)
59 my_real ,
INTENT(IN) :: dgapload,drad
60
61
62
63 INTEGER I,K_STOK,I_STOK,N,NE,J
64 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ)
65
67 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
68 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
69 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
70 . xi(mvsiz), yi(mvsiz), zi(mvsiz),
71 . nnx1(mvsiz), nnx2(mvsiz), nnx3(mvsiz), nnx4(mvsiz),
72 . nny1(mvsiz), nny2(mvsiz), nny3(mvsiz), nny4(mvsiz),
73 . nnz1(mvsiz), nnz2(mvsiz), nnz3(mvsiz), nnz4(mvsiz),
74 . pene(mvsiz), gapv(mvsiz)
75
76 CALL i21cor3t(j_stok ,xloc ,irect ,prov_e ,
77 1 prov_n ,igap ,gap ,x1 ,x2 ,
78 2 x3 ,x4 ,y1 ,y2 ,y3 ,
79 3 y4 ,z1 ,z2 ,z3 ,z4 ,
80 4 xi ,yi ,zi ,ix1 ,ix2 ,
81 5 ix3 ,ix4 ,nsn ,gap_s ,gapv ,
82 6 gapmax ,gapmin,curv_max,depth ,xm0 ,
83 7 drad ,dgapload)
84
85 CALL i7pen3(j_stok ,marge ,x1 ,x2 ,x3 ,
86 . x4 ,y1 ,y2 ,y3 ,y4 ,
87 . z1 ,z2
88 . yi ,zi ,pene
89 .
90
91 k_stok = 0
92 DO i=1,j_stok
93 IF(pene(i)/=zero) k_stok = k_stok + 1
94 ENDDO
95 IF(k_stok==0)RETURN
96
97#include "lockon.inc"
98 i_stok = ii_stok
99 IF(i_stok+k_stok>mulnsn) THEN
100 i_mem = 2
101#include "lockoff.inc"
102 RETURN
103 ENDIF
104 ii_stok = i_stok + k_stok
105#include "lockoff.inc"
106
107 DO i=1,j_stok
108 IF(pene(i)/=zero)THEN
109 i_stok = i_stok + 1
110 cand_n(i_stok) = prov_n(i)
111 cand_e(i_stok) = prov_e(i)+eshift
112 ENDIF
113 ENDDO
114
115 RETURN
subroutine i21cor3t(x, irect, nsv, cand_e, cand_n, gapv, igap, gap, gap_s, gapmin, gapmax, xm0, depth, drad, ix1, ix2, ix3, ix4, nsvg, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, dgapload)
subroutine i7pen3(marge, gapv, n1, n2, n3, pene, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, last)