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)
63   do k=1,nz
84            else if (k == nz) then
106   k = nz/2+1
107   f(i,j,k:nz) = f(i,j,k:nz) + 300.0