FIXME : Add documentation
subroutine allocate_arrays
use theta_grid, only: ntgrid
use kt_grids, only: ntheta0, naky
implicit none
if (.not. allocated(D_res)) then
allocate (D_res(ntheta0, naky))
endif
if (.not. allocated(hypervisc_filter)) then
allocate (hypervisc_filter(-ntgrid:ntgrid,ntheta0,naky)) ; hypervisc_filter = 1.0
end if
D_res = 0.
end subroutine allocate_arrays