e_smagorinsky_filter.c File Reference

Subroutines for filtering variables needed for the Smagorinsky model. More...

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include "../../util/source/useful.h"
#include "../../util/source/compiler.h"
#include "../../util/source/logfile.h"
#include "../../gas_models/source/gas.h"
#include "../../flux_calc/source/flux_calc.h"
#include "../../mb_cns/source/cns_case_id.h"
#include "../../mb_cns/source/cns_work_vector.h"
#include "../../mb_cns/source/cns_debug_level.h"
#include "../../mb_cns/source/cns_cell.h"
#include "../../mb_cns/source/cns_echo.h"
#include "../../mb_cns/source/cns_bc_defs.h"
#include "../../mb_cns/source/mb_cfg.h"
#include "./e_block.h"
#include "../../mb_cns/source/interpolate.h"
#include "./e_stats.h"
#include "./e_smagorinsky_filter.h"

Defines
#define	COPY_TO_FILTER(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define	COPY_BACK_TO_BLOCK(WORK, QUANTITY, ILAST, JLAST, KLAST)
Functions
int	initialise_smagorinsky_filter_memory (int nnx, int nny, int nnz)
int	top_hat_filter (struct block_data_3D *A)
int	top_hat_filter_center (struct block_data_3D *A)
int	filter_x_top_hat (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int	filter_y_top_hat (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int	filter_z_top_hat (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int	filter_coefficients_ratio_2 (struct block_data_3D *A)
Variables
int	smagorinsky_filter_memory_alloc = 0
double ***	work_array
double ***	work_array_2

---

Detailed Description

Subroutines for filtering variables needed for the Smagorinsky model.

Author:

AD

Version:

September 2005

---

Define Documentation

#define COPY_BACK_TO_BLOCK (  WORK, QUANTITY, ILAST, JLAST, KLAST   )

Value: for ( i = 0; i <= ILAST; ++i ) { \ for ( j = 0; j <= JLAST; ++j ) { \ for ( k = 0; k <= KLAST; ++k) { \ A->ctr[i][j][k].QUANTITY = WORK[i][j][k] ; \ } \ } \ }

#define COPY_TO_FILTER (  WORK, QUANTITY, ILAST, JLAST, KLAST   )

Value: for ( i = 0; i <= ILAST; ++i ) { \ for ( j = 0; j <= JLAST; ++j ) { \ for ( k = 0; k <= KLAST; ++k) { \ WORK[i][j][k] = A->ctr[i][j][k].QUANTITY ; \ } \ } \ }

---