#include "Alloc_Save_Load_Field_Thread.h"
#include "fdtd_3d_alloc.h"
#include <sys/stat.h> //for mkdir()
#include "Matrix.h"

#include <iostream>
using std::cout;
using std::endl;

////////////////////////////////////////////////////////////////////
//Collect the information - Ex field
////////////////////////////////////////////////////////////////////
void Init_Alloc_FxGxEx(Allocate_E_comp *Allocate_FGE_x, double ***&Fx_1,  double ***&Gx_1,
				  double ***&Fx_2,  double ***&Gx_2, double ***&Fx_3,  double ***&Gx_3,  
				  double ***&Fx_4,  double ***&Fx_5, double ***&Fx_6,  double ***&Fx_7,  
				  double ***&Gx_7,  double ***&Fx_8,  double ***&Gx_8,  double ***&Fx_9,
				  double ***&Gx_9,  double ***&Fx_10, double ***&Gx_10, double ***&Fx_11, 
				  double ***&Fx_12, double ***&Gx_12, double ***&Fx_13, double ***&Fx_15, 
				  double ***&Fx_16, double ***&Gx_16, double ***&Fx_17, double ***&Fx_18, 
				  double ***&Gx_18, double ***&Fx_19, double ***&Gx_19, double ***&Fx_20,
				  double ***&Gx_20, double ***&Fx_21, double ***&Gx_21, double ***&Fx_22, 
				  double ***&Fx_23, double ***&Fx_24, double ***&Fx_25, double ***&Gx_25, 
				  double ***&Fx_26, double ***&Gx_26, double ***&Fx_27, double ***&Gx_27, 
				  double ***&Ex, long nPML_x_1, long nPML_y_1, long nPML_z_1, long nPML_x_2,
				  long nPML_y_2, long nPML_z_2, long nx, long ny, long nz)
{
	long nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2 = nx - 1 - nPML_x_1 - nPML_x_2;
	long ny_MIN_nPML_y_1_MIN_nPML_y_2 = ny - nPML_y_1 - nPML_y_2;
	long nz_MIN_nPML_z_1_MIN_nPML_z_2 = nz - nPML_z_1 - nPML_z_2;

	Allocate_FGE_x->A[0]  = &Fx_1;  Allocate_FGE_x->A[1]  = &Gx_1;  Allocate_FGE_x->A[2]  = &Fx_2;  
	Allocate_FGE_x->A[3]  = &Gx_2;  Allocate_FGE_x->A[4]  = &Fx_3;  Allocate_FGE_x->A[5]  = &Gx_3;
	Allocate_FGE_x->A[6]  = &Fx_4;  Allocate_FGE_x->A[7]  = &Fx_5;  Allocate_FGE_x->A[8]  = &Fx_6;
	Allocate_FGE_x->A[9]  = &Fx_7;  Allocate_FGE_x->A[10] = &Gx_7;  Allocate_FGE_x->A[11] = &Fx_8; 
	Allocate_FGE_x->A[12] = &Gx_8;  Allocate_FGE_x->A[13] = &Fx_9;  Allocate_FGE_x->A[14] = &Gx_9;  
	Allocate_FGE_x->A[15] = &Fx_10; Allocate_FGE_x->A[16] = &Gx_10; Allocate_FGE_x->A[17] = &Fx_11; 
	Allocate_FGE_x->A[18] = &Fx_12; Allocate_FGE_x->A[19] = &Gx_12; Allocate_FGE_x->A[20] = &Fx_13; 
	Allocate_FGE_x->A[21] = &Fx_15; Allocate_FGE_x->A[22] = &Fx_16; Allocate_FGE_x->A[23] = &Gx_16;
	Allocate_FGE_x->A[24] = &Fx_17; Allocate_FGE_x->A[25] = &Fx_18; Allocate_FGE_x->A[26] = &Gx_18; 
	Allocate_FGE_x->A[27] = &Fx_19; Allocate_FGE_x->A[28] = &Gx_19; Allocate_FGE_x->A[29] = &Fx_20; 
	Allocate_FGE_x->A[30] = &Gx_20; Allocate_FGE_x->A[31] = &Fx_21; Allocate_FGE_x->A[32] = &Gx_21; 
	Allocate_FGE_x->A[33] = &Fx_22; Allocate_FGE_x->A[34] = &Fx_23; Allocate_FGE_x->A[35] = &Fx_24;
	Allocate_FGE_x->A[36] = &Fx_25; Allocate_FGE_x->A[37] = &Gx_25; Allocate_FGE_x->A[38] = &Fx_26; 
	Allocate_FGE_x->A[39] = &Gx_26; Allocate_FGE_x->A[40] = &Fx_27; Allocate_FGE_x->A[41] = &Gx_27; 
	Allocate_FGE_x->A[42] = &Ex;    
	
	Allocate_FGE_x->er = 0;

	strcpy(&Allocate_FGE_x->name[0][0],"Fx_1");   strcpy(&Allocate_FGE_x->name[1][0],"Gx_1");  
	strcpy(&Allocate_FGE_x->name[2][0],"Fx_2");   strcpy(&Allocate_FGE_x->name[3][0],"Gx_2");  
	strcpy(&Allocate_FGE_x->name[4][0],"Fx_3");   strcpy(&Allocate_FGE_x->name[5][0],"Gx_3");
	strcpy(&Allocate_FGE_x->name[6][0],"Fx_4");   strcpy(&Allocate_FGE_x->name[7][0],"Fx_5");  
	strcpy(&Allocate_FGE_x->name[8][0],"Fx_6");   strcpy(&Allocate_FGE_x->name[9][0],"Fx_7");  
	strcpy(&Allocate_FGE_x->name[10][0],"Gx_7");  strcpy(&Allocate_FGE_x->name[11][0],"Fx_8"); 
	strcpy(&Allocate_FGE_x->name[12][0],"Gx_8");  strcpy(&Allocate_FGE_x->name[13][0],"Fx_9");  
	strcpy(&Allocate_FGE_x->name[14][0],"Gx_9");  strcpy(&Allocate_FGE_x->name[15][0],"Fx_10"); 
	strcpy(&Allocate_FGE_x->name[16][0],"Gx_10"); strcpy(&Allocate_FGE_x->name[17][0],"Fx_11"); 
	strcpy(&Allocate_FGE_x->name[18][0],"Fx_12"); strcpy(&Allocate_FGE_x->name[19][0],"Gx_12"); 
	strcpy(&Allocate_FGE_x->name[20][0],"Fx_13"); strcpy(&Allocate_FGE_x->name[21][0],"Fx_15"); 
	strcpy(&Allocate_FGE_x->name[22][0],"Fx_16"); strcpy(&Allocate_FGE_x->name[23][0],"Gx_16");
	strcpy(&Allocate_FGE_x->name[24][0],"Fx_17"); strcpy(&Allocate_FGE_x->name[25][0],"Fx_18");
	strcpy(&Allocate_FGE_x->name[26][0],"Gx_18"); strcpy(&Allocate_FGE_x->name[27][0],"Fx_19"); 
	strcpy(&Allocate_FGE_x->name[28][0],"Gx_19"); strcpy(&Allocate_FGE_x->name[29][0],"Fx_20"); 
	strcpy(&Allocate_FGE_x->name[30][0],"Gx_20"); strcpy(&Allocate_FGE_x->name[31][0],"Fx_21"); 
	strcpy(&Allocate_FGE_x->name[32][0],"Gx_21"); strcpy(&Allocate_FGE_x->name[33][0],"Fx_22"); 
	strcpy(&Allocate_FGE_x->name[34][0],"Fx_23"); strcpy(&Allocate_FGE_x->name[35][0],"Fx_24");
	strcpy(&Allocate_FGE_x->name[36][0],"Fx_25"); strcpy(&Allocate_FGE_x->name[37][0],"Gx_25");
	strcpy(&Allocate_FGE_x->name[38][0],"Fx_26"); strcpy(&Allocate_FGE_x->name[39][0],"Gx_26");
	strcpy(&Allocate_FGE_x->name[40][0],"Fx_27"); strcpy(&Allocate_FGE_x->name[41][0],"Gx_27"); 
	strcpy(&Allocate_FGE_x->name[42][0],"Ex");  
	
	//Fx_1
	Allocate_FGE_x->size[0][0] = nPML_x_1; 
	Allocate_FGE_x->size[0][1] = nPML_y_1; 
	Allocate_FGE_x->size[0][2] = nPML_z_1;
	//Gx_1
	Allocate_FGE_x->size[1][0] = nPML_x_1; 
	Allocate_FGE_x->size[1][1] = nPML_y_1; 
	Allocate_FGE_x->size[1][2] = nPML_z_1;
	//Fx_2
	Allocate_FGE_x->size[2][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[2][1] = nPML_y_1; 
	Allocate_FGE_x->size[2][2] = nPML_z_1;
	//Gx_2
	Allocate_FGE_x->size[3][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[3][1] = nPML_y_1; 
	Allocate_FGE_x->size[3][2] = nPML_z_1;
	//Fx_3
	Allocate_FGE_x->size[4][0] = nPML_x_2; 
	Allocate_FGE_x->size[4][1] = nPML_y_1; 
	Allocate_FGE_x->size[4][2] = nPML_z_1;
	//Gx_3
	Allocate_FGE_x->size[5][0] = nPML_x_2; 
	Allocate_FGE_x->size[5][1] = nPML_y_1; 
	Allocate_FGE_x->size[5][2] = nPML_z_1;
	//Fx_4
	Allocate_FGE_x->size[6][0] = nPML_x_1;
	Allocate_FGE_x->size[6][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2; 
	Allocate_FGE_x->size[6][2] = nPML_z_1;
	//Fx_5
	Allocate_FGE_x->size[7][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2; 
	Allocate_FGE_x->size[7][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2; 
	Allocate_FGE_x->size[7][2] = nPML_z_1;
	//Fx_6
	Allocate_FGE_x->size[8][0] = nPML_x_2; 
	Allocate_FGE_x->size[8][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2; 
	Allocate_FGE_x->size[8][2] = nPML_z_1;
	//Fx_7
	Allocate_FGE_x->size[9][0] = nPML_x_1;
	Allocate_FGE_x->size[9][1] = nPML_y_2; 
	Allocate_FGE_x->size[9][2] = nPML_z_1;
	//Gx_7
	Allocate_FGE_x->size[10][0] = nPML_x_1; 
	Allocate_FGE_x->size[10][1] = nPML_y_2; 
	Allocate_FGE_x->size[10][2] = nPML_z_1;
	//Fx_8
	Allocate_FGE_x->size[11][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2; 
	Allocate_FGE_x->size[11][1] = nPML_y_2; 
	Allocate_FGE_x->size[11][2] = nPML_z_1;
	//Gx_8
	Allocate_FGE_x->size[12][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[12][1] = nPML_y_2; 
	Allocate_FGE_x->size[12][2] = nPML_z_1;
	//Fx_9
	Allocate_FGE_x->size[13][0] = nPML_x_2;
	Allocate_FGE_x->size[13][1] = nPML_y_2; 
	Allocate_FGE_x->size[13][2] = nPML_z_1;
	//Gx_9
	Allocate_FGE_x->size[14][0] = nPML_x_2;
	Allocate_FGE_x->size[14][1] = nPML_y_2;
	Allocate_FGE_x->size[14][2] = nPML_z_1;
	//Fx_10
	Allocate_FGE_x->size[15][0] = nPML_x_1; 
	Allocate_FGE_x->size[15][1] = nPML_y_1;
	Allocate_FGE_x->size[15][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Gx_10
	Allocate_FGE_x->size[16][0] = nPML_x_1; 
	Allocate_FGE_x->size[16][1] = nPML_y_1; 
	Allocate_FGE_x->size[16][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_11
	Allocate_FGE_x->size[17][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[17][1] = nPML_y_1; 
	Allocate_FGE_x->size[17][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_12
	Allocate_FGE_x->size[18][0] = nPML_x_2; 
	Allocate_FGE_x->size[18][1] = nPML_y_1; 
	Allocate_FGE_x->size[18][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Gx_12
	Allocate_FGE_x->size[19][0] = nPML_x_2; 
	Allocate_FGE_x->size[19][1] = nPML_y_1; 
	Allocate_FGE_x->size[19][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_13
	Allocate_FGE_x->size[20][0] = nPML_x_1; 
	Allocate_FGE_x->size[20][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2;
	Allocate_FGE_x->size[20][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_15
	Allocate_FGE_x->size[21][0] = nPML_x_2;
	Allocate_FGE_x->size[21][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2; 
	Allocate_FGE_x->size[21][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_16
	Allocate_FGE_x->size[22][0] = nPML_x_1;
	Allocate_FGE_x->size[22][1] = nPML_y_2; 
	Allocate_FGE_x->size[22][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Gx_16
	Allocate_FGE_x->size[23][0] = nPML_x_1; 
	Allocate_FGE_x->size[23][1] = nPML_y_2; 
	Allocate_FGE_x->size[23][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_17
	Allocate_FGE_x->size[24][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[24][1] = nPML_y_2; 
	Allocate_FGE_x->size[24][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_18
	Allocate_FGE_x->size[25][0] = nPML_x_2;
	Allocate_FGE_x->size[25][1] = nPML_y_2; 
	Allocate_FGE_x->size[25][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Gx_18
	Allocate_FGE_x->size[26][0] = nPML_x_2;
	Allocate_FGE_x->size[26][1] = nPML_y_2;
	Allocate_FGE_x->size[26][2] = nz_MIN_nPML_z_1_MIN_nPML_z_2;
	//Fx_19
	Allocate_FGE_x->size[27][0] = nPML_x_1;
	Allocate_FGE_x->size[27][1] = nPML_y_1;
	Allocate_FGE_x->size[27][2] = nPML_z_2;
	//Gx_19
	Allocate_FGE_x->size[28][0] = nPML_x_1; 
	Allocate_FGE_x->size[28][1] = nPML_y_1; 
	Allocate_FGE_x->size[28][2] = nPML_z_2;
	//Fx_20
	Allocate_FGE_x->size[29][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2; 
	Allocate_FGE_x->size[29][1] = nPML_y_1;
	Allocate_FGE_x->size[29][2] = nPML_z_2;
	//Gx_20
	Allocate_FGE_x->size[30][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[30][1] = nPML_y_1; 
	Allocate_FGE_x->size[30][2] = nPML_z_2;
	//Fx_21
	Allocate_FGE_x->size[31][0] = nPML_x_2;
	Allocate_FGE_x->size[31][1] = nPML_y_1; 
	Allocate_FGE_x->size[31][2] = nPML_z_2;
	//Gx_21
	Allocate_FGE_x->size[32][0] = nPML_x_2; 
	Allocate_FGE_x->size[32][1] = nPML_y_1; 
	Allocate_FGE_x->size[32][2] = nPML_z_2;
	//Fx_22
	Allocate_FGE_x->size[33][0] = nPML_x_1; 
	Allocate_FGE_x->size[33][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2; 
	Allocate_FGE_x->size[33][2] = nPML_z_2;
	//Fx_23
	Allocate_FGE_x->size[34][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2; 
	Allocate_FGE_x->size[34][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2; 
	Allocate_FGE_x->size[34][2] = nPML_z_2;
	//Fx_24
	Allocate_FGE_x->size[35][0] = nPML_x_2;
	Allocate_FGE_x->size[35][1] = ny_MIN_nPML_y_1_MIN_nPML_y_2;
	Allocate_FGE_x->size[35][2] = nPML_z_2;
	//Fx_25
	Allocate_FGE_x->size[36][0] = nPML_x_1; 
	Allocate_FGE_x->size[36][1] = nPML_y_2; 
	Allocate_FGE_x->size[36][2] = nPML_z_2;
	//Gx_25
	Allocate_FGE_x->size[37][0] = nPML_x_1;
	Allocate_FGE_x->size[37][1] = nPML_y_2;
	Allocate_FGE_x->size[37][2] = nPML_z_2;
	//Fx_26
	Allocate_FGE_x->size[38][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2; 
	Allocate_FGE_x->size[38][1] = nPML_y_2; 
	Allocate_FGE_x->size[38][2] = nPML_z_2;
	//Gx_26
	Allocate_FGE_x->size[39][0] = nx_MIN_1_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_x->size[39][1] = nPML_y_2; 
	Allocate_FGE_x->size[39][2] = nPML_z_2;
	//Fx_27
	Allocate_FGE_x->size[40][0] = nPML_x_2; 
	Allocate_FGE_x->size[40][1] = nPML_y_2;
	Allocate_FGE_x->size[40][2] = nPML_z_2;
	//Gx_27
	Allocate_FGE_x->size[41][0] = nPML_x_2; 
	Allocate_FGE_x->size[41][1] = nPML_y_2;
	Allocate_FGE_x->size[41][2] = nPML_z_2;
	//Ex
	Allocate_FGE_x->size[42][0] = nx - 1; 
	Allocate_FGE_x->size[42][1] = ny;
	Allocate_FGE_x->size[42][2] = nz;
}

////////////////////////////////////////////////////////////////////
//Collect the information - Ey field
////////////////////////////////////////////////////////////////////
void Init_Alloc_FyGyEy(Allocate_E_comp *Allocate_FGE_y, double ***&Fy_1, double ***&Gy_1, 
				  double ***&Fy_2,  double ***&Gy_2,  double ***&Fy_3,  double ***&Gy_3, 
				  double ***&Fy_4,  double ***&Gy_4,  double ***&Fy_5,  double ***&Fy_6,  
				  double ***&Gy_6,  double ***&Fy_7,  double ***&Gy_7,  double ***&Fy_8, 
				  double ***&Gy_8,  double ***&Fy_9,  double ***&Gy_9,  double ***&Fy_10, 
				  double ***&Fy_11, double ***&Fy_12, double ***&Fy_13, double ***&Fy_15, 
				  double ***&Fy_16, double ***&Fy_17, double ***&Fy_18, double ***&Fy_19, 
				  double ***&Gy_19, double ***&Fy_20, double ***&Gy_20, double ***&Fy_21, 
				  double ***&Gy_21, double ***&Fy_22, double ***&Gy_22, double ***&Fy_23, 
				  double ***&Fy_24, double ***&Gy_24, double ***&Fy_25, double ***&Gy_25, 
				  double ***&Fy_26, double ***&Gy_26, double ***&Fy_27, double ***&Gy_27, 
				  double ***&Ey, long nPML_x_1, long nPML_y_1, long nPML_z_1, long nPML_x_2,
				  long nPML_y_2, long nPML_z_2, long nx, long ny, long nz)
{
	long nx_MIN_nPML_x_1_MIN_nPML_x_2 = nx - nPML_x_1 - nPML_x_2;
	long ny_MIN_1_MIN_nPML_y_1_MIN_nPML_y_2 = ny - 1 - nPML_y_1 - nPML_y_2;
	long nz_MIN_nPML_z_1_MIN_nPML_z_2 = nz - nPML_z_1 - nPML_z_2;

	Allocate_FGE_y->A[0]  = &Fy_1;  Allocate_FGE_y->A[1]  = &Gy_1;  Allocate_FGE_y->A[2]  = &Fy_2;  
	Allocate_FGE_y->A[3]  = &Gy_2;  Allocate_FGE_y->A[4]  = &Fy_3;  Allocate_FGE_y->A[5]  = &Gy_3;
	Allocate_FGE_y->A[6]  = &Fy_4;  Allocate_FGE_y->A[7]  = &Gy_4;  Allocate_FGE_y->A[8]  = &Fy_5;
	Allocate_FGE_y->A[9]  = &Fy_6;  Allocate_FGE_y->A[10] = &Gy_6;  Allocate_FGE_y->A[11] = &Fy_7;
	Allocate_FGE_y->A[12] = &Gy_7;  Allocate_FGE_y->A[13] = &Fy_8;  Allocate_FGE_y->A[14] = &Gy_8;
	Allocate_FGE_y->A[15] = &Fy_9;  Allocate_FGE_y->A[16] = &Gy_9;  Allocate_FGE_y->A[17] = &Fy_10;
	Allocate_FGE_y->A[18] = &Fy_11; Allocate_FGE_y->A[19] = &Fy_12; Allocate_FGE_y->A[20] = &Fy_13; 
	Allocate_FGE_y->A[21] = &Fy_15; Allocate_FGE_y->A[22] = &Fy_16; Allocate_FGE_y->A[23] = &Fy_17; 
	Allocate_FGE_y->A[24] = &Fy_18; Allocate_FGE_y->A[25] = &Fy_19; Allocate_FGE_y->A[26] = &Gy_19;
	Allocate_FGE_y->A[27] = &Fy_20; Allocate_FGE_y->A[28] = &Gy_20; Allocate_FGE_y->A[29] = &Fy_21;
	Allocate_FGE_y->A[30] = &Gy_21; Allocate_FGE_y->A[31] = &Fy_22; Allocate_FGE_y->A[32] = &Gy_22;
	Allocate_FGE_y->A[33] = &Fy_23; Allocate_FGE_y->A[34] = &Fy_24; Allocate_FGE_y->A[35] = &Gy_24;
	Allocate_FGE_y->A[36] = &Fy_25; Allocate_FGE_y->A[37] = &Gy_25; Allocate_FGE_y->A[38] = &Fy_26;
	Allocate_FGE_y->A[39] = &Gy_26; Allocate_FGE_y->A[40] = &Fy_27; Allocate_FGE_y->A[41] = &Gy_27;
	Allocate_FGE_y->A[42] = &Ey;

	Allocate_FGE_y->er = 0;

	strcpy(&Allocate_FGE_y->name[0][0],"Fy_1");   strcpy(&Allocate_FGE_y->name[1][0],"Gy_1"); 
	strcpy(&Allocate_FGE_y->name[2][0],"Fy_2");   strcpy(&Allocate_FGE_y->name[3][0],"Gy_2");  
	strcpy(&Allocate_FGE_y->name[4][0],"Fy_3");   strcpy(&Allocate_FGE_y->name[5][0],"Gy_3");
	strcpy(&Allocate_FGE_y->name[6][0],"Fy_4");   strcpy(&Allocate_FGE_y->name[7][0],"Gy_4");
	strcpy(&Allocate_FGE_y->name[8][0],"Fy_5");   strcpy(&Allocate_FGE_y->name[9][0],"Fy_6");  
	strcpy(&Allocate_FGE_y->name[10][0],"Gy_6");  strcpy(&Allocate_FGE_y->name[11][0],"Fy_7");
	strcpy(&Allocate_FGE_y->name[12][0],"Gy_7");  strcpy(&Allocate_FGE_y->name[13][0],"Fy_8");  
	strcpy(&Allocate_FGE_y->name[14][0],"Gy_8");  strcpy(&Allocate_FGE_y->name[15][0],"Fy_9"); 
	strcpy(&Allocate_FGE_y->name[16][0],"Gy_9");  strcpy(&Allocate_FGE_y->name[17][0],"Fy_10");
	strcpy(&Allocate_FGE_y->name[18][0],"Fy_11"); strcpy(&Allocate_FGE_y->name[19][0],"Fy_12");
	strcpy(&Allocate_FGE_y->name[20][0],"Fy_13"); strcpy(&Allocate_FGE_y->name[21][0],"Fy_15");
	strcpy(&Allocate_FGE_y->name[22][0],"Fy_16"); strcpy(&Allocate_FGE_y->name[23][0],"Fy_17"); 
	strcpy(&Allocate_FGE_y->name[24][0],"Fy_18"); strcpy(&Allocate_FGE_y->name[25][0],"Fy_19"); 
	strcpy(&Allocate_FGE_y->name[26][0],"Gy_19"); strcpy(&Allocate_FGE_y->name[27][0],"Fy_20");
	strcpy(&Allocate_FGE_y->name[28][0],"Gy_20"); strcpy(&Allocate_FGE_y->name[29][0],"Fy_21");
	strcpy(&Allocate_FGE_y->name[30][0],"Gy_21"); strcpy(&Allocate_FGE_y->name[31][0],"Fy_22");
	strcpy(&Allocate_FGE_y->name[32][0],"Gy_22"); strcpy(&Allocate_FGE_y->name[33][0],"Fy_23");
	strcpy(&Allocate_FGE_y->name[34][0],"Fy_24"); strcpy(&Allocate_FGE_y->name[35][0],"Gy_24");
	strcpy(&Allocate_FGE_y->name[36][0],"Fy_25"); strcpy(&Allocate_FGE_y->name[37][0],"Gy_25");
	strcpy(&Allocate_FGE_y->name[38][0],"Fy_26"); strcpy(&Allocate_FGE_y->name[39][0],"Gy_26");
	strcpy(&Allocate_FGE_y->name[40][0],"Fy_27"); strcpy(&Allocate_FGE_y->name[41][0],"Gy_27");
	strcpy(&Allocate_FGE_y->name[42][0],"Ey"); 
	
	//Fy_1
	Allocate_FGE_y->size[0][0] = nPML_x_1;
	Allocate_FGE_y->size[0][1] = nPML_y_1;
	Allocate_FGE_y->size[0][2] = nPML_z_1;
	//Gy_1
	Allocate_FGE_y->size[1][0] = nPML_x_1;
	Allocate_FGE_y->size[1][1] = nPML_y_1;
	Allocate_FGE_y->size[1][2] = nPML_z_1;
	//Fy_2
	Allocate_FGE_y->size[2][0] = nx_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_y->size[2][1] = nPML_y_1;
	Allocate_FGE_y->size[2][2] = nPML_z_1;
	//Gy_2
	Allocate_FGE_y->size[3][0] = nx_MIN_nPML_x_1_MIN_nPML_x_2;
	Allocate_FGE_y->size[3][1] = nPML_y_1;
	Allocate_FGE_y->size[3][2] = nPML_z_1;
	//Fy_3
	Allocate_FGE_y->size[4][0] = nPML_x_2;
	Allocate_FGE_y->size[4][1] = nPML_y_1;
	Allocate_FGE_y->size[4][2] = nPML_z_1;
	//Gy_3
	Allocate_FGE_y->size[5][0] = nPML_x_2;
	Allocate_FGE_y->size[5][1] = nPML_y_1;
	Allocate_FGE_y->size[5][2] = nPML_z_1;
	//Fy_4
	Allocate_FGE_y->size[6][0] = nPML_x_1;
	Allocate_FGE_y->size[6][1] = ny_MIN_1_MIN_nPML_y_1_MIN_nPML_y_2;
	Allocate_FGE_y->size[6][2] = nPML_z_1;
	//Gy_4
	Allocate_FGE_y->size[7][0] = nPML_x_1;