1*55a74a43SLisandro Dalcin! -*- f90 -*- 2*55a74a43SLisandro Dalcin 3*55a74a43SLisandro Dalcinpython module App 4*55a74a43SLisandro Dalcininterface 5*55a74a43SLisandro Dalcin 6*55a74a43SLisandro Dalcin subroutine formFunction(nx, ny, nz, h, t, x, xdot, f) 7*55a74a43SLisandro Dalcin intent(c) formFunction 8*55a74a43SLisandro Dalcin integer, intent(in), intent(hide) :: nx = shape(x,0) 9*55a74a43SLisandro Dalcin integer, intent(in), intent(hide) :: ny = shape(x,1) 10*55a74a43SLisandro Dalcin integer, intent(in), intent(hide) :: nz = shape(x,2) 11*55a74a43SLisandro Dalcin real(kind=8), intent(in) :: h(3), t 12*55a74a43SLisandro Dalcin real(kind=8), intent(in) :: x(nx,ny,nz), xdot(nx,ny,nz) 13*55a74a43SLisandro Dalcin real(kind=8), intent(inout) :: f(nx,ny,nz) 14*55a74a43SLisandro Dalcin end subroutine formFunction 15*55a74a43SLisandro Dalcin 16*55a74a43SLisandro Dalcin subroutine formInitial(nx, ny, nz, h, t, x) 17*55a74a43SLisandro Dalcin intent(c) formInitial 18*55a74a43SLisandro Dalcin integer, intent(in), intent(hide) :: nx = shape(x,0) 19*55a74a43SLisandro Dalcin integer, intent(in), intent(hide) :: ny = shape(x,1) 20*55a74a43SLisandro Dalcin integer, intent(in), intent(hide) :: nz = shape(x,2) 21*55a74a43SLisandro Dalcin real(kind=8), intent(in) :: h(3), t 22*55a74a43SLisandro Dalcin real(kind=8), intent(inout) :: x(nx,ny,nz) 23*55a74a43SLisandro Dalcin end subroutine formFunction 24*55a74a43SLisandro Dalcin 25*55a74a43SLisandro Dalcinend python module App 26