/*
 * @(#)Kul.java	 
 *  PYL septembre 2006!!!!!!!!!!!
 *
 * This software is not designed or intended for use in on-line control of
 * aircraft, air traffic, aircraft navigation or aircraft communications; or in
 * the design, construction, operation or maintenance of any nuclear
 * facility. Licensee represents and warrants that it will not use or
 * redistribute the Software for such purposes.
 */
//javac guiKnudt.java Kul.java kolor.java  ; java guiKnudt
//javac guiKnudt.java Kul.java kolor.java  ; open index.html       


import java.util.*;
import java.awt.*;
import java.applet.Applet;
import java.awt.event.*;
import javax.swing.*;
import java.io.*;
 
 
class Kul extends Thread    
 {
 public double t;
 public double phase;
 private double dt=0;
 public double dx;
 public double visc;
 public double frein;
 
 public double p[];
 public double QPp[];
 public double QPm[];
 public double QUp[];
 public double QUm[];
 public double po[];
 public double u[];
 public double uo[];
 public double uoo[];
 public double ustet[];
 public double p0=0;
 public double p1=0;
 public double ps=0;
 public double u0=0;
 public double us=0;
 public double u1=0;
 public double q=0;
 boolean duq  = false;
 boolean libre= false;
 boolean ferme= false;
 boolean uimp =true;

 public double omega;
 public int n;
 public int nstet;
 public int cours=0;
 boolean pleaseWait = true;
 
private int property;
public int getProperty  (         ){ return property; }
public void setProperty (int value){ property = value; }
public void setdt (double value){ dt = value;System.out.println("------------>>>dt"+dt); }
public double getdt  (         ){/*System.out.println("-------------<<<dt"+dt);*/return dt; }

public void Kul( )
 {    
 }
 
 public void zzz()
 {
 p0=0;
 p1=0;
 u0=0;
 u1=0;
 us=0;
 q=0;
 libre = false;
 ferme = false;
 duq =false;
 uimp = true;
 }
 
 public void cl00( )
 {    
 zzz();
 }

 public void cl10( )
 {
 zzz();    
 u0=1; 
 uimp = true;
 }
 
public void cl01( )
 {    
 zzz();
 u1=1;
 uimp = true;
 } 
public void cl11( )
 {    
 zzz();
 u0=1;
 u1=1;
 uimp = true;
 } 
 
 public void cll( )
 {  
 zzz(); 
 us=0.05; 
 libre=true;
 uimp = true;
 } 
 
 public void clf( )
 {    
 zzz(); 
 us=0.05; 
 ferme=true;
 uimp = true;
 } 
 
  public void clpl( )
 {  
 zzz(); 
 ps=0.05; 
 libre=true;
 uimp = false;
 System.out.println("clpl"); 
 } 
 
 public void clpf( )
 {    
 zzz(); 
 ps=0.05; 
 ferme=true;
 uimp =false;
 } 
 
 public void clz( )
 {   
 zzz(); 
 us=0.05;  
 duq = true;
 uimp = true;
  } 
 
 
  
 
 
  
 public void prems()
 { int i;
    //pleaseWait = true;
    t=0;
	phase=0;
	n=250;
	dx=.004;
	omega=0;
	nstet=300;
   System.out.println(" initialize "); 
   po = new double[n+1];
   p  = new double[n+1];
   uo = new double[n+1];
   uoo= new double[n+1];
   QPp= new double[n+1];
   QPm= new double[n+1];
   QUp= new double[n+1];
   QUm= new double[n+1];
   u  = new double[n+1];
   ustet= new double[nstet+1];
   for(i=0;i<=n;i++)
		{po[i]=0;
		 p[i]=po[i];
		 u[i]=0;uo[i]=0;uoo[i]=0;}
    }
 
   public void run() {
   int i,it=0,itt=0;
   double lax;
   
   System.out.println("params " + dt + " " +dx +" dt/dx="+(dt/dx)+ " " + visc);
  
	while (!interrupted()) {
	//System.out.println(pleaseWait);// visc=dx/10.*1.00;
	     
	      lax=dt/2;
	      visc=dx/5;;
	 if(dt>0){ 
								 	
	   for(it=1;it<50000;it++){
		phase=phase+omega*dt;
		
		if(it%5000==0){
		  ustet[nstet]= u[n/2];
          for(itt=1;itt<nstet;itt++){ustet[itt]=ustet[itt+1];}}

		t=t+dt; 
		if(uimp) {
		u[0]=u0+us*Math.sin(phase); } //System.out.println("params " + t + " " +T[0]);
		else{
		p[0]=p0+ps*Math.sin(phase);	
			}
		//uo[0]=u[0];
		  
		for(i=1;i<n;i++)
		 {
		  //if (Math.abs(uo[i+1]-uo[i-1]) >.1 ) {visc=.005;}else{ visc=.000;}	
		  QPp[i]= (po[i+1]+po[i]  )/2 - (lax+visc)*(uo[i+1]-uo[i]  )/dx ;
 	 	  QPm[i]= (po[i]  +po[i-1])/2 - (lax+visc)*(uo[i]  -uo[i-1])/dx ;}
        
        for(i=1;i<n;i++)
		 {QUp[i]= (uo[i+1]+uo[i]  )/2 - lax*(po[i+1]-po[i]  )/dx ;
 	 	  QUm[i]= (uo[i]  +uo[i-1])/2 - lax*(po[i]  -po[i-1])/dx ;}
 

		 
		 for(i=1;i<n;i++)
		{u[i]= uo[i] + dt * (- (QPp[i]-QPm[i])/dx - frein*uo[i] ) ;
		 p[i]= po[i] + dt * (- (QUp[i]-QUm[i])/dx);} 
		
		// CONDITIONS AUX LIMITES
		if(uimp) {
		p[0]= po[0] - dt * (uo[1]-uo[0]  )/dx ;}
		else{
		u[0]= uo[0] - dt * ((po[1]-po[0]  )/dx - frein*uo[0]);	}
		p[n]= po[n] - dt * (uo[n]-uo[n-1])/dx;
		u[n]=u1;
										
		if(libre) {
		p[n]=0;
		u[n]= uo[n] + dt * (-frein*uo[n]  - (po[n]-po[n-1])/dx ) ;}
		
		if(ferme){
		p[n]= po[n] - dt * (uo[n]-uo[n-1])/dx;
		u[n]=0;}

		if(duq){
		 u[n]= uo[n] + dt * (  - (uo[n]-uo[n-1])/dx ) ;
		 p[n]= po[n] + dt * (  - (po[n]-po[n-1])/dx ) ;} 	
		
        for(i=0;i<=n;i++)
		 {uo[i]= u[i] ;
		  po[i]= p[i] ;}
		}
 }
		try {
		Thread.sleep(10);
	    } catch (InterruptedException e) {}
				// Do work
                // Check if should wait
                synchronized (this) {
                    while (pleaseWait) {
                        try {
                            wait();
                        } catch (Exception e) {}
                    }
                }
    }
                // Do work http://javaalmanac.com/egs/java.lang/PauseThread.html?l=rel
}
		
 }