da_balance_geoterm_adj.inc
References to this file elsewhere.
1 subroutine da_balance_geoterm_adj( cori, rho, term_x, term_y, u, v)
2
3 !---------------------------------------------------------------------------
4 ! Purpose: Adjoint of da_balance_geoterm.
5 !---------------------------------------------------------------------------
6
7 implicit none
8
9 real, dimension(ims:ime,jms:jme),intent(in):: cori, & ! Coriolis factor.
10 rho, & ! Density
11 term_x,& ! x component of term.
12 term_y ! y component of term.
13
14 real, dimension(ims:ime,jms:jme), intent(inout) :: u, & ! u wind increment
15 v ! v wind increment
16
17 !---------------------------------------------------------------------------
18 ! [2.0] Calculate term_y = f rho u~:
19 !---------------------------------------------------------------------------
20
21 u(its:ite,jts:jte) = u(its:ite,jts:jte) + &
22 rho(its:ite,jts:jte) * cori(its:ite,jts:jte) * &
23 term_y(its:ite,jts:jte)
24
25 !---------------------------------------------------------------------------
26 ! [1.0] Calculate term_x = -f rho v~:
27 !---------------------------------------------------------------------------
28
29 v(its:ite,jts:jte) = v(its:ite,jts:jte) - &
30 rho(its:ite,jts:jte) * cori(its:ite,jts:jte) * &
31 term_x(its:ite,jts:jte)
32
33 end subroutine da_balance_geoterm_adj
34
35