Subroutine to smooth a MS_FE_FILM field, the pressure for instance.
By default, the smoothing kernel is a 5x5 Gaussian filter
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(MS_FE_FILM), | intent(inout) | :: | ms_fe_f | |||
integer(kind=I4), | intent(in) | :: | code | saved information like P_N |
||
logical(kind=I4), | intent(in) | :: | nodal | if false : cell value, if true : nodal value |
subroutine smooth_ms_fe_f(ms_fe_f, code, nodal)
implicit none
type(MS_FE_FILM), intent(inout) :: ms_fe_f
integer(kind=I4), intent(in) :: code !! *saved information like P_N*
logical(kind=I4), intent(in) :: nodal !! *if false : cell value, if true : nodal value*
integer(kind=I4) :: nx, ny
real(kind=R8), allocatable, dimension(:,:) :: mat
call ms_fe_f_2_mat(ms_fe_f, code, nodal, mat)
nx = ubound(mat, 1) ; ny = ubound(mat, 2)
call smooth_mat(mat, nx, ny, s=5)
call mat_2_ms_fe_f(ms_fe_f, code, nodal, mat)
deallocate(mat)
return
endsubroutine smooth_ms_fe_f