1
Calling fluids from mainfl
2
Checking momentum discretisation options
3
Finished checking momentum discretisation options
4
Options sanity check successful
6
Using the new makeeelist
8
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
9
Halo nodes are consistent with trailing receive ordering
10
Creating universal numbering for a trailing receives halo
11
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
12
Halo nodes are consistent with trailing receive ordering
13
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
14
Halo nodes are consistent with trailing receive ordering
15
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
16
Halo nodes are consistent with trailing receive ordering
17
Checking nodes in halo CoordinateMeshLevel2Halo for consistency with trailing receive ordering
18
Halo nodes are consistent with trailing receive ordering
19
Creating universal numbering for a trailing receives halo
20
Checking nodes in halo CoordinateMeshLevel2Halo for consistency with trailing receive ordering
21
Halo nodes are consistent with trailing receive ordering
22
Checking nodes in halo CoordinateMeshLevel2Halo for consistency with trailing receive ordering
23
Halo nodes are consistent with trailing receive ordering
24
Checking nodes in halo CoordinateMeshLevel2Halo for consistency with trailing receive ordering
25
Halo nodes are consistent with trailing receive ordering
26
In derive_maximal_element_halo
27
Checking nodes in halo CoordinateMeshLevel2Halo for consistency with trailing receive ordering
28
Halo nodes are consistent with trailing receive ordering
29
Checking nodes in halo CoordinateMeshMaximalElementHalo for consistency with trailing receive ordering
30
Halo nodes are consistent with trailing receive ordering
31
Creating universal numbering for a trailing receives halo
32
Checking nodes in halo CoordinateMeshMaximalElementHalo for consistency with trailing receive ordering
33
Halo nodes are consistent with trailing receive ordering
34
Checking nodes in halo CoordinateMeshMaximalElementHalo for consistency with trailing receive ordering
35
Halo nodes are consistent with trailing receive ordering
36
Checking nodes in halo CoordinateMeshMaximalElementHalo for consistency with trailing receive ordering
37
Halo nodes are consistent with trailing receive ordering
38
Exiting derive_maximal_element_halo
39
Exiting read_halos_mesh
40
Creating universal numbering for a trailing receives halo
41
Checking nodes in halo CoordinateMeshMaximalElementHaloDG for consistency with trailing receive ordering
42
Halo nodes are consistent with trailing receive ordering
43
Checking nodes in halo CoordinateMeshMaximalElementHaloDG for consistency with trailing receive ordering
44
Halo nodes are consistent with trailing receive ordering
45
Checking nodes in halo CoordinateMeshMaximalElementHaloDG for consistency with trailing receive ordering
46
Halo nodes are consistent with trailing receive ordering
47
Creating universal numbering for a trailing receives halo
48
Checking nodes in halo CoordinateMeshMaximalElementHaloDG for consistency with trailing receive ordering
49
Halo nodes are consistent with trailing receive ordering
50
Checking nodes in halo CoordinateMeshMaximalElementHaloDG for consistency with trailing receive ordering
51
Halo nodes are consistent with trailing receive ordering
52
Checking nodes in halo CoordinateMeshMaximalElementHaloDG for consistency with trailing receive ordering
53
Halo nodes are consistent with trailing receive ordering
54
Using the new makeeelist
56
In derive_maximal_surface_element_halo
57
Checking nodes in halo CoordinateMeshMaximalElementHalo for consistency with trailing receive ordering
58
Halo nodes are consistent with trailing receive ordering
59
Exiting derive_maximal_surface_element_halo
60
Exiting merge_surface_ids
61
In allocate_and_insert_scalar_field, field is: Density
62
Is field prognostic? F
63
Is field prescribed? F
65
Is field diagnostic? T
66
In allocate_and_insert_scalar_field, field is: Temperature
67
Is field prognostic? T
68
Is field prescribed? F
70
Is field diagnostic? F
71
In allocate_and_insert_scalar_field, field is: Salinity
72
Is field prognostic? T
73
Is field prescribed? F
75
Is field diagnostic? F
76
In allocate_and_insert_scalar_field, field is: PerturbationDensity
77
Is field prognostic? F
78
Is field prescribed? F
80
Is field diagnostic? T
81
In allocate_and_insert_scalar_field, field is: CFLNumber
82
Is field prognostic? F
83
Is field prescribed? F
85
Is field diagnostic? T
86
In allocate_and_insert_scalar_field, field is: DistanceToTop
87
Is field prognostic? F
88
Is field prescribed? T
90
Is field diagnostic? F
91
In allocate_and_insert_scalar_field, field is: DistanceToBottom
92
Is field prognostic? F
93
Is field prescribed? T
95
Is field diagnostic? F
96
In allocate_and_insert_vector_field, field is: Velocity
97
Is field prognostic? T
98
Is field prescribed? F
100
Is field diagnostic? F
101
In allocate_and_insert_vector_field, field is: VelocityAbsorption
102
Is field prognostic? F
103
Is field prescribed? T
105
Is field diagnostic? F
106
In allocate_and_insert_tensor_field, field is: Viscosity
107
Is field prescribed? F
108
Is field diagnostic? T
110
In allocate_and_insert_vector_field, field is: GravityDirection
111
Is field prognostic? F
112
Is field prescribed? T
114
Is field diagnostic? F
115
In allocate_and_insert_scalar_field, field is: GLSTurbulentKineticEnergy
116
Is field prognostic? T
117
Is field prescribed? F
119
Is field diagnostic? F
120
In allocate_and_insert_scalar_field, field is: GLSTurbulentKineticEnergySource
121
Is field prognostic? F
122
Is field prescribed? F
124
Is field diagnostic? T
125
In allocate_and_insert_scalar_field, field is: GLSTurbulentKineticEnergyAbsorption
126
Is field prognostic? F
127
Is field prescribed? F
129
Is field diagnostic? T
130
In allocate_and_insert_tensor_field, field is: GLSTurbulentKineticEnergyDiffusivity
131
Is field prescribed? F
132
Is field diagnostic? T
134
In allocate_and_insert_scalar_field, field is: GLSGenericSecondQuantity
135
Is field prognostic? T
136
Is field prescribed? F
138
Is field diagnostic? F
139
In allocate_and_insert_scalar_field, field is: GLSGenericSecondQuantitySource
140
Is field prognostic? F
141
Is field prescribed? F
143
Is field diagnostic? T
144
In allocate_and_insert_scalar_field, field is: GLSGenericSecondQuantityAbsorption
145
Is field prognostic? F
146
Is field prescribed? F
148
Is field diagnostic? T
149
In allocate_and_insert_tensor_field, field is: GLSGenericSecondQuantityDiffusivity
150
Is field prescribed? F
151
Is field diagnostic? T
153
In allocate_and_insert_scalar_field, field is: GLSLengthScale
154
Is field prognostic? F
155
Is field prescribed? F
157
Is field diagnostic? T
158
In allocate_and_insert_scalar_field, field is: GLSBuoyancyFrequency
159
Is field prognostic? F
160
Is field prescribed? F
162
Is field diagnostic? T
163
In allocate_and_insert_scalar_field, field is: GLSVelocityShear
164
Is field prognostic? F
165
Is field prescribed? F
167
Is field diagnostic? T
168
In allocate_and_insert_scalar_field, field is: GLSShearProduction
169
Is field prognostic? F
170
Is field prescribed? F
172
Is field diagnostic? T
173
In allocate_and_insert_scalar_field, field is: GLSBuoyancyProduction
174
Is field prognostic? F
175
Is field prescribed? F
177
Is field diagnostic? T
178
In allocate_and_insert_scalar_field, field is: GLSDissipationEpsilon
179
Is field prognostic? F
180
Is field prescribed? F
182
Is field diagnostic? T
183
In allocate_and_insert_scalar_field, field is: GLSStabilityFunctionSM
184
Is field prognostic? F
185
Is field prescribed? F
187
Is field diagnostic? T
188
In allocate_and_insert_scalar_field, field is: GLSStabilityFunctionSH
189
Is field prognostic? F
190
Is field prescribed? F
192
Is field diagnostic? T
193
In allocate_and_insert_scalar_field, field is: GLSSource1
194
Is field prognostic? F
195
Is field prescribed? F
197
Is field diagnostic? T
198
In allocate_and_insert_scalar_field, field is: GLSSource2
199
Is field prognostic? F
200
Is field prescribed? F
202
Is field diagnostic? T
203
In allocate_and_insert_scalar_field, field is: GLSAbsorption1
204
Is field prognostic? F
205
Is field prescribed? F
207
Is field diagnostic? T
208
In allocate_and_insert_scalar_field, field is: GLSAbsorption2
209
Is field prognostic? F
210
Is field prescribed? F
212
Is field diagnostic? T
213
In allocate_and_insert_scalar_field, field is: GLSWallFunction
214
Is field prognostic? F
215
Is field prescribed? F
217
Is field diagnostic? T
218
In allocate_and_insert_scalar_field, field is: GLSVerticalViscosity
219
Is field prognostic? F
220
Is field prescribed? F
222
Is field diagnostic? T
223
In allocate_and_insert_scalar_field, field is: GLSVerticalDiffusivity
224
Is field prognostic? F
225
Is field prescribed? F
227
Is field diagnostic? T
228
In allocate_and_insert_tensor_field, field is: GLSBackgroundViscosity
229
Is field prescribed? T
230
Is field diagnostic? F
232
In allocate_and_insert_tensor_field, field is: GLSBackgroundDiffusivity
233
Is field prescribed? T
234
Is field diagnostic? F
236
In allocate_and_insert_tensor_field, field is: GLSEddyViscosityKM
237
Is field prescribed? F
238
Is field diagnostic? T
240
In allocate_and_insert_tensor_field, field is: GLSEddyDiffusivityKH
241
Is field prescribed? F
242
Is field diagnostic? T
244
In initialise_prognostic_fields
245
In set_prescribed_field_values
246
In populate_boundary_conditions
247
Initialising GLS stable boundaries
248
In set_boundary_conditions
249
In set_dirichlet_consistent
250
Exiting populate_state
251
In allocate_and_insert_auxilliary_fields
252
In enforce_discrete_properties
253
Exiting enforce_discrete_properties
254
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
255
% Some quantities associated with the initial set-up of this problem. %
256
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
258
The time step (DT) is: 15.000000000000000
259
The end time (LTIME) is set to: 45.000000000000000
260
This corresponds to this many time steps in simulation: 3.0000000000000000
262
The maximum volume of element in the mesh is: 15000.000000000187
263
The minimum volume of element in the mesh is: 14999.999999999809
264
The ratio of max to min volumes is: 1.0000000000000253
265
The following material phases and their contents are in existence
274
+Density (Density) on VelocityMesh
275
+Temperature (Temperature) on VelocityMesh
276
+Salinity (Salinity) on VelocityMesh
277
+PerturbationDensity (PerturbationDensity) on VelocityMesh
278
+CFLNumber (CFLNumber) on VelocityMesh
279
+DistanceToTop (DistanceToTop) on VelocityMesh
280
+DistanceToBottom (DistanceToBottom) on VelocityMesh
281
+GLSTurbulentKineticEnergy (GLSTurbulentKineticEnergy) on VelocityMesh
282
+GLSTurbulentKineticEnergySource (GLSTurbulentKineticEnergySource) on VelocityMesh
283
+GLSTurbulentKineticEnergyAbsorption (GLSTurbulentKineticEnergyAbsorption) on VelocityMesh
284
+GLSGenericSecondQuantity (GLSGenericSecondQuantity) on VelocityMesh
285
+GLSGenericSecondQuantitySource (GLSGenericSecondQuantitySource) on VelocityMesh
286
+GLSGenericSecondQuantityAbsorption (GLSGenericSecondQuantityAbsorption) on VelocityMesh
287
+GLSLengthScale (GLSLengthScale) on VelocityMesh
288
+GLSBuoyancyFrequency (GLSBuoyancyFrequency) on VelocityMesh
289
+GLSVelocityShear (GLSVelocityShear) on VelocityMesh
290
+GLSShearProduction (GLSShearProduction) on VelocityMesh
291
+GLSBuoyancyProduction (GLSBuoyancyProduction) on VelocityMesh
292
+GLSDissipationEpsilon (GLSDissipationEpsilon) on VelocityMesh
293
+GLSStabilityFunctionSM (GLSStabilityFunctionSM) on VelocityMesh
294
+GLSStabilityFunctionSH (GLSStabilityFunctionSH) on VelocityMesh
295
+GLSSource1 (GLSSource1) on VelocityMesh
296
+GLSSource2 (GLSSource2) on VelocityMesh
297
+GLSAbsorption1 (GLSAbsorption1) on VelocityMesh
298
+GLSAbsorption2 (GLSAbsorption2) on VelocityMesh
299
+GLSWallFunction (GLSWallFunction) on VelocityMesh
300
+GLSVerticalViscosity (GLSVerticalViscosity) on VelocityMesh
301
+GLSVerticalDiffusivity (GLSVerticalDiffusivity) on VelocityMesh
302
+OldDensity (OldDensity) on VelocityMesh
303
+IteratedDensity (IteratedDensity) on VelocityMesh
304
+OldTemperature (OldTemperature) on VelocityMesh
305
+IteratedTemperature (IteratedTemperature) on VelocityMesh
306
+OldSalinity (OldSalinity) on VelocityMesh
307
+IteratedSalinity (IteratedSalinity) on VelocityMesh
308
+OldPerturbationDensity (OldPerturbationDensity) on VelocityMesh
309
+IteratedPerturbationDensity (IteratedPerturbationDensity) on VelocityMesh
310
+OldCFLNumber (OldCFLNumber) on VelocityMesh
311
+IteratedCFLNumber (IteratedCFLNumber) on VelocityMesh
312
+OldDistanceToTop (OldDistanceToTop) on VelocityMesh
313
+IteratedDistanceToTop (IteratedDistanceToTop) on VelocityMesh
314
+OldDistanceToBottom (OldDistanceToBottom) on VelocityMesh
315
+IteratedDistanceToBottom (IteratedDistanceToBottom) on VelocityMesh
316
+OldGLSTurbulentKineticEnergy (OldGLSTurbulentKineticEnergy) on VelocityMesh
317
+IteratedGLSTurbulentKineticEnergy (IteratedGLSTurbulentKineticEnergy) on VelocityMesh
318
+OldGLSTurbulentKineticEnergySource (OldGLSTurbulentKineticEnergySource) on VelocityMesh
319
+IteratedGLSTurbulentKineticEnergySource (IteratedGLSTurbulentKineticEnergySource) on VelocityMesh
320
+OldGLSTurbulentKineticEnergyAbsorption (OldGLSTurbulentKineticEnergyAbsorption) on VelocityMesh
321
+IteratedGLSTurbulentKineticEnergyAbsorption (IteratedGLSTurbulentKineticEnergyAbsorption) on VelocityMesh
322
+OldGLSGenericSecondQuantity (OldGLSGenericSecondQuantity) on VelocityMesh
323
+IteratedGLSGenericSecondQuantity (IteratedGLSGenericSecondQuantity) on VelocityMesh
324
+OldGLSGenericSecondQuantitySource (OldGLSGenericSecondQuantitySource) on VelocityMesh
325
+IteratedGLSGenericSecondQuantitySource (IteratedGLSGenericSecondQuantitySource) on VelocityMesh
326
+OldGLSGenericSecondQuantityAbsorption (OldGLSGenericSecondQuantityAbsorption) on VelocityMesh
327
+IteratedGLSGenericSecondQuantityAbsorption (IteratedGLSGenericSecondQuantityAbsorption) on VelocityMesh
328
+OldGLSLengthScale (OldGLSLengthScale) on VelocityMesh
329
+IteratedGLSLengthScale (IteratedGLSLengthScale) on VelocityMesh
330
+OldGLSBuoyancyFrequency (OldGLSBuoyancyFrequency) on VelocityMesh
331
+IteratedGLSBuoyancyFrequency (IteratedGLSBuoyancyFrequency) on VelocityMesh
332
+OldGLSVelocityShear (OldGLSVelocityShear) on VelocityMesh
333
+IteratedGLSVelocityShear (IteratedGLSVelocityShear) on VelocityMesh
334
+OldGLSShearProduction (OldGLSShearProduction) on VelocityMesh
335
+IteratedGLSShearProduction (IteratedGLSShearProduction) on VelocityMesh
336
+OldGLSBuoyancyProduction (OldGLSBuoyancyProduction) on VelocityMesh
337
+IteratedGLSBuoyancyProduction (IteratedGLSBuoyancyProduction) on VelocityMesh
338
+OldGLSDissipationEpsilon (OldGLSDissipationEpsilon) on VelocityMesh
339
+IteratedGLSDissipationEpsilon (IteratedGLSDissipationEpsilon) on VelocityMesh
340
+OldGLSStabilityFunctionSM (OldGLSStabilityFunctionSM) on VelocityMesh
341
+IteratedGLSStabilityFunctionSM (IteratedGLSStabilityFunctionSM) on VelocityMesh
342
+OldGLSStabilityFunctionSH (OldGLSStabilityFunctionSH) on VelocityMesh
343
+IteratedGLSStabilityFunctionSH (IteratedGLSStabilityFunctionSH) on VelocityMesh
344
+OldGLSSource1 (OldGLSSource1) on VelocityMesh
345
+IteratedGLSSource1 (IteratedGLSSource1) on VelocityMesh
346
+OldGLSSource2 (OldGLSSource2) on VelocityMesh
347
+IteratedGLSSource2 (IteratedGLSSource2) on VelocityMesh
348
+OldGLSAbsorption1 (OldGLSAbsorption1) on VelocityMesh
349
+IteratedGLSAbsorption1 (IteratedGLSAbsorption1) on VelocityMesh
350
+OldGLSAbsorption2 (OldGLSAbsorption2) on VelocityMesh
351
+IteratedGLSAbsorption2 (IteratedGLSAbsorption2) on VelocityMesh
352
+OldGLSWallFunction (OldGLSWallFunction) on VelocityMesh
353
+IteratedGLSWallFunction (IteratedGLSWallFunction) on VelocityMesh
354
+OldGLSVerticalViscosity (OldGLSVerticalViscosity) on VelocityMesh
355
+IteratedGLSVerticalViscosity (IteratedGLSVerticalViscosity) on VelocityMesh
356
+OldGLSVerticalDiffusivity (OldGLSVerticalDiffusivity) on VelocityMesh
357
+IteratedGLSVerticalDiffusivity (IteratedGLSVerticalDiffusivity) on VelocityMesh
358
+VelocityBuoyancyDensity (VelocityBuoyancyDensity) on VelocityMesh
359
+Time (Time) on CoordinateMesh
361
+Coordinate (Coordinate) on CoordinateMesh
362
+Velocity (Velocity) on VelocityMesh
363
+VelocityAbsorption (VelocityAbsorption) on VelocityMesh
364
+GravityDirection (GravityDirection) on CoordinateMesh
365
+OldVelocity (OldVelocity) on VelocityMesh
366
+IteratedVelocity (IteratedVelocity) on VelocityMesh
367
+NonlinearVelocity (NonlinearVelocity) on VelocityMesh
368
+OldVelocityAbsorption (OldVelocityAbsorption) on VelocityMesh
369
+IteratedVelocityAbsorption (IteratedVelocityAbsorption) on VelocityMesh
370
+OldGravityDirection (OldGravityDirection) on CoordinateMesh
371
+IteratedGravityDirection (IteratedGravityDirection) on CoordinateMesh
372
+OriginalCoordinate (OriginalCoordinate) on CoordinateMesh
373
+OldCoordinate (OldCoordinate) on CoordinateMesh
374
+IteratedCoordinate (IteratedCoordinate) on CoordinateMesh
375
+GridVelocity (GridVelocity) on CoordinateMesh
377
+Viscosity (Viscosity) on VelocityMesh
378
+GLSTurbulentKineticEnergyDiffusivity (GLSTurbulentKineticEnergyDiffusivity) on VelocityMesh
379
+GLSGenericSecondQuantityDiffusivity (GLSGenericSecondQuantityDiffusivity) on VelocityMesh
380
+GLSBackgroundViscosity (GLSBackgroundViscosity) on VelocityMesh
381
+GLSBackgroundDiffusivity (GLSBackgroundDiffusivity) on VelocityMesh
382
+GLSEddyViscosityKM (GLSEddyViscosityKM) on VelocityMesh
383
+GLSEddyDiffusivityKH (GLSEddyDiffusivityKH) on VelocityMesh
384
+TemperatureDiffusivity (TemperatureDiffusivity) on VelocityMesh
385
+SalinityDiffusivity (SalinityDiffusivity) on VelocityMesh
392
In calculate_diagnostic_variables
393
In calculate_perturbation_density
394
In calculate_perturbation_density
395
Exiting calculate_diagnostic_variables
396
In calculate_diagnostic_variables
397
Exiting calculate_diagnostic_variables
398
In initialise_diagnostics
399
Exiting initialise_diagnostics
400
do_checkpoint returning .false.
402
In vtk_write_state_new_options
403
Exiting vtk_write_state_new_options
406
--------------------------------------------
407
cm0: 0.52646469697902409
408
kappa: 0.36474544142360915
409
p: 3.0000000000000000
410
m: 1.5000000000000000
411
n: -1.0000000000000000
412
sigma_k: 1.3000000000000000
413
sigma_psi: 1.0000000000000000
415
Using wall function: none
418
--------------------------------------------
419
Initialising the GLS surfaces required for BCs
420
In compute_lumped_mass
422
Calculating GLS stability functions
427
Current simulation time: 0.0000000000000000
429
Timestep size (dt): 15.000000000000000
430
In set_boundary_conditions
431
In set_prescribed_field_values
432
In enforce_discrete_properties
433
Exiting enforce_discrete_properties
435
Start of another nonlinear iteration; ITS,nonlinear_iterations= 1 2
437
Considering scalar field Density in state Fluid
438
Finished field Density in state Fluid
439
Considering scalar field Temperature in state Fluid
440
In solve_field_equation_cg
441
In assemble_advection_diffusion_cg
442
Exiting assemble_advection_diffusion_cg
444
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
445
Halo nodes are consistent with trailing receive ordering
446
Matrix assembly completed.
447
Inside setup_ksp_from_options
449
RHS assembly completed.
450
Assembling initial guess.
451
Initial guess assembly completed.
454
Exiting solve_field_equation_cg
455
Finished field Temperature in state Fluid
456
Considering scalar field Salinity in state Fluid
457
In solve_field_equation_cg
458
In assemble_advection_diffusion_cg
459
Exiting assemble_advection_diffusion_cg
461
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
462
Halo nodes are consistent with trailing receive ordering
463
Matrix assembly completed.
464
Inside setup_ksp_from_options
466
RHS assembly completed.
467
Assembling initial guess.
468
Initial guess assembly completed.
471
Exiting solve_field_equation_cg
472
Finished field Salinity in state Fluid
473
Considering scalar field PerturbationDensity in state Fluid
474
Finished field PerturbationDensity in state Fluid
475
Considering scalar field CFLNumber in state Fluid
476
Finished field CFLNumber in state Fluid
477
Considering scalar field DistanceToTop in state Fluid
478
Finished field DistanceToTop in state Fluid
479
Considering scalar field DistanceToBottom in state Fluid
480
Finished field DistanceToBottom in state Fluid
481
Considering scalar field GLSTurbulentKineticEnergy in state Fluid
483
Calculating GLS stability functions
485
In assemble_advectiondiffusion_m_cv
486
In find_upwind_values
487
In calculate_all_upwind_project
488
Exiting calculate_all_upwind_project
489
Exiting find_upwind_values
490
Exiting assemble_advectiondiffusion_m_cv
491
In assemble_field_eqn_cv
492
Exiting assemble_field_eqn_cv
494
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
495
Halo nodes are consistent with trailing receive ordering
496
Matrix assembly completed.
497
Inside setup_ksp_from_options
499
RHS assembly completed.
500
Assembling initial guess.
501
Initial guess assembly completed.
504
In assemble_advectiondiffusion_m_cv
505
In find_upwind_values
506
In calculate_all_upwind_project
507
Exiting calculate_all_upwind_project
508
Exiting find_upwind_values
509
Exiting assemble_advectiondiffusion_m_cv
510
In assemble_field_eqn_cv
511
Exiting assemble_field_eqn_cv
513
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
514
Halo nodes are consistent with trailing receive ordering
515
Matrix assembly completed.
516
Inside setup_ksp_from_options
518
RHS assembly completed.
519
Assembling initial guess.
520
Initial guess assembly completed.
523
In assemble_advectiondiffusion_m_cv
524
In find_upwind_values
525
In calculate_all_upwind_project
526
Exiting calculate_all_upwind_project
527
Exiting find_upwind_values
528
Exiting assemble_advectiondiffusion_m_cv
529
In assemble_field_eqn_cv
530
Exiting assemble_field_eqn_cv
532
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
533
Halo nodes are consistent with trailing receive ordering
534
Matrix assembly completed.
535
Inside setup_ksp_from_options
537
RHS assembly completed.
538
Assembling initial guess.
539
Initial guess assembly completed.
542
Finished field GLSTurbulentKineticEnergy in state Fluid
543
Considering scalar field GLSGenericSecondQuantity in state Fluid
545
In assemble_advectiondiffusion_m_cv
546
In find_upwind_values
547
In calculate_all_upwind_project
548
Exiting calculate_all_upwind_project
549
Exiting find_upwind_values
550
Exiting assemble_advectiondiffusion_m_cv
551
In assemble_field_eqn_cv
552
Exiting assemble_field_eqn_cv
554
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
555
Halo nodes are consistent with trailing receive ordering
556
Matrix assembly completed.
557
Inside setup_ksp_from_options
559
RHS assembly completed.
560
Assembling initial guess.
561
Initial guess assembly completed.
564
In assemble_advectiondiffusion_m_cv
565
In find_upwind_values
566
In calculate_all_upwind_project
567
Exiting calculate_all_upwind_project
568
Exiting find_upwind_values
569
Exiting assemble_advectiondiffusion_m_cv
570
In assemble_field_eqn_cv
571
Exiting assemble_field_eqn_cv
573
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
574
Halo nodes are consistent with trailing receive ordering
575
Matrix assembly completed.
576
Inside setup_ksp_from_options
578
RHS assembly completed.
579
Assembling initial guess.
580
Initial guess assembly completed.
583
In assemble_advectiondiffusion_m_cv
584
In find_upwind_values
585
In calculate_all_upwind_project
586
Exiting calculate_all_upwind_project
587
Exiting find_upwind_values
588
Exiting assemble_advectiondiffusion_m_cv
589
In assemble_field_eqn_cv
590
Exiting assemble_field_eqn_cv
592
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
593
Halo nodes are consistent with trailing receive ordering
594
Matrix assembly completed.
595
Inside setup_ksp_from_options
597
RHS assembly completed.
598
Assembling initial guess.
599
Initial guess assembly completed.
602
Finished field GLSGenericSecondQuantity in state Fluid
604
Entering solve_momentum
605
Entering get_pressure_options
606
Entering get_phase_submaterials
607
Number of sub-materials = 1
608
Exiting get_phase_submaterials
609
Entering calculate_momentum_diagnostics
610
Entering calculate_diagnostic_phase_volume_fraction
611
Exiting calculate_diagnostic_phase_volume_fraction
612
In calculate_perturbation_density
613
Exiting calculate_momentum_diagnostics
614
Entering get_velocity_options
615
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
616
Halo nodes are consistent with trailing receive ordering
617
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
618
Halo nodes are consistent with trailing receive ordering
619
Entering construct_momentum_cg
620
Initialising Coriolis
622
Entering advance_velocity
623
Inside setup_ksp_from_options
625
RHS assembly completed.
626
Assembling initial guess.
627
Initial guess assembly completed.
630
Entering finalise_state
632
Start of another nonlinear iteration; ITS,nonlinear_iterations= 2 2
634
Considering scalar field Density in state Fluid
635
Finished field Density in state Fluid
636
Considering scalar field Temperature in state Fluid
637
In solve_field_equation_cg
638
In assemble_advection_diffusion_cg
639
Exiting assemble_advection_diffusion_cg
641
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
642
Halo nodes are consistent with trailing receive ordering
643
Matrix assembly completed.
644
Inside setup_ksp_from_options
646
RHS assembly completed.
647
Assembling initial guess.
648
Initial guess assembly completed.
651
Exiting solve_field_equation_cg
652
Finished field Temperature in state Fluid
653
Considering scalar field Salinity in state Fluid
654
In solve_field_equation_cg
655
In assemble_advection_diffusion_cg
656
Exiting assemble_advection_diffusion_cg
658
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
659
Halo nodes are consistent with trailing receive ordering
660
Matrix assembly completed.
661
Inside setup_ksp_from_options
663
RHS assembly completed.
664
Assembling initial guess.
665
Initial guess assembly completed.
668
Exiting solve_field_equation_cg
669
Finished field Salinity in state Fluid
670
Considering scalar field PerturbationDensity in state Fluid
671
Finished field PerturbationDensity in state Fluid
672
Considering scalar field CFLNumber in state Fluid
673
Finished field CFLNumber in state Fluid
674
Considering scalar field DistanceToTop in state Fluid
675
Finished field DistanceToTop in state Fluid
676
Considering scalar field DistanceToBottom in state Fluid
677
Finished field DistanceToBottom in state Fluid
678
Considering scalar field GLSTurbulentKineticEnergy in state Fluid
680
Calculating GLS stability functions
682
In assemble_advectiondiffusion_m_cv
683
In find_upwind_values
684
In calculate_all_upwind_project
685
Exiting calculate_all_upwind_project
686
Exiting find_upwind_values
687
Exiting assemble_advectiondiffusion_m_cv
688
In assemble_field_eqn_cv
689
Exiting assemble_field_eqn_cv
691
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
692
Halo nodes are consistent with trailing receive ordering
693
Matrix assembly completed.
694
Inside setup_ksp_from_options
696
RHS assembly completed.
697
Assembling initial guess.
698
Initial guess assembly completed.
701
In assemble_advectiondiffusion_m_cv
702
In find_upwind_values
703
In calculate_all_upwind_project
704
Exiting calculate_all_upwind_project
705
Exiting find_upwind_values
706
Exiting assemble_advectiondiffusion_m_cv
707
In assemble_field_eqn_cv
708
Exiting assemble_field_eqn_cv
710
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
711
Halo nodes are consistent with trailing receive ordering
712
Matrix assembly completed.
713
Inside setup_ksp_from_options
715
RHS assembly completed.
716
Assembling initial guess.
717
Initial guess assembly completed.
720
In assemble_advectiondiffusion_m_cv
721
In find_upwind_values
722
In calculate_all_upwind_project
723
Exiting calculate_all_upwind_project
724
Exiting find_upwind_values
725
Exiting assemble_advectiondiffusion_m_cv
726
In assemble_field_eqn_cv
727
Exiting assemble_field_eqn_cv
729
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
730
Halo nodes are consistent with trailing receive ordering
731
Matrix assembly completed.
732
Inside setup_ksp_from_options
734
RHS assembly completed.
735
Assembling initial guess.
736
Initial guess assembly completed.
739
Finished field GLSTurbulentKineticEnergy in state Fluid
740
Considering scalar field GLSGenericSecondQuantity in state Fluid
742
In assemble_advectiondiffusion_m_cv
743
In find_upwind_values
744
In calculate_all_upwind_project
745
Exiting calculate_all_upwind_project
746
Exiting find_upwind_values
747
Exiting assemble_advectiondiffusion_m_cv
748
In assemble_field_eqn_cv
749
Exiting assemble_field_eqn_cv
751
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
752
Halo nodes are consistent with trailing receive ordering
753
Matrix assembly completed.
754
Inside setup_ksp_from_options
756
RHS assembly completed.
757
Assembling initial guess.
758
Initial guess assembly completed.
761
In assemble_advectiondiffusion_m_cv
762
In find_upwind_values
763
In calculate_all_upwind_project
764
Exiting calculate_all_upwind_project
765
Exiting find_upwind_values
766
Exiting assemble_advectiondiffusion_m_cv
767
In assemble_field_eqn_cv
768
Exiting assemble_field_eqn_cv
770
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
771
Halo nodes are consistent with trailing receive ordering
772
Matrix assembly completed.
773
Inside setup_ksp_from_options
775
RHS assembly completed.
776
Assembling initial guess.
777
Initial guess assembly completed.
780
In assemble_advectiondiffusion_m_cv
781
In find_upwind_values
782
In calculate_all_upwind_project
783
Exiting calculate_all_upwind_project
784
Exiting find_upwind_values
785
Exiting assemble_advectiondiffusion_m_cv
786
In assemble_field_eqn_cv
787
Exiting assemble_field_eqn_cv
789
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
790
Halo nodes are consistent with trailing receive ordering
791
Matrix assembly completed.
792
Inside setup_ksp_from_options
794
RHS assembly completed.
795
Assembling initial guess.
796
Initial guess assembly completed.
799
Finished field GLSGenericSecondQuantity in state Fluid
801
Entering solve_momentum
802
Entering get_pressure_options
803
Entering get_phase_submaterials
804
Number of sub-materials = 1
805
Exiting get_phase_submaterials
806
Entering calculate_momentum_diagnostics
807
Entering calculate_diagnostic_phase_volume_fraction
808
Exiting calculate_diagnostic_phase_volume_fraction
809
In calculate_perturbation_density
810
Exiting calculate_momentum_diagnostics
811
Entering get_velocity_options
812
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
813
Halo nodes are consistent with trailing receive ordering
814
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
815
Halo nodes are consistent with trailing receive ordering
816
Entering construct_momentum_cg
818
Entering advance_velocity
819
Inside setup_ksp_from_options
821
RHS assembly completed.
822
Assembling initial guess.
823
Initial guess assembly completed.
826
Entering finalise_state
827
In calculate_diagnostic_variables
828
In calculate_perturbation_density
829
In calculate_perturbation_density
830
Exiting calculate_diagnostic_variables
831
In calculate_diagnostic_variables
832
Exiting calculate_diagnostic_variables
834
Registered diagnostics:
836
In set_dirichlet_consistent
840
Current simulation time: 15.000000000000000
842
Timestep size (dt): 15.000000000000000
843
In set_boundary_conditions
844
In set_prescribed_field_values
845
In enforce_discrete_properties
846
Exiting enforce_discrete_properties
848
Start of another nonlinear iteration; ITS,nonlinear_iterations= 1 2
850
Considering scalar field Density in state Fluid
851
Finished field Density in state Fluid
852
Considering scalar field Temperature in state Fluid
853
In solve_field_equation_cg
854
In assemble_advection_diffusion_cg
855
Exiting assemble_advection_diffusion_cg
857
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
858
Halo nodes are consistent with trailing receive ordering
859
Matrix assembly completed.
860
Inside setup_ksp_from_options
862
RHS assembly completed.
863
Assembling initial guess.
864
Initial guess assembly completed.
867
Exiting solve_field_equation_cg
868
Finished field Temperature in state Fluid
869
Considering scalar field Salinity in state Fluid
870
In solve_field_equation_cg
871
In assemble_advection_diffusion_cg
872
Exiting assemble_advection_diffusion_cg
874
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
875
Halo nodes are consistent with trailing receive ordering
876
Matrix assembly completed.
877
Inside setup_ksp_from_options
879
RHS assembly completed.
880
Assembling initial guess.
881
Initial guess assembly completed.
884
Exiting solve_field_equation_cg
885
Finished field Salinity in state Fluid
886
Considering scalar field PerturbationDensity in state Fluid
887
Finished field PerturbationDensity in state Fluid
888
Considering scalar field CFLNumber in state Fluid
889
Finished field CFLNumber in state Fluid
890
Considering scalar field DistanceToTop in state Fluid
891
Finished field DistanceToTop in state Fluid
892
Considering scalar field DistanceToBottom in state Fluid
893
Finished field DistanceToBottom in state Fluid
894
Considering scalar field GLSTurbulentKineticEnergy in state Fluid
896
Calculating GLS stability functions
898
In assemble_advectiondiffusion_m_cv
899
In find_upwind_values
900
In calculate_all_upwind_project
901
Exiting calculate_all_upwind_project
902
Exiting find_upwind_values
903
Exiting assemble_advectiondiffusion_m_cv
904
In assemble_field_eqn_cv
905
Exiting assemble_field_eqn_cv
907
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
908
Halo nodes are consistent with trailing receive ordering
909
Matrix assembly completed.
910
Inside setup_ksp_from_options
912
RHS assembly completed.
913
Assembling initial guess.
914
Initial guess assembly completed.
917
In assemble_advectiondiffusion_m_cv
918
In find_upwind_values
919
In calculate_all_upwind_project
920
Exiting calculate_all_upwind_project
921
Exiting find_upwind_values
922
Exiting assemble_advectiondiffusion_m_cv
923
In assemble_field_eqn_cv
924
Exiting assemble_field_eqn_cv
926
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
927
Halo nodes are consistent with trailing receive ordering
928
Matrix assembly completed.
929
Inside setup_ksp_from_options
931
RHS assembly completed.
932
Assembling initial guess.
933
Initial guess assembly completed.
936
In assemble_advectiondiffusion_m_cv
937
In find_upwind_values
938
In calculate_all_upwind_project
939
Exiting calculate_all_upwind_project
940
Exiting find_upwind_values
941
Exiting assemble_advectiondiffusion_m_cv
942
In assemble_field_eqn_cv
943
Exiting assemble_field_eqn_cv
945
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
946
Halo nodes are consistent with trailing receive ordering
947
Matrix assembly completed.
948
Inside setup_ksp_from_options
950
RHS assembly completed.
951
Assembling initial guess.
952
Initial guess assembly completed.
955
Finished field GLSTurbulentKineticEnergy in state Fluid
956
Considering scalar field GLSGenericSecondQuantity in state Fluid
958
In assemble_advectiondiffusion_m_cv
959
In find_upwind_values
960
In calculate_all_upwind_project
961
Exiting calculate_all_upwind_project
962
Exiting find_upwind_values
963
Exiting assemble_advectiondiffusion_m_cv
964
In assemble_field_eqn_cv
965
Exiting assemble_field_eqn_cv
967
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
968
Halo nodes are consistent with trailing receive ordering
969
Matrix assembly completed.
970
Inside setup_ksp_from_options
972
RHS assembly completed.
973
Assembling initial guess.
974
Initial guess assembly completed.
977
In assemble_advectiondiffusion_m_cv
978
In find_upwind_values
979
In calculate_all_upwind_project
980
Exiting calculate_all_upwind_project
981
Exiting find_upwind_values
982
Exiting assemble_advectiondiffusion_m_cv
983
In assemble_field_eqn_cv
984
Exiting assemble_field_eqn_cv
986
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
987
Halo nodes are consistent with trailing receive ordering
988
Matrix assembly completed.
989
Inside setup_ksp_from_options
991
RHS assembly completed.
992
Assembling initial guess.
993
Initial guess assembly completed.
996
In assemble_advectiondiffusion_m_cv
997
In find_upwind_values
998
In calculate_all_upwind_project
999
Exiting calculate_all_upwind_project
1000
Exiting find_upwind_values
1001
Exiting assemble_advectiondiffusion_m_cv
1002
In assemble_field_eqn_cv
1003
Exiting assemble_field_eqn_cv
1005
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1006
Halo nodes are consistent with trailing receive ordering
1007
Matrix assembly completed.
1008
Inside setup_ksp_from_options
1010
RHS assembly completed.
1011
Assembling initial guess.
1012
Initial guess assembly completed.
1015
Finished field GLSGenericSecondQuantity in state Fluid
1017
Entering solve_momentum
1018
Entering get_pressure_options
1019
Entering get_phase_submaterials
1020
Number of sub-materials = 1
1021
Exiting get_phase_submaterials
1022
Entering calculate_momentum_diagnostics
1023
Entering calculate_diagnostic_phase_volume_fraction
1024
Exiting calculate_diagnostic_phase_volume_fraction
1025
In calculate_perturbation_density
1026
Exiting calculate_momentum_diagnostics
1027
Entering get_velocity_options
1028
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1029
Halo nodes are consistent with trailing receive ordering
1030
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1031
Halo nodes are consistent with trailing receive ordering
1032
Entering construct_momentum_cg
1034
Entering advance_velocity
1035
Inside setup_ksp_from_options
1037
RHS assembly completed.
1038
Assembling initial guess.
1039
Initial guess assembly completed.
1042
Entering finalise_state
1044
Start of another nonlinear iteration; ITS,nonlinear_iterations= 2 2
1046
Considering scalar field Density in state Fluid
1047
Finished field Density in state Fluid
1048
Considering scalar field Temperature in state Fluid
1049
In solve_field_equation_cg
1050
In assemble_advection_diffusion_cg
1051
Exiting assemble_advection_diffusion_cg
1053
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1054
Halo nodes are consistent with trailing receive ordering
1055
Matrix assembly completed.
1056
Inside setup_ksp_from_options
1058
RHS assembly completed.
1059
Assembling initial guess.
1060
Initial guess assembly completed.
1063
Exiting solve_field_equation_cg
1064
Finished field Temperature in state Fluid
1065
Considering scalar field Salinity in state Fluid
1066
In solve_field_equation_cg
1067
In assemble_advection_diffusion_cg
1068
Exiting assemble_advection_diffusion_cg
1070
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1071
Halo nodes are consistent with trailing receive ordering
1072
Matrix assembly completed.
1073
Inside setup_ksp_from_options
1075
RHS assembly completed.
1076
Assembling initial guess.
1077
Initial guess assembly completed.
1080
Exiting solve_field_equation_cg
1081
Finished field Salinity in state Fluid
1082
Considering scalar field PerturbationDensity in state Fluid
1083
Finished field PerturbationDensity in state Fluid
1084
Considering scalar field CFLNumber in state Fluid
1085
Finished field CFLNumber in state Fluid
1086
Considering scalar field DistanceToTop in state Fluid
1087
Finished field DistanceToTop in state Fluid
1088
Considering scalar field DistanceToBottom in state Fluid
1089
Finished field DistanceToBottom in state Fluid
1090
Considering scalar field GLSTurbulentKineticEnergy in state Fluid
1092
Calculating GLS stability functions
1094
In assemble_advectiondiffusion_m_cv
1095
In find_upwind_values
1096
In calculate_all_upwind_project
1097
Exiting calculate_all_upwind_project
1098
Exiting find_upwind_values
1099
Exiting assemble_advectiondiffusion_m_cv
1100
In assemble_field_eqn_cv
1101
Exiting assemble_field_eqn_cv
1103
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1104
Halo nodes are consistent with trailing receive ordering
1105
Matrix assembly completed.
1106
Inside setup_ksp_from_options
1108
RHS assembly completed.
1109
Assembling initial guess.
1110
Initial guess assembly completed.
1113
In assemble_advectiondiffusion_m_cv
1114
In find_upwind_values
1115
In calculate_all_upwind_project
1116
Exiting calculate_all_upwind_project
1117
Exiting find_upwind_values
1118
Exiting assemble_advectiondiffusion_m_cv
1119
In assemble_field_eqn_cv
1120
Exiting assemble_field_eqn_cv
1122
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1123
Halo nodes are consistent with trailing receive ordering
1124
Matrix assembly completed.
1125
Inside setup_ksp_from_options
1127
RHS assembly completed.
1128
Assembling initial guess.
1129
Initial guess assembly completed.
1132
In assemble_advectiondiffusion_m_cv
1133
In find_upwind_values
1134
In calculate_all_upwind_project
1135
Exiting calculate_all_upwind_project
1136
Exiting find_upwind_values
1137
Exiting assemble_advectiondiffusion_m_cv
1138
In assemble_field_eqn_cv
1139
Exiting assemble_field_eqn_cv
1141
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1142
Halo nodes are consistent with trailing receive ordering
1143
Matrix assembly completed.
1144
Inside setup_ksp_from_options
1146
RHS assembly completed.
1147
Assembling initial guess.
1148
Initial guess assembly completed.
1151
Finished field GLSTurbulentKineticEnergy in state Fluid
1152
Considering scalar field GLSGenericSecondQuantity in state Fluid
1154
In assemble_advectiondiffusion_m_cv
1155
In find_upwind_values
1156
In calculate_all_upwind_project
1157
Exiting calculate_all_upwind_project
1158
Exiting find_upwind_values
1159
Exiting assemble_advectiondiffusion_m_cv
1160
In assemble_field_eqn_cv
1161
Exiting assemble_field_eqn_cv
1163
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1164
Halo nodes are consistent with trailing receive ordering
1165
Matrix assembly completed.
1166
Inside setup_ksp_from_options
1168
RHS assembly completed.
1169
Assembling initial guess.
1170
Initial guess assembly completed.
1173
In assemble_advectiondiffusion_m_cv
1174
In find_upwind_values
1175
In calculate_all_upwind_project
1176
Exiting calculate_all_upwind_project
1177
Exiting find_upwind_values
1178
Exiting assemble_advectiondiffusion_m_cv
1179
In assemble_field_eqn_cv
1180
Exiting assemble_field_eqn_cv
1182
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1183
Halo nodes are consistent with trailing receive ordering
1184
Matrix assembly completed.
1185
Inside setup_ksp_from_options
1187
RHS assembly completed.
1188
Assembling initial guess.
1189
Initial guess assembly completed.
1192
In assemble_advectiondiffusion_m_cv
1193
In find_upwind_values
1194
In calculate_all_upwind_project
1195
Exiting calculate_all_upwind_project
1196
Exiting find_upwind_values
1197
Exiting assemble_advectiondiffusion_m_cv
1198
In assemble_field_eqn_cv
1199
Exiting assemble_field_eqn_cv
1201
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1202
Halo nodes are consistent with trailing receive ordering
1203
Matrix assembly completed.
1204
Inside setup_ksp_from_options
1206
RHS assembly completed.
1207
Assembling initial guess.
1208
Initial guess assembly completed.
1211
Finished field GLSGenericSecondQuantity in state Fluid
1213
Entering solve_momentum
1214
Entering get_pressure_options
1215
Entering get_phase_submaterials
1216
Number of sub-materials = 1
1217
Exiting get_phase_submaterials
1218
Entering calculate_momentum_diagnostics
1219
Entering calculate_diagnostic_phase_volume_fraction
1220
Exiting calculate_diagnostic_phase_volume_fraction
1221
In calculate_perturbation_density
1222
Exiting calculate_momentum_diagnostics
1223
Entering get_velocity_options
1224
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1225
Halo nodes are consistent with trailing receive ordering
1226
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1227
Halo nodes are consistent with trailing receive ordering
1228
Entering construct_momentum_cg
1230
Entering advance_velocity
1231
Inside setup_ksp_from_options
1233
RHS assembly completed.
1234
Assembling initial guess.
1235
Initial guess assembly completed.
1238
Entering finalise_state
1239
In calculate_diagnostic_variables
1240
In calculate_perturbation_density
1241
In calculate_perturbation_density
1242
Exiting calculate_diagnostic_variables
1243
In calculate_diagnostic_variables
1244
Exiting calculate_diagnostic_variables
1245
In write_diagnostics
1246
Registered diagnostics:
1248
In set_dirichlet_consistent
1249
********************
1250
*** NEW TIMESTEP ***
1251
********************
1252
Current simulation time: 30.000000000000000
1254
Timestep size (dt): 15.000000000000000
1255
In set_boundary_conditions
1256
In set_prescribed_field_values
1257
In enforce_discrete_properties
1258
Exiting enforce_discrete_properties
1260
Start of another nonlinear iteration; ITS,nonlinear_iterations= 1 2
1262
Considering scalar field Density in state Fluid
1263
Finished field Density in state Fluid
1264
Considering scalar field Temperature in state Fluid
1265
In solve_field_equation_cg
1266
In assemble_advection_diffusion_cg
1267
Exiting assemble_advection_diffusion_cg
1269
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1270
Halo nodes are consistent with trailing receive ordering
1271
Matrix assembly completed.
1272
Inside setup_ksp_from_options
1274
RHS assembly completed.
1275
Assembling initial guess.
1276
Initial guess assembly completed.
1279
Exiting solve_field_equation_cg
1280
Finished field Temperature in state Fluid
1281
Considering scalar field Salinity in state Fluid
1282
In solve_field_equation_cg
1283
In assemble_advection_diffusion_cg
1284
Exiting assemble_advection_diffusion_cg
1286
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1287
Halo nodes are consistent with trailing receive ordering
1288
Matrix assembly completed.
1289
Inside setup_ksp_from_options
1291
RHS assembly completed.
1292
Assembling initial guess.
1293
Initial guess assembly completed.
1296
Exiting solve_field_equation_cg
1297
Finished field Salinity in state Fluid
1298
Considering scalar field PerturbationDensity in state Fluid
1299
Finished field PerturbationDensity in state Fluid
1300
Considering scalar field CFLNumber in state Fluid
1301
Finished field CFLNumber in state Fluid
1302
Considering scalar field DistanceToTop in state Fluid
1303
Finished field DistanceToTop in state Fluid
1304
Considering scalar field DistanceToBottom in state Fluid
1305
Finished field DistanceToBottom in state Fluid
1306
Considering scalar field GLSTurbulentKineticEnergy in state Fluid
1308
Calculating GLS stability functions
1310
In assemble_advectiondiffusion_m_cv
1311
In find_upwind_values
1312
In calculate_all_upwind_project
1313
Exiting calculate_all_upwind_project
1314
Exiting find_upwind_values
1315
Exiting assemble_advectiondiffusion_m_cv
1316
In assemble_field_eqn_cv
1317
Exiting assemble_field_eqn_cv
1319
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1320
Halo nodes are consistent with trailing receive ordering
1321
Matrix assembly completed.
1322
Inside setup_ksp_from_options
1324
RHS assembly completed.
1325
Assembling initial guess.
1326
Initial guess assembly completed.
1329
In assemble_advectiondiffusion_m_cv
1330
In find_upwind_values
1331
In calculate_all_upwind_project
1332
Exiting calculate_all_upwind_project
1333
Exiting find_upwind_values
1334
Exiting assemble_advectiondiffusion_m_cv
1335
In assemble_field_eqn_cv
1336
Exiting assemble_field_eqn_cv
1338
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1339
Halo nodes are consistent with trailing receive ordering
1340
Matrix assembly completed.
1341
Inside setup_ksp_from_options
1343
RHS assembly completed.
1344
Assembling initial guess.
1345
Initial guess assembly completed.
1348
In assemble_advectiondiffusion_m_cv
1349
In find_upwind_values
1350
In calculate_all_upwind_project
1351
Exiting calculate_all_upwind_project
1352
Exiting find_upwind_values
1353
Exiting assemble_advectiondiffusion_m_cv
1354
In assemble_field_eqn_cv
1355
Exiting assemble_field_eqn_cv
1357
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1358
Halo nodes are consistent with trailing receive ordering
1359
Matrix assembly completed.
1360
Inside setup_ksp_from_options
1362
RHS assembly completed.
1363
Assembling initial guess.
1364
Initial guess assembly completed.
1367
Finished field GLSTurbulentKineticEnergy in state Fluid
1368
Considering scalar field GLSGenericSecondQuantity in state Fluid
1370
In assemble_advectiondiffusion_m_cv
1371
In find_upwind_values
1372
In calculate_all_upwind_project
1373
Exiting calculate_all_upwind_project
1374
Exiting find_upwind_values
1375
Exiting assemble_advectiondiffusion_m_cv
1376
In assemble_field_eqn_cv
1377
Exiting assemble_field_eqn_cv
1379
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1380
Halo nodes are consistent with trailing receive ordering
1381
Matrix assembly completed.
1382
Inside setup_ksp_from_options
1384
RHS assembly completed.
1385
Assembling initial guess.
1386
Initial guess assembly completed.
1389
In assemble_advectiondiffusion_m_cv
1390
In find_upwind_values
1391
In calculate_all_upwind_project
1392
Exiting calculate_all_upwind_project
1393
Exiting find_upwind_values
1394
Exiting assemble_advectiondiffusion_m_cv
1395
In assemble_field_eqn_cv
1396
Exiting assemble_field_eqn_cv
1398
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1399
Halo nodes are consistent with trailing receive ordering
1400
Matrix assembly completed.
1401
Inside setup_ksp_from_options
1403
RHS assembly completed.
1404
Assembling initial guess.
1405
Initial guess assembly completed.
1408
In assemble_advectiondiffusion_m_cv
1409
In find_upwind_values
1410
In calculate_all_upwind_project
1411
Exiting calculate_all_upwind_project
1412
Exiting find_upwind_values
1413
Exiting assemble_advectiondiffusion_m_cv
1414
In assemble_field_eqn_cv
1415
Exiting assemble_field_eqn_cv
1417
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1418
Halo nodes are consistent with trailing receive ordering
1419
Matrix assembly completed.
1420
Inside setup_ksp_from_options
1422
RHS assembly completed.
1423
Assembling initial guess.
1424
Initial guess assembly completed.
1427
Finished field GLSGenericSecondQuantity in state Fluid
1429
Entering solve_momentum
1430
Entering get_pressure_options
1431
Entering get_phase_submaterials
1432
Number of sub-materials = 1
1433
Exiting get_phase_submaterials
1434
Entering calculate_momentum_diagnostics
1435
Entering calculate_diagnostic_phase_volume_fraction
1436
Exiting calculate_diagnostic_phase_volume_fraction
1437
In calculate_perturbation_density
1438
Exiting calculate_momentum_diagnostics
1439
Entering get_velocity_options
1440
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1441
Halo nodes are consistent with trailing receive ordering
1442
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1443
Halo nodes are consistent with trailing receive ordering
1444
Entering construct_momentum_cg
1446
Entering advance_velocity
1447
Inside setup_ksp_from_options
1449
RHS assembly completed.
1450
Assembling initial guess.
1451
Initial guess assembly completed.
1454
Entering finalise_state
1456
Start of another nonlinear iteration; ITS,nonlinear_iterations= 2 2
1458
Considering scalar field Density in state Fluid
1459
Finished field Density in state Fluid
1460
Considering scalar field Temperature in state Fluid
1461
In solve_field_equation_cg
1462
In assemble_advection_diffusion_cg
1463
Exiting assemble_advection_diffusion_cg
1465
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1466
Halo nodes are consistent with trailing receive ordering
1467
Matrix assembly completed.
1468
Inside setup_ksp_from_options
1470
RHS assembly completed.
1471
Assembling initial guess.
1472
Initial guess assembly completed.
1475
Exiting solve_field_equation_cg
1476
Finished field Temperature in state Fluid
1477
Considering scalar field Salinity in state Fluid
1478
In solve_field_equation_cg
1479
In assemble_advection_diffusion_cg
1480
Exiting assemble_advection_diffusion_cg
1482
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1483
Halo nodes are consistent with trailing receive ordering
1484
Matrix assembly completed.
1485
Inside setup_ksp_from_options
1487
RHS assembly completed.
1488
Assembling initial guess.
1489
Initial guess assembly completed.
1492
Exiting solve_field_equation_cg
1493
Finished field Salinity in state Fluid
1494
Considering scalar field PerturbationDensity in state Fluid
1495
Finished field PerturbationDensity in state Fluid
1496
Considering scalar field CFLNumber in state Fluid
1497
Finished field CFLNumber in state Fluid
1498
Considering scalar field DistanceToTop in state Fluid
1499
Finished field DistanceToTop in state Fluid
1500
Considering scalar field DistanceToBottom in state Fluid
1501
Finished field DistanceToBottom in state Fluid
1502
Considering scalar field GLSTurbulentKineticEnergy in state Fluid
1504
Calculating GLS stability functions
1506
In assemble_advectiondiffusion_m_cv
1507
In find_upwind_values
1508
In calculate_all_upwind_project
1509
Exiting calculate_all_upwind_project
1510
Exiting find_upwind_values
1511
Exiting assemble_advectiondiffusion_m_cv
1512
In assemble_field_eqn_cv
1513
Exiting assemble_field_eqn_cv
1515
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1516
Halo nodes are consistent with trailing receive ordering
1517
Matrix assembly completed.
1518
Inside setup_ksp_from_options
1520
RHS assembly completed.
1521
Assembling initial guess.
1522
Initial guess assembly completed.
1525
In assemble_advectiondiffusion_m_cv
1526
In find_upwind_values
1527
In calculate_all_upwind_project
1528
Exiting calculate_all_upwind_project
1529
Exiting find_upwind_values
1530
Exiting assemble_advectiondiffusion_m_cv
1531
In assemble_field_eqn_cv
1532
Exiting assemble_field_eqn_cv
1534
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1535
Halo nodes are consistent with trailing receive ordering
1536
Matrix assembly completed.
1537
Inside setup_ksp_from_options
1539
RHS assembly completed.
1540
Assembling initial guess.
1541
Initial guess assembly completed.
1544
In assemble_advectiondiffusion_m_cv
1545
In find_upwind_values
1546
In calculate_all_upwind_project
1547
Exiting calculate_all_upwind_project
1548
Exiting find_upwind_values
1549
Exiting assemble_advectiondiffusion_m_cv
1550
In assemble_field_eqn_cv
1551
Exiting assemble_field_eqn_cv
1553
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1554
Halo nodes are consistent with trailing receive ordering
1555
Matrix assembly completed.
1556
Inside setup_ksp_from_options
1558
RHS assembly completed.
1559
Assembling initial guess.
1560
Initial guess assembly completed.
1563
Finished field GLSTurbulentKineticEnergy in state Fluid
1564
Considering scalar field GLSGenericSecondQuantity in state Fluid
1566
In assemble_advectiondiffusion_m_cv
1567
In find_upwind_values
1568
In calculate_all_upwind_project
1569
Exiting calculate_all_upwind_project
1570
Exiting find_upwind_values
1571
Exiting assemble_advectiondiffusion_m_cv
1572
In assemble_field_eqn_cv
1573
Exiting assemble_field_eqn_cv
1575
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1576
Halo nodes are consistent with trailing receive ordering
1577
Matrix assembly completed.
1578
Inside setup_ksp_from_options
1580
RHS assembly completed.
1581
Assembling initial guess.
1582
Initial guess assembly completed.
1585
In assemble_advectiondiffusion_m_cv
1586
In find_upwind_values
1587
In calculate_all_upwind_project
1588
Exiting calculate_all_upwind_project
1589
Exiting find_upwind_values
1590
Exiting assemble_advectiondiffusion_m_cv
1591
In assemble_field_eqn_cv
1592
Exiting assemble_field_eqn_cv
1594
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1595
Halo nodes are consistent with trailing receive ordering
1596
Matrix assembly completed.
1597
Inside setup_ksp_from_options
1599
RHS assembly completed.
1600
Assembling initial guess.
1601
Initial guess assembly completed.
1604
In assemble_advectiondiffusion_m_cv
1605
In find_upwind_values
1606
In calculate_all_upwind_project
1607
Exiting calculate_all_upwind_project
1608
Exiting find_upwind_values
1609
Exiting assemble_advectiondiffusion_m_cv
1610
In assemble_field_eqn_cv
1611
Exiting assemble_field_eqn_cv
1613
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1614
Halo nodes are consistent with trailing receive ordering
1615
Matrix assembly completed.
1616
Inside setup_ksp_from_options
1618
RHS assembly completed.
1619
Assembling initial guess.
1620
Initial guess assembly completed.
1623
Finished field GLSGenericSecondQuantity in state Fluid
1625
Entering solve_momentum
1626
Entering get_pressure_options
1627
Entering get_phase_submaterials
1628
Number of sub-materials = 1
1629
Exiting get_phase_submaterials
1630
Entering calculate_momentum_diagnostics
1631
Entering calculate_diagnostic_phase_volume_fraction
1632
Exiting calculate_diagnostic_phase_volume_fraction
1633
In calculate_perturbation_density
1634
Exiting calculate_momentum_diagnostics
1635
Entering get_velocity_options
1636
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1637
Halo nodes are consistent with trailing receive ordering
1638
Checking nodes in halo CoordinateMeshLevel1Halo for consistency with trailing receive ordering
1639
Halo nodes are consistent with trailing receive ordering
1640
Entering construct_momentum_cg
1642
Entering advance_velocity
1643
Inside setup_ksp_from_options
1645
RHS assembly completed.
1646
Assembling initial guess.
1647
Initial guess assembly completed.
1650
Entering finalise_state
1651
In calculate_diagnostic_variables
1652
In calculate_perturbation_density
1653
In calculate_perturbation_density
1654
Exiting calculate_diagnostic_variables
1655
In calculate_diagnostic_variables
1656
Exiting calculate_diagnostic_variables
1657
In write_diagnostics
1658
Registered diagnostics:
1660
In set_dirichlet_consistent
1663
In vtk_write_state_new_options
1664
Exiting vtk_write_state_new_options
1666
Cleaning up GLS variables
1667
Cleaning up GLS surface variables
1668
Finished gls_cleanup
1669
Printing references before final deallocation
1670
Printing out all currently allocated references:
1671
mesh_type VelocityMeshDGSurfaceMesh has reference count 1
1673
csr_sparsity VelocityMeshVelocityMeshSparsity has reference count 1
1675
scalar_field VelocityMeshLumpedMass has reference count 1
1677
scalar_field Time has reference count 1
1679
vector_field GridVelocity has reference count 1
1681
scalar_field VelocityBuoyancyDensity has reference count 1
1683
vector_field NonlinearVelocity has reference count 1
1685
vector_field IteratedVelocity has reference count 1
1687
vector_field OldVelocity has reference count 1
1689
scalar_field IteratedGLSVerticalDiffusivity has reference count 1
1691
scalar_field OldGLSVerticalDiffusivity has reference count 1
1693
scalar_field IteratedGLSVerticalViscosity has reference count 1
1695
scalar_field OldGLSVerticalViscosity has reference count 1
1697
scalar_field IteratedGLSWallFunction has reference count 1
1699
scalar_field OldGLSWallFunction has reference count 1
1701
scalar_field IteratedGLSAbsorption2 has reference count 1
1703
scalar_field OldGLSAbsorption2 has reference count 1
1705
scalar_field IteratedGLSAbsorption1 has reference count 1
1707
scalar_field OldGLSAbsorption1 has reference count 1
1709
scalar_field IteratedGLSSource2 has reference count 1
1711
scalar_field OldGLSSource2 has reference count 1
1713
scalar_field IteratedGLSSource1 has reference count 1
1715
scalar_field OldGLSSource1 has reference count 1
1717
scalar_field IteratedGLSStabilityFunctionSH has reference count 1
1719
scalar_field OldGLSStabilityFunctionSH has reference count 1
1721
scalar_field IteratedGLSStabilityFunctionSM has reference count 1
1723
scalar_field OldGLSStabilityFunctionSM has reference count 1
1725
scalar_field IteratedGLSDissipationEpsilon has reference count 1
1727
scalar_field OldGLSDissipationEpsilon has reference count 1
1729
scalar_field IteratedGLSBuoyancyProduction has reference count 1
1731
scalar_field OldGLSBuoyancyProduction has reference count 1
1733
scalar_field IteratedGLSShearProduction has reference count 1
1735
scalar_field OldGLSShearProduction has reference count 1
1737
scalar_field IteratedGLSVelocityShear has reference count 1
1739
scalar_field OldGLSVelocityShear has reference count 1
1741
scalar_field IteratedGLSBuoyancyFrequency has reference count 1
1743
scalar_field OldGLSBuoyancyFrequency has reference count 1
1745
scalar_field IteratedGLSLengthScale has reference count 1
1747
scalar_field OldGLSLengthScale has reference count 1
1749
scalar_field OldGLSGenericSecondQuantityAbsorption has reference count 1
1751
scalar_field OldGLSGenericSecondQuantitySource has reference count 1
1753
scalar_field IteratedGLSGenericSecondQuantity has reference count 1
1755
scalar_field OldGLSGenericSecondQuantity has reference count 1
1757
scalar_field OldGLSTurbulentKineticEnergyAbsorption has reference count 1
1759
scalar_field OldGLSTurbulentKineticEnergySource has reference count 1
1761
scalar_field IteratedGLSTurbulentKineticEnergy has reference count 1
1763
scalar_field OldGLSTurbulentKineticEnergy has reference count 1
1765
scalar_field IteratedCFLNumber has reference count 1
1767
scalar_field OldCFLNumber has reference count 1
1769
scalar_field IteratedPerturbationDensity has reference count 1
1771
scalar_field OldPerturbationDensity has reference count 1
1773
scalar_field IteratedSalinity has reference count 1
1775
scalar_field OldSalinity has reference count 1
1777
scalar_field IteratedTemperature has reference count 1
1779
scalar_field OldTemperature has reference count 1
1781
scalar_field IteratedDensity has reference count 1
1783
scalar_field OldDensity has reference count 1
1785
scalar_field value has reference count 1
1787
scalar_field value has reference count 1
1789
mesh_type psi_bottom_boundaryMesh has reference count 2
1791
mesh_type tke_bottom_boundaryMesh has reference count 2
1793
scalar_field value has reference count 1
1795
scalar_field value has reference count 1
1797
mesh_type psi_top_boundaryMesh has reference count 2
1799
mesh_type tke_top_boundaryMesh has reference count 2
1801
vector_field WindSurfaceField has reference count 1
1803
mesh_type BulkForcingVelocityMesh has reference count 2
1805
vector_field value has reference count 1
1807
mesh_type topbottomMesh has reference count 2
1809
vector_field value has reference count 1
1811
mesh_type lateralMesh has reference count 2
1813
scalar_field value has reference count 1
1815
mesh_type BulkForcingSalinityMesh has reference count 2
1817
scalar_field value has reference count 1
1819
mesh_type BulkForcingTemperatureMesh has reference count 2
1821
tensor_field GLSEddyDiffusivityKH has reference count 3
1823
tensor_field GLSEddyViscosityKM has reference count 1
1825
tensor_field GLSBackgroundDiffusivity has reference count 1
1827
tensor_field GLSBackgroundViscosity has reference count 1
1829
scalar_field GLSVerticalDiffusivity has reference count 1
1831
scalar_field GLSVerticalViscosity has reference count 1
1833
scalar_field GLSWallFunction has reference count 1
1835
scalar_field GLSAbsorption2 has reference count 1
1837
scalar_field GLSAbsorption1 has reference count 1
1839
scalar_field GLSSource2 has reference count 1
1841
scalar_field GLSSource1 has reference count 1
1843
scalar_field GLSStabilityFunctionSH has reference count 1
1845
scalar_field GLSStabilityFunctionSM has reference count 1
1847
scalar_field GLSDissipationEpsilon has reference count 1
1849
scalar_field GLSBuoyancyProduction has reference count 1
1851
scalar_field GLSShearProduction has reference count 1
1853
scalar_field GLSVelocityShear has reference count 1
1855
scalar_field GLSBuoyancyFrequency has reference count 1
1857
scalar_field GLSLengthScale has reference count 1
1859
tensor_field GLSGenericSecondQuantityDiffusivity has reference count 1
1861
scalar_field GLSGenericSecondQuantityAbsorption has reference count 2
1863
scalar_field GLSGenericSecondQuantitySource has reference count 2
1865
scalar_field GLSGenericSecondQuantity has reference count 1
1867
tensor_field GLSTurbulentKineticEnergyDiffusivity has reference count 1
1869
scalar_field GLSTurbulentKineticEnergyAbsorption has reference count 2
1871
scalar_field GLSTurbulentKineticEnergySource has reference count 2
1873
scalar_field GLSTurbulentKineticEnergy has reference count 1
1875
vector_field GravityDirection has reference count 3
1877
tensor_field Viscosity has reference count 1
1879
vector_field VelocityAbsorption has reference count 3
1881
vector_field Velocity has reference count 1
1883
scalar_field DistanceToBottom has reference count 3
1885
scalar_field DistanceToTop has reference count 3
1887
scalar_field CFLNumber has reference count 1
1889
scalar_field PerturbationDensity has reference count 1
1891
scalar_field Salinity has reference count 1
1893
scalar_field Temperature has reference count 1
1895
scalar_field Density has reference count 1
1897
csr_sparsity EEListSparsity has reference count 2
1899
csr_sparsity NEListSparsity has reference count 1
1901
halo_type CoordinateMeshMaximalElementHalo has reference count 2
1903
halo_type CoordinateMeshLevel2Halo has reference count 1
1905
halo_type CoordinateMeshLevel1Halo has reference count 3
1907
mesh_type SurfaceCoordinateMesh has reference count 1
1909
element_type has reference count 12
1911
quadrature_type has reference count 2
1913
csr_matrix CoordinateMeshFaceList has reference count 1
1915
csr_sparsity EEListSparsity has reference count 2
1917
csr_sparsity NEListSparsity has reference count 1
1919
vector_field Coordinate has reference count 4
1921
mesh_type CoordinateMesh has reference count 97
1923
element_type has reference count 1
1925
quadrature_type has reference count 1
1927
Printing references after final deallocation
1928
Printing out all currently allocated references:
1929
There are no references left.
1930
No registered memory in use.