weno_face_conc_y Function

private pure function weno_face_conc_y(nx, ny, nz, tr, wet_T, ii, cc, kk, d, avail_up, avail_down, rung_max, cfl) result(conc)

Meridional analogue of weno_face_conc_x. cc is the donor cell (j-index); the stencil steps along j with d = +1 (u>0) / -1 (u<0).

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nx
integer, intent(in) :: ny
integer, intent(in) :: nz
real(kind=wp), intent(in) :: tr(nx,ny,nz)
real(kind=wp), intent(in) :: wet_T(nx,ny)
integer, intent(in) :: ii
integer, intent(in) :: cc
integer, intent(in) :: kk
integer, intent(in) :: d
integer, intent(in) :: avail_up
integer, intent(in) :: avail_down
integer, intent(in) :: rung_max
real(kind=wp), intent(in) :: cfl

Return Value real(kind=wp)


Calls

proc~~weno_face_conc_y~~CallsGraph proc~weno_face_conc_y weno_face_conc_y proc~plm_face_swept plm_face_swept proc~weno_face_conc_y->proc~plm_face_swept proc~ppm_mirror_h ppm_mirror_h proc~weno_face_conc_y->proc~ppm_mirror_h proc~recon_rung_for_face recon_rung_for_face proc~weno_face_conc_y->proc~recon_rung_for_face proc~weno5_face_swept weno5_face_swept proc~weno_face_conc_y->proc~weno5_face_swept proc~weno7_face_swept weno7_face_swept proc~weno_face_conc_y->proc~weno7_face_swept proc~weno9_face_swept weno9_face_swept proc~weno_face_conc_y->proc~weno9_face_swept

Called by

proc~~weno_face_conc_y~~CalledByGraph proc~weno_face_conc_y weno_face_conc_y proc~drain_swept_flux_y_weno drain_swept_flux_y_weno proc~drain_swept_flux_y_weno->proc~weno_face_conc_y proc~continuity_tracer_drain continuity_tracer_drain proc~continuity_tracer_drain->proc~drain_swept_flux_y_weno proc~ocean_dyn_flush_tracer_window ocean_dyn_flush_tracer_window proc~ocean_dyn_flush_tracer_window->proc~continuity_tracer_drain proc~ocean_dyn_step ocean_dyn_step proc~ocean_dyn_step->proc~continuity_tracer_drain proc~ocean_dyn_step_split ocean_dyn_step_split proc~ocean_dyn_step_split->proc~continuity_tracer_drain proc~driver_run_ocean driver_run_ocean proc~driver_run_ocean->proc~ocean_dyn_flush_tracer_window proc~engine_step engine_step proc~driver_run_ocean->proc~engine_step proc~engine_step->proc~ocean_dyn_step proc~engine_step->proc~ocean_dyn_step_split proc~rdb_ocean_set_tracer rdb_ocean_set_tracer proc~rdb_ocean_set_tracer->proc~ocean_dyn_flush_tracer_window proc~driver_run driver_run proc~driver_run->proc~driver_run_ocean proc~rdb_ocean_step rdb_ocean_step proc~rdb_ocean_step->proc~engine_step

Variables

Type Visibility Attributes Name Initial
integer, private :: jm
integer, private :: jp
real(kind=wp), private :: q0
real(kind=wp), private :: qm1
real(kind=wp), private :: qm2
real(kind=wp), private :: qm3
real(kind=wp), private :: qm4
real(kind=wp), private :: qp1
real(kind=wp), private :: qp2
real(kind=wp), private :: qp3
real(kind=wp), private :: qp4
integer, private :: rung

Source Code

   pure function weno_face_conc_y(nx, ny, nz, tr, wet_T, ii, cc, kk, d, &
                                  avail_up, avail_down, rung_max, cfl) result(conc)
      !! Meridional analogue of weno_face_conc_x.  `cc` is the donor cell
      !! (j-index); the stencil steps along j with `d = +1` (u>0) / `-1` (u<0).
      !$acc routine seq
      integer, intent(in) :: nx, ny, nz, ii, cc, kk, d, avail_up, avail_down, rung_max
      real(wp), intent(in) :: tr(nx, ny, nz), wet_T(nx, ny)
      real(wp), intent(in) :: cfl
      real(wp) :: conc
      integer :: rung, jp, jm
      real(wp) :: q0, qp1, qp2, qp3, qp4, qm1, qm2, qm3, qm4
      rung = recon_rung_for_face(avail_up, avail_down, rung_max)
      q0 = tr(ii, cc, kk)
      if (rung <= 0) then
         conc = q0
         return
      end if
      jp = min(max(cc + d, 1), ny)
      jm = min(max(cc - d, 1), ny)
      qp1 = ppm_mirror_h(tr(ii, jp, kk), q0, wet_T(ii, jp))
      qm1 = ppm_mirror_h(tr(ii, jm, kk), q0, wet_T(ii, jm))
      if (rung == 1) then
         conc = plm_face_swept(qm1, q0, qp1, cfl)
         return
      end if
      jp = min(max(cc + 2*d, 1), ny)
      jm = min(max(cc - 2*d, 1), ny)
      qp2 = ppm_mirror_h(tr(ii, jp, kk), qp1, wet_T(ii, jp))
      qm2 = ppm_mirror_h(tr(ii, jm, kk), qm1, wet_T(ii, jm))
      if (rung == 2) then
         conc = weno5_face_swept(qm2, qm1, q0, qp1, qp2, cfl)
         return
      end if
      jp = min(max(cc + 3*d, 1), ny)
      jm = min(max(cc - 3*d, 1), ny)
      qp3 = ppm_mirror_h(tr(ii, jp, kk), qp2, wet_T(ii, jp))
      qm3 = ppm_mirror_h(tr(ii, jm, kk), qm2, wet_T(ii, jm))
      if (rung == 3) then
         conc = weno7_face_swept(qm3, qm2, qm1, q0, qp1, qp2, qp3, cfl)
         return
      end if
      jp = min(max(cc + 4*d, 1), ny)
      jm = min(max(cc - 4*d, 1), ny)
      qp4 = ppm_mirror_h(tr(ii, jp, kk), qp3, wet_T(ii, jp))
      qm4 = ppm_mirror_h(tr(ii, jm, kk), qm3, wet_T(ii, jm))
      conc = weno9_face_swept(qm4, qm3, qm2, qm1, q0, qp1, qp2, qp3, qp4, cfl)
   end function weno_face_conc_y