.. _iriclib_output_result_gridedge:

格子エッジごとに値を持つ計算結果
=================================

格子エッジごとに値を持つ計算結果を出力する場合、
:numref:`table_iriclib_output_gridedge_functions` に示す関数を
使用します。

出力するプログラムの例は、:numref:`example_output_calc_result_gridedge`
を参照してください。

.. _table_iriclib_output_gridedge_functions:

.. list-table:: 格子エッジごとに値を持つ計算結果の出力に利用する関数
   :header-rows: 1

   * - 関数
     - 備考
   * - cg_iric_write_sol_iface_integer_f
     - 整数のI方向の格子エッジごとに値を持つ計算結果を出力する
   * - cg_iric_write_sol_iface_real_f
     - 倍精度実数のI方向の格子エッジごとに値を持つ計算結果を出力する
   * - cg_iric_write_sol_jface_integer_f
     - 整数のJ方向の格子エッジごとに値を持つ計算結果を出力する
   * - cg_iric_write_sol_jface_real_f
     - 倍精度実数のJ方向の格子エッジごとに値を持つ計算結果を出力する

.. code-block:: fortran
   :caption: サンプルプログラム (格子エッジごとに値を持つ計算結果)
   :name: example_output_calc_result_gridedge
   :linenos:

   program SampleProgram
     implicit none
     include 'cgnslib_f.h'

     integer:: fin, ier, isize, jsize
     integer:: canceled
     integer:: locked
     double precision:: time
     double precision, dimension(:,:), allocatable::grid_x, grid_y
     double precision, dimension(:,:), allocatable:: fluxi, fluxj
     character(len=20):: condFile

     condFile = 'test.cgn'

     ! CGNS ファイルのオープン
     call cg_open_f(condFile, CG_MODE_MODIFY, fin, ier)
     if (ier /=0) STOP "*** Open error of CGNS file ***"

     ! 内部変数の初期化
     call cg_iric_init_f(fin, ier)
     if (ier /=0) STOP "*** Initialize error of CGNS file ***"

     ! 格子のサイズを調べる
     call cg_iric_gotogridcoord2d_f(isize, jsize, ier)
     ! 格子を読み込むためのメモリを確保
     allocate(grid_x(isize, jsize), grid_y(isize, jsize))
     ! 計算結果を保持するメモリも確保
     allocate(fluxi(isize, jsize - 1), fluxj(isize - 1, jsize))
     ! 格子を読み込む
     call cg_iric_getgridcoord2d_f (grid_x, grid_y, ier)

     ! 初期状態の情報を出力
     time = 0
     convergence = 0.1
     call cg_iric_write_sol_time_f(time, ier)
     call cg_iric_write_sol_iface_real_f('FluxI', fluxi, ier)
     call cg_iric_write_sol_jface_real_f('FluxJ', fluxj, ier)
     do
       time = time + 10.0

       ! (ここで計算を実行)

       call iric_check_cancel_f(canceled)
       if (canceled == 1) exit

       ! 計算結果を出力
       call iric_write_sol_start_f(condFile, ier)
       call cg_iric_write_sol_time_f(time, ier)
       call cg_iric_write_sol_iface_real_f('FluxI', fluxi, ier)
       call cg_iric_write_sol_jface_real_f('FluxJ', fluxj, ier)
       call cg_iric_flush_f(condFile, fin, ier)
       call iric_write_sol_end_f(condFile, ier)

       if (time > 1000) exit
     end do

     ! CGNS ファイルのクローズ
     call cg_close_f(fin, ier)
     stop
   end program SampleProgram