write_global_field Subroutine

private subroutine write_global_field(global_field, state, filename, is_ok, message, xc, yc)

Write the assembled global field to filename in the standard ASCII format.

i (x) is fastest; a blank line separates each constant-j row so gnuplot/ParaView parse the output as a structured 2D grid. When xc/yc are present (curvilinear), those cell centroids are written instead of the uniform x_left+(i-0.5)*dx formula. The uniform branch (no optional args) is byte-identical to the prior implementation.

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: global_field(:,:,:)
type(solver_state_2d_t), intent(in) :: state
character(len=*), intent(in) :: filename
logical, intent(out) :: is_ok
character(len=*), intent(out) :: message
real(kind=wp), intent(in), optional :: xc(:,:)
real(kind=wp), intent(in), optional :: yc(:,:)