implicit real*8(a-h,o-z)
dimension x(315500),y(315500),z(315500),e66(6,6),b(3,3)&
&,meshmoku(300000,8),meshkou(300000,8),meshkou2(300000,8),meshhogou(300000,8)
b1=8.4d-1
b2=0.9d-2
b3=0.7d-2!陬懷央
b4=6.8d-2!陬懷央
h2=2.55d-1
h3=0.5d-2!陬懷央
h4=1.13d-1!陬懷央
h5=0.7d-2!陬懷央
h1=h4+h5
ell=3.500d0
p=1.000d-1!闕ｷ驥�
ellhaji=1.500d-1!蜃ｺ縺ｦ縺�ｋ驛ｨ蛻��髟ｷ縺�
ellall=ell*2+ellhaji*2!ell*2
!ellall=ell+ellhaji*2!ell*1
!xyz譁ｹ蜷代�隕∫ｴ�蛻�牡謨ｰ
nx1=8!譛ｨ譚�
nx2=7!驪ｼ譚ｿ(nx2=1縺ｮ譎る蕎譚ｿ縺ｮ荳企Κ縺ｮ荳驛ｨ蛻�′縺ｪ縺上↑繧�)
nx3=5!陬懷央
nx4=5!陬懷央
ny1=5!荳企Κ蛻�→荳矩Κ蛻��鬮倥＆譁ｹ蜷代�蛻�牡謨ｰ
ny2=10!髢薙�驪ｼ譚ｿ縺ｮ鬮倥＆譁ｹ蜷代�蛻�牡謨ｰ(3莉･荳�)
ny3=2!陬懷央
ny4=3!陬懷央
ny5=3!陬懷央
nx=nx1+2*(nx2+nx3+nx4)+1
nznaka=50
nzhaji=10!蜃ｺ縺ｦ縺�ｋ驛ｨ蛻��蛻�牡謨ｰ
nz=nzhaji*2+nznaka*2
!繝､繝ｳ繧ｰ邇��蜈･蜉�
ezz=7.694d3!z霆ｸ譁ｹ蜷代�繝､繝ｳ繧ｰ邇�
exx=ezz/25.0d0!x霆ｸ譁ｹ蜷代�繝､繝ｳ繧ｰ邇�
eyy=ezz/25.0d0!y霆ｸ譁ｹ蜷代�繝､繝ｳ繧ｰ邇�
ezz2=206.0d3!z霆ｸ譁ｹ蜷代�繝､繝ｳ繧ｰ邇�
!繝昴い繧ｽ繝ｳ豈斐�蜈･蜉�
poixy=0.016d0
poixz=poixy
poiyx=poixy
poiyz=poixy
poizx=0.4d0
poizy=poizx
!縺帙ｓ譁ｭ蠑ｾ諤ｧ菫よ焚縺ｮ蜈･蜉�
gxy=ezz/15.0d0
gxz=gxy
gyz=gxy
!1縺､縺ゅ◆繧翫�隕∫ｴ�縺ｮ髟ｷ縺�?
xn1=b1/nx1
xn2=b2/nx2
xn3=b3/nx3
xn4=b4/nx4
yn1=h1/ny1
yn2=h2/ny2
yn3=h3/ny3
yn4=h4/ny4
yn5=h5/ny5
zn=ell/nznaka
znhaji=ellhaji/nzhaji
! 遽轤ｹ逡ｪ蜿ｷ縺斐→縺ｮ遽轤ｹ蠎ｧ讓�!i~m縺ｧ蟋九∪繧九�縺慧o縺ｧ菴ｿ縺医ｋ

