Reading calculation conditions

Reads calculation conditions from the CGNS file.

Table 61 Subroutines to use

Subroutine

Remarks

cg_iric_read_integer_f

Reads an integer calculation-condition value

cg_iric_read_real_f

Reads a double-precision real calculation-condition value

cg_iric_read_realsingle_f

Reads a single-precision real calculation-condition value

cg_iric_read_string_f

Reads a string calculation-condition value

cg_iric_read_functionalsize_f

Checks the size of a functional-type calculation condition

cg_iric_read_functional_f

Reads functional calculation condition data in double-precision real type

cg_iric_read_functional_realsingle_f

Reads functional calculation condition data in single-precision real type

cg_iric_read_functionalwithname_f

Reads functional calculation condition data (with multiple values)

For reading calculation condition data other than in functional type, a subroutine reads a single calculation condition. An example of reading an integer calculation condition value is shown in List 93.

List 93 Example of source code to read calculation conditions
 1program Sample1
 2  implicit none
 3  include 'cgnslib_f.h'
 4
 5  integer:: fin, ier, i_flow
 6
 7  ! Open CGNS file
 8  call cg_open_f('test.cgn', CG_MODE_MODIFY, fin, ier)
 9  if (ier /=0) STOP "*** Open error of CGNS file ***"
10
11  ! Initialize iRIClib
12  call cg_iric_init_f(fin, ier)
13  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
14
15  call cg_iric_read_integer_f('i_flow', i_flow, ier)
16  print *, i_flow;
17
18  ! Close CGNS file
19  call cg_close_f(fin, ier)
20  stop
21end program Sample1

In contrast, for getting functional-type calculation conditions, it is necessary to use two subroutines: cg_iric_read_functionalsize_f and cg_iric_read_functional_f. An example of getting functional-type calculation condition data is shown in List 94.

List 94 Example of source code to read functional-type calculation conditions
 1program Sample2
 2  implicit none
 3  include 'cgnslib_f.h'
 4
 5  integer:: fin, ier, discharge_size, i
 6  double precision, dimension(:), allocatable:: discharge_time, discharge_value ! Array for storing discharge time and discharge value
 7
 8  ! Open CGNS file
 9  call cg_open_f('test.cgn', CG_MODE_MODIFY, fin, ier)
10  if (ier /=0) STOP "*** Open error of CGNS file ***"
11
12  ! Initialize iRIClib
13  call cg_iric_init_f(fin, ier)
14  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
15
16  ! First, check the size of the functional-type input conditions
17  call cg_iric_read_functionalsize_f('discharge', discharge_size, ier)
18  ! Allocate memory
19  allocate(discharge_time(discharge_size), discharge_value(discharge_size))
20  ! Read values into the allocated memory
21  call cg_iric_read_functional_f('discharge', discharge_time, discharge_value, ier)
22
23  ! (Output)
24  if (ier ==0) then
25    print *, 'discharge: discharge_size=', discharge_size
26    do i = 1, min(discharge_size, 5)
27      print *, ' i,time,value:', i, discharge_time(i), discharge_value(i)
28    end do
29  end if
30
31  ! Deallocate memory that has been allocated
32  deallocate(discharge_time, discharge_value)
33
34  ! Close CGNS file
35  call cg_close_f(fin, ier)
36  stop
37end program Sample2

Refer to Examples of calculation conditions, boundary conditions, and grid generating condition for examples of codes to load calculation conditions (or grid generating conditions).