<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_THICKNESS_TL'><A href='../../html_code/physics/da_thickness_tl.inc.html#DA_THICKNESS_TL' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_thickness_tl(pre_ma,tv_ma,ks,ke,tv1,tv2,layer1,layer2,pre1,pre2, & 1,2
TGL_pre_ma,TGL_tv_ma,TGL_tv1,TGL_tv2,TGL_thk)
!-----------------------------------------------------------------------
! Purpose: tangent-linear routine for thickness
!-----------------------------------------------------------------------
implicit none
integer, intent(in) :: layer1,layer2
integer, intent(in) :: ks,ke
real, intent(in) :: pre_ma(ks-1:ke+1)
real, intent(in) :: tv_ma(ks-1:ke+1)
real, intent(in) :: TGL_pre_ma(ks-1:ke+1)
real, intent(in) :: TGL_tv_ma(ks-1:ke+1)
real, intent(in) :: tv1,tv2
real, intent(in) :: TGL_tv1,TGL_tv2
real, intent(in) :: pre1,pre2
real, intent(out) :: TGL_thk
integer :: k
real :: p_tmp(ks-1:ke+1)
real :: TGL_p_tmp(ks-1:ke+1)
if (trace_use_dull) call da_trace_entry
("da_thickness_tl")
TGL_thk = TGL_tv1*log(pre_ma(layer1-1)/pre1) + &
TGL_pre_ma(layer1-1)*tv1/pre_ma(layer1-1) + &
TGL_tv2*log(pre2/pre_ma(layer2)) - &
TGL_pre_ma(layer2)*tv2/pre_ma(layer2)
TGL_thk = 0.5 * gas_constant/gravity * TGL_thk
TGL_p_tmp(layer1) = 0.0
p_tmp(layer1) = pre1
TGL_p_tmp(layer2-1) = 0.0
p_tmp(layer2-1) = pre2
do k=layer2,layer1-1
TGL_p_tmp(k) = TGL_pre_ma(k)
p_tmp(k) = pre_ma(k)
end do
! Vertical integration of the virtual temperature
do k=layer2,layer1-1
TGL_thk = TGL_thk + 0.5 * gas_constant/gravity * ( &
TGL_tv_ma(k)*log(p_tmp(k-1)/p_tmp(k+1)) + &
TGL_p_tmp(k-1)*tv_ma(k)/p_tmp(k-1) - &
TGL_p_tmp(k+1)*tv_ma(k)/p_tmp(k+1) )
end do
if (trace_use_dull) call da_trace_exit("da_thickness_tl")
end subroutine da_thickness_tl