var/da/da_physics/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