5   integer(kind=C_INT), intent(in)    :: nx, ny, nz  local
7   real(kind=C_DOUBLE), intent(in)    :: x(nx,ny,nz), xdot(nx,ny,nz)
8   real(kind=C_DOUBLE), intent(inout) :: f(nx,ny,nz)
9   call formfunction_f(nx, ny, nz, h, t, x, xdot, f)
16   integer(kind=C_INT), intent(in)    :: nx, ny, nz  local
18   real(kind=C_DOUBLE), intent(inout) :: x(nx,ny,nz)
19   call forminitial_f(nx, ny, nz, h, t, x)
45   integer, intent(in)         :: nx, ny, nz  local
47   real(kind=8), intent(inout) :: x(nx,ny,nz)
54   integer, intent(in)         :: nx, ny, nz  local
56   real(kind=8), intent(in)    :: x(nx,ny,nz), xdot(nx,ny,nz)
57   real(kind=8), intent(inout) :: f(nx,ny,nz)
65         do i=1,nx
70            else if (i == nx) then
104   i = nx/4+1