1 subroutine da_transform_xtoy_buoy_adj(grid, iv, jo_grad_y, jo_grad_x)
3 !--------------------------------------------------------------------------
5 ! Updated for Analysis on Arakawa-C grid
6 ! Author: Syed RH Rizvi, MMM/ESSL/NCAR, Date: 10/22/2008
7 !--------------------------------------------------------------------------
11 type (domain), intent(in) :: grid
12 type (iv_type), intent(in) :: iv ! obs. inc vector (o-b).
13 type (y_type) , intent(inout) :: jo_grad_y ! grad_y(jo)
14 type (x_type) , intent(inout) :: jo_grad_x ! grad_x(jo)
16 integer :: n ! Loop counter.
18 real, allocatable :: model_u(:,:)
19 real, allocatable :: model_v(:,:)
20 real, allocatable :: model_t(:,:)
21 real, allocatable :: model_q(:,:)
22 real, allocatable :: model_psfc(:)
23 real, allocatable :: ub(:,:)
24 real, allocatable :: vb(:,:)
26 if (trace_use_dull) call da_trace_entry("da_transform_xtoy_buoy_adj")
28 if (sfc_assi_options == sfc_assi_options_1) then
29 allocate (model_u(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
30 allocate (model_v(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
31 allocate (model_t(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
32 allocate (model_q(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
33 allocate (model_psfc(iv%info(buoy)%n1:iv%info(buoy)%n2))
35 allocate (ub(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
36 allocate (vb(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
38 call da_interp_lin_3d (grid%xb%u, iv%info(buoy), ub)
39 call da_interp_lin_3d (grid%xb%v, iv%info(buoy), vb)
41 ! [1.2] Interpolate horizontally:
42 do n=iv%info(buoy)%n1,iv%info(buoy)%n2
44 call da_uv_to_sd_adj(jo_grad_y%buoy(n)%u, &
45 jo_grad_y%buoy(n)%v, model_u(1,n), model_v(1,n), ub(1,n), vb(1,n))
47 model_u(1,n) = jo_grad_y%buoy(n)%u
48 model_v(1,n) = jo_grad_y%buoy(n)%v
51 model_t(1,n) = jo_grad_y%buoy(n)%t
52 model_q(1,n) = jo_grad_y%buoy(n)%q
53 model_psfc(n) = jo_grad_y%buoy(n)%p
56 call da_interp_lin_3d_adj (jo_grad_x%u, iv%info(buoy), model_u,'u')
57 call da_interp_lin_3d_adj (jo_grad_x%v, iv%info(buoy), model_v,'v')
59 call da_interp_lin_3d_adj (jo_grad_x%u, iv%info(buoy), model_u)
60 call da_interp_lin_3d_adj (jo_grad_x%v, iv%info(buoy), model_v)
62 call da_interp_lin_3d_adj (jo_grad_x%t, iv%info(buoy), model_t)
63 call da_interp_lin_3d_adj (jo_grad_x%q, iv%info(buoy), model_q)
65 call da_interp_lin_2d_adj (jo_grad_x%psfc, iv%info(buoy), 1, model_psfc)
70 deallocate (model_psfc)
74 else if (sfc_assi_options == sfc_assi_options_2) then
75 call da_transform_xtopsfc_adj(grid,iv, buoy, iv%buoy(:), jo_grad_y%buoy(:), jo_grad_x)
78 if (trace_use_dull) call da_trace_exit("da_transform_xtoy_buoy_adj")
80 end subroutine da_transform_xtoy_buoy_adj