!荳逡ｪ荳矩Κ蛻��驪ｼ譚舌�遽轤ｹ
do k=1,nz+1
do j=0,ny1
do l=0,1
do i=0,nx2-1
nset=1+i+l*(nx2+nx1+1)+j*(nx1+2*nx2+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn2+l*(b1+b2+xn2)+b4+b3
y(nset)=j*yn1
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳逡ｪ荳矩Κ蛻��遽轤ｹ 譛ｨ譚�
do k=1,nz+1
do j=0,ny1
do i=0,nx1
nset=nx2+1+i+j*(nx2*2+nx1+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn1+b2+b4+b3!i縺�0~xn2縺ｾ縺ｧ縺�▲縺ｦi=0,l=1縺ｫ縺ｪ繧�
y(nset)=j*yn1
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do

!荳九°繧�2逡ｪ逶ｮ縺ｮ驪ｼ譚舌�遽轤ｹ
do k=1,nz+1
do j=0,ny2-2
do l=0,1
do i=0,nx2
nset=(nx2*2+nx1+1)*(ny1+1)+i+1+l*(nx2+1)+j*2*(nx2+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn2+l*(b1+b2)+b4+b3
y(nset)=j*yn2+h1+yn2
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�3逡ｪ逶ｮ縺ｮ陬懷央譚宣Κ蛻��遽轤ｹ
do k=1,nz+1
do j=0,ny3-1
do l=0,1
do i=0,nx3-1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+1+i+l*(nx3+nx2+1+nx2+1)&
&+j*2*(nx3+nx2+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn3+l*(b3+b2+b1+b2+xn3)+b4
y(nset)=j*yn3+h1+h2
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�3逡ｪ逶ｮ縺ｮ驪ｼ譚宣Κ蛻��遽轤ｹ
do k=1,nz+1
do j=0,ny3-1
do l=0,1
do i=0,nx2
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2+nx3&
&+1+i+l*(nx2+1)&
&+j*2*(nx3+nx2+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn2+l*(b1+b2)+b4+b3
y(nset)=j*yn3+h1+h2
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ陬懷央譚宣Κ蛻��遽轤ｹ!縺薙％縺ｾ縺ｧ
do k=1,nz+1
do j=0,ny4-1
do l=0,1
do i=0,nx3-1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*ny3*2&
&+1+i+l*(nx3+nx2*2+nx1+1)&
&+j*(nx3*2+nx2*2+nx1+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn3+l*(b3+b2+b1+b2+xn3)+b4
y(nset)=j*yn4+h1+h2+h3
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ驪ｼ譚宣Κ蛻��遽轤ｹ
do k=1,nz+1
do j=0,ny4-1
do l=0,1
do i=0,nx2-1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+nx3&
&+1+i+l*(nx2+nx1+1)&
&+j*(nx3*2+(nx2+1)*2+nx1-1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn2+l*(b1+b2+xn2)+b3+b4
y(nset)=j*yn4+h1+h2+h3
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ譛ｨ譚宣Κ蛻��遽轤ｹ
do k=1,nz+1
do j=0,ny4-1
do i=0,nx1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+nx3+nx2&
&+1+i&
&+j*(nx3*2+(nx2+1)*2+nx1-1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn1+b4+b3+b2
y(nset)=j*yn4+h1+h2+h3
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
!荳九°繧�5逡ｪ逶ｮ縺ｮ陬懷央譚宣Κ蛻��遽轤ｹ縺昴�1�育ｪ√″蜃ｺ縺ｦ縺�ｋ驛ｨ蛻�ｼ�
do k=1,nz+1
do j=0,ny5
do l=0,1
do i=0,nx4-1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+1+i+l*(nx4+2*nx3+2*nx2+nx1+1)&
&+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn4+l*(b4+b3*2+b2*2+b1+xn4)
y(nset)=j*yn5+h1+h2+h3+h4
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�5逡ｪ逶ｮ縺ｮ陬懷央譚宣Κ蛻��遽轤ｹ縺昴�2
do k=1,nz+1
do j=0,ny5
do l=0,1
do i=0,nx3-1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+nx4&
&+1+i+l*(nx3+2*nx2+nx1+1)&
&+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn3+l*(b3+b2*2+b1+xn3)+b4
y(nset)=j*yn5+h1+h2+h3+h4
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do
!荳九°繧�5逡ｪ逶ｮ縺ｮ驪ｼ譚宣Κ蛻��遽轤ｹ
do k=1,nz+1
do j=0,ny5
do l=0,1
do i=0,nx2-1
nset=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+nx4+nx3&
&+1+i+l*(nx2+nx1+1)&
&+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn2+l*(b2+b1+xn2)+b4+b3
y(nset)=j*yn5+h1+h2+h3+h4
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
end do

!荳九°繧�5逡ｪ逶ｮ縺ｮ譛ｨ譚宣Κ蛻��遽轤ｹ
do k=1,nz+1
do j=0,ny5
do i=0,nx1
nset=(1+nx1+2*nx2)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+nx4+nx3+nx2&
&+1+i&
&+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+(k-1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
x(nset)=i*xn1+b4+b3+b2
y(nset)=j*yn5+h1+h2+h3+h4
!print'(i6,a,f9.6,a,f9.6,a,f9.6)',&
!&nset,',',x(nset),',',y(nset),',',z(nset)
if (k>=0.and.k<=nzhaji) then
z(nset)=znhaji*(k-1)
else if (k>nzhaji.and.k<=(nznaka*2+1+nzhaji)) then
z(nset)=zn*(k-nzhaji-1)+ellhaji
else
z(nset)=znhaji*(k-nzhaji-nznaka*2-1)+ellhaji+ell*2.0
end if
end do
end do
end do
!i5~6縺梧紛謨ｰ,a縺梧枚蟄�,f9,6縺悟ｰ乗焚.1pe縺ｨ縺九�E縺ｮ�贋ｹ�
!謗･轤ｹ逡ｪ蜿ｷ縺ｨ蠎ｧ讓吶�蜃ｺ蜉�
nsettensuu=(nz+1)*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
print'("*NODE, NSET=Nall")'
do i=1,nsettensuu
print'(i6,a,1pe13.5,a,1pe13.5,a,1pe13.5)',&
&i,',',x(i),',',y(i),',',z(i)
end do

!縺薙％縺九ｉ
!隕∫ｴ�逡ｪ蜿ｷ縺ｨ縺昴�隕∫ｴ�縺ｮ8謗･轤ｹ縺ｮ謗･轤ｹ逡ｪ蜿ｷ
!譛ｨ譚舌�隕∫ｴ�縲荳逡ｪ荳九°繧蛾��分縺ｫ
nyou=0
do k=0,nz-1
do j=0,ny1-1
do i=1,nx1
nyou=nyou+1
meshmoku(nyou,1)=nx2+i+j*(2*nx2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,2)=meshmoku(nyou,1)+(2*nx2+nx1+1)
meshmoku(nyou,3)=meshmoku(nyou,2)&
&+(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1)
meshmoku(nyou,4)=meshmoku(nyou,1)&
&+(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1)
meshmoku(nyou,5)=meshmoku(nyou,1)+1
meshmoku(nyou,6)=meshmoku(nyou,2)+1
meshmoku(nyou,7)=meshmoku(nyou,3)+1
meshmoku(nyou,8)=meshmoku(nyou,4)+1
end do
end do
end do
!荳九°繧�4谿ｵ逶ｮ縺ｮ譛ｨ譚舌�隕∫ｴ�(荳逡ｪ荳九�谿ｵ縲應ｸ逡ｪ荳翫�谿ｵ縺ｮ1縺､蜑阪�谿ｵ縺ｾ縺ｧ)
do k=0,nz-1
do j=0,ny4-2
do i=1,nx1
nyou=nyou+1
meshmoku(nyou,1)=(2*nx2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)+2*(nx3+nx2+1)*ny3&
&+nx3+nx2+i+j*((nx3+nx2)*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,2)=meshmoku(nyou,1)+((nx3+nx2)*2+nx1+1)
meshmoku(nyou,3)=meshmoku(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,4)=meshmoku(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,5)=meshmoku(nyou,1)+1
meshmoku(nyou,6)=meshmoku(nyou,2)+1
meshmoku(nyou,7)=meshmoku(nyou,3)+1
meshmoku(nyou,8)=meshmoku(nyou,4)+1
end do
end do
end do
!荳九°繧�4谿ｵ逶ｮ縺ｮ譛ｨ譚舌�隕∫ｴ�(荳逡ｪ荳翫�谿ｵ)
do k=0,nz-1
do i=1,nx1
nyou=nyou+1
meshmoku(nyou,1)=(2*nx2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)+2*(nx3+nx2+1)*ny3&
&+(nx3*2+nx2*2+nx1+1)*(ny4-1)+nx3+nx2+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,2)=(2*nx2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)+2*(nx3+nx2+1)*ny3&
&+(nx3*2+nx2*2+nx1+1)*(ny4)+nx4+nx3+nx2+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,3)=meshmoku(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,4)=meshmoku(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,5)=meshmoku(nyou,1)+1
meshmoku(nyou,6)=meshmoku(nyou,2)+1
meshmoku(nyou,7)=meshmoku(nyou,3)+1
meshmoku(nyou,8)=meshmoku(nyou,4)+1
end do
end do
!荳九°繧�5谿ｵ逶ｮ縺ｮ譛ｨ譚舌�隕∫ｴ�
do k=0,nz-1
do j=0,ny5-1
do i=1,nx1
nyou=nyou+1
meshmoku(nyou,1)=(2*nx2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+nx4+nx3+nx2+i+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,2)=meshmoku(nyou,1)+nx4*2+nx3*2+nx2*2+nx1+1
meshmoku(nyou,3)=meshmoku(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,4)=meshmoku(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshmoku(nyou,5)=meshmoku(nyou,1)+1
meshmoku(nyou,6)=meshmoku(nyou,2)+1
meshmoku(nyou,7)=meshmoku(nyou,3)+1
meshmoku(nyou,8)=meshmoku(nyou,4)+1
end do
end do
end do
!荳逡ｪ荳九�谿ｵ縺ｮ驪ｼ譚ｿ
nyou=0
do m=0,1
do k=0,nz-1
do j=0,ny1-1
do i=1,nx2
nyou=nyou+1
meshkou(nyou,1)=i+j*(nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou(nyou,2)=meshkou(nyou,1)+nx2*2+nx1+1
meshkou(nyou,3)=meshkou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,4)=meshkou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,5)=meshkou(nyou,1)+1
meshkou(nyou,6)=meshkou(nyou,2)+1
meshkou(nyou,7)=meshkou(nyou,3)+1
meshkou(nyou,8)=meshkou(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�2谿ｵ逶ｮ縺ｮ驪ｼ譚ｿ(譛蛻昴�荳谿ｵ逶ｮ)
do m=0,1
do k=0,nz-1
do i=1,nx2
nyou=nyou+1
meshkou(nyou,1)=(nx2*2+nx1+1)*(ny1)+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou(nyou,2)=(2*nx2+nx1+1)*(ny1+1)+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+1)
meshkou(nyou,3)=meshkou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,4)=meshkou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,5)=meshkou(nyou,1)+1
meshkou(nyou,6)=meshkou(nyou,2)+1
meshkou(nyou,7)=meshkou(nyou,3)+1
meshkou(nyou,8)=meshkou(nyou,4)+1
end do
end do
end do
!荳九°繧�2谿ｵ逶ｮ縺ｮ驪ｼ譚ｿ(2谿ｵ逶ｮ縺九ｉ譛蠕後�蜑阪∪縺ｧ)
do m=0,1
do k=0,nz-1
do j=0,ny2-3!ny2>3
do i=1,nx2
nyou=nyou+1
meshkou(nyou,1)=(2*nx2+nx1+1)*(ny1+1)+i&
&+j*2*(nx2+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+1)
meshkou(nyou,2)=meshkou(nyou,1)+(nx2+1)*2
meshkou(nyou,3)=meshkou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,4)=meshkou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,5)=meshkou(nyou,1)+1
meshkou(nyou,6)=meshkou(nyou,2)+1
meshkou(nyou,7)=meshkou(nyou,3)+1
meshkou(nyou,8)=meshkou(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�2谿ｵ逶ｮ縺ｮ驪ｼ譚ｿ(譛蠕後�谿ｵ)
do m=0,1
do k=0,nz-1
do i=1,nx2
nyou=nyou+1
meshkou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-2)*2+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+1)
meshkou(nyou,2)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2+nx3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+1)
meshkou(nyou,3)=meshkou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,4)=meshkou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,5)=meshkou(nyou,1)+1
meshkou(nyou,6)=meshkou(nyou,2)+1
meshkou(nyou,7)=meshkou(nyou,3)+1
meshkou(nyou,8)=meshkou(nyou,4)+1
end do
end do
end do
!荳九°繧�3逡ｪ逶ｮ縺ｮ驪ｼ譚舌�隕∫ｴ�(荳逡ｪ荳九應ｸ逡ｪ荳翫�谿ｵ縺ｮ1縺､荳九�谿ｵ縺ｾ縺ｧ)
do m=0,1
do k=0,nz-1
do j=0,ny3-2
do i=1,nx2
nyou=nyou+1
meshkou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)+nx3+i&
&+j*2*(nx3+nx2+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+1)
meshkou(nyou,2)=meshkou(nyou,1)+2*(nx3+nx2+1)
meshkou(nyou,3)=meshkou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,4)=meshkou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,5)=meshkou(nyou,1)+1
meshkou(nyou,6)=meshkou(nyou,2)+1
meshkou(nyou,7)=meshkou(nyou,3)+1
meshkou(nyou,8)=meshkou(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�3逡ｪ逶ｮ縺ｮ驪ｼ譚舌�隕∫ｴ�(1逡ｪ荳翫�谿ｵ)
do m=0,1
do k=0,nz-1
do i=1,nx2
nyou=nyou+1
meshkou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3-1)+nx3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+1)
meshkou(nyou,2)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+i+2*(nx3+nx2+1)*ny3+nx3&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)!m=1縺ｮ譎ゅ↓+1
meshkou(nyou,3)=meshkou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,4)=meshkou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou(nyou,5)=meshkou(nyou,1)+1
meshkou(nyou,6)=meshkou(nyou,2)+1
meshkou(nyou,7)=meshkou(nyou,3)+1
meshkou(nyou,8)=meshkou(nyou,4)+1
end do
end do
end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ驪ｼ譚舌�隕∫ｴ�縺昴�1�井ｸ逡ｪ荳九�谿ｵ縺九ｉ荳逡ｪ荳翫�1縺､荳九�谿ｵ縺ｾ縺ｧ)
!縺薙％縺九ｉmeshkou2縺ｫ縺励※荳九°繧�4谿ｵ逶ｮ縺ｨ5谿ｵ逶ｮ縺ｮ驪ｼ譚舌′縺ｪ縺��ｴ蜷医ｒ
!nx2=1縺ｮ縺ｨ縺咲ｰ｡蜊倥↓豎ゅａ繧峨ｌ繧九ｈ縺�↓縺吶ｋ
nyou=0
do m=0,1
do k=0,nz-1
do j=0,ny4-2
do i=1,nx2
nyou=nyou+1
meshkou2(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3)+nx3+i+j*(nx3*2+nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou2(nyou,2)=meshkou2(nyou,1)+nx3*2+nx2*2+nx1+1
meshkou2(nyou,3)=meshkou2(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,4)=meshkou2(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,5)=meshkou2(nyou,1)+1
meshkou2(nyou,6)=meshkou2(nyou,2)+1
meshkou2(nyou,7)=meshkou2(nyou,3)+1
meshkou2(nyou,8)=meshkou2(nyou,4)+1
end do
end do
end do
end do
!do i=1,2*nx2*(ny4+ny5)*nz
!print'(i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6)',&
!&i,',',&
!&meshkou2(i,1),',',&
!&meshkou2(i,2),',',&
!&meshkou2(i,3),',',&
!&meshkou2(i,4),',',&
!&meshkou2(i,5),',',&
!&meshkou2(i,6),',',&
!&meshkou2(i,7),',',&
!&meshkou2(i,8)
!end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ驪ｼ譚舌�隕∫ｴ�(荳逡ｪ荳翫�1縺､荳九�谿ｵ縺九ｉ荳逡ｪ荳翫�谿ｵ縺ｾ縺ｧ)
do m=0,1
do k=0,nz-1
do i=1,nx2
nyou=nyou+1
meshkou2(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3)+(nx3*2+nx2*2+nx1+1)*(ny4-1)+nx3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou2(nyou,2)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3)+(nx3*2+nx2*2+nx1+1)*(ny4)+nx4+nx3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou2(nyou,3)=meshkou2(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,4)=meshkou2(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,5)=meshkou2(nyou,1)+1
meshkou2(nyou,6)=meshkou2(nyou,2)+1
meshkou2(nyou,7)=meshkou2(nyou,3)+1
meshkou2(nyou,8)=meshkou2(nyou,4)+1
end do
end do
end do

