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