da_interp_lin_2d_new.inc

1 subroutine da_interp_lin_2d_new(fm2d, is, ie, js, je, &
2                          i, j, dx, dy, dxm, dym, &
3                          fo2d, num)
4 
5    !-----------------------------------------------------------------------
6    ! Purpose: TBD
7    !-----------------------------------------------------------------------
8 
9    implicit none
10 
11    integer,                intent(in)  :: num
12    integer,                intent(in)  :: is, ie, js, je
13    integer,                intent(in)  :: i(num), j(num)
14    real   ,                intent(in)  :: dx(num), dxm(num), dy(num), dym(num)
15    real   ,                intent(in)  :: fm2d(is:ie,js:je) ! Input variable
16    real   ,                intent(out) :: fo2d(num)  ! Output variable 
17    
18    integer :: point
19 
20    if (trace_use) call da_trace_entry("da_interp_lin_2d_new")
21 
22    do point =1,num
23       fo2d(point) = dym(point)*(dxm(point)*fm2d(i(point),j(point)) + dx(point)*fm2d(i(point)+1,j(point))) &
24          + dy(point) *(dxm(point)*fm2d(i(point),j(point)+1) + dx(point)*fm2d(i(point)+1,j(point)+1))
25    end do
26 
27    if (trace_use) call da_trace_exit("da_interp_lin_2d_new")
28 
29 end subroutine da_interp_lin_2d_new
30 
31