// shspe2.java  plot 3d data  Utah Graphics  to shape2.dat
//                 javac -cp . shape2.java
//             _   java  -cp . shape2 > shape2_java.out
// |\____/|   / \
// | ____ |  |   |
// |/    \|   \_/

import java.io.*;

public class shape2
{
  boolean debug = true;
  int pts, pt, k;
  datread DR; // = new datread(filename); // defines data structures
  int new_points = 0; // for each x,y,z point kept
  int new_polys = 0; // changed for saved new_points z >= 0.0
  double xyz[][] = new double[10000][3]; // x, y, z for output
  int newpts[] = new int[40000]; // 3, xyz1, xyz2, xyz3 for output
  String fileout = "shape2.dat";
  double ang, dang;
  int n = 32;
  int nz = 19;
  int np = 0;
  double Pi =  Math.PI;
  double x, y, z, dz;
  double rxy[] = new double[19];     

  public shape2()
  {
    System.out.println("shape2.java running");
    new_polys = 0;
    new_points = 0;

    rxy[0]=2.0;     // radius at z = 0.0
    rxy[1]=1.7369;
    rxy[2]=1.4517;
    rxy[3]=1.2022;
    rxy[4]=0.9894;
    rxy[5]=0.8141;
    rxy[6]=0.6771;
    rxy[7]=0.5788;
    rxy[8]=0.5197;
    rxy[9]=0.5;     // radius at z = 2.5
    rxy[10]=0.5197;
    rxy[11]=0.5788;
    rxy[12]=0.6771;
    rxy[13]=0.8141;
    rxy[14]=0.9894;
    rxy[15]=1.2022;
    rxy[16]=1.4517;
    rxy[17]=1.7369;
    rxy[18]=2.0;    // radius at z = 5.0

    System.out.println(" ");
    System.out.println("generate new points ");

    z = 0.0;
    dz = 5.0/(double)(nz-1);
    dang = Pi/(double)(n/2);
    for(int j=0; j<nz; j++)
    {
      ang = 0.0;
      dang = Pi/(double)(n/2);	
      for(int i=0; i<n; i++) // top
      {
        x = rxy[j]*Math.cos(ang);
        y = rxy[j]*Math.sin(ang);
        xyz[new_points][0] = x;
        xyz[new_points][1] = y;
        xyz[new_points][2] = z;
        new_points = new_points +1;
        ang = ang + dang;
      } // end i
      z = z + dz;
      System.out.println("j="+j+" xyz="+xyz[new_points-1][0]+" "+
			 xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end j
    System.out.println("new_points= "+new_points);

    System.out.println(" ");
    System.out.println("generate new polys ");
    
    k = 0;
    np = 0;
    for(int j=0; j<nz-1; j++)
    {
      for(int i=0; i<n; i++) // top and 1
      {
        if(i==n-1) // close
        {
          newpts[k++] = 3; // tri 1
          newpts[k++] = np+i+1;
          newpts[k++] = np+i+1+n;
          newpts[k++] = np+i+2; // back
          newpts[k++] = 3; // tri 2 square 1
          newpts[k++] = np+i+2; // back
          newpts[k++] = np+i+2-n; // back
          newpts[k++] = np+i+1;
        }
        else
        {
          newpts[k++] = 3; // tri 1
          newpts[k++] = np+i+1;
          newpts[k++] = np+i+1+n;
          newpts[k++] = np+i+2+n;
          newpts[k++] = 3; // tri 2 square 1
          newpts[k++] = np+i+2+n;
          newpts[k++] = np+i+2;
          newpts[k++] = np+i+1;
        }
        System.out.println("j="+j+" newpts[k-8]="+newpts[k-8]+" "+newpts[k-7]+
	  		            " "+newpts[k-6]+" "+newpts[k-5]); 
        System.out.println("j="+j+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
			            " "+newpts[k-2]+" "+newpts[k-1]); 
      } // end i
      np = np + n;
    } // end j
    
    System.out.println("new_points="+new_points+"  np="+np);
    new_polys = k/4;
    System.out.println("new_polys="+new_polys);

    xyz[new_points][0] = 0.0;
    xyz[new_points][1] = 0.0;
    xyz[new_points][2] = 0.0;
    new_points = new_points +1;
    for(int i=0; i<n; i++) // fill botom circle
    {
      if(i==n-1) // close
      {
        newpts[k++] = 3;
        newpts[k++] = i+1;
        newpts[k++] = i+2-n;
        newpts[k++] = new_points;
      }
      else
      {
        newpts[k++] = 3; 
        newpts[k++] = i+1;
        newpts[k++] = i+2;
        newpts[k++] = new_points;
      }
      System.out.println("bot newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
			            " "+newpts[k-2]+" "+newpts[k-1]); 
    } // end i	


    xyz[new_points][0] = 0.0;
    xyz[new_points][1] = 0.0;
    xyz[new_points][2] = 5.0;
    new_points = new_points +1;
    np = n*(nz-1);
    System.out.println("np bot="+np+" n="+n+" nz="+nz+" new_points="+new_points); 
    
    for(int i=0; i<n; i++) // fill botom circle
    {
      if(i==n-1) // close
      {
        newpts[k++] = 3;
        newpts[k++] = np+i+1;
        newpts[k++] = np+i+2-n;
        newpts[k++] = new_points;
      }
      else
      {
        newpts[k++] = 3; 
        newpts[k++] = np+i+1;
        newpts[k++] = np+i+2;
        newpts[k++] = new_points;
      }
      System.out.println("bot newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
			            " "+newpts[k-2]+" "+newpts[k-1]); 
    } // end i	
    new_polys = k/4;


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

    new datwrite(fileout, xyz, new_points, newpts, new_polys);
    System.out.println("shape2.java finished");

  } // shape2 constructor

  public static void main(String[] args)
  {
    new shape2(); // writes shape2.dat
  }
} // end class  shape2.java