!c Description: 
!c   Èùʬ±é»»
!c
!c Current Code Owner: 
!c   sugiyama@gfd-dennou.org
!c
!c Histry: 
!c   Version    Date          Comment
!c   -------    ----------    --------
!c   1.0        2004-08-12    ¿ù»³¹Ì°ìϯ 
!c
!c Copyright (C) SUGIYAMA Ko-ichiro, 2004, All rights reserved

function u_d2x_u(var)
  use gridset
  use if_boundary

  implicit none

  real(8), intent(in) :: var(imin:imax, kmin:kmax)
  real(8)             :: u_d2x_u(imin:imax, kmin:kmax)
  
  u_d2x_u(imin+1:imax-1,kmin:kmax) =           &
       &  (                                       &
       &             var(imin+2:imax,kmin:kmax)   &
       &   - 2.0d0 * var(imin+1:imax-1,kmin:kmax) &
       &   +         var(imin:imax-2,kmin:kmax)   &
       &  ) / dx
  
  call boundary(bdc_u, u_d2x_u)    
end function u_d2x_u


function u_d2z_u(var) 
  use gridset
  use if_boundary

  implicit none

  real(8), intent(in) :: var(imin:imax, kmin:kmax)
  real(8)             :: u_d2z_u(imin:imax, kmin:kmax)
  
  u_d2z_u(imin:imax,kmin+1:kmax-1) =           & 
       &  (                                       &
       &             var(imin:imax,kmin+2:kmax)   &
       &   - 2.0d0 * var(imin:imax,kmin+1:kmax-1) &
       &   +         var(imin:imax,kmin:kmax-2)   &
       &  ) / dz
  
  call boundary(bdc_u, u_d2z_u)    
end function u_d2z_u


function w_d2x_w(var) 
  use gridset
  use if_boundary

  implicit none

  real(8), intent(in) :: var(imin:imax, kmin:kmax)
  real(8)             :: w_d2x_w(imin:imax, kmin:kmax)
  
  w_d2x_w(imin+1:imax-1,kmin:kmax) =           &
       &  (                                       &
       &             var(imin+2:imax,kmin:kmax)   &
       &   - 2.0d0 * var(imin+1:imax-1,kmin:kmax) &
       &   +         var(imin:imax-2,kmin:kmax)   &
       &  ) / dx
  
  call boundary(bdc_w, w_d2x_w)    
end function w_d2x_w


function w_d2z_w(var) 
  use gridset
  use if_boundary

  implicit none

  real(8), intent(in) :: var(imin:imax, kmin:kmax)
  real(8)             :: w_d2z_w(imin:imax, kmin:kmax)
    
  w_d2z_w(imin:imax,kmin+1:kmax-1) =           & 
       &  (                                       &
       &             var(imin:imax,kmin+2:kmax)   &
       &   - 2.0d0 * var(imin:imax,kmin+1:kmax-1) &
       &   +         var(imin:imax,kmin:kmax-2)   &
       &  ) / dz
  
  call boundary(bdc_w, w_d2z_w)    
end function w_d2z_w