// stella.java  develope stellaoctamgula 8 points from a cube

import java.io.*;

public class stella
{
  boolean debug = true;
  int pts, pt, np, j;
  int k = 0; // index into newpts k++
  datread DR; // = new datread(filename); // defines data structures
  int new_points = 0; // for each x,y,z point kept
  int new_polys = 0; // k/4 for datwrite
  double xyz[][] = new double[10000][3]; // x, y, z for output
  int newpts[] = new int[40000]; // 3, xyz1, xyz2, xyz3 for output
  String fileout = "stella_java.dat";
  String fileout1 = "stella1_java.dat";
  String fileout2 = "stella2_java.dat";
  double x, y, z, dz;
  int box = 1;
    
  public  stella()
  {
    System.out.println("stella.java running");
    new_points = 0;

    System.out.println(" ");
    System.out.println("generate vertex for a cube");
    xyz[new_points][0] =-1.0; // 1
    xyz[new_points][1] =-1.0;
    xyz[new_points][2] =-1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] =-1.0; // 2
    xyz[new_points][1] =-1.0;
    xyz[new_points][2] = 1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] =-1.0; // 3
    xyz[new_points][1] = 1.0;
    xyz[new_points][2] =-1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] =-1.0; // 4
    xyz[new_points][1] = 1.0;
    xyz[new_points][2] = 1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 1.0; // 5
    xyz[new_points][1] =-1.0;
    xyz[new_points][2] =-1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 1.0; // 6
    xyz[new_points][1] = 1.0;
    xyz[new_points][2] =-1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 1.0; // 7
    xyz[new_points][1] =-1.0;
    xyz[new_points][2] = 1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 1.0; // 8
    xyz[new_points][1] = 1.0;
    xyz[new_points][2] = 1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 0.0; // 9  centers, 6
    xyz[new_points][1] = 0.0;
    xyz[new_points][2] =-1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 1.0; // 10
    xyz[new_points][1] = 0.0;
    xyz[new_points][2] = 0.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 0.0; // 11
    xyz[new_points][1] = 0.0;
    xyz[new_points][2] = 1.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] =-1.0; // 12
    xyz[new_points][1] = 0.0;
    xyz[new_points][2] = 0.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 0.0; // 13
    xyz[new_points][1] = 1.0;
    xyz[new_points][2] = 0.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;
    xyz[new_points][0] = 0.0; // 14
    xyz[new_points][1] =-1.0;
    xyz[new_points][2] = 0.0;
    System.out.println("pt "+(new_points+1)+" x="+xyz[new_points][0]+" y="+xyz[new_points][1]+" z="+xyz[new_points][2]);
    new_points = new_points +1;

    k = 0; // all triangles
    new_polys = 0;
    System.out.println("tetra 1");
    newpts[k++] = 3;
    newpts[k++] = 5;
    newpts[k++] = 9;
    newpts[k++] = 10;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 1b
    newpts[k++] = 5;
    newpts[k++] = 9;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 1c
    newpts[k++] = 5;
    newpts[k++] = 10;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 1d
    newpts[k++] = 9;
    newpts[k++] = 10;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 2");
    newpts[k++] = 3;
    newpts[k++] = 7;
    newpts[k++] = 10;
    newpts[k++] = 11;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 2b
    newpts[k++] = 7;
    newpts[k++] = 10;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 2c
    newpts[k++] = 7;
    newpts[k++] = 11;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 2d
    newpts[k++] = 10;
    newpts[k++] = 11;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 3");
    newpts[k++] = 3;
    newpts[k++] = 8;
    newpts[k++] = 10;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 3b
    newpts[k++] = 8;
    newpts[k++] = 10;
    newpts[k++] = 11;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 3c
    newpts[k++] = 8;
    newpts[k++] = 11;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 3d
    newpts[k++] = 10;
    newpts[k++] = 11;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 4");
    newpts[k++] = 3;
    newpts[k++] = 3;
    newpts[k++] = 12;
    newpts[k++] = 9;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 4b
    newpts[k++] = 3;
    newpts[k++] = 12;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 4c
    newpts[k++] = 3;
    newpts[k++] = 9;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 4d
    newpts[k++] = 9;
    newpts[k++] = 12;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 5 ");
    newpts[k++] = 3;
    newpts[k++] = 6;
    newpts[k++] = 10;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 5b
    newpts[k++] = 6;
    newpts[k++] = 10;
    newpts[k++] = 9;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 5c
    newpts[k++] = 6;
    newpts[k++] = 9;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 5d
    newpts[k++] = 9;
    newpts[k++] = 10;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 6 ");
    newpts[k++] = 3;
    newpts[k++] = 1;
    newpts[k++] = 9;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 6b
    newpts[k++] = 1;
    newpts[k++] = 9;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 6c
    newpts[k++] = 1;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 6d
    newpts[k++] = 9;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 7 ");
    newpts[k++] = 3;
    newpts[k++] = 4;
    newpts[k++] = 13;
    newpts[k++] = 11;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 7b
    newpts[k++] = 4;
    newpts[k++] = 13;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 7c
    newpts[k++] = 4;
    newpts[k++] = 11;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 7d
    newpts[k++] = 11;
    newpts[k++] = 12;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 8 ");
    newpts[k++] = 3;
    newpts[k++] = 2;
    newpts[k++] = 11;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 8b
    newpts[k++] = 2;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 8c
    newpts[k++] = 2;
    newpts[k++] = 11;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 8d
    newpts[k++] = 11;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("new_points="+new_points);
    System.out.println("new_polys="+new_polys);

    new datwrite(fileout, xyz, new_points, newpts, new_polys);
    System.out.println("file "+fileout+" written");


    k = 0; // only external triangles
    new_polys = 0;
    System.out.println("tetra 1");
    newpts[k++] = 3;
    newpts[k++] = 5;
    newpts[k++] = 9;
    newpts[k++] = 10;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 1b
    newpts[k++] = 5;
    newpts[k++] = 9;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 1c
    newpts[k++] = 5;
    newpts[k++] = 10;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 2");
    newpts[k++] = 3;
    newpts[k++] = 7;
    newpts[k++] = 10;
    newpts[k++] = 11;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 2b
    newpts[k++] = 7;
    newpts[k++] = 10;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 2c
    newpts[k++] = 7;
    newpts[k++] = 11;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 3");
    newpts[k++] = 3;
    newpts[k++] = 8;
    newpts[k++] = 10;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 3b
    newpts[k++] = 8;
    newpts[k++] = 10;
    newpts[k++] = 11;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 3c
    newpts[k++] = 8;
    newpts[k++] = 11;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 4");
    newpts[k++] = 3;
    newpts[k++] = 3;
    newpts[k++] = 12;
    newpts[k++] = 9;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 4b
    newpts[k++] = 3;
    newpts[k++] = 12;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 4c
    newpts[k++] = 3;
    newpts[k++] = 9;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 5 ");
    newpts[k++] = 3;
    newpts[k++] = 6;
    newpts[k++] = 10;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 5b
    newpts[k++] = 6;
    newpts[k++] = 10;
    newpts[k++] = 9;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 5c
    newpts[k++] = 6;
    newpts[k++] = 9;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 6 ");
    newpts[k++] = 3;
    newpts[k++] = 1;
    newpts[k++] = 9;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 6b
    newpts[k++] = 1;
    newpts[k++] = 9;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 6c
    newpts[k++] = 1;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 7 ");
    newpts[k++] = 3;
    newpts[k++] = 4;
    newpts[k++] = 13;
    newpts[k++] = 11;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 7b
    newpts[k++] = 4;
    newpts[k++] = 13;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 7c
    newpts[k++] = 4;
    newpts[k++] = 11;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 8 ");
    newpts[k++] = 3;
    newpts[k++] = 2;
    newpts[k++] = 11;
    newpts[k++] = 12;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3;  // 8b
    newpts[k++] = 2;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);
    newpts[k++] = 3; // 8c
    newpts[k++] = 2;
    newpts[k++] = 11;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("new_points="+new_points);
    System.out.println("new_polys="+new_polys);

    new datwrite(fileout1, xyz, new_points, newpts, new_polys);
    System.out.println("file "+fileout+" written");


    k = 0; // only internal triangles
    new_polys = 0;
    System.out.println("tetra 1");
    newpts[k++] = 3; // 1d
    newpts[k++] = 9;
    newpts[k++] = 10;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 2");
    newpts[k++] = 3; // 2d
    newpts[k++] = 10;
    newpts[k++] = 11;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 3");
    newpts[k++] = 3; // 3d
    newpts[k++] = 10;
    newpts[k++] = 11;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 4");
    newpts[k++] = 3; // 4d
    newpts[k++] = 9;
    newpts[k++] = 12;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 5 ");
    newpts[k++] = 3; // 5d
    newpts[k++] = 9;
    newpts[k++] = 10;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 6 ");
    newpts[k++] = 3; // 6d
    newpts[k++] = 9;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 7 ");
    newpts[k++] = 3; // 7d
    newpts[k++] = 11;
    newpts[k++] = 12;
    newpts[k++] = 13;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("tetra 8 ");
    newpts[k++] = 3; // 8d
    newpts[k++] = 11;
    newpts[k++] = 12;
    newpts[k++] = 14;
    new_polys = new_polys + 1;
    System.out.println("tetra "+(new_polys)+"   "+newpts[k-4]+", "+newpts[k-3]+" ,"+newpts[k-2]+" ,"+newpts[k-1]);

    System.out.println("new_points="+new_points);
    System.out.println("new_polys="+new_polys);

    new datwrite(fileout2, xyz, new_points, newpts, new_polys);
    System.out.println("file "+fileout+" written");
    
    System.out.println("stella.java finished");
    
  } // end stella
  
  public static void main(String[] args)
  {
    new stella(); // writes stella.dat
  }
} // end class  stella.java