Skip to content

Latest commit

 

History

History
109 lines (84 loc) · 3.84 KB

12_05_output_result_particle.rst

File metadata and controls

109 lines (84 loc) · 3.84 KB
Raw

Value defined at particles

When you output value defined at particles, please use the functions in :numref:`table_iriclib_output_particle_functions`.

:numref:`example_output_calc_result_particle` shows an example of the process to output value defined at grid nodes.

Note

These functions are deprecated

When you want to output value defined at particles, we recommend that you use functions described at :ref:`iriclib_output_result_particlegroup`. Using those new functions, you can output particles with multiple groups, and you can visualize particles with different settings for color and size for each group.

Subroutines to use for outputting result defined at particles
Subroutine Remarks
cg_iric_write_sol_particle_pos2d Outputs particle positions (two-dimensions)
cg_iric_write_sol_particle_pos3d Outputs particle positions (three-dimensions)
cg_iric_write_sol_particle_integer Outputs integer-type calculation results, giving a value for each particle.
cg_iric_write_sol_particle_real Outputs double-precision real-type calculation results, giving a value for each particle. 
program SampleProgram
  use iric
  implicit none

  integer:: fin, ier, isize, jsize
  integer:: canceled
  integer:: locked
  double precision:: time
  double precision, dimension(:,:), allocatable::grid_x, grid_y
  integer:: numparticles = 10
  double precision, dimension(:), allocatable:: particle_x, particle_y
  double precision, dimension(:), allocatable:: velocity_x, velocity_y, temperature
  character(len=20):: condFile

  condFile = 'test.cgn'

  ! Open CGNS file
  call cg_iric_open(condFile, IRIC_MODE_MODIFY, fin, ier)
  if (ier /=0) STOP "*** Open error of CGNS file ***"

  ! Check the grid size
  call cg_iric_read_grid2d_str_size(fin, isize, jsize, ier)
  ! Allocate memory for loading the grid
  allocate(grid_x(isize, jsize), grid_y(isize, jsize))
  ! Allocate memory for calculation result
  allocate(particle_x(numparticles), particle_y(numparticles))
  allocate(velocity_x(numparticles), velocity_y(numparticles), temperature(numparticles))

  ! Read the grid into memory
  call cg_iric_read_grid2d_coords(fin, grid_x, grid_y, ier)

  ! Output the initial state information.
  time = 0
  call cg_iric_write_sol_start(fin, ier)
  call cg_iric_write_sol_time(fin, time, ier)
  call cg_iric_write_sol_particle_pos2d(fin, numparticles, particle_x, particle_y, ier)
  call cg_iric_write_sol_particle_real(fin, 'VelocityX', velocity_x, ier)
  call cg_iric_write_sol_particle_real(fin, 'VelocityY', velocity_y, ier)
  call cg_iric_write_sol_particle_real(fin, 'Temperature', temperature, ier)
  call cg_iric_write_sol_end(fin, ier)
  do
    time = time + 10.0

    ! (Perform calculation here)

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

    ! Output calculation results
    call cg_iric_write_sol_start(fin, ier)
    call cg_iric_write_sol_time(fin, time, ier)
    call cg_iric_write_sol_particle_pos2d(fin, numparticles, particle_x, particle_y, ier)
    call cg_iric_write_sol_particle_real(fin, 'VelocityX', velocity_x, ier)
    call cg_iric_write_sol_particle_real(fin, 'VelocityY', velocity_y, ier)
    call cg_iric_write_sol_particle_real(fin, 'Temperature', temperature, ier)
    call cg_iric_write_sol_end(fin, ier)

    if (time > 1000) exit
  end do

  ! Close CGNS file
  call cg_iric_close(fin, ier)
  stop
end program SampleProgram