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

public class plot_s12 extends Frame{

	static int ND=128;					//組織１辺の分割数
	static int nd=ND, ndm=ND-1;
	static int width;					// Window全体の幅
	static int height;					// Window全体の高さ
	static int xwidth;					// 描画領域の幅
	static int yheight;					// 描画領域の高さ
	static int insetx;					// Windowの枠の幅（左右および下）
	static int insety;					// Windowの枠の幅（上）
	static double [][] s1h=new double[ND][ND];		//組織内の秩序変数デ−タs1配列
	static double [][] s2h=new double[ND][ND];		//組織内の秩序変数デ−タs2配列

	//グローバル変数化
	static double time1;	//計算時間（カウント数）
	Image buff;
	static Graphics bg, g;

	//*************** コンストラクタ ****************************
	public plot_s12(){
		xwidth=400; yheight=400; 			//描画画面の横と縦の長さ（ピクセル単位）
		insetx=4; insety=30;					//描画画面のふちの長さ
		width=xwidth+insetx*2;  			//描画Window全体の横の長さ
		height=yheight+insetx+insety; 	//描画Window全体の縦の長さ
		setSize(width, height);				//描画Windowのセット
		setBackground(Color.white); 		//描画Windowの描画部の色を白に設定
		setVisible(true); 						//描画Windowを見えるようにする
		addWindowListener(new WindowAdapter(){ 
			public void windowClosing(WindowEvent e){ System.exit(0); }
												//Windowを閉じる時の操作（Windowの右上×の設定）
		});
	}

//************** main *****************************
	public static void main(String[] args) throws Exception{//例外処理は行わない

		plot_s12 prog=new plot_s12();

		int  i, j, k, l; 						//整数
		double sumc;
    String str_data, str_time1;

		BufferedReader infile=new BufferedReader(new FileReader("test.dat"));
		while( (str_time1=infile.readLine())!=null ){
			time1=new Double(str_time1).doubleValue();
			for(i=0;i<=ndm;i++){
				for(j=0;j<=ndm;j++){
					str_data=infile.readLine();
					s1h[i][j]=new Double(str_data).doubleValue();
					str_data=infile.readLine();
					s2h[i][j]=new Double(str_data).doubleValue();
				}
			}
			prog.update_draw(g);
			if(time1==1800.0){ prog.datsave(); }		//カウント数を指定して場を保存
		}

		infile.close();

//----------------------------------------------------------------
		System.out.printf("\n 終了しました。グラフの右上×をクリックして終了してください。\n");

	}//main

// *******************************************************************************
	public void update_draw( Graphics g ){ g=getGraphics(); paint(g); }

// **** [phase fieldの描画] *******************************************************
	public void paint(Graphics g){
		//g.clearRect(0, 0, width, height);//Windowをクリア

		int i, j, ii, jj;
		int icol, icol_r, icol_g, icol_b;
		double c_r, c_g, c_b;
		int ixmin=0, iymin=0, igx, igy, irad0;
		int ixmax=xwidth, iymax=yheight;
		double c, x, xmax, xmin, y, ymax, ymin, rad0;

		xmin=0.; xmax=1.; 						//横軸の最小値、最大値
		ymin=0.; ymax=1.; 						//縦軸の最小値、最大値
		rad0=1.0/(double)nd/2.0; 				//差分ブロックの長さの半分
		irad0=1+(int)( ((double)ixmax-(double)ixmin)/(xmax-xmin)*rad0 ); 	//rad0のピクセル化

		System.out.printf("%f \n", time1); 	//計算の繰返し回数を標準入出力に表示

		for(i=0;i<=nd;i++){
			for(j=0;j<=nd;j++){
				//濃度場の位置座標（実際の値）
				x=1.0/(double)nd*(double)i+rad0;
				y=1.0/(double)nd*(double)j+rad0;
				//濃度場の位置座標（スクリーン座標に変換）
				igx=(int)( ((double)ixmax-(double)ixmin)*(x-xmin)/(xmax-xmin)+(double)ixmin );
				igy=(int)( ((double)iymax-(double)iymin)*(y-ymin)/(ymax-ymin)+(double)iymin );

				//個々の差分ブロックの濃度値
				ii=i; jj=j;
				if(i==nd){ii=0;} if(j==nd){jj=0;}			//周期的境界条件

				c_r= s1h[ii][jj];							//s1を赤
				c_g= s2h[ii][jj]; 							//s2を赤
				c_b=1.0-c_r-c_g; 						//変態前の相を青
				if(c_r>1.0){c_r=1.0;}  if(c_r<0.0){c_r=0.0;}
				if(c_g>1.0){c_g=1.0;}  if(c_g<0.0){c_g=0.0;}
				if(c_b>1.0){c_b=1.0;}  if(c_b<0.0){c_b=0.0;}

				icol_r=(int)(255.*c_r);  	icol_g=(int)(255.*c_g);  icol_b=(int)(255.*c_b);		//256階層に変換
				g.setColor(new Color(icol_r, icol_g, icol_b)); 									//差分ブロックの色を設定
				g.fillRect(insetx+igx-irad0,insety+igy-irad0, irad0*2, irad0*2);					//個々の差分ブロック描画
			}
		}
	}

//**** [デ−タの保存] ************************************
	private void datsave() throws Exception{
		int	i, j;

		//保存ファイル名をtest.datとする。
		PrintWriter outfile= new PrintWriter(							//ファイルのオープン
			new BufferedWriter(new FileWriter("test_selected.dat", true)) );	//追記

		outfile.println(time1);											//カウントの書き込み
		for(i=0;i<=ndm;i++){
			for(j=0;j<=ndm;j++){
				//outfile.printf("%e  %e  ",s1h[i][j], s2h[i][j]);									//phase fieldの書き込み
				outfile.println(s1h[i][j]);									//phase fieldの書き込み
				outfile.println(s2h[i][j]);									//phase fieldの書き込み
			}
		}
		outfile.close();													//ファイルのクローズ
	}

//****************************************************************
}