<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_TV_PROFILE_TL'><A href='../../html_code/physics/da_tv_profile_tl.inc.html#DA_TV_PROFILE_TL' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_tv_profile_tl(grid,i,j,dx,dxm,dy,dym,pre_ma,tv_ma,TGL_pre_ma,TGL_tv_ma) 1,2
!--------------------------------------------------------------------------
! Purpose: tangent-linear routine for da_tv_profile
!--------------------------------------------------------------------------
implicit none
type (domain), intent(in) :: grid
integer, intent(in) :: i, j ! OBS location
real, intent(in) :: dx, dxm ! interpolation weights.
real, intent(in) :: dy, dym ! interpolation weights.
real, intent(out) :: TGL_pre_ma(kts-1:kte+1)
real, intent(out) :: TGL_tv_ma(kts-1:kte+1)
real, intent(out) :: pre_ma(kts-1:kte+1)
real, intent(out) :: tv_ma(kts-1:kte+1)
integer :: ii,jj
real :: tv_m(2,2,kts:kte)
real :: TGL_tv_m(2,2,kts:kte)
if (trace_use_dull) call da_trace_entry
("da_tv_profile_tl")
do ii=i,i+1
do jj=j,j+1
TGL_tv_m(ii-i+1,jj-j+1,kts:kte) = grid%xa%t(ii,jj,kts:kte)*(1.0+0.61*grid%xb%q(ii,jj,kts:kte)) + &
0.61*grid%xb%t(ii,jj,kts:kte)*grid%xa%q(ii,jj,kts:kte)
tv_m(ii-i+1,jj-j+1,kts:kte) = grid%xb%t(ii,jj,kts:kte)*(1.0+0.61*grid%xb%q(ii,jj,kts:kte))
end do
end do
TGL_pre_ma(kts:kte) = dym* ( dxm * grid%xa%p(i,j,kts:kte) + dx * grid%xa%p(i+1,j,kts:kte) ) + &
dy * ( dxm * grid%xa%p(i,j+1,kts:kte) + dx * grid%xa%p(i+1,j+1,kts:kte) )
TGL_tv_ma (kts:kte) = dym* ( dxm * TGL_tv_m(1,1,kts:kte) + dx * TGL_tv_m(2,1,kts:kte) ) + &
dy * ( dxm * TGL_tv_m(1,2,kts:kte) + dx * TGL_tv_m(2,2,kts:kte) )
pre_ma(kts:kte) = dym* ( dxm * grid%xb%p(i,j,kts:kte) + dx * grid%xb%p(i+1,j,kts:kte) ) + &
dy * ( dxm * grid%xb%p(i,j+1,kts:kte) + dx * grid%xb%p(i+1,j+1,kts:kte) )
tv_ma (kts:kte) = dym* ( dxm * tv_m (1,1,kts:kte) + dx * tv_m (2,1,kts:kte) ) + &
dy * ( dxm * tv_m (1,2,kts:kte) + dx * tv_m (2,2,kts:kte) )
if (trace_use_dull) call da_trace_exit("da_tv_profile_tl")
end subroutine da_tv_profile_tl