<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_TV_PROFILE_ADJ'><A href='../../html_code/physics/da_tv_profile_adj.inc.html#DA_TV_PROFILE_ADJ' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_tv_profile_adj(grid,jo_grad_x,info, n,pre_ma,tv_ma,ADJ_pre_ma, ADJ_tv_ma) 1,2
!-----------------------------------------------------------------------
! Purpose: adjoint routine for da_tv_profile
!-----------------------------------------------------------------------
implicit none
type (x_type), intent(inout) :: jo_grad_x ! grad_x(jo)
type (infa_type), intent(in) :: info
integer, intent(in) :: n
type (domain), intent(in) :: grid
real, intent(in) :: pre_ma(kts-1:kte+1)
real, intent(in) :: tv_ma(kts-1:kte+1)
real, intent(inout) :: ADJ_pre_ma(kts-1:kte+1)
real, intent(inout) :: ADJ_tv_ma(kts-1:kte+1)
integer :: ii,jj
real :: ADJ_tv_m(2,2,kts:kte)
integer :: i, j ! OBS location
real :: dx, dxm ! interpolation weights.
real :: dy, dym ! interpolation weights.
if (trace_use_dull) call da_trace_entry
("da_tv_profile_adj")
i = info%i(1,n)
j = info%j(1,n)
dx = info%dx(1,n)
dy = info%dy(1,n)
dxm = info%dxm(1,n)
dym = info%dym(1,n)
ADJ_tv_m(1,1,kts:kte) = dym*dxm * ADJ_tv_ma (kts:kte)
ADJ_tv_m(2,1,kts:kte) = dym*dx * ADJ_tv_ma (kts:kte)
ADJ_tv_m(1,2,kts:kte) = dy*dxm* ADJ_tv_ma (kts:kte)
ADJ_tv_m(2,2,kts:kte) = dy*dx* ADJ_tv_ma (kts:kte)
jo_grad_x%p(i,j,kts:kte) = jo_grad_x%p(i,j,kts:kte) + dym*dxm * ADJ_pre_ma(kts:kte)
jo_grad_x%p(i+1,j,kts:kte) = jo_grad_x%p(i+1,j,kts:kte) + dym*dx * ADJ_pre_ma(kts:kte)
jo_grad_x%p(i,j+1,kts:kte) = jo_grad_x%p(i,j+1,kts:kte) + dy*dxm * ADJ_pre_ma(kts:kte)
jo_grad_x%p(i+1,j+1,kts:kte) = jo_grad_x%p(i+1,j+1,kts:kte) + dy*dx* ADJ_pre_ma(kts:kte)
ADJ_tv_ma (kts:kte) = 0.0
ADJ_pre_ma(kts:kte) = 0.0
do ii=i,i+1
do jj=j,j+1
jo_grad_x%t(ii,jj,kts:kte) = jo_grad_x%t(ii,jj,kts:kte) + &
ADJ_tv_m(ii-i+1,jj-j+1,kts:kte)*(1.0+0.61*grid%xb%q(ii,jj,kts:kte))
jo_grad_x%q(ii,jj,kts:kte) = jo_grad_x%q(ii,jj,kts:kte) + &
0.61*grid%xb%t(ii,jj,kts:kte)*ADJ_tv_m(ii-i+1,jj-j+1,kts:kte)
end do
end do
if (trace_use_dull) call da_trace_exit("da_tv_profile_adj")
end subroutine da_tv_profile_adj