c
c This file should not be run by itself
c It is INCLUDED in the file 'filter'
c
update

C LOWER FACE OF HOLE
sd %nofsrf plan
   %xcent [0.25*%width+(%i-1)*%width] %zcent
   0 [cos(%xang)] [sin(%xang)]
sfi 2 3;[-2*%i];;sd %nofsrf
para nofsrf [%nofsrf+1];

C UPPER FACE OF HOLE
sd %nofsrf plan
   %xcent [0.75*%width+(%i-1)*%width] %zcent
   0 [cos(%xang)] [sin(%xang)]
sfi 2 3;[-(2*%i+1)];;sd %nofsrf
para nofsrf [%nofsrf+1];

C LEFT FACE OF HOLE
sd %nofsrf plan
   [0.25*%width] %ycent %zcent
   [cos(%yang)] 0 [-sin(%yang)]
sfi -2;[2*%i] [2*%i+1];; sd %nofsrf
para nofsrf [%nofsrf+1];

C RIGHT FACE OF HOLE
sd %nofsrf plan
   [0.75*%width] %ycent %zcent
   [cos(%yang)] 0 [-sin(%yang)]
sfi -3;[2*%i] [2*%i+1];; sd %nofsrf
para nofsrf [%nofsrf+1];

C FREEZE THE MESH
update

C PROJECT HOLES TO CYLINDERS
sd %nofsrf cyli 
      %xcent [0.5*%width+(%i-1)*%width] %zcent
      [sin(%yang)*cos(%xang)] [-sin(%xang)] [cos(%yang)*cos(%xang)]
      [0.25*%width]
sfi -2 -3;[-2*%i] [-(2*%i+1)];; sd %nofsrf
para %nofsrf [%nofsrf+1];

C PUT THE EDGES OF THE HOLES BACK ON THE OUTER CYLINDER
C (MEMORY IS LOST AFTER AN UPDATE)
sfi 2 3;[-2*%i];-1; sd 1
sfi 2 3;[-(2*%i+1)];-1; sd 1
sfi 2 3;[-2*%i];-2; sd 1
sfi 2 3;[-(2*%i+1)];-2; sd 1
sfi -2;[2*%i] [2*%i+1];-1;sd 1
sfi -2;[2*%i] [2*%i+1];-2;sd 1
sfi -3;[2*%i] [2*%i+1];-1;sd 1
sfi -3;[2*%i] [2*%i+1];-2;sd 1

C INTERPOLATE THE FACES AND EDGES OF THE HOLE
C (AFTER AN UPDATE, INTERPOLATIONS ARE NO LONGER AUTOMATIC)
lini -2 0 -3;[-2*%i] 0 [-(2*%i+1)];;
lini -2 -3;[-2*%i] [-(2*%i+1)];;
