Lines Matching refs:i
115 do i=1,nshl
116 u1 = u1 + shp(i,intp)*yl(:,i,2)
117 u2 = u2 + shp(i,intp)*yl(:,i,3)
118 u3 = u3 + shp(i,intp)*yl(:,i,4)
120 S11 = S11 + shg(:,i,1)*yl(:,i,2)
121 S22 = S22 + shg(:,i,2)*yl(:,i,3)
122 S33 = S33 + shg(:,i,3)*yl(:,i,4)
124 S12 = S12 + shg(:,i,2)*yl(:,i,2)
125 & +shg(:,i,1)*yl(:,i,3)
126 S13 = S13 + shg(:,i,3)*yl(:,i,2)
127 & +shg(:,i,1)*yl(:,i,4)
128 S23 = S23 + shg(:,i,3)*yl(:,i,3)
129 & +shg(:,i,2)*yl(:,i,4)
157 do i = 1, 16 ! Add contribution of each quad. point
158 fresl(:,i) = fresl(:,i) + fresli(:,i)
303 do i = 1, nenl
304 fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1) ! hat{u1}
305 fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2) ! hat{u2}
306 fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3) ! hat{u3}
308 fresli(:,4) = fresli(:,4) + shp(i,intp)*hfresl(:,i,1)*
309 & hfresl(:,i,1) ! hat{u1}*hat{u1}
310 fresli(:,5) = fresli(:,5) + shp(i,intp)*hfresl(:,i,2)*
311 & hfresl(:,i,2) ! hat{u2}*hat{u2}
312 fresli(:,6) = fresli(:,6) + shp(i,intp)*hfresl(:,i,3)*
313 & hfresl(:,i,3) ! hat{u3}*hat{u3}
314 fresli(:,7) = fresli(:,7) + shp(i,intp)*hfresl(:,i,1)*
315 & hfresl(:,i,2) ! hat{u1}*hat{u2}
316 fresli(:,8) = fresli(:,8) + shp(i,intp)*hfresl(:,i,1)*
317 & hfresl(:,i,3) ! hat{u1}*hat{u3}
318 fresli(:,9) = fresli(:,9) + shp(i,intp)*hfresl(:,i,2)*
319 & hfresl(:,i,3) ! hat{u2}*hat{u3}
321 fresli(:,10) = fresli(:,10) + shp(i,intp)*hfresl(:,i,10) ! hat{S11}
322 fresli(:,11) = fresli(:,11) + shp(i,intp)*hfresl(:,i,11) ! hat{S22}
323 fresli(:,12) = fresli(:,12) + shp(i,intp)*hfresl(:,i,12) ! hat{S33}
324 fresli(:,13) = fresli(:,13) + shp(i,intp)*hfresl(:,i,13) ! hat{S12}
325 fresli(:,14) = fresli(:,14) + shp(i,intp)*hfresl(:,i,14) ! hat{S13}
326 fresli(:,15) = fresli(:,15) + shp(i,intp)*hfresl(:,i,15) ! hat{S23}
328 S11p = S11p + shg(:,i,1)*hfresl(:,i,1)
329 S22p = S22p + shg(:,i,2)*hfresl(:,i,2)
330 S33p = S33p + shg(:,i,3)*hfresl(:,i,3)
332 S12p = S12p + shg(:,i,2)*hfresl(:,i,1) +
333 & shg(:,i,1)*hfresl(:,i,2)
334 S13p = S13p + shg(:,i,3)*hfresl(:,i,1) +
335 & shg(:,i,1)*hfresl(:,i,3)
336 S23p = S23p + shg(:,i,3)*hfresl(:,i,2) +
337 & shg(:,i,2)*hfresl(:,i,3)
359 do i = 1, nenl
360 fresli(:,16) =fresli(:,16)+shp(i,intp)*SS11(:,i)! |hat{S}|*hat{S11}
361 fresli(:,17) =fresli(:,17)+shp(i,intp)*SS22(:,i)! |hat{S}|*hat{S22}
362 fresli(:,18) =fresli(:,18)+shp(i,intp)*SS33(:,i)! |hat{S}|*hat{S33}
363 fresli(:,19) =fresli(:,19)+shp(i,intp)*SS12(:,i)! |hat{S}|*hat{S12}
364 fresli(:,20) =fresli(:,20)+shp(i,intp)*SS13(:,i)! |hat{S}|*hat{S13}
365 fresli(:,21) =fresli(:,21)+shp(i,intp)*SS23(:,i)! |hat{S}|*hat{S23}
383 do i = 1, 22
384 fresl(:,i) = fresl(:,i) + fresli(:,i)
515 do i=1,nshl
516 u1 = u1 + shp(i,intp)*yl(:,i,2)
517 u2 = u2 + shp(i,intp)*yl(:,i,3)
518 u3 = u3 + shp(i,intp)*yl(:,i,4)
520 S11 = S11 + shg(:,i,1)*yl(:,i,2)
521 S22 = S22 + shg(:,i,2)*yl(:,i,3)
522 S33 = S33 + shg(:,i,3)*yl(:,i,4)
524 S12 = S12 + shg(:,i,2)*yl(:,i,2)
525 & +shg(:,i,1)*yl(:,i,3)
526 S13 = S13 + shg(:,i,3)*yl(:,i,2)
527 & +shg(:,i,1)*yl(:,i,4)
528 S23 = S23 + shg(:,i,3)*yl(:,i,3)
529 & +shg(:,i,2)*yl(:,i,4)
572 do i = 1, 22 ! Add contribution of each quad. point
573 fresl(:,i) = fresl(:,i) + fresli(:,i)
716 do i = 1, nenl
717 fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1) ! hat{u1}
718 fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2) ! hat{u2}
719 fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3) ! hat{u3}
721 fresli(:,4) = fresli(:,4) + shp(i,intp)*hfresl(:,i,4) ! hat{u1*u1}
722 fresli(:,5) = fresli(:,5) + shp(i,intp)*hfresl(:,i,5) ! hat{u2*u2}
723 fresli(:,6) = fresli(:,6) + shp(i,intp)*hfresl(:,i,6) ! hat{u3*u3}
724 fresli(:,7) = fresli(:,7) + shp(i,intp)*hfresl(:,i,7) ! hat{u1*u2}
725 fresli(:,8) = fresli(:,8) + shp(i,intp)*hfresl(:,i,8) ! hat{u1*u3}
726 fresli(:,9) = fresli(:,9) + shp(i,intp)*hfresl(:,i,9) ! hat{u2*u3}
728 fresli(:,10) = fresli(:,10) + shp(i,intp)*hfresl(:,i,10) ! hat{S11}
729 fresli(:,11) = fresli(:,11) + shp(i,intp)*hfresl(:,i,11) ! hat{S22}
730 fresli(:,12) = fresli(:,12) + shp(i,intp)*hfresl(:,i,12) ! hat{S33}
731 fresli(:,13) = fresli(:,13) + shp(i,intp)*hfresl(:,i,13) ! hat{S12}
732 fresli(:,14) = fresli(:,14) + shp(i,intp)*hfresl(:,i,14) ! hat{S13}
733 fresli(:,15) = fresli(:,15) + shp(i,intp)*hfresl(:,i,15) ! hat{S23}
735 fresli(:,16) = fresli(:,16) +shp(i,intp)*hfresl(:,i,17)! hat{S11*|S|}
736 fresli(:,17) = fresli(:,17) +shp(i,intp)*hfresl(:,i,18)! hat{S22*|S|}
737 fresli(:,18) = fresli(:,18) +shp(i,intp)*hfresl(:,i,19)! hat{S33*|S|}
738 fresli(:,19) = fresli(:,19) +shp(i,intp)*hfresl(:,i,20)! hat{S12*|S|}
739 fresli(:,20) = fresli(:,20) +shp(i,intp)*hfresl(:,i,21)! hat{S13*|S|}
740 fresli(:,21) = fresli(:,21) +shp(i,intp)*hfresl(:,i,22)! hat{S23*|S|}
760 do i = 1, 22
761 fresl(:,i) = fresl(:,i) + fresli(:,i)
894 do i=1,nshl
895 u1 = u1 + shp(i,intp)*yl(:,i,2)
896 u2 = u2 + shp(i,intp)*yl(:,i,3)
897 u3 = u3 + shp(i,intp)*yl(:,i,4)
899 S11 = S11 + shg(:,i,1)*yl(:,i,2)
900 S22 = S22 + shg(:,i,2)*yl(:,i,3)
901 S33 = S33 + shg(:,i,3)*yl(:,i,4)
903 S12 = S12 + shg(:,i,2)*yl(:,i,2)
904 & +shg(:,i,1)*yl(:,i,3)
905 S13 = S13 + shg(:,i,3)*yl(:,i,2)
906 & +shg(:,i,1)*yl(:,i,4)
907 S23 = S23 + shg(:,i,3)*yl(:,i,3)
908 & +shg(:,i,2)*yl(:,i,4)
924 do i = 1, nenl ! Loop over element nodes
926 fresli(:,i,1) = WdetJ * u1 * shp(i,intp) !G * u1 * WdetJ
927 fresli(:,i,2) = WdetJ * u2 * shp(i,intp) !G * u2 * WdetJ
928 fresli(:,i,3) = WdetJ * u3 * shp(i,intp) !G * u3 * WdetJ
930 fresli(:,i,4) = WdetJ * u1 * u1 * shp(i,intp) ! G*u1*u1*WdetJ
931 fresli(:,i,5) = WdetJ * u2 * u2 * shp(i,intp) ! G*u2*u2*WdetJ
932 fresli(:,i,6) = WdetJ * u3 * u3 * shp(i,intp) ! G*u3*u3*WdetJ
933 fresli(:,i,7) = WdetJ * u1 * u2 * shp(i,intp) ! G*u1*u2*WdetJ
934 fresli(:,i,8) = WdetJ * u1 * u3 * shp(i,intp) ! G*u1*u3*WdetJ
935 fresli(:,i,9) = WdetJ * u2 * u3 * shp(i,intp) ! G*u2*u3*WdetJ
937 fresli(:,i,10) = S11 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
938 fresli(:,i,11) = S22 * shp(i,intp) * WdetJ ! G*S_{2,2}*WdetJ
939 fresli(:,i,12) = S33 * shp(i,intp) * WdetJ ! G*S_{3,3}*WdetJ
940 fresli(:,i,13) = S12 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
941 fresli(:,i,14) = S13 * shp(i,intp) * WdetJ ! G*S_{1,3}*WdetJ
942 fresli(:,i,15) = S23 * shp(i,intp) * WdetJ ! G*S_{2,3}*WdetJ
944 fresli(:,i,22) = WdetJ*shp(i,intp) !Integral of filter kernel, G,
949 fresli(:,i,16) = S11 * strnrmi * shp(i,intp) * WdetJ
950 fresli(:,i,17) = S22 * strnrmi * shp(i,intp) * WdetJ
951 fresli(:,i,18) = S33 * strnrmi * shp(i,intp) * WdetJ
952 fresli(:,i,19) = S12 * strnrmi * shp(i,intp) * WdetJ
953 fresli(:,i,20) = S13 * strnrmi * shp(i,intp) * WdetJ
954 fresli(:,i,21) = S23 * strnrmi * shp(i,intp) * WdetJ
958 fresl(:,i,j) = fresl(:,i,j) + fresli(:,i,j)
1103 i = iper(j)
1104 if (i .ne. j) then
1105 stress(i,:) = stress(i,:) + stress(j,:)
1106 diss(i,:) = diss(i,:) + diss(j,:)
1107 vol(i) = vol(i) + vol(j)
1112 i = iper(j)
1113 if (i .ne. j) then
1114 stress(j,:) = stress(i,:)
1115 diss(j,:) = diss(i,:)
1116 vol(j) = vol(i)
1125 do i = 1, 9
1126 stress(:,i) = stress(:,i)*vol(:)
1128 do i = 1, 3
1129 diss(:,i) = diss(:,i)*vol(:)
1134 do i = 1, 9
1135 xave(:,i) = stress(:,i)
1143 i = iper(j)
1144 if (i .ne. j) then
1177 do i = 1,nshg
1178 xaveg(ifath(i),:) = xaveg(ifath(i),:) + xave(i,:)
1215 do i = 1, nfath
1216 write(376,*)xavegt(i,1),xavegt(i,2),xavegt(i,3)
1217 write(377,*)xavegt(i,4),xavegt(i,5),xavegt(i,6)
1218 write(378,*)xavegt(i,7),xavegt(i,8),xavegt(i,9)
1219 write(379,*)xavegt(i,10),xavegt(i,11),xavegt(i,12)
1348 do i=1,nshl
1349 u1 = u1 + shp(i,intp)*yl(:,i,2)
1350 u2 = u2 + shp(i,intp)*yl(:,i,3)
1351 u3 = u3 + shp(i,intp)*yl(:,i,4)
1353 S11 = S11 + shg(:,i,1)*yl(:,i,2)
1354 S22 = S22 + shg(:,i,2)*yl(:,i,3)
1355 S33 = S33 + shg(:,i,3)*yl(:,i,4)
1357 S12 = S12 + shg(:,i,2)*yl(:,i,2)
1358 & +shg(:,i,1)*yl(:,i,3)
1359 S13 = S13 + shg(:,i,3)*yl(:,i,2)
1360 & +shg(:,i,1)*yl(:,i,4)
1361 S23 = S23 + shg(:,i,3)*yl(:,i,3)
1362 & +shg(:,i,2)*yl(:,i,4)
1378 do i = 1, nenl ! Loop over element nodes
1380 fresli(:,i,1) = WdetJ * u1 * shp(i,intp) !G * u1 * WdetJ
1381 fresli(:,i,2) = WdetJ * u2 * shp(i,intp) !G * u2 * WdetJ
1382 fresli(:,i,3) = WdetJ * u3 * shp(i,intp) !G * u3 * WdetJ
1384 fresli(:,i,4) = WdetJ * u1 * u1 * shp(i,intp) ! G*u1*u1*WdetJ
1385 fresli(:,i,5) = WdetJ * u2 * u2 * shp(i,intp) ! G*u2*u2*WdetJ
1386 fresli(:,i,6) = WdetJ * u3 * u3 * shp(i,intp) ! G*u3*u3*WdetJ
1387 fresli(:,i,7) = WdetJ * u1 * u2 * shp(i,intp) ! G*u1*u2*WdetJ
1388 fresli(:,i,8) = WdetJ * u1 * u3 * shp(i,intp) ! G*u1*u3*WdetJ
1389 fresli(:,i,9) = WdetJ * u2 * u3 * shp(i,intp) ! G*u2*u3*WdetJ
1391 fresli(:,i,10) = S11 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
1392 fresli(:,i,11) = S22 * shp(i,intp) * WdetJ ! G*S_{2,2}*WdetJ
1393 fresli(:,i,12) = S33 * shp(i,intp) * WdetJ ! G*S_{3,3}*WdetJ
1394 fresli(:,i,13) = S12 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
1395 fresli(:,i,14) = S13 * shp(i,intp) * WdetJ ! G*S_{1,3}*WdetJ
1396 fresli(:,i,15) = S23 * shp(i,intp) * WdetJ ! G*S_{2,3}*WdetJ
1399 fresli(:,i,16) = strnrmi*strnrmi*strnrmi*shp(i,intp)*WdetJ
1401 fresli(:,i,22) = WdetJ*shp(i,intp) !Integral of filter kernel, G,
1407 fresli(:,i,17) = S22 * strnrmi * shp(i,intp) * WdetJ
1408 fresli(:,i,18) = S33 * strnrmi * shp(i,intp) * WdetJ
1409 fresli(:,i,19) = S12 * strnrmi * shp(i,intp) * WdetJ
1410 fresli(:,i,20) = S13 * strnrmi * shp(i,intp) * WdetJ
1411 fresli(:,i,21) = S23 * strnrmi * shp(i,intp) * WdetJ
1415 fresl(:,i,j) = fresl(:,i,j) + fresli(:,i,j)
1784 i = iper(j)
1785 if (i .ne. j) then
1786 fres(i,:) = fres(i,:) + fres(j,:)
1790 i = iper(j)
1791 if (i .ne. j) then
1792 fres(j,:) = fres(i,:)
1872 i = iper(j)
1873 if (i .ne. j) then
1938 do i = 1,nshg
1939 xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
1940 xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
1942 ynude(ifath(i),1) = ynude(ifath(i),1) + xfac(i,1)
1943 ynude(ifath(i),2) = ynude(ifath(i),2) + xfac(i,2)
1944 ynude(ifath(i),3) = ynude(ifath(i),3) + xfac(i,3)
1945 ynude(ifath(i),4) = ynude(ifath(i),4) + xfac(i,4)
1946 ynude(ifath(i),5) = ynude(ifath(i),5) + xfac(i,5)
1948 xm(ifath(i),1) = xm(ifath(i),1) + xmij(i,1)
1949 xm(ifath(i),2) = xm(ifath(i),2) + xmij(i,2)
1950 xm(ifath(i),3) = xm(ifath(i),3) + xmij(i,3)
1951 xm(ifath(i),4) = xm(ifath(i),4) + xmij(i,4)
1952 xm(ifath(i),5) = xm(ifath(i),5) + xmij(i,5)
1953 xm(ifath(i),6) = xm(ifath(i),6) + xmij(i,6)
1955 xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
1956 xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
1957 xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
1958 xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
1959 xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
1960 xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)
1962 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
1963 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
1964 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
1965 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
1966 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
1967 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
1969 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
1970 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
1971 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
1972 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
1973 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
1974 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
1976 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
1977 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
1978 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
1980 snorm(ifath(i)) = snorm(ifath(i)) + strnrm(i)
1998 do i = 1, nfath
1999 ynuder(i,6) = ( ynuder(i,4) - fwr1*ynuder(i,2) ) /
2000 & ( two*ynuder(i,5) - four*fwr1*ynuder(i,3)
2001 & + two*fwr1*fwr1*ynuder(i,1) )
2100 do i = 1, nfath
2101 ynude(i,6) = ( ynude(i,4) - fwr1*ynude(i,2) ) /
2102 & ( two*ynude(i,5) - four*fwr1*ynude(i,3)
2103 & + fwr1*fwr1*ynude(i,1) )
2154 do i = 1, nfath
2155 write(365,*)xmodcomp(i,1),xmodcomp(i,2),xmodcomp(i,3),
2156 & xmodcomp(i,4),xmodcomp(i,5)
2158 write(366,*)xmcomp(i,1),xmcomp(i,2),xmcomp(i,3)
2159 write(367,*)xmcomp(i,4),xmcomp(i,5),xmcomp(i,6)
2161 write(368,*)xlcomp(i,1),xlcomp(i,2),xlcomp(i,3)
2162 write(369,*)xlcomp(i,4),xlcomp(i,5),xlcomp(i,6)
2164 write(370,*)xl1comp(i,1),xl1comp(i,2),xl1comp(i,3)
2165 write(371,*)xl1comp(i,4),xl1comp(i,5),xl1comp(i,6)
2167 write(372,*)xl2comp(i,1),xl2comp(i,2),xl2comp(i,3)
2168 write(373,*)xl2comp(i,4),xl2comp(i,5),xl2comp(i,6)
2170 write(374,*)xnd(i,1),xnd(i,2),scomp(i)
2171 write(375,*)ucomp(i,1),ucomp(i,2),ucomp(i,3)
2287 do i = 1,numnp
2288 xnude(ifath(i),2) = xnude(ifath(i),2) + strnrm(i)
2304 do i = 1,numnp
2305 xnude(ifath(i),1) = xnude(ifath(i),1) + x(i,1)
2306 xnude(ifath(i),2) = xnude(ifath(i),2) + x(i,2)
2317 do i=1,nfath ! cdelsq * |S|
2318 write(444,*) xnuder(i,1),xnuder(i,2),xnutf(i)
2444 i = iper(j)
2445 if (i .ne. j) then
2446 hfres(i,:) = hfres(i,:) + hfres(j,:)
2450 i = iper(j)
2451 if (i .ne. j) then
2452 hfres(j,:) = hfres(i,:)
2503 i = iper(j)
2504 if (i .ne. j) then
2505 fres(i,:) = fres(i,:) + fres(j,:)
2510 i = iper(j)
2511 if (i .ne. j) then
2512 fres(j,:) = fres(i,:)
2613 i = iper(j)
2614 if (i .ne. j) then
2679 do i = 1,nshg
2680 xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
2681 xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
2683 ynude(ifath(i),1) = ynude(ifath(i),1) + xfac(i,1)
2684 ynude(ifath(i),2) = ynude(ifath(i),2) + xfac(i,2)
2685 ynude(ifath(i),3) = ynude(ifath(i),3) + xfac(i,3)
2686 ynude(ifath(i),4) = ynude(ifath(i),4) + xfac(i,4)
2687 ynude(ifath(i),5) = ynude(ifath(i),5) + xfac(i,5)
2689 xm(ifath(i),1) = xm(ifath(i),1) + xmij(i,1)
2690 xm(ifath(i),2) = xm(ifath(i),2) + xmij(i,2)
2691 xm(ifath(i),3) = xm(ifath(i),3) + xmij(i,3)
2692 xm(ifath(i),4) = xm(ifath(i),4) + xmij(i,4)
2693 xm(ifath(i),5) = xm(ifath(i),5) + xmij(i,5)
2694 xm(ifath(i),6) = xm(ifath(i),6) + xmij(i,6)
2696 xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
2697 xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
2698 xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
2699 xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
2700 xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
2701 xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)
2703 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
2704 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
2705 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
2706 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
2707 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
2708 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
2710 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
2711 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
2712 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
2713 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
2714 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
2715 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
2717 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
2718 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
2719 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
2721 snorm(ifath(i)) = snorm(ifath(i)) + strnrm(i)
2739 do i = 1, nfath
2740 ynuder(i,6) = ( ynuder(i,4) - fwr1*ynuder(i,2) ) /
2741 & ( two*ynuder(i,5) - four*fwr1*ynuder(i,3)
2742 & + two*fwr1*fwr1*ynuder(i,1) )
2841 do i = 1, nfath
2842 ynude(i,6) = ( ynude(i,4) - fwr1*ynude(i,2) ) /
2843 & ( two*ynude(i,5) - four*fwr1*ynude(i,3)
2844 & + fwr1*fwr1*ynude(i,1) )
2893 do i = 1, nfath
2894 write(365,*)xmodcomp(i,1),xmodcomp(i,2),xmodcomp(i,3),
2895 & xmodcomp(i,4),xmodcomp(i,5)
2897 write(366,*)xmcomp(i,1),xmcomp(i,2),xmcomp(i,3)
2898 write(367,*)xmcomp(i,4),xmcomp(i,5),xmcomp(i,6)
2900 write(368,*)xlcomp(i,1),xlcomp(i,2),xlcomp(i,3)
2901 write(369,*)xlcomp(i,4),xlcomp(i,5),xlcomp(i,6)
2903 write(370,*)xl1comp(i,1),xl1comp(i,2),xl1comp(i,3)
2904 write(371,*)xl1comp(i,4),xl1comp(i,5),xl1comp(i,6)
2906 write(372,*)xl2comp(i,1),xl2comp(i,2),xl2comp(i,3)
2907 write(373,*)xl2comp(i,4),xl2comp(i,5),xl2comp(i,6)
2909 write(374,*)xnd(i,1),xnd(i,2),scomp(i)
2910 write(375,*)ucomp(i,1),ucomp(i,2),ucomp(i,3)
3039 do i = 1,numnp
3040 xnude(ifath(i),2) = xnude(ifath(i),2) + strnrm(i)
3056 do i = 1,numnp
3057 xnude(ifath(i),1) = xnude(ifath(i),1) + x(i,1)
3058 xnude(ifath(i),2) = xnude(ifath(i),2) + x(i,2)
3069 do i=1,nfath ! cdelsq * |S|
3070 write(444,*) xnuder(i,1),xnuder(i,2),xnutf(i)