!荳九°繧�5逡ｪ逶ｮ縺ｮ驪ｼ譚舌�隕∫ｴ�(荳逡ｪ荳九°繧我ｸ逡ｪ荳翫�谿ｵ縺ｮ1縺､荳九�谿ｵ縺ｾ縺ｧ��
do m=0,1
do k=0,nz-1
do j=0,ny5-2
do i=1,nx2
nyou=nyou+1
meshkou2(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3)+(nx3*2+nx2*2+nx1+1)*(ny4)+nx4+nx3&
&+i+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou2(nyou,2)=meshkou2(nyou,1)+nx3*2+nx2*2+nx1+1+nx4*2
meshkou2(nyou,3)=meshkou2(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,4)=meshkou2(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,5)=meshkou2(nyou,1)+1
meshkou2(nyou,6)=meshkou2(nyou,2)+1
meshkou2(nyou,7)=meshkou2(nyou,3)+1
meshkou2(nyou,8)=meshkou2(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�5逡ｪ逶ｮ縺ｮ驪ｼ譚舌�隕∫ｴ�(荳逡ｪ荳翫�谿ｵ縺ｮ1縺､荳九�谿ｵ縺九ｉ荳逡ｪ荳翫�谿ｵ縺ｾ縺ｧ��
do m=0,1
do k=0,nz-1
do i=1,nx2
nyou=nyou+1
meshkou2(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3)+(nx3*2+nx2*2+nx1+1)*(ny4)&
&+(nx4*2+nx3*2+nx2*2+nx1+1)*(ny5-1)+nx4+nx3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou2(nyou,2)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+2*(nx3+nx2+1)*(ny3)+(nx3*2+nx2*2+nx1+1)*(ny4)&
&+(nx4*2+nx3*2+nx2*2+nx1+1)*(ny5)+nx4+nx3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx2+nx1)
meshkou2(nyou,3)=meshkou2(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,4)=meshkou2(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshkou2(nyou,5)=meshkou2(nyou,1)+1
meshkou2(nyou,6)=meshkou2(nyou,2)+1
meshkou2(nyou,7)=meshkou2(nyou,3)+1
meshkou2(nyou,8)=meshkou2(nyou,4)+1
end do
end do
end do
!縺薙％縺ｾ縺ｧmeshkou2
!荳九°繧�3逡ｪ逶ｮ縺ｮ陬懷央譚舌�隕∫ｴ�(荳逡ｪ荳救荳逡ｪ荳翫�谿ｵ縺ｮ1縺､荳九�谿ｵ縺ｾ縺ｧ)
nyou=0
do m=0,1
do k=0,nz-1
do j=0,ny3-2
do i=1,nx3
nyou=nyou+1
meshhogou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+i+j*(nx3+nx2+1)*2&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2*2+1)
meshhogou(nyou,2)=meshhogou(nyou,1)+(nx3+nx2+1)*2
meshhogou(nyou,3)=meshhogou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,4)=meshhogou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,5)=meshhogou(nyou,1)+1
meshhogou(nyou,6)=meshhogou(nyou,2)+1
meshhogou(nyou,7)=meshhogou(nyou,3)+1
meshhogou(nyou,8)=meshhogou(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�3逡ｪ逶ｮ縺ｮ陬懷央譚舌�隕∫ｴ�(荳逡ｪ荳翫�谿ｵ)
do m=0,1
do k=0,nz-1
do i=1,nx3
nyou=nyou+1
meshhogou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+2*(nx2+1)*(ny2-1)&
&+(nx3+nx2+1)*2*(ny3-1)+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2*2+1)
meshhogou(nyou,2)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2*2+nx1)
meshhogou(nyou,3)=meshhogou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,4)=meshhogou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,5)=meshhogou(nyou,1)+1
meshhogou(nyou,6)=meshhogou(nyou,2)+1
meshhogou(nyou,7)=meshhogou(nyou,3)+1
meshhogou(nyou,8)=meshhogou(nyou,4)+1
end do
end do
end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ陬懷央譚舌�隕∫ｴ�(荳逡ｪ荳救荳逡ｪ荳翫�1縺､荳九�谿ｵ縺ｾ縺ｧ)
do m=0,1
do k=0,nz-1
do j=0,ny4-2
do i=1,nx3
nyou=nyou+1
meshhogou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2+(nx3+nx2+1)*2*ny3&
&+i+j*(nx3*2+nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2+nx1+nx2)
meshhogou(nyou,2)=meshhogou(nyou,1)+(nx3*2+nx2*2+nx1+1)
meshhogou(nyou,3)=meshhogou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,4)=meshhogou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,5)=meshhogou(nyou,1)+1
meshhogou(nyou,6)=meshhogou(nyou,2)+1
meshhogou(nyou,7)=meshhogou(nyou,3)+1
meshhogou(nyou,8)=meshhogou(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�4逡ｪ逶ｮ縺ｮ陬懷央譚舌�隕∫ｴ�(荳逡ｪ荳翫�谿ｵ)
do m=0,1
do k=0,nz-1
do i=1,nx3
nyou=nyou+1
meshhogou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2+(nx3+nx2+1)*2*ny3&
&+(nx3*2+nx2*2+nx1+1)*(ny4-1)+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2+nx1+nx2)
meshhogou(nyou,2)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+(nx3*2+nx2*2+nx1+1)*ny4+nx4+i&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2+nx1+nx2)
meshhogou(nyou,3)=meshhogou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,4)=meshhogou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,5)=meshhogou(nyou,1)+1
meshhogou(nyou,6)=meshhogou(nyou,2)+1
meshhogou(nyou,7)=meshhogou(nyou,3)+1
meshhogou(nyou,8)=meshhogou(nyou,4)+1
end do
end do
end do
!荳九°繧�5逡ｪ逶ｮ縺ｮ陬懷央譚舌�隕∫ｴ�縺昴�1�井ｸ逡ｪ螟門�縺ｮ遯√″蜃ｺ縺励※縺�ｋ驛ｨ蛻�)
do m=0,1
do k=0,nz-1
do j=0,ny5-1
do i=1,nx4
nyou=nyou+1
meshhogou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2+(nx3+nx2+1)*2*ny3&
&+(nx3*2+nx2*2+nx1+1)*ny4+i+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx4+nx3*2+nx2*2+nx1)
meshhogou(nyou,2)=meshhogou(nyou,1)+(nx4*2+nx3*2+nx2*2+nx1+1)
meshhogou(nyou,3)=meshhogou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,4)=meshhogou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,5)=meshhogou(nyou,1)+1
meshhogou(nyou,6)=meshhogou(nyou,2)+1
meshhogou(nyou,7)=meshhogou(nyou,3)+1
meshhogou(nyou,8)=meshhogou(nyou,4)+1
end do
end do
end do
end do
!荳九°繧�5逡ｪ逶ｮ縺ｮ陬懷央譚舌�隕∫ｴ�縺昴�2�亥�蛛ｴ縺ｮ驛ｨ蛻�)
do m=0,1
do k=0,nz-1
do j=0,ny5-1
do i=1,nx3
nyou=nyou+1
meshhogou(nyou,1)=(nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2+(nx3+nx2+1)*2*ny3&
&+(nx3*2+nx2*2+nx1+1)*ny4+nx4+i+j*(nx4*2+nx3*2+nx2*2+nx1+1)&
&+k*((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))&
&+m*(nx3+nx2*2+nx1)
meshhogou(nyou,2)=meshhogou(nyou,1)+(nx4*2+nx3*2+nx2*2+nx1+1)
meshhogou(nyou,3)=meshhogou(nyou,2)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,4)=meshhogou(nyou,1)&
&+((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
meshhogou(nyou,5)=meshhogou(nyou,1)+1
meshhogou(nyou,6)=meshhogou(nyou,2)+1
meshhogou(nyou,7)=meshhogou(nyou,3)+1
meshhogou(nyou,8)=meshhogou(nyou,4)+1
end do
end do
end do
end do
!險ｭ螳�

print'("*ELEMENT, TYPE=C3D8, ELSET=Ekou")'
do i=1,2*nx2*(ny1+ny2+ny3)*nz
print'(i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6)',&
&i,',',&
&meshkou(i,1),',',&
&meshkou(i,2),',',&
&meshkou(i,3),',',&
&meshkou(i,4),',',&
&meshkou(i,5),',',&
&meshkou(i,6),',',&
&meshkou(i,7),',',&
&meshkou(i,8)
end do
nyousosuukou=2*nx2*(ny1+ny2+ny3)*nz
print'("*ELEMENT, TYPE=C3D8, ELSET=Ekou2")'
do i=1,2*nx2*(ny4+ny5)*nz
print'(i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6)',&
&i+nyousosuukou,',',&
&meshkou2(i,1),',',&
&meshkou2(i,2),',',&
&meshkou2(i,3),',',&
&meshkou2(i,4),',',&
&meshkou2(i,5),',',&
&meshkou2(i,6),',',&
&meshkou2(i,7),',',&
&meshkou2(i,8)
end do
nyousosuukou2=2*nx2*(ny4+ny5)*nz
print'("*ELEMENT, TYPE=C3D8, ELSET=Ehogou")'
do i=1,2*(nx4*ny5+nx3*(ny3+ny4+ny5))*nz
print'(i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6)',&
&i+nyousosuukou+nyousosuukou2,',',&
&meshhogou(i,1),',',&
&meshhogou(i,2),',',&
&meshhogou(i,3),',',&
&meshhogou(i,4),',',&
&meshhogou(i,5),',',&
&meshhogou(i,6),',',&
&meshhogou(i,7),',',&
&meshhogou(i,8)
end do
nyousosuuhogou=2*(nx4*ny5+nx3*(ny3+ny4+ny5))*nz
print'("*ELEMENT, TYPE=C3D8, ELSET=Emoku")'
do i=1,nx1*(ny1+ny4+ny5)*nz
print'(i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6,a,i6)',&
&i+nyousosuukou+nyousosuuhogou+nyousosuukou2,',',&
&meshmoku(i,1),',',&
&meshmoku(i,2),',',&
&meshmoku(i,3),',',&
&meshmoku(i,4),',',&
&meshmoku(i,5),',',&
&meshmoku(i,6),',',&
&meshmoku(i,7),',',&
&meshmoku(i,8)
end do

!霈芽差謗･轤ｹ縲∵鋸譚滓磁轤ｹ縺ｮ險ｭ螳�
nitimen=((nx2*2+nx1+1)*(ny1+1)+(nx2+1)*(ny2-1)*2&
&+(nx3+nx2+1)*2*ny3+((nx3+nx2)*2+nx1+1)*ny4&
&+((nx4+nx3+nx2)*2+nx1+1)*(ny5+1))
print'("*NSET,NSET=kousoku")'

do i=1,nitimen*(nz+1)
if (z(i)>ellhaji-0.000001d0.and.z(i)<ellhaji+0.000001d0&
&.and.y(i)<0.000001d0) then
nkousoku=i
print'(i6,a)',nkousoku,','
end if
end do
print'("*NSET,NSET=kousoku2")'
do i=1,nitimen*(nz+1)
if (z(i)>ellall-ellhaji-0.000001d0.and.z(i)<ellall-ellhaji+0.000001d0&
&.and.y(i)<0.000001d0) then
nkousoku=i
print'(i6,a)',nkousoku,','
end if
end do
print'("*NSET,NSET=saika")'
do i=1,nitimen*(nz+1)
if (z(i)>ellall/2-0.000001d0.and.z(i)<ellall/2+0.000001d0&
&.and.y(i)>h1+h2+h3+h4+h5-0.0000001d0) then
nsaika=i
print'(i6,a)',nsaika,','
end if
end do
! 蠅�阜譚｡莉ｶ
print'("*BOUNDARY")'
print'("kousoku,1,3")'
print'("kousoku2,1,2")'
!print'("kousoku2,2")'
! 譛ｨ譚舌�險ｭ螳�
print'("*MATERIAL,NAME=ELmoku")'
print'("*ELASTIC,TYPE=ORTHOROPIC")'
call ortho(e66,exx,eyy,ezz,poixy,poixz,poiyz,gyz,&
& gxz,gxy,poiyx,poizx,poizy)
print'(f12.4,a,f12.4,a,f12.4,a,&
&f12.4,a,f12.4,a,f12.4,a,f12.4,&
& a ,f12.4,a)',&
&e66(1,1),',',e66(1,2),',',e66(2,2),',',e66(1,3),','&
&,e66(2,3),',',e66(3,3),',',e66(6,6),',',e66(5,5),','
print'(f12.4,a,a)',e66(4,4),',','293'
print'("*MATERIAL,NAME=ELkou")'
print'("*ELASTIC")'
print'("206.d3, 0.3")'
print'("*SOLID SECTION,ELSET=Emoku,MATERIAL=ELmoku")'
print'("*SOLID SECTION,ELSET=Ekou,MATERIAL=ELkou")'
print'("*SOLID SECTION,ELSET=Ekou2,MATERIAL=ELkou")'
print'("*SOLID SECTION,ELSET=Ehogou,MATERIAL=ELkou")'
! 霈芽差譚｡莉ｶ
print'("*STEP")'
print'("*STATIC,SOLVER=SPOOLES")'
print'("*CLOAD")'
print'(A,a,i6,a,1pe13.5)','saika',',',2,',',p/nx
!!p/((nx+1)*(ny+1))

! dat縺ｮ繧ｪ繝励す繝ｧ繝ｳ
print'("*NODE PRINT,NSET=saika")'
print'("U")'
print'("*NODE FILE")'
print'("U")'
print'("*END STEP")'
!
!
!
end
!######## 逶ｴ莠､逡ｰ譁ｹ諤ｧ縺ｮ螳壽焚繧貞�縺吶�繝ｭ繧ｰ繝ｩ繝� ########
subroutine ortho(e66,exx,eyy,ezz,poixy,poixz,poiyz,gyz,&
&gxz,gxy,poiyx,poizx,poizy)
implicit real*8(a-h,o-z)
dimension e66(6,6),b(3,3)
!蠢懷鴨�阪�縺壹∩陦悟��啼66
do i=1,6
do j=1,6
e66(i,j)=0.d0
end do
end do
!
e66(1,1)=1.d0/exx
e66(1,2)=-poixy/exx
e66(1,3)=-poixz/exx
e66(2,1)=-poiyx/eyy
e66(2,2)=1.d0/eyy
e66(2,3)=-poiyz/eyy
e66(3,1)=-poizx/ezz
e66(3,2)=-poizy/ezz
e66(3,3)=1.d0/ezz
!
!****** 騾�｡悟�繧貞叙縺｣縺ｦ縲∝ｿ懷鴨�阪�縺壹∩陦悟�縺ｫ螟画鋤縺吶ｋ******
det=e66(1,1)*e66(2,2)*e66(3,3)+e66(2,1)*e66(3,2)*e66(1,3)&
&+e66(3,1)*e66(1,2)*e66(2,3)-e66(2,1)*e66(1,2)*e66(3,3)&
&-e66(3,1)*e66(2,2)*e66(1,3)-e66(1,1)*e66(3,2)*e66(2,3)
b(1,1)=(-1.d0)**(1+1)*(e66(2,2)*e66(3,3)-e66(3,2)*e66(2,3))
b(1,2)=(-1.d0)**(1+2)*(e66(2,1)*e66(3,3)-e66(3,1)*e66(2,3))
b(1,3)=(-1.d0)**(1+3)*(e66(2,1)*e66(3,2)-e66(3,1)*e66(2,2))
b(2,1)=(-1.d0)**(2+1)*(e66(1,2)*e66(3,3)-e66(3,2)*e66(1,3))
b(2,2)=(-1.d0)**(2+2)*(e66(1,1)*e66(3,3)-e66(3,1)*e66(1,3))
b(2,3)=(-1.d0)**(2+3)*(e66(1,1)*e66(3,2)-e66(3,1)*e66(1,2))
b(3,1)=(-1.d0)**(3+1)*(e66(1,2)*e66(2,3)-e66(2,2)*e66(1,3))
b(3,2)=(-1.d0)**(3+2)*(e66(1,1)*e66(2,3)-e66(2,1)*e66(1,3))
b(3,3)=(-1.d0)**(3+3)*(e66(1,1)*e66(2,2)-e66(2,1)*e66(1,2))
do i=1,3
do j=1,3
e66(i,j)=b(j,i)/det
end do
end do
!
!****** 縺帙ｓ譁ｭ陬懈ｭ｣菫よ焚繧貞�蜉� ******
!print*,'xy髱｢縺ｮ縺帙ｓ譁ｭ蠑ｾ諤ｧ菫よ焚繧貞�繧後※荳九＆縺�'
!read*, gxy
e66(4,4)=gyz
!print*,'xz髱｢縺ｮ縺帙ｓ譁ｭ蠑ｾ諤ｧ菫よ焚繧貞�繧後※荳九＆縺�'
!read*, gxz
e66(5,5)=gxz
!print*,'yz髱｢縺ｮ縺帙ｓ譁ｭ蠑ｾ諤ｧ菫よ焚繧貞�繧後※荳九＆縺�'
!read*, gyz
e66(6,6)=gxy
!
!
return
end