var/da/da_vtox_transforms/da_transform_vpatox_adj.inc

   1 subroutine da_transform_vpatox_adj(grid, nens, ep, vp, nobwin)
   2
   3    !----------------------------------------------------------------------
   4    ! Purpose:  Does adjoint transform of alpha control variable (vp%alpha) to x
   5    !----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type(ep_type),              intent(in)    :: ep     ! Ensemble perturbation structure.
  10    type(vp_type),              intent(inout) :: vp     ! Grdipt/level CV.
  11    type(domain),               intent(inout) :: grid
  12    integer, intent(in)                       :: nens   ! number of ensembles
  13    integer, intent(in), optional             :: nobwin
  14
  15    integer :: iobwin
  16
  17    if (nens <= 0 .or. alphacv_method /= alphacv_method_xa) return
  18
  19    if (trace_use) call da_trace_entry("da_transform_vpatox_adj")
  20
  21    iobwin = 1
  22    if ( present(nobwin) ) then
  23       iobwin = nobwin
  24    end if
  25
  26    if (nens > 0 .and. alphacv_method == alphacv_method_xa) then
  27       !this is for getting vp%alpha only
  28       vp%alpha = 0.0
  29       if ( anal_type_hybrid_dual_res ) then
  30          call da_calc_flow_dependence_xa_adj_dual_res(nens, ep, grid, grid%xa, vp, iobwin)
  31       else
  32          call da_calc_flow_dependence_xa_adj(nens, ep, grid%xa, vp, iobwin)
  33       end if
  34    end if
  35
  36    if (trace_use) call da_trace_exit("da_transform_vpatox_adj")
  37
  38 end subroutine da_transform_vpatox_adj
  39
  40