mof indent100.inp
para r0 1000;
para r1 2000;
para r2 3000;
para e0 100;
para e1 200;
para h 460;
para i1 1;
para i2 [%i1+int(%r0/%e0/2+0.5)];
para i3 [%i2+int(%r0/%e0/2+0.5)];
para i4 [%i3+int((%r1-%r0)/%e0+0.5)];
para i5 [%i4+int((%r2-%r1)/%e1+0.5)];
para j1 1;
para j2 [%j1+int(%r0/%e0/2+0.5)];
para j3 [%j2+int(%r0/%e0/2+0.5)];
para j4 [%j3+int((%r1-%r0)/%e0+0.5)];
para j5 [%j4+int((%r2-%r1)/%e1+0.5)];
para k1 1;
para k2 [%k1+int(%h/%e0+0.5)];
  
block %i1 %i2 %i3 %i4 %i5; %j1 %j2 %j3 %j4 %j5; %k1 %k2; 0 [%r0/2] %r0 %r1 %r2; 0 [%r0/2] %r0 %r1 %r2; 0 %h;
dei 2 5; 2 5;;
mbi -5; -2;; xy [(%r2-%r0/2)*cos(45)-(%r2-%r0/2)] [(%r2-%r0/2)*sin(45)]
mbi -4; -2;; xy [(%r1-%r0/2)*cos(45)-(%r1-%r0/2)] [(%r1-%r0/2)*sin(45)]
mbi -3; -2;; xy [(%r0-%r0/2)*cos(45)-(%r0-%r0/2)] [(%r0-%r0/2)*sin(45)]
bb 2 2 1 5 2 2 1;
bb 2 2 1 2 5 2 1;
sd 1 cy 0 0 0 0 0 1 %r0 ;
sd 2 cy 0 0 0 0 0 1 %r1;
sd 3 cy 0 0 0 0 0 1 %r2;
sfi -5;;;sd 3
sfi -4;;;sd 2
sfi -3;;;sd 1
sfi ; -5;;sd 3
sfi ; -4;;sd 2
sfi ; -3;;sd 1
res 3 1 1 4 5 2 i 1.1
res 1 3 1 5 4 2 j 1.1
res 2 1 1 3 5 2 i .8
res 1 2 1 5 3 2 j .8
res 1 1 1 5 5 2 k .9
eset 1 1 1 5 5 2 = eall
nset 1 1 1 5 5 2 = nall
nset 1 1 1 1 5 2 = nxm
fset 1 1 1 1 5 2 = fxm
nset 1 1 1 5 1 2 = nym
fset 1 1 1 5 1 2 = fym
nset 1 1 1 5 5 1 = nzm
fset 1 1 1 5 5 1 = fzm
nset 1 1 2 3 3 2 = nind
fset 1 1 2 3 3 2 = find
nset 3 1 2 5 5 2 = nzp
fset 3 1 2 5 5 2 = fzp
nset 1 3 2 2 5 2 or nzp
fset 1 3 2 2 5 2 or fzp
nset 5 1 1 5 5 2 = nperi
fset 5 1 1 5 5 2 = fperi
nset 1 5 1 5 5 2 or nperi
fset 1 5 1 5 5 2 or fperi
block %i1 %i2 %i3 [%i3+int(100/%e0+0.5)] [%i3+2*int(100/%e0+0.5)];
%j1 %j2 %j3 [%j3+int(100/%e0+0.5)] [%j3+2*int(100/%e0+0.5)];
1 [1+int(100/%e0+0.5)] [1+2*int(100/%e0+0.5)] [2+2*int(100/%e0+0.5)];
0 [%r0/2] [%r0-100] [%r0-50] %r0;
0 [%r0/2] [%r0-100] [%r0-50] %r0;
[%h+1] [%h+51] [%h+101] [2*%h];
dei 2 5; 2 5;;
mbi -3; -2;; xy [(%r0-100-%r0/2)*cos(45)-(%r0-100-%r0/2)] [(%r0-100-%r0/2)*sin(45)]
mbi -4; -2;; xy [(%r0-50-%r0/2)*cos(45)-(%r0-50-%r0/2)] [(%r0-50-%r0/2)*sin(45)]
mbi -5; -2;; xy [(%r0-%r0/2)*cos(45)-(%r0-%r0/2)] [(%r0-%r0/2)*sin(45)]
dei 4 5; 1 2; 1 2;
dei 1 2; 4 5; 1 2;
bb 2 2 1 4 2 4 2;
bb 4 2 2 5 2 4 3;
bb 2 2 1 2 4 4 2;
bb 2 4 2 2 5 4 3;
mbi -5; 1 2; -2; xz [50*cos(45)-50] [50*sin(45)-50]
mbi -1; -5; -2; yz [50*cos(45)-50] [50*sin(45)-50]
bb 4 1 2 5 2 2 4;
bb 4 1 1 4 2 2 4;
bb 1 4 2 2 5 2 5;
bb 1 4 1 2 4 2 5;
sfi -5; 1 2; 3 4;sd 1
sfi 1 2; -5; 3 4;sd 1
sd 8 pr 0 0 0 0 0 1 [%r0-100] [%h+1] [%r0-100+100*cos(45)] [%h+1+100*sin(45)] %r0 [%h+101];