write_netcdf_helmholtz Subroutine

public subroutine write_netcdf_helmholtz(fgid, mesh, bnd_type_core, bnd_type_wall, bnd_type_dome, bnd_type_out, co, lambda, xi, rhs, guess, sol, hcsr_write_on)

Writes a Helmholtz problem to NETCDF file

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: fgid

NETCDF file or group id

type(mesh_cart_t), intent(in) :: mesh

Mesh

integer, intent(in) :: bnd_type_core

Boundary descriptor for core boundary

integer, intent(in) :: bnd_type_wall

Boundary descriptor for wall boundary

integer, intent(in) :: bnd_type_dome

Boundary descriptor for dome boundary

integer, intent(in) :: bnd_type_out

Boundary descriptor for outer(mask) boundary

real(kind=FP), intent(in), dimension(mesh%get_n_points()) :: co

Coefficient within Helmholtz operator

real(kind=FP), intent(in), dimension(mesh%get_n_points_inner()) :: lambda

Coefficient within Helmholtz operator

real(kind=FP), intent(in), dimension(mesh%get_n_points_inner()) :: xi

Coefficient within Helmholtz operator

real(kind=FP), intent(in), dimension(mesh%get_n_points()) :: rhs

Right hand side

real(kind=FP), intent(in), dimension(mesh%get_n_points()) :: guess

Initial guess

real(kind=FP), intent(in), dimension(mesh%get_n_points()) :: sol

Solution

logical, intent(in) :: hcsr_write_on

Switch to write also Helmoltz matrix in CSR format Useful for debugging, benchmarking wwith external tools