1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
#include "cctk.h"
#include "cctk_Functions.h"
#include "cctk_Parameters.h"
subroutine deriv2_mixed ( cctkGH, dir1, dir2, var, ni, nj, nk, &
dvar2, table_handle )
implicit none
DECLARE_CCTK_FUNCTIONS
DECLARE_CCTK_PARAMETERS
CCTK_POINTER, intent(IN) :: cctkGH
CCTK_INT, intent(IN) :: dir1, dir2
CCTK_REAL, dimension(ni,nj,nk), intent(IN) :: var
CCTK_INT, intent(IN) :: ni, nj, nk
CCTK_REAL, dimension(ni,nj,nk), intent(OUT) :: dvar2
CCTK_INT, intent(IN) :: table_handle
CCTK_REAL, dimension(:,:,:), allocatable :: tmp
allocate ( tmp(ni,nj,nk) )
call Diff_gv ( cctkGH, dir1, var, tmp, table_handle)
call Diff_gv ( cctkGH, dir2, tmp, dvar2, table_handle )
deallocate ( tmp )
end subroutine deriv2_mixed
|