da_prho_to_t_lin.inc

1 subroutine da_prho_to_t_lin(grid) 
2 
3    !---------------------------------------------------------------------------
4    ! Purpose: Calculates temperature increments from pressure/density increments
5    !
6    ! Method:  Linearitsed ideal gas law: T~ = T (p"/p - rho"/rho)
7    !
8    ! Assumptions: 1) Model level stored top down.
9    !---------------------------------------------------------------------------
10 
11    implicit none
12    
13    type (domain), intent(inout)    :: grid
14 
15    if (trace_use) call da_trace_entry("da_prho_to_t_lin")
16 
17    grid%xa % t(its:ite,jts:jte,kts:kte) = grid%xb % t(its:ite,jts:jte,kts:kte) * &
18                             (grid%xa % p(its:ite,jts:jte,kts:kte) / &
19                               grid%xb % p(its:ite,jts:jte,kts:kte) - &
20                               grid%xa % rho(its:ite,jts:jte,kts:kte) / &
21                               grid%xb % rho(its:ite,jts:jte,kts:kte)) 
22 
23    if (trace_use) call da_trace_exit("da_prho_to_t_lin")    
24 
25 end subroutine da_prho_to_t_lin
26 
27