var/da/da_physics/da_tp_to_qs_adj1.inc

   1 subroutine da_tp_to_qs_adj1( t, p, es, t_prime, p_prime, qs_prime_over_qs )
   2
   3    !---------------------------------------------------------------------------
   4    !  Purpose: Adjoint of da_tp_to_qs_lin.
   5    !---------------------------------------------------------------------------
   6
   7    implicit none
   8
   9    real, intent(in)    :: t                ! Temperature.
  10    real, intent(in)    :: p                ! Pressure.
  11    real, intent(in)    :: es               ! Sat. vapour pressure.
  12    real, intent(inout) :: t_prime          ! Temperature increment.
  13    real, intent(inout) :: p_prime          ! Pressure increment.
  14    real, intent(in)    :: qs_prime_over_qs ! qs~/qs.
  15
  16    real :: temp             ! Temporary storage.
  17    real :: es_prime_over_es ! es~/es
  18
  19    if (trace_use) call da_trace_entry("da_tp_to_qs_adj1")
  20
  21    !---------------------------------------------------------------------------
  22    ! [3.0] Calculate saturation specific humidity increment:
  23    !---------------------------------------------------------------------------
  24
  25    temp = qs_prime_over_qs / ( p - rd_over_rv1 * es )
  26
  27    es_prime_over_es = temp * p
  28
  29    p_prime = p_prime - temp
  30
  31    !---------------------------------------------------------------------------
  32    ! [2.0] Calculate saturation vapour pressure increment:
  33    !---------------------------------------------------------------------------
  34
  35    temp = t + es_gammakelvin
  36
  37    t_prime = t_prime + es_gammabeta * es_prime_over_es / ( temp * temp )
  38
  39    if (trace_use) call da_trace_exit("da_tp_to_qs_adj1")
  40
  41 end subroutine da_tp_to_qs_adj1
  42
  43