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
|
program write_cgns_1
implicit none
! author: Diane Poirier (diane@icemcfd.com)
! last revised on March 8 2000
! This example test the complete SIDS for multi-block data.
! It creates a dummy mesh composed of 2 structured blocks in 3D.
#ifdef WINNT
include 'cgnswin_f.h'
#endif
include 'cgnslib_f.h'
integer Ndim
parameter (Ndim = 3)
integer index_dim, cell_dim, phys_dim
integer base_no, zone_no, coord_no, sol_no, discr_no, conn_no
integer hole_no, boco_no, field_no, dset_no
integer num, size(Ndim*3), npnts, NormalIndex(Ndim)
integer cg, ier, zone, coord, i, j, k, n, pos, sol, field
integer pnts(Ndim,120), donor_pnts(Ndim,120)
integer transform(Ndim)
real*4 data(120), normals(360)
double precision Dxyz(120), values(120)
character*32 zonename, solname, fieldname
character*32 coordname(Ndim)
character*32 donorname
coordname(1) = 'CoordinateX'
coordname(2) = 'CoordinateY'
coordname(3) = 'CoordinateZ'
! *** initialize
ier = 0
index_dim=Ndim
cell_dim=Ndim
phys_dim=Ndim
! *** open CGNS file for writing
call cg_open_f('cgtest.cgns', CG_MODE_WRITE, cg, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** base
call cg_base_write_f(cg, 'Basename', cell_dim, phys_dim,
& base_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** zone
do zone=1, 2
write(zonename,'(a5,i1)') 'zone#',zone
num = 1
do i=1,index_dim ! zone#1: 3*4*5, zone#2: 4*5*6
size(i) = i+zone+1 ! nr of nodes in i,j,k
size(i+Ndim) = size(i)-1 ! nr of elements in i,j,k
size(i+2*Ndim) = 0 ! nr of bnd nodes if ordered
num = num * size(i) ! nr of nodes
enddo
!234567890!234567890!234567890!234567890!234567890!234567890!23456789012
call cg_zone_write_f(cg, base_no, zonename, size,
& Structured, zone_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** coordinate
do coord=1, phys_dim
do k=1, size(3)
do j=1, size(2)
do i=1, size(1)
pos = i + (j-1)*size(1) + (k-1)*size(1)*size(2)
! * make up some dummy coordinates just for the test:
if (coord.eq.1) Dxyz(pos) = i
if (coord.eq.2) Dxyz(pos) = j
if (coord.eq.3) Dxyz(pos) = k
enddo
enddo
enddo
call cg_coord_write_f(cg, base_no, zone_no, RealDouble,
& coordname(coord), Dxyz, coord_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
enddo
! *** solution
do sol=1, 2
write(solname,'(a5,i1,a5,i1)') 'Zone#',zone,' sol#',sol
call cg_sol_write_f(cg, base_no, zone_no, solname,
& Vertex, sol_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** solution field
do field=1, 2
! make up some dummy solution values
do i=1, num
values(i) = i*field*sol
enddo
write(fieldname,'(a6,i1)') 'Field#',field
call cg_field_write_f(cg, base_no, zone_no, sol_no,
& RealDouble, fieldname, values, field_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
enddo ! field loop
enddo ! solution loop
! *** discrete data
call cg_discrete_write_f(cg, base_no, zone_no, 'discrete#1',
& discr_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** discrete data arrays, defined on vertices:
call cg_goto_f(cg, base_no, ier, 'Zone_t', zone,
& 'DiscreteData_t', discr_no, 'end')
if (ier .eq. ERROR) call cg_error_exit_f
do 123 k=1, size(3)
do 123 j=1, size(2)
do 123 i=1, size(1)
pos = i + (j-1)*size(1) + (k-1)*size(1)*size(2)
data(pos) = pos ! * make up some dummy data
123 continue
call cg_array_write_f('arrayname', RealSingle, index_dim,
& size, data, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** discrete data arrays attribute: GOTO DataArray node
call cg_goto_f(cg, base_no, ier, 'Zone_t', zone,
& 'DiscreteData_t', discr_no, 'DataArray_t', 1, 'end')
if (ier .eq. ERROR) call cg_error_exit_f
call cg_units_write_f(Kilogram, Meter, Second, Kelvin,
& Radian, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** overset holes
! create dummy data
do i=1,3
! Define 2 separate PointRange, for 2 patches in the hole
pnts(i,1)=1
pnts(i,2)=size(i)
! second PointRange of hole
pnts(i,3)=2
pnts(i,4)=size(i)
enddo
! Hole defined with 2 point set type PointRange, so 4 points:
call cg_hole_write_f(cg, base_no, zone_no, 'hole#1', Vertex,
& PointRange, 2, 4, pnts, hole_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** general connectivity
do 100 n=1, 5
do 100 i=1,3
pnts(i,n)=i ! * dummy data
donor_pnts(i,n)=i*2
100 continue
! create a point matching connectivity
call cg_conn_write_f(cg, base_no, zone_no, 'Connect#1',
& Vertex, Abutting1to1, PointList, 5, pnts, 'zone#2',
& Structured, PointListDonor, Integer, 5, donor_pnts,
& conn_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** connectivity 1to1
! generate data
do i=1,3
!**make up some dummy data:
pnts(i,1)=1
pnts(i,2)=size(i)
donor_pnts(i,1)=1
donor_pnts(i,2)=size(i)
transform(i)=i*(-1)
enddo
if (zone .eq. 1) then
donorname='zone#2'
else if (zone .eq. 2) then
donorname='zone#1'
endif
call cg_1to1_write_f(cg, base_no, zone_no, '1to1_#1',
& donorname, pnts, donor_pnts, transform, conn_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** ZoneGridConnectivity attributes: GOTO ZoneGridConnectivity_t node
call cg_goto_f(cg, base_no, ier, 'Zone_t', zone,
& 'ZoneGridConnectivity_t', 1, 'end')
if (ier .eq. ERROR) call cg_error_exit_f
! *** ZoneGridConnectivity attributes: Descriptor_t
!234567890!234567890!234567890!234567890!234567890!234567890!23456789012
call cg_descriptor_write_f('DescriptorName',
& 'Zone Connectivity', ier)
! *** bocos
call cg_boco_write_f(cg, base_no, zone_no, 'boco#1',
& BCInflow, PointRange, 2, pnts, boco_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! *** boco normal
npnts = 1
do i=1,Ndim
NormalIndex(i)=0
! compute nr of points on bc patch:
npnts = npnts * (pnts(i,2)-pnts(i,1)+1)
enddo
NormalIndex(1)=1
do i=1,phys_dim*npnts
normals(i)=i
enddo
call cg_boco_normal_write_f(cg, base_no, zone_no, boco_no,
& NormalIndex, 1, RealSingle, normals, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! ** boundary condition attributes: GOTO BC_t node
call cg_goto_f(cg, base_no, ier, 'Zone_t', zone, 'ZoneBC_t',
& 1, 'BC_t', boco_no, 'end')
if (ier .eq. ERROR) call cg_error_exit_f
! ** boundary condition attributes: GridLocation_t
call cg_gridlocation_write_f(Vertex, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! ** boundary condition dataset
call cg_dataset_write_f(cg, base_no, zone,
& boco_no, 'DataSetName', BCInflow, dset_no, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! ** boundary condition data:
call cg_bcdata_write_f(cg, base_no, zone,
& boco_no, dset_no, Neumann, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! ** boundary condition data arrays: GOTO BCData_t node
call cg_goto_f(cg, base_no, ier, 'Zone_t', zone_no,
& 'ZoneBC_t', 1, 'BC_t', boco_no, 'BCDataSet_t', dset_no,
& 'BCData_t', Neumann, 'end')
if (ier .eq. ERROR) call cg_error_exit_f
do i=1, npnts
data(i) = i
enddo
call cg_array_write_f('dataset_arrayname', RealSingle,
& 1, npnts, data, ier)
if (ier .eq. ERROR) call cg_error_exit_f
! ** boundary condition data attributes:
call cg_dataclass_write_f(NormalizedByDimensional, ier)
if (ier .eq. ERROR) call cg_error_exit_f
enddo ! zone loop
! *** close CGNS file
call cg_close_f(cg, ier)
if (ier .eq. ERROR) call cg_error_exit_f
end
|
