// humpc.java       compile    javac -cp . humpc.java
//                  execute    java  -cp . humpc # writes humpc2.dat

public class humpc
{
  int nr = 8;  // not including 0.0 pt 1 for poly
  double rmin = 0.0;
  double rmax = 1.0;
  double r = 0.0;
  double dr = (rmax-rmin)/(double)nr;
  int na = 16;  // angles
  double Pi = Math.PI;
  double ang=0.0;
  double dang = 2.0*Pi/(double)(na);
  double xyz[][] = new double[4000][3];
  int npt = 0;       // xyz index and count
  int poly[] = new int[12000];  // for lines
  int poly2[] = new int[12000]; // for triangles
  int npy = 0;       // poly count
  int k = 0;         // poly index
  double zmin = 0.0; // all values plus this
  double zmax = 1.0; // max at center
  double x, y, z, rp;
  int ii, jj, ptb;
  String fileout = "humpc.dat";
  String fileout2 = "humpc2.dat";
    
  public humpc()
  {
    System.out.println("humpc.java running");
    System.out.println("rmin ="+rmin+"  rmax="+rmax);
    System.out.println("radius nr="+nr+"  dr="+dr);
    System.out.println("ang 0 ="+0.0+"  angmax="+(2.0*Pi));
    System.out.println("angle na="+na+"  dang="+dang);
    System.out.println(" ");
    System.out.println("generated xyz");
    // unique center top
    x = 0.0;
    y = 0.0;
    z = zmax;
    xyz[npt][0] = x;
    xyz[npt][1] = y;
    xyz[npt][2] = z;
    npt = npt+1;
    System.out.println("npt="+npt+"  x="+x+"  y="+y+"  z="+z);
    r = dr;
    for(int i=0; i<nr; i++)
    {
      ang = 0.0;
      System.out.println("r="+r+"  ang="+ang);
      for(int j=0; j<na; j++)
      {
	x = r*Math.cos(ang);
	y = r*Math.sin(ang);
	z = humpc2(r);
        xyz[npt][0] = x;
        xyz[npt][1] = y;
        xyz[npt][2] = z;
        npt = npt+1;
        System.out.println("npt="+npt+"  x="+x+"  y="+y+"  z="+z);
	ang = ang + dang;
      } // end j
      r = r + dr;
    } // end i

    System.out.println("generate poly, lines");
    for(int j=0; j<na; j++)
    {
      poly2[k++] = 2; // out from center
      poly2[k++] = 1;
      poly2[k++] = j+2;
    } // end j

    ii = 2;
    for(int i=0; i<nr-1; i++)
    {
      for(int j=0; j<na; j++)
      {
        poly2[k++] = 2; // out from center
        poly2[k++] = ii+j;
        poly2[k++] = ii+j+na;
      } // end j
      ii = ii + na;
    } // end i

    ii = 2;
    for(int i=0; i<nr; i++)
    {
      for(int j=0; j<na; j++)
      {
        poly2[k++] = 2; // around
        poly2[k++] = ii+j;
        poly2[k++] = ii+j+1;
      } // end j
      poly2[k-1] = ii;
      ii = ii + na;
    } // end i

    npy = k/3;
    System.out.println("npt xyz="+npt+",  npy poly2="+npy);
    new datwrite(fileout2, xyz, npt, poly2, npy);
    System.out.println("datwrite wrote "+fileout2);

    System.out.println("generate poly, triangles, later");
    // unique center bottom
    x = 0.0;
    y = 0.0;
    z = zmin;
    xyz[npt][0] = x;
    xyz[npt][1] = y;
    xyz[npt][2] = z;
    npt = npt+1;
    System.out.println("npt="+npt+"  x="+x+"  y="+y+"  z="+z);
    ptb = npt;
    k = 0;
    
    for(int j=0; j<na; j++)
    {
      poly[k++] = 3; // tri from top center
      poly[k++] = 1;
      poly[k++] = j+2;
      poly[k++] = j+3;
    } // end j
    poly[k-1] = 2;

    jj = (nr-1)*na;
    for(int j=0; j<na; j++)
    {
      poly[k++] = 3; // tri from bottom center
      poly[k++] = ptb;
      poly[k++] = jj+j+2;
      poly[k++] = jj+j+3;
    } // end j
    poly[k-1] = jj+2;

    ii = 2;
    for(int i=0; i<nr-1; i++)
    {
      for(int j=0; j<na; j++)
      {
        poly[k++] = 3; // tri 
        poly[k++] = ii+j;
        poly[k++] = ii+j+na;
        poly[k++] = ii+j+1+na; // k-5
        poly[k++] = 3;  
        poly[k++] = ii+j+1+na; // k-3
        poly[k++] = ii+j+1;    // k-2
        poly[k++] = ii+j;      // k-1
      } // end j
      poly[k-5] = ii+na;
      poly[k-3] = ii+na;
      poly[k-2] = ii;
      poly[k-1] = ii+na-1;
      System.out.println("poly k-8 "+poly[k-8]+" "+poly[k-7]+"  "+
			             poly[k-6]+" "+poly[k-5]);
      System.out.println("poly k-4 "+poly[k-4]+" "+poly[k-3]+"  "+
			             poly[k-2]+" "+poly[k-1]);
      ii = ii + na;
    } // end i
    
    npy = k/4;
    System.out.println("npt xyz="+npt+",  npy poly="+npy);
    new datwrite(fileout, xyz, npt, poly, npy);
    System.out.println("datwrite wrote "+fileout);

    
    System.out.println("humpc.java finished");
  }

  double humpc2(double rr)
  {
    double zz = 1.0-rr*rr;
	// (rmax-rr)*(rmax-rr)+zmin;
    return zz;
  }  // end humpc2
    
  public static void main (String[] args)
  {
    new humpc();
  }
} // end humpc.java