1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
|
/*============================================================================
* Define postprocessing output.
*============================================================================*/
/* VERS */
/*
This file is part of Code_Saturne, a general-purpose CFD tool.
Copyright (C) 1998-2019 EDF S.A.
This program is free software; you can redistribute it and/or modify it under
the terms of the GNU General Public License as published by the Free Software
Foundation; either version 2 of the License, or (at your option) any later
version.
This program is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
details.
You should have received a copy of the GNU General Public License along with
this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
/*----------------------------------------------------------------------------*/
#include "cs_defs.h"
/*----------------------------------------------------------------------------
* Standard C library headers
*----------------------------------------------------------------------------*/
#include "stdlib.h"
#include "string.h"
/*----------------------------------------------------------------------------
* Local headers
*----------------------------------------------------------------------------*/
#include "bft_mem.h"
#include "bft_error.h"
#include "cs_base.h"
#include "cs_field.h"
#include "cs_geom.h"
#include "cs_interpolate.h"
#include "cs_mesh.h"
#include "cs_selector.h"
#include "cs_parall.h"
#include "cs_post.h"
#include "cs_post_util.h"
#include "cs_probe.h"
#include "cs_time_plot.h"
#include "cs_field_pointer.h"
#include "cs_notebook.h"
#include "cs_parameters.h"
#include "cs_physical_constants.h"
#include "cs_turbulence_model.h"
/*----------------------------------------------------------------------------
* Header for the current file
*----------------------------------------------------------------------------*/
#include "cs_prototypes.h"
/*----------------------------------------------------------------------------*/
BEGIN_C_DECLS
/*============================================================================
* Local (user defined) function definitions
*============================================================================*/
/*----------------------------------------------------------------------------
* Example function for selection of cells with scalar field values above
* a certain threshold.
*
* In this case, the selection is base on the value of the temperaturefield
* being above 21 degrees.
*
* parameters:
* input <-> pointer to input (unused here)
* n_cells --> number of selected cells
* cell_ids --> array of selected cell ids (0 to n-1 numbering)
*----------------------------------------------------------------------------*/
static void
_t_lt_21_select(void *input,
cs_lnum_t *n_cells,
cs_lnum_t **cell_ids)
{
cs_lnum_t i;
cs_lnum_t _n_cells = 0;
cs_lnum_t *_cell_ids = NULL;
const cs_mesh_t *m = cs_glob_mesh;
cs_field_t *f = cs_field_by_name("temperature"); /* Get access to field */
if (f == NULL)
bft_error(__FILE__, __LINE__, 0,
"No field with name \"temperature\" defined");
/* Before time loop, field is defined, but has no values yet,
so ignore that case (postprocessing mesh will be initially empty) */
if (f->val != NULL) {
BFT_MALLOC(_cell_ids, m->n_cells, cs_lnum_t); /* Allocate selection list */
for (i = 0; i < m->n_cells; i++) {
if (f->val[i] < 21) {
_cell_ids[_n_cells] = i;
_n_cells += 1;
}
}
BFT_REALLOC(_cell_ids, _n_cells, cs_lnum_t); /* Adjust size (good practice,
but not required) */
}
/* Set return values */
*n_cells = _n_cells;
*cell_ids = _cell_ids;
}
/*============================================================================
* User function definitions
*============================================================================*/
/*----------------------------------------------------------------------------*/
/*!
* \brief Define post-processing writers.
*
* The default output format and frequency may be configured, and additional
* post-processing writers allowing outputs in different formats or with
* different format options and output frequency than the main writer may
* be defined.
*/
/*----------------------------------------------------------------------------*/
void
cs_user_postprocess_writers(void)
{
/* Define additional writers */
cs_post_define_writer(1, /* writer_id */
"user", /* writer name */
"postprocessing", /* directory name */
"EnSight Gold", /* format_name */
"", /* format_options */
FVM_WRITER_TRANSIENT_CONNECT, /* time dependency */
false, /* output_at_start */
true, /* output_at_end */
5, /* frequency_n */
-1); /* frequency_t */
}
/*----------------------------------------------------------------------------*/
/*!
* \brief Define post-processing meshes.
*
* The main post-processing meshes may be configured, and additional
* post-processing meshes may be defined as a subset of the main mesh's
* cells or faces (both interior and boundary).
*/
/*----------------------------------------------------------------------------*/
void
cs_user_postprocess_meshes(void)
{
/* Build a (time varying) volume mesh containing cells
with values of field named "temperature" > < 21 */
const int n_writers = 1;
const int writer_ids[] = {1}; /* Associate to writer 1 */
/* Define postprocessing mesh */
cs_post_define_volume_mesh_by_func(1, /* mesh id */
"T_lt_21",
_t_lt_21_select,
NULL, /* _t_lt_21_select_input */
true, /* time varying */
false, /* add_groups */
false, /* auto_variables */
n_writers,
writer_ids);
}
/*----------------------------------------------------------------------------*/
/*!
* \brief User function for output of values on a post-processing mesh.
*
* \param[in] mesh_name name of the output mesh for the current call
* \param[in] mesh_id id of the output mesh for the current call
* \param[in] cat_id category id of the output mesh for the
* current call
* \param[in] probes pointer to associated probe set structure if
* the mesh is a probe set, NULL otherwise
* \param[in] n_cells local number of cells of post_mesh
* \param[in] n_i_faces local number of interior faces of post_mesh
* \param[in] n_b_faces local number of boundary faces of post_mesh
* \param[in] n_vertices local number of vertices faces of post_mesh
* \param[in] cell_list list of cells (0 to n-1) of post-processing
* mesh
* \param[in] i_face_list list of interior faces (0 to n-1) of
* post-processing mesh
* \param[in] b_face_list list of boundary faces (0 to n-1) of
* post-processing mesh
* \param[in] vertex_list list of vertices (0 to n-1) of
* post-processing mesh
* \param[in] ts time step status structure, or NULL
*/
/*----------------------------------------------------------------------------*/
void
cs_user_postprocess_values(const char *mesh_name,
int mesh_id,
int cat_id,
cs_probe_set_t *probes,
cs_lnum_t n_cells,
cs_lnum_t n_i_faces,
cs_lnum_t n_b_faces,
cs_lnum_t n_vertices,
const cs_lnum_t cell_list[],
const cs_lnum_t i_face_list[],
const cs_lnum_t b_face_list[],
const cs_lnum_t vertex_list[],
const cs_time_step_t *ts)
{
/* Output of the temperature on output mesh 1 */
if (mesh_id == 1) {
cs_field_t *f = cs_field_by_name("temperature"); /* Get access to field */
cs_post_write_var(mesh_id,
CS_POST_WRITER_ALL_ASSOCIATED, /* writer id filter */
cs_field_get_label(f), /* var_name */
1, /* var_dim */
true, /* interlace, */
true, /* use_parent */
CS_POST_TYPE_cs_real_t, /* var_type */
f->val, /* cel_vals */
NULL, /* i_face_vals */
NULL, /* b_face_vals */
ts);
}
}
/*----------------------------------------------------------------------------*/
END_C_DECLS
|
