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)
64      do j=1,ny
77            else if (j == ny) then
105   j = ny/4+1