~reducedmodelling/fluidity/ReducedModel

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
*realformat "%16.6f"
*intformat "%7i"
*Set Cond 01 *elems *canrepeat
*Add Cond 02 *elems *canrepeat
*Add Cond 03 *elems *canrepeat
*Add Cond 04 *elems *canrepeat
*Add Cond 05 *elems *canrepeat
*Add Cond 06 *elems *canrepeat
*Add Cond 07 *elems *canrepeat
*Add Cond 08 *elems *canrepeat
*Add Cond 09 *elems *canrepeat
*Add Cond 10 *elems *canrepeat
*Add Cond 11 *elems *canrepeat
*Add Cond 12 *elems *canrepeat
*Add Cond 13 *elems *canrepeat
*Add Cond 14 *elems *canrepeat
*Add Cond 15 *elems *canrepeat
*Add Cond 16 *elems *canrepeat
*Add Cond 17 *elems *canrepeat
*Add Cond 18 *elems *canrepeat
*Add Cond 19 *elems *canrepeat
*Add Cond 20 *elems *canrepeat
*set var nsurfaces(int)=CondNumEntities(int) 
$MeshFormat
2 0 8
$EndMeshFormat
$Nodes
*format "%i"
*npoin
*loop nodes
*format "%i%f %f %f"
*NodesNum *NodesCoord
*end nodes
$EndNodes
$Elements
*format "%i"
*operation(nsurfaces+nelem)
*loop elems *OnlyinCond
*format "%i %i %i %i %i %i"
*LoopVar 2 2 *Cond(1,int) *Cond(1,int) *GlobalNodes
*end elems
*Set Cond region_ids *elems 
*loop elems
*operation(nsurfaces+loopvar) *ElemsType 2 1 0 *ElemsConec
*end elems
$EndElements