scalar_residual_sx.f90

Go to the documentation of this file.

 
module scalar_residual_sx
  use scalar_residual, only : scalar_residual_t
  use ax_product, only : ax_t
  use field, only : field_t
  use coefs, only : coef_t
  use space, only : space_t
  use mesh, only : mesh_t
  use num_types, only : rp
  use math, only : copy
  use field, only : field_t
  use mesh, only : mesh_t
  use ax_product, only : ax_t
  use space, only : space_t
  use coefs, only : coef_t
  implicit none
  private
 
  type, public, extends(scalar_residual_t) :: scalar_residual_sx_t
   contains
     procedure, nopass :: compute => scalar_residual_sx_compute
  end type scalar_residual_sx_t
 
contains
 
  subroutine scalar_residual_sx_compute(Ax, s, s_res, f_Xh, c_Xh, msh, Xh, &
       lambda, rho_cp, bd, dt, n)
    class(ax_t), intent(in) :: Ax
    type(mesh_t), intent(inout) :: msh
    type(space_t), intent(inout) :: Xh
    type(field_t), intent(inout) :: s
    type(field_t), intent(inout) :: s_res
    type(field_t), intent(in) :: f_Xh
    type(coef_t), intent(inout) :: c_Xh
    type(field_t), intent(in) :: lambda, rho_cp
    real(kind=rp), intent(in) :: bd
    real(kind=rp), intent(in) :: dt
    integer, intent(in) :: n
    integer :: i
 
    call copy(c_xh%h1, lambda%x, n)
    do i = 1, n
       c_xh%h2(i,1,1,1) = rho_cp%x(i,1,1,1) * bd / dt
    end do
    c_xh%ifh2 = .true.
 
    call ax%compute(s_res%x, s%x, c_xh, msh, xh)
 
    do i = 1, n
       s_res%x(i,1,1,1) = (-s_res%x(i,1,1,1)) + f_xh%x(i,1,1,1)
    end do
 
  end subroutine scalar_residual_sx_compute
 
end module scalar_residual_sx