// w0015vert.java  wing in wind tunnel, non uniform grid
//                 data read from w0015a22.grid produced by w0015grid.java
//                 order=5 e.g for size computation
//                 wing has a chord of 1 meter,

import java.awt.*;
import java.awt.event.*;
import java.io.*;

// needs  read_grid.java
// and a data file  w0015a22.grid or equavelent

public class w0015vert extends Frame
{
  double xmin = -1.0;
  double xmax = 2.0;
  double xoff = 340.0;
  double yoff = 500.0;
  double ymin = -1.25;
  double ymax = 1.15;
  double scale = 320.0; // for plotting
  double xminw; // from data read
  double xmaxw;
  double yminw;
  double ymaxw;
  double xadj = -0.5; // approx center for output 
  double grid = 0.01; // grid spacing, may be computed 
  int width = 1000;
  int height = 1000; // upper part for title 
  int ixbound[] = new int[100]; // output boundary
  int jybound[] = new int[100];
  double xg[] = new double[100]; // output grid
  double yg[] = new double[100];
  int nxg;
  int nyg;
  int nbound;
  int ndel = 10; // grids outside wing
  String fname = "w0015a22.grid";
  int neqn; // number of equations
  int nDOF; // number of degrees of freedom, one independent variable
  int ntot; // independent plus squared variables, max

  int order = 5; // fifth order nonuniform solution
 
  w0015vert(String gname, int tndel)
  {
    System.out.println("w0015vert.java running");
    fname = gname;
    ndel = tndel;
    System.out.println("grid file "+fname+"  steps "+ndel);
    build_vert();  

    neqn = 3*nDOF; // also number of linear unknowns
    ntot = neqn + order*order*nDOF; // maximum columns with order u*v terms
    System.out.println("neqn="+neqn+", ntot="+ntot+", A is "+
                       (neqn*ntot*8)+" bytes");

    setTitle("w0015vert");
    setSize(width,height);
    setBackground(Color.white);
    setForeground(Color.black);
    addWindowListener(new WindowAdapter()
    {
      public void windowClosing(WindowEvent e)
      {
        System.exit(0);
      }
    });
    setVisible(true);
    this.addMouseListener (new mousePressHandler());

  } // end constructor w0015vert

  void build_vert()
  {
    double x, y;

    // read in a wing
    read_grid G = new read_grid(fname);
    nxg = G.read_xgrid(xg);
    nyg = G.read_ygrid(yg);
    nbound = G.read_bound(ixbound, jybound);
    for(int k=0; k<nbound; k++) System.out.println("ixb["+k+"]="+ixbound[k]+
                                                   ", jyb="+jybound[k]);

    // offset and build surrounding wind tunnel grid
    if(nxg>0)
    {
      xminw = xg[0];
      xmaxw = xg[nxg-1];
      yminw = yg[0];
      ymaxw = yg[nyg-1];
    }
    else
    {
      xminw = 0.0;
      xmaxw = 0.0;
      yminw = 0.0;
      ymaxw = 0.0;
    }
    System.out.println("xminw="+xminw+", xmaxw="+xmaxw+
                       ", yminw="+yminw+", ymaxw="+ymaxw);
    // xmin, xmax, ymin, ymax set manually for now
    System.out.println("xmin="+xmin+", xmax="+xmax+
                       ", ymin="+ymin+", ymax="+ymax);
    double xmind = -(xmin-xminw)/(double)(ndel);
    double xmaxd = (xmax-xmaxw)/(double)(ndel);
    double ymind = -(ymin-yminw)/(double)(ndel);
    double ymaxd = (ymax-ymaxw)/(double)(ndel);
    System.out.println("xmind="+xmind+", xmaxd="+xmaxd+
                       ", ymind="+ymind+", ymaxd="+ymaxd);
    // move wing and insert additional grid
    for(int i=nxg-1; i>=0; i--) xg[i+ndel] = xg[i]; 
    for(int i=nyg-1; i>=0; i--) yg[i+ndel] = yg[i]; 
    for(int i=0; i<nbound; i++) ixbound[i] = ixbound[i]+ndel;
    for(int i=0; i<nbound; i++) jybound[i] = jybound[i]+ndel;
    x = xmin;
    for(int i=0; i<ndel; i++) {xg[i]=x; x+=xmind;};
    x = xmaxw+xmaxd;
    for(int i=ndel+nxg; i<ndel+ndel+nxg; i++) {xg[i]=x; x+=xmaxd;};
    nxg = nxg+ndel+ndel;
    for(int i=0; i<nxg; i++) System.out.println("xg["+i+"]="+xg[i]);

    y = ymin;
    for(int i=0; i<ndel; i++) {yg[i]=y; y+=ymind;};
    y = ymaxw+ymaxd;
    for(int i=ndel+nyg; i<ndel+ndel+nyg; i++) {yg[i]=y; y+=ymaxd;};
    nyg = nyg+ndel+ndel;
    for(int j=0; j<nyg; j++) System.out.println("yg["+j+"]="+yg[j]);
  } // end build_vert
    
  int bound(int i, int j) // 0 for DOF, 1 for external, 2 for wing, 3 inside
  {
    int jval;

    if(i==0 || i==nxg-1) return 1;
    if(j==0 || j==nyg-1) return 1;
    if(i<ndel || j<ndel) return 0;
    if(i>nxg-ndel || j>nyg-ndel) return 0;
    for(int k=0; k<nbound; k++)
      if(i==ixbound[k] && j==jybound[k]) return 2;
    jval=0;
    for(int k=0; k<nbound; k++)
    {
      if(i==ixbound[k])
      {
	 if(jval==0) jval = jybound[k];
         else
	 {
	   // may not be any
           if((j<jval&&j>jybound[k]) || (j<jybound[k]&&j>jval)) return 3;
           jval = 0;
	 }
      }
    }
    return 0;
  } // end bound

  class mousePressHandler extends MouseAdapter
  {
    public void mousePressed (MouseEvent e)
    {
      int x, y, b;
      x = e.getX();
      y = e.getY();
      b = e.getButton();
      if(b==3)
      {
      }
      requestFocus();
      System.out.println("at x="+x+"   y="+y+"   b="+b); // debug print
      repaint();
    } // end mousePressed
  } // end class mousePressedHandler

  public void paint(Graphics g)
  {
    int ix,jy;
    g.setColor(Color.black);
    g.drawString(fname, 20, 50);
    g.setColor(Color.red);
    for(int i=0; i<nxg; i++)
    {
      for(int j=0; j<nyg; j++)
      {
        if(i==0 && j== 0)
           System.out.println("in paint, x="+((int)(scale*xg[i]+xoff))+
                              ", y="+(height-(int)(scale*yg[j]+yoff)));
        if(bound(i,j)==0) g.setColor(Color.red);
        if(bound(i,j)==1) g.setColor(Color.black);
        if(bound(i,j)==2) g.setColor(Color.blue);
        if(bound(i,j)==3) g.setColor(Color.yellow);
 
        g.drawOval((int)(scale*xg[i]+xoff), 
	  	   height-(int)(scale*yg[j]+yoff),
                   1, 1);
      }
    }
  } // end paint

  public static void main(String args[])
  {
      // read in name of .grid file
      String gname = "w0015a22.grid";
      if(args.length>0)
	  gname = args[0];
      int tndel = 10;
      if(args.length>1)
	  tndel = Integer.parseInt(args[1]);

      new w0015vert(gname, tndel);
    System.out.println("w0015vert.java ends");
  }
} // end class w0015vert