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(:)
24 if (trace_use_dull) call da_trace_entry("da_transform_xtoy_buoy_adj")
26 if (sfc_assi_options == sfc_assi_options_1) then
27 allocate (model_u(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
28 allocate (model_v(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
29 allocate (model_t(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
30 allocate (model_q(1,iv%info(buoy)%n1:iv%info(buoy)%n2))
31 allocate (model_psfc(iv%info(buoy)%n1:iv%info(buoy)%n2))
32 ! [1.2] Interpolate horizontally:
33 do n=iv%info(buoy)%n1,iv%info(buoy)%n2
34 model_u(1,n) = jo_grad_y%buoy(n)%u
35 model_v(1,n) = jo_grad_y%buoy(n)%v
36 model_t(1,n) = jo_grad_y%buoy(n)%t
37 model_q(1,n) = jo_grad_y%buoy(n)%q
38 model_psfc(n) = jo_grad_y%buoy(n)%p
41 call da_interp_lin_3d_adj (jo_grad_x%u, iv%info(buoy), model_u,'u')
42 call da_interp_lin_3d_adj (jo_grad_x%v, iv%info(buoy), model_v,'v')
44 call da_interp_lin_3d_adj (jo_grad_x%u, iv%info(buoy), model_u)
45 call da_interp_lin_3d_adj (jo_grad_x%v, iv%info(buoy), model_v)
47 call da_interp_lin_3d_adj (jo_grad_x%t, iv%info(buoy), model_t)
48 call da_interp_lin_3d_adj (jo_grad_x%q, iv%info(buoy), model_q)
50 call da_interp_lin_2d_adj (jo_grad_x%psfc, iv%info(buoy), 1, model_psfc)
55 deallocate (model_psfc)
56 else if (sfc_assi_options == sfc_assi_options_2) then
57 call da_transform_xtopsfc_adj(grid,iv, buoy, iv%buoy(:), jo_grad_y%buoy(:), jo_grad_x)
60 if (trace_use_dull) call da_trace_exit("da_transform_xtoy_buoy_adj")
62 end subroutine da_transform_xtoy_buoy_adj