da_transform_xtoy_sound.inc

1 subroutine da_transform_xtoy_sound (grid, iv, y)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    type (domain),     intent(in)    :: grid
10    type (iv_type),    intent(in)    :: iv       ! Innovation vector (O-B).
11    type (y_type),     intent(inout) :: y        ! y = h (grid%xa) (linear)
12 
13    real, allocatable :: u(:,:)
14    real, allocatable :: v(:,:)
15    real, allocatable :: t(:,:)
16    real, allocatable :: q(:,:)
17 
18    integer :: n
19 
20    if (trace_use_dull) call da_trace_entry("da_transform_xtoy_sound")
21 
22    allocate (u(iv%info(sound)%max_lev,iv%info(sound)%n1:iv%info(sound)%n2))
23    allocate (v(iv%info(sound)%max_lev,iv%info(sound)%n1:iv%info(sound)%n2))
24    allocate (t(iv%info(sound)%max_lev,iv%info(sound)%n1:iv%info(sound)%n2))
25    allocate (q(iv%info(sound)%max_lev,iv%info(sound)%n1:iv%info(sound)%n2))
26   
27    call da_interp_lin_3d (grid%xa%u, iv%info(sound), u)
28    call da_interp_lin_3d (grid%xa%v, iv%info(sound), v)
29    call da_interp_lin_3d (grid%xa%t, iv%info(sound), t)
30    call da_interp_lin_3d (grid%xa%q, iv%info(sound), q)
31 
32    do n=iv%info(sound)%n1,iv%info(sound)%n2
33       y%sound(n)%u(:) = u(1:size(y%sound(n)%u),n)
34       y%sound(n)%v(:) = v(1:size(y%sound(n)%v),n)
35       y%sound(n)%t(:) = t(1:size(y%sound(n)%t),n)
36       y%sound(n)%q(:) = q(1:size(y%sound(n)%q),n)
37    end do
38 
39    deallocate (u)
40    deallocate (v)
41    deallocate (t)
42    deallocate (q)
43 
44    if (trace_use_dull) call da_trace_exit("da_transform_xtoy_sound")
45 
46 end subroutine da_transform_xtoy_sound
47 
48