// ==================================================
// 波モデルの頂点座標を更新する処理
// ==================================================

void updateModelVertex() {
	
	for( int yi=0; yi<N; yi++ ) {
		for( int xi=0; xi<N; xi++ ) {
			
			// 以下、メッシュモデル内の [yi][xi] 番目の頂点（格子点）に注目している
			
			double x = X_MIN + xi * DX; // 注目頂点の x 座標値
			double y = Y_MIN + yi * DY; // 注目頂点の y 座標値
			
			double rA = sqrt( (x-xA)*(x-xA) + (y-yA)*(y-yA) ); // 波源Aからの距離
			double rB = sqrt( (x-xB)*(x-xB) + (y-yB)*(y-yB) ); // 波源Bからの距離
			
			// 円形波の式から、波源Aからの波（波A）と波源Bからの波（波B）のZ値をそれぞれに求める
			double zA = amplitudeA * sin( 2.0 * PI * ( t/periodA - rA/wavelengthA ) ); // 波AのZ値
			double zB = amplitudeB * sin( 2.0 * PI * ( t/periodB - rB/wavelengthB ) ); // 波BのZ値
			
			// 頂点座標値配列に値を格納
			waveVertex[yi][xi][X] = x;
			waveVertex[yi][xi][Y] = y;
			waveVertex[yi][xi][Z] = zA + zB; // 合成波のZ値は、重ね合わせの原理にり波Aと波Bの和
		}
	}
	
	// モデルに頂点座標値配列を渡して形状を更新
	setModelVertex(waveModel, waveVertex, QUADRANGLE_GRID);
}