 
 	return integral;
 	*/
-	return (def.P_atm / def.ro0_w) + (T - T_0) * (c0_w - (T - T_0) * (C_w / 2 + C_w2 * (T - T_0) / 3));
+	return (gpu_def->P_atm / gpu_def->ro0_w) + (T - T_0) * (c0_w - (T - T_0) * (C_w / 2 + C_w2 * (T - T_0) / 3));
 }
 
 __device__ double assign_H_n (double T)
 {
-	return (def.P_atm / def.ro0_n) + (T - T_0) * (c0_n + C_n * (T - T_0) / 2);
+	return (gpu_def->P_atm / gpu_def->ro0_n) + (T - T_0) * (c0_n + C_n * (T - T_0) / 2);
 }
 
 __device__ double assign_H_g (double T)
 {
-	return (def.P_atm / def.ro0_g) + (T - T_0) * (c0_g + C_g * (T - T_0) / 2);
+	return (gpu_def->P_atm / gpu_def->ro0_g) + (T - T_0) * (c0_g + C_g * (T - T_0) / 2);
 }
 
 __device__ double assign_H_r (double T)
 {
-	return (def.P_atm / ro_r) + (T - T_0) * (c0_r + C_r * (T - T_0) / 2);
+	return (gpu_def->P_atm / ro_r) + (T - T_0) * (c0_r + C_r * (T - T_0) / 2);
 }
 
 __global__ void assign_H (ptr_Arrays DevArraysPtr, int local)

	case 'w':
		if (var == 'P')
		{
			result_d_ro = def.ro0_w * (def.beta_w);
		} 
		else if (var == 'T')
		{
			result_d_ro = (-1) * def.ro0_w * alfa_w;
		}
		else 
		{

__device__ double assign_E_flow (ptr_Arrays DevArraysPtr, i, j, k)
{
	
	if ((((i != 0) && (i != (gpu_def->locNx) - 1)) || ((gpu_def->locNx) < 2)) && (j != 0) && (j != (gpu_def->locNy) - 1) && (((k != 0) && (k != (gpu_def->locNz) - 1)) || ((gpu_def->locNz) < 2)))
	{
		double E_flow = 0;
		int local=i + j * (gpu_def->locNx) + k * (gpu_def->locNx) * (gpu_def->locNy);

		if ((gpu_def->locNx) > 2)
		{
			E_flow += (DevArraysPtr.ro_w[local + 1] * DevArraysPtr.H_w[local + 1] * DevArraysPtr.ux_w[local + 1]
				- DevArraysPtr.ro_w[local - 1] * DevArraysPtr.H_w[local - 1] * DevArraysPtr.ux_w[local - 1]
				+ DevArraysPtr.ro_n[local + 1] * DevArraysPtr.H_n[local + 1] * DevArraysPtr.ux_n[local + 1]
				- DevArraysPtr.ro_n[local - 1] * DevArraysPtr.H_n[local - 1] * DevArraysPtr.ux_n[local - 1]
				+ DevArraysPtr.ro_g[local + 1] * DevArraysPtr.H_g[local + 1] * DevArraysPtr.ux_g[local + 1]
				- DevArraysPtr.ro_g[local - 1] * DevArraysPtr.H_g[local - 1] * DevArraysPtr.ux_g[local - 1]
				) / (2. * (gpu_def->hx));
		}
		if ((gpu_def->locNy) > 2)
		{
			E_flow += (DevArraysPtr.ro_w[local + gpu_def->locNx] * DevArraysPtr.H_w[local + gpu_def->locNx] * DevArraysPtr.uy_w[local + gpu_def->locNx]
			- DevArraysPtr.ro_w[local - gpu_def->locNx] * DevArraysPtr.H_w[local - gpu_def->locNx] * DevArraysPtr.uy_w[local - gpu_def->locNx]
			+ DevArraysPtr.ro_n[local + gpu_def->locNx] * DevArraysPtr.H_n[local + gpu_def->locNx] * DevArraysPtr.uy_n[local + gpu_def->locNx]
			- DevArraysPtr.ro_n[local - gpu_def->locNx] * DevArraysPtr.H_n[local - gpu_def->locNx] * DevArraysPtr.uy_n[local - gpu_def->locNx]
			+ DevArraysPtr.ro_g[local + gpu_def->locNx] * DevArraysPtr.H_g[local + gpu_def->locNx] * DevArraysPtr.uy_g[local + gpu_def->locNx]
			- DevArraysPtr.ro_g[local - gpu_def->locNx] * DevArraysPtr.H_g[local - gpu_def->locNx] * DevArraysPtr.uy_g[local - gpu_def->locNx]
			)/ (2. * (gpu_def->hy));	
		}

		if ((gpu_def->locNz) > 2)
		{
			E_flow += (DevArraysPtr.ro_w[local + (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.H_w[local + (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.uy_w[local + (gpu_def->locNx) * (gpu_def->locNy)]
			- DevArraysPtr.ro_w[local - (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.H_w[local - (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.uy_w[local - (gpu_def->locNx) * (gpu_def->locNy)]
			+ DevArraysPtr.ro_n[local + (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.H_n[local + (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.uy_n[local + (gpu_def->locNx) * (gpu_def->locNy)]
			- DevArraysPtr.ro_n[local - (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.H_n[local - (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.uy_n[local - (gpu_def->locNx) * (gpu_def->locNy)]
			+ DevArraysPtr.ro_g[local + (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.H_g[local + (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.uy_g[local + (gpu_def->locNx) * (gpu_def->locNy)]
			- DevArraysPtr.ro_g[local - (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.H_g[local - (gpu_def->locNx) * (gpu_def->locNy)] * DevArraysPtr.uy_g[local - (gpu_def->locNx) * (gpu_def->locNy)]
			)/ (2. * (gpu_def->hz));	
		}

		device_test_u(E_flow, __FILE__, __LINE__);
		return E_flow;
	}
	else
