mom_coriolisadv Module Reference

Detailed Description

Accelerations due to the Coriolis force and momentum advection.

This file contains the subroutine that calculates the time derivatives of the velocities due to Coriolis acceleration and momentum advection. This subroutine uses either a vorticity advection scheme from Arakawa and Hsu, Mon. Wea. Rev. 1990, or Sadourny's (JAS 1975) energy conserving scheme. Both have been modified to use general orthogonal coordinates as described in Arakawa and Lamb, Mon. Wea. Rev. 1981. Both schemes are second order accurate, and allow for vanishingly small layer thicknesses. The Arakawa and Hsu scheme globally conserves both total energy and potential enstrophy in the limit of nondivergent flow. Sadourny's energy conserving scheme conserves energy if the flow is nondivergent or centered difference thickness fluxes are used.

Two sets of boundary conditions have been coded in the definition of relative vorticity. These are written as: NOSLIP defined (in spherical coordinates): relvort = dv/dx (east & west), with v = 0. relvort = -sec(Q) * d(u cos(Q))/dy (north & south), with u = 0.

NOSLIP not defined (free slip): relvort = 0 (all boundaries)

with Q temporarily defined as latitude. The free slip boundary condition is much more natural on a C-grid.

A small fragment of the grid is shown below:

    j+1  x ^ x ^ x   At x:  q, CoriolisBu
    j+1  > o > o >   At ^:  v, CAv, vh
    j    x ^ x ^ x   At >:  u, CAu, uh, a, b, c, d
    j    > o > o >   At o:  h, KE
    j-1  x ^ x ^ x
        i-1  i  i+1  At x & ^:
           i  i+1    At > & o:

The boundaries always run through q grid points (x).

Data Types
type	coriolisadv_cs
	Control structure for mom_coriolisadv. More...

Functions/Subroutines
subroutine, public	coradcalc (u, v, h, uh, vh, CAu, CAv, OBC, AD, G, GV, CS)
	Calculates the Coriolis and momentum advection contributions to the acceleration. More...
subroutine	gradke (u, v, h, KE, KEx, KEy, k, OBC, G, CS)
	Calculates the acceleration due to the gradient of kinetic energy. More...
subroutine, public	coriolisadv_init (Time, G, param_file, diag, AD, CS)
	Initializes the control structure for coriolisadv_cs. More...
subroutine, public	coriolisadv_end (CS)
	Destructor for coriolisadv_cs. More...

Variables
integer, parameter	sadourny75_energy = 1
integer, parameter	arakawa_hsu90 = 2
integer, parameter	robust_enstro = 3
integer, parameter	sadourny75_enstro = 4
integer, parameter	arakawa_lamb81 = 5
integer, parameter	al_blend = 6
character *(20), parameter	sadourny75_energy_string = "SADOURNY75_ENERGY"
character *(20), parameter	arakawa_hsu_string = "ARAKAWA_HSU90"
character *(20), parameter	robust_enstro_string = "ROBUST_ENSTRO"
character *(20), parameter	sadourny75_enstro_string = "SADOURNY75_ENSTRO"
character *(20), parameter	arakawa_lamb_string = "ARAKAWA_LAMB81"
character *(20), parameter	al_blend_string = "ARAKAWA_LAMB_BLEND"
integer, parameter	ke_arakawa = 10
integer, parameter	ke_simple_gudonov = 11
integer, parameter	ke_gudonov = 12
character *(20), parameter	ke_arakawa_string = "KE_ARAKAWA"
character *(20), parameter	ke_simple_gudonov_string = "KE_SIMPLE_GUDONOV"
character *(20), parameter	ke_gudonov_string = "KE_GUDONOV"
integer, parameter	pv_adv_centered = 21
integer, parameter	pv_adv_upwind1 = 22
character *(20), parameter	pv_adv_centered_string = "PV_ADV_CENTERED"
character *(20), parameter	pv_adv_upwind1_string = "PV_ADV_UPWIND1"

Function/Subroutine Documentation

coradcalc()

subroutine, public mom_coriolisadv::coradcalc	(	real, dimension(szib_(g),szj_(g),szk_(g)), intent(in)	u,
		real, dimension(szi_(g),szjb_(g),szk_(g)), intent(in)	v,
		real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	h,
		real, dimension(szib_(g),szj_(g),szk_(g)), intent(in)	uh,
		real, dimension(szi_(g),szjb_(g),szk_(g)), intent(in)	vh,
		real, dimension(szib_(g),szj_(g),szk_(g)), intent(out)	CAu,
		real, dimension(szi_(g),szjb_(g),szk_(g)), intent(out)	CAv,
		type(ocean_obc_type), pointer	OBC,
		type(accel_diag_ptrs), intent(inout)	AD,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(coriolisadv_cs), pointer	CS
	)

Calculates the Coriolis and momentum advection contributions to the acceleration.

Parameters

[in]	g	Ocen grid structure
[in]	gv	Vertical grid structure
[in]	u	Zonal velocity (m/s)
[in]	v	Meridional velocity (m/s)
[in]	h	Layer thickness (m or kg/m2)
[in]	uh	Zonal transport u*h*dy (m3/s or kg/s)
[in]	vh	Meridional transport v*h*dx (m3/s or kg/s)
[out]	cau	Zonal acceleration due to Coriolis and momentum advection, in m/s2.
[out]	cav	Meridional acceleration due to Coriolis and momentum advection, in m/s2.
	obc	Open boundary control structure
[in,out]	ad	Storage for acceleration diagnostics
	cs	Control structure for MOM_CoriolisAdv

Definition at line 107 of file MOM_CoriolisAdv.F90.

References al_blend, arakawa_hsu90, arakawa_lamb81, gradke(), mom_error_handler::mom_error(), mom_open_boundary::obc_direction_n, pv_adv_centered, pv_adv_upwind1, robust_enstro, sadourny75_energy, and sadourny75_enstro.

Referenced by mom_dynamics_legacy_split::step_mom_dyn_legacy_split(), mom_dynamics_split_rk2::step_mom_dyn_split_rk2(), mom_dynamics_unsplit::step_mom_dyn_unsplit(), and mom_dynamics_unsplit_rk2::step_mom_dyn_unsplit_rk2().

  type(ocean_grid_type),                     intent(in)    :: g !< Ocen grid structure
  type(verticalgrid_type),                   intent(in)    :: gv !< Vertical grid structure
  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)), intent(in)    :: u !< Zonal velocity (m/s)
  real, dimension(SZI_(G),SZJB_(G),SZK_(G)), intent(in)    :: v !< Meridional velocity (m/s)
  real, dimension(SZI_(G),SZJ_(G),SZK_(G)),  intent(in)    :: h !< Layer thickness (m or kg/m2)
  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)), intent(in)    :: uh !< Zonal transport u*h*dy (m3/s or kg/s)
  real, dimension(SZI_(G),SZJB_(G),SZK_(G)), intent(in)    :: vh !< Meridional transport v*h*dx (m3/s or kg/s)
  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)), intent(out)   :: cau !< Zonal acceleration due to Coriolis
                                                                  !! and momentum advection, in m/s2.
  real, dimension(SZI_(G),SZJB_(G),SZK_(G)), intent(out)   :: cav !< Meridional acceleration due to Coriolis
                                                                  !! and momentum advection, in m/s2.
  type(ocean_obc_type),                      pointer       :: obc !< Open boundary control structure
  type(accel_diag_ptrs),                     intent(inout) :: ad !< Storage for acceleration diagnostics
  type(coriolisadv_cs),                      pointer       :: cs !< Control structure for MOM_CoriolisAdv
  ! Local variables

  real, dimension(SZIB_(G),SZJB_(G)) :: &
    q, &        ! Layer potential vorticity, in m-1 s-1.
    ih_q, &     ! The inverse of thickness interpolated to q pointes, in
                ! units of m-1 or m2 kg-1.
    area_q      ! The sum of the ocean areas at the 4 adjacent thickness
                ! points, in m2.

  real, dimension(SZIB_(G),SZJ_(G)) :: &
    a, b, c, d  ! a, b, c, & d are combinations of the potential vorticities
                ! surrounding an h grid point.  At small scales, a = q/4,
                ! b = q/4, etc.  All are in units of m-1 s-1 or m2 kg-1 s-1,
                ! and use the indexing of the corresponding u point.

  real, dimension(SZI_(G),SZJ_(G)) :: &
    area_h, &   ! The ocean area at h points, in m2.  Area_h is used to find the
                ! average thickness in the denominator of q.  0 for land points.
    ke          ! Kinetic energy per unit mass, KE = (u^2 + v^2)/2, in m2 s-2.
  real, dimension(SZIB_(G),SZJ_(G)) :: &
    harea_u, &  ! The cell area weighted thickness interpolated to u points
                ! times the effective areas, in H m2.
    kex, &      ! The zonal gradient of Kinetic energy per unit mass,
                ! KEx = d/dx KE, in m s-2.
    uh_center   ! centered u times h at u-points
  real, dimension(SZI_(G),SZJB_(G)) :: &
    harea_v, &  ! The cell area weighted thickness interpolated to v points
                ! times the effective areas, in H m2.
    key, &      ! The meridonal gradient of Kinetic energy per unit mass,
                ! KEy = d/dy KE, in m s-2.
    vh_center   ! centered v times h at v-points
  real, dimension(SZI_(G),SZJ_(G)) :: &
    uh_min, uh_max, &   ! The smallest and largest estimates of the volume
    vh_min, vh_max, &   ! fluxes through the faces (i.e. u*h*dy & v*h*dx),
                        ! in m3 s-1 or kg s-1.
    ep_u, ep_v  ! Additional pseudo-Coriolis terms in the Arakawa and Lamb
                ! discretization, in m-1 s-1 or m2 kg-1 s-1.
  real, dimension(SZIB_(G),SZJB_(G)) :: &
    dvdx,dudy, &! Contributions to the circulation around q-points (m2 s-1)
    abs_vort, & ! Absolute vorticity at q-points, in s-1.
    q2, &       ! Relative vorticity over thickness.
    max_fvq, &  ! The maximum or minimum of the
    min_fvq, &  ! adjacent values of (-u) or v times
    max_fuq, &  ! the absolute vorticity, in m s-2.
    min_fuq     ! All are defined at q points.
  real, dimension(SZIB_(G),SZJB_(G),SZK_(G)) :: &
    pv, &       ! A diagnostic array of the potential vorticities, in m-1 s-1.
    rv          ! A diagnostic array of the relative vorticities, in s-1.
  real :: fv1, fv2, fu1, fu2   ! (f+rv)*v or (f+rv)*u in m s-2.
  real :: max_fv, max_fu       ! The maximum or minimum of the neighbor-
  real :: min_fv, min_fu       ! max(min)_fu(v)q, in m s-2.

  real, parameter :: c1_12=1.0/12.0 ! C1_12 = 1/12
  real, parameter :: c1_24=1.0/24.0 ! C1_24 = 1/24
  real :: absolute_vorticity     ! Absolute vorticity, in s-1.
  real :: relative_vorticity     ! Relative vorticity, in s-1.
  real :: ih                     ! Inverse of thickness, m-1 or m2 kg-1.
  real :: max_ihq, min_ihq       ! The maximum and minimum of the nearby Ihq.
  real :: harea_q                ! The sum of area times thickness of the cells
                                 ! surrounding a q point, in m3 or kg.
  real :: h_neglect              ! A thickness that is so small it is usually
                                 ! lost in roundoff and can be neglected, in m.
  real :: temp1, temp2           ! Temporary variables, in m2 s-2.
  real, parameter :: eps_vel=1.0e-10 ! A tiny, positive velocity, in m s-1.

  real :: uhc, vhc               ! Centered estimates of uh and vh in m3 s-1 or kg s-1.
  real :: uhm, vhm               ! The input estimates of uh and vh in m3 s-1 or kg s-1.
  real :: c1, c2, c3, slope      ! Nondimensional parameters for the Coriolis
                                 ! limiter scheme.

  real :: fe_m2         ! Nondimensional temporary variables asssociated with
  real :: rat_lin       ! the ARAKAWA_LAMB_BLEND scheme.
  real :: rat_m1        ! The ratio of the maximum neighboring inverse thickness
                        ! to the minimum inverse thickness minus 1. rat_m1 >= 0.
  real :: al_wt         ! The relative weight of the Arakawa & Lamb scheme to the
                        ! Arakawa & Hsu scheme, nondimensional between 0 and 1.
  real :: sad_wt        ! The relative weight of the Sadourny energy scheme to
                        ! the other two with the ARAKAWA_LAMB_BLEND scheme,
                        ! nondimensional between 0 and 1.

  real :: heff1, heff2  ! Temporary effective H at U or V points in m or kg m-2.
  real :: heff3, heff4  ! Temporary effective H at U or V points in m or kg m-2.
  real :: h_tiny        ! A very small thickness, in m or kg m-2.
  real :: uheff, vheff  ! More temporary variables, in m3 s-1 or kg s-1.
  real :: quheff,qvheff ! More temporary variables, in m3 s-2 or kg s-2.
  integer :: i, j, k, n, is, ie, js, je, isq, ieq, jsq, jeq, nz

! To work, the following fields must be set outside of the usual
! is to ie range before this subroutine is called:
!  v[is-1,ie+1,ie+2], u[is-1,ie+1], vh[ie+1], uh[is-1], and
!  h[is-1,ie+1,ie+2].
! In the y-direction, the following fields must be set:
!  v[js-1,je+1], u[js-1,je+1,je+2], vh[js-1], uh[je+1], and
!  h[js-1,je+1,je+2].

  if (.not.associated(cs)) call mom_error(fatal, &
         "MOM_CoriolisAdv: Module must be initialized before it is used.")
  is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec
  isq = g%IscB ; ieq = g%IecB ; jsq = g%JscB ; jeq = g%JecB ; nz = g%ke
  h_neglect = gv%H_subroundoff
  h_tiny = gv%Angstrom  ! Perhaps this should be set to h_neglect instead.

  !$OMP parallel do default(private) shared(Isq,Ieq,Jsq,Jeq,G,Area_h)
  do j=jsq-1,jeq+2 ; do i=isq-1,ieq+2
    area_h(i,j) = g%mask2dT(i,j) * g%areaT(i,j)
  enddo ; enddo
  !$OMP parallel do default(private) shared(Isq,Ieq,Jsq,Jeq,G,Area_h,Area_q)
  do j=jsq-1,jeq+1 ; do i=isq-1,ieq+1
    area_q(i,j) = (area_h(i,j) + area_h(i+1,j+1)) + &
                  (area_h(i+1,j) + area_h(i,j+1))
  enddo ; enddo

  !$OMP parallel do default(private) shared(u,v,h,uh,vh,CAu,CAv,G,CS,AD,Area_h,Area_q,&
  !$OMP                        RV,PV,is,ie,js,je,Isq,Ieq,Jsq,Jeq,nz,h_neglect,h_tiny,OBC)
  do k=1,nz
    ! Here the second order accurate layer potential vorticities, q,
    ! are calculated.  hq is  second order accurate in space.  Relative
    ! vorticity is second order accurate everywhere with free slip b.c.s,
    ! but only first order accurate at boundaries with no slip b.c.s.
    ! First calculate the contributions to the circulation around the q-point.
    do j=jsq-1,jeq+1 ; do i=isq-1,ieq+1
      dvdx(i,j) = v(i+1,j,k)*g%dyCv(i+1,j) - v(i,j,k)*g%dyCv(i,j)
      dudy(i,j) = u(i,j+1,k)*g%dxCu(i,j+1) - u(i,j,k)*g%dxCu(i,j)
    enddo ; enddo
    do j=jsq-1,jeq+1 ; do i=isq-1,ieq+2
      harea_v(i,j) = 0.5*(area_h(i,j) * h(i,j,k) + area_h(i,j+1) * h(i,j+1,k))
    enddo ; enddo
    do j=jsq-1,jeq+2 ; do i=isq-1,ieq+1
      harea_u(i,j) = 0.5*(area_h(i,j) * h(i,j,k) + area_h(i+1,j) * h(i+1,j,k))
    enddo ; enddo
    if (cs%Coriolis_En_Dis) then
      do j=jsq,jeq+1 ; do i=is-1,ie
        uh_center(i,j) = 0.5 * (g%dy_Cu(i,j) * u(i,j,k)) * (h(i,j,k) + h(i+1,j,k))
      enddo ; enddo
      do j=js-1,je ; do i=isq,ieq+1
        vh_center(i,j) = 0.5 * (g%dx_Cv(i,j) * v(i,j,k)) * (h(i,j,k) + h(i,j+1,k))
      enddo ; enddo
    endif

    ! Adjust circulation components to relative vorticity and thickness projected onto
    ! velocity points on open boundaries.
    if (associated(obc)) then ; do n=1,obc%number_of_segments
      if (.not. obc%segment(n)%on_pe) cycle
      i = obc%segment(n)%HI%IsdB ; j = obc%segment(n)%HI%JsdB
      if (obc%segment(n)%is_N_or_S .and. (j >= jsq-1) .and. (j <= jeq+1)) then
        if (obc%zero_vorticity) then ; do i=obc%segment(n)%HI%IsdB,obc%segment(n)%HI%IedB
          dvdx(i,j) = 0. ; dudy(i,j) = 0.
        enddo ; endif
        if (obc%freeslip_vorticity) then ; do i=obc%segment(n)%HI%IsdB,obc%segment(n)%HI%IedB
          dudy(i,j) = 0.
        enddo ; endif

        ! Project thicknesses across OBC points with a no-gradient condition.
        do i = max(isq-1,obc%segment(n)%HI%isd), min(ieq+2,obc%segment(n)%HI%ied)
          if (obc%segment(n)%direction == obc_direction_n) then
            harea_v(i,j) = 0.5 * (area_h(i,j) + area_h(i,j+1)) * h(i,j,k)
          else ! (OBC%segment(n)%direction == OBC_DIRECTION_S)
            harea_v(i,j) = 0.5 * (area_h(i,j) + area_h(i,j+1)) * h(i,j+1,k)
          endif
        enddo

        if (cs%Coriolis_En_Dis) then
          do i = max(isq-1,obc%segment(n)%HI%isd), min(ieq+2,obc%segment(n)%HI%ied)
            if (obc%segment(n)%direction == obc_direction_n) then
              vh_center(i,j) = g%dx_Cv(i,j) * v(i,j,k) * h(i,j,k)
            else ! (OBC%segment(n)%direction == OBC_DIRECTION_S)
              vh_center(i,j) = g%dx_Cv(i,j) * v(i,j,k) * h(i,j+1,k)
            endif
          enddo
        endif
      elseif (obc%segment(n)%is_E_or_W .and. (i >= isq-1) .and. (i <= ieq+1)) then
        if (obc%zero_vorticity) then ; do j=obc%segment(n)%HI%JsdB,obc%segment(n)%HI%JedB
          dvdx(i,j) = 0. ; dudy(i,j) = 0.
        enddo ; endif
        if (obc%freeslip_vorticity) then ; do j=obc%segment(n)%HI%JsdB,obc%segment(n)%HI%JedB
          dvdx(i,j) = 0.
        enddo ; endif

        ! Project thicknesses across OBC points with a no-gradient condition.
        do j = max(jsq-1,obc%segment(n)%HI%jsd), min(jeq+2,obc%segment(n)%HI%jed)
          if (obc%segment(n)%direction == obc_direction_e) then
            harea_u(i,j) = 0.5*(area_h(i,j) + area_h(i+1,j)) * h(i,j,k)
          else ! (OBC%segment(n)%direction == OBC_DIRECTION_W)
            harea_u(i,j) = 0.5*(area_h(i,j) + area_h(i+1,j)) * h(i+1,j,k)
          endif
        enddo
        if (cs%Coriolis_En_Dis) then
          do j = max(jsq-1,obc%segment(n)%HI%jsd), min(jeq+2,obc%segment(n)%HI%jed)
            if (obc%segment(n)%direction == obc_direction_e) then
              uh_center(i,j) = g%dy_Cu(i,j) * u(i,j,k) * h(i,j,k)
            else ! (OBC%segment(n)%direction == OBC_DIRECTION_W)
              uh_center(i,j) = g%dy_Cu(i,j) * u(i,j,k) * h(i+1,j,k)
            endif
          enddo
        endif
      endif
    enddo ; endif

    if (associated(obc)) then ; do n=1,obc%number_of_segments
      if (.not. obc%segment(n)%on_pe) cycle
      ! Now project thicknesses across cell-corner points in the OBCs.  The two
      ! projections have to occur in sequence and can not be combined easily.
      i = obc%segment(n)%HI%IsdB ; j = obc%segment(n)%HI%JsdB
      if (obc%segment(n)%is_N_or_S .and. (j >= jsq-1) .and. (j <= jeq+1)) then
        do i = max(isq-1,obc%segment(n)%HI%IsdB), min(ieq+1,obc%segment(n)%HI%IedB)
          if (obc%segment(n)%direction == obc_direction_n) then
            if (area_h(i,j) + area_h(i+1,j) > 0.0) then
              harea_u(i,j+1) = harea_u(i,j) * ((area_h(i,j+1) + area_h(i+1,j+1)) / &
                                               (area_h(i,j) + area_h(i+1,j)))
            else ; harea_u(i,j+1) = 0.0 ; endif
          else ! (OBC%segment(n)%direction == OBC_DIRECTION_S)
            if (area_h(i,j+1) + area_h(i+1,j+1) > 0.0) then
              harea_u(i,j) = harea_u(i,j+1) * ((area_h(i,j) + area_h(i+1,j)) / &
                                               (area_h(i,j+1) + area_h(i+1,j+1)))
            else ; harea_u(i,j) = 0.0 ; endif
          endif
        enddo
      elseif (obc%segment(n)%is_E_or_W .and. (i >= isq-1) .and. (i <= ieq+1)) then
        do j = max(jsq-1,obc%segment(n)%HI%JsdB), min(jeq+1,obc%segment(n)%HI%JedB)
          if (obc%segment(n)%direction == obc_direction_e) then
            if (area_h(i,j) + area_h(i,j+1) > 0.0) then
              harea_v(i+1,j) = harea_v(i,j) * ((area_h(i+1,j) + area_h(i+1,j+1)) / &
                                               (area_h(i,j) + area_h(i,j+1)))
            else ; harea_v(i+1,j) = 0.0 ; endif
          else ! (OBC%segment(n)%direction == OBC_DIRECTION_W)
            harea_v(i,j) = 0.5 * (area_h(i,j) + area_h(i,j+1)) * h(i,j+1,k)
            if (area_h(i+1,j) + area_h(i+1,j+1) > 0.0) then
              harea_v(i,j) = harea_v(i+1,j) * ((area_h(i,j) + area_h(i,j+1)) / &
                                               (area_h(i+1,j) + area_h(i+1,j+1)))
            else ; harea_v(i,j) = 0.0 ; endif
          endif
        enddo
      endif
    enddo ; endif

    do j=jsq-1,jeq+1 ; do i=isq-1,ieq+1
      if (cs%no_slip ) then
        relative_vorticity = (2.0-g%mask2dBu(i,j)) * (dvdx(i,j) - dudy(i,j)) * &
                             g%IareaBu(i,j)
      else
        relative_vorticity = g%mask2dBu(i,j) * (dvdx(i,j) - dudy(i,j)) * &
                             g%IareaBu(i,j)
      endif
      absolute_vorticity = g%CoriolisBu(i,j) + relative_vorticity
      ih = 0.0
      if (area_q(i,j) > 0.0) then
        harea_q = (harea_u(i,j) + harea_u(i,j+1)) + (harea_v(i,j) + harea_v(i+1,j))
        ih = area_q(i,j) / (harea_q + h_neglect*area_q(i,j))
      endif
      q(i,j) = absolute_vorticity * ih
      abs_vort(i,j) = absolute_vorticity
      ih_q(i,j) = ih

      if (cs%bound_Coriolis) then
        fv1 = absolute_vorticity*v(i+1,j,k)
        fv2 = absolute_vorticity*v(i,j,k)
        fu1 = -absolute_vorticity*u(i,j+1,k)
        fu2 = -absolute_vorticity*u(i,j,k)
        if (fv1 > fv2) then
          max_fvq(i,j) = fv1 ; min_fvq(i,j) = fv2
        else
          max_fvq(i,j) = fv2 ; min_fvq(i,j) = fv1
        endif
        if (fu1 > fu2) then
          max_fuq(i,j) = fu1 ; min_fuq(i,j) = fu2
        else
          max_fuq(i,j) = fu2 ; min_fuq(i,j) = fu1
        endif
      endif

      if (cs%id_rv > 0) rv(i,j,k) = relative_vorticity
      if (cs%id_PV > 0) pv(i,j,k) = q(i,j)
      if (ASSOCIATED(ad%rv_x_v) .or. ASSOCIATED(ad%rv_x_u)) &
        q2(i,j) = relative_vorticity * ih
    enddo ; enddo

    !   a, b, c, and d are combinations of neighboring potential
    ! vorticities which form the Arakawa and Hsu vorticity advection
    ! scheme.  All are defined at u grid points.

    if (cs%Coriolis_Scheme == arakawa_hsu90) then
      do j=jsq,jeq+1
        do i=is-1,ieq
          a(i,j) = (q(i,j) + (q(i+1,j) + q(i,j-1))) * c1_12
          d(i,j) = ((q(i,j) + q(i+1,j-1)) + q(i,j-1)) * c1_12
        enddo
        do i=isq,ieq
          b(i,j) = (q(i,j) + (q(i-1,j) + q(i,j-1))) * c1_12
          c(i,j) = ((q(i,j) + q(i-1,j-1)) + q(i,j-1)) * c1_12
        enddo
      enddo
    elseif (cs%Coriolis_Scheme == arakawa_lamb81) then
      do j=jsq,jeq+1 ; do i=isq,ieq+1
        a(i-1,j) = (2.0*(q(i,j) + q(i-1,j-1)) + (q(i-1,j) + q(i,j-1))) * c1_24
        d(i-1,j) = ((q(i,j) + q(i-1,j-1)) + 2.0*(q(i-1,j) + q(i,j-1))) * c1_24
        b(i,j) =   ((q(i,j) + q(i-1,j-1)) + 2.0*(q(i-1,j) + q(i,j-1))) * c1_24
        c(i,j) =   (2.0*(q(i,j) + q(i-1,j-1)) + (q(i-1,j) + q(i,j-1))) * c1_24
        ep_u(i,j) = ((q(i,j) - q(i-1,j-1)) + (q(i-1,j) - q(i,j-1))) * c1_24
        ep_v(i,j) = (-(q(i,j) - q(i-1,j-1)) + (q(i-1,j) - q(i,j-1))) * c1_24
      enddo ; enddo
    elseif (cs%Coriolis_Scheme == al_blend) then
      fe_m2 = cs%F_eff_max_blend - 2.0
      rat_lin = 1.5 * fe_m2 / max(cs%wt_lin_blend, 1.0e-16)

      ! This allows the code to always give Sadourny Energy
      if (cs%F_eff_max_blend <= 2.0) then ; fe_m2 = -1. ; rat_lin = -1.0 ; endif

      do j=jsq,jeq+1 ; do i=isq,ieq+1
        min_ihq = min(ih_q(i-1,j-1), ih_q(i,j-1), ih_q(i-1,j), ih_q(i,j))
        max_ihq = max(ih_q(i-1,j-1), ih_q(i,j-1), ih_q(i-1,j), ih_q(i,j))
        rat_m1 = 1.0e15
        if (max_ihq < 1.0e15*min_ihq) rat_m1 = max_ihq / min_ihq - 1.0
        ! The weights used here are designed to keep the effective Coriolis
        ! acceleration from any one point on its neighbors within a factor
        ! of F_eff_max.  The minimum permitted value is 2 (the factor for
        ! Sadourny's energy conserving scheme).

        ! Determine the relative weights of Arakawa & Lamb vs. Arakawa and Hsu.
        if (rat_m1 <= fe_m2) then ; al_wt = 1.0
        elseif (rat_m1 < 1.5*fe_m2) then ; al_wt = 3.0*fe_m2 / rat_m1 - 2.0
        else ; al_wt = 0.0 ; endif

        ! Determine the relative weights of Sadourny Energy vs. the other two.
        if (rat_m1 <= 1.5*fe_m2) then ; sad_wt = 0.0
        elseif (rat_m1 <= rat_lin) then
          sad_wt = 1.0 - (1.5*fe_m2) / rat_m1
        elseif (rat_m1 < 2.0*rat_lin) then
          sad_wt = 1.0 - (cs%wt_lin_blend / rat_lin) * (rat_m1 - 2.0*rat_lin)
        else ; sad_wt = 1.0 ; endif

        a(i-1,j) = sad_wt * 0.25 * q(i-1,j) + (1.0 - sad_wt) * &
                   ( ((2.0-al_wt)* q(i-1,j) + al_wt*q(i,j-1)) + &
                      2.0 * (q(i,j) + q(i-1,j-1)) ) * c1_24
        d(i-1,j) = sad_wt * 0.25 * q(i-1,j-1) + (1.0 - sad_wt) * &
                   ( ((2.0-al_wt)* q(i-1,j-1) + al_wt*q(i,j)) + &
                      2.0 * (q(i-1,j) + q(i,j-1)) ) * c1_24
        b(i,j) =   sad_wt * 0.25 * q(i,j) + (1.0 - sad_wt) * &
                   ( ((2.0-al_wt)* q(i,j) + al_wt*q(i-1,j-1)) + &
                      2.0 * (q(i-1,j) + q(i,j-1)) ) * c1_24
        c(i,j) =   sad_wt * 0.25 * q(i,j-1) + (1.0 - sad_wt) * &
                   ( ((2.0-al_wt)* q(i,j-1) + al_wt*q(i-1,j)) + &
                      2.0 * (q(i,j) + q(i-1,j-1)) ) * c1_24
        ep_u(i,j) = al_wt  * ((q(i,j) - q(i-1,j-1)) + (q(i-1,j) - q(i,j-1))) * c1_24
        ep_v(i,j) = al_wt * (-(q(i,j) - q(i-1,j-1)) + (q(i-1,j) - q(i,j-1))) * c1_24
      enddo ; enddo
    endif

    if (cs%Coriolis_En_Dis) then
    !  c1 = 1.0-1.5*RANGE ; c2 = 1.0-RANGE ; c3 = 2.0 ; slope = 0.5
      c1 = 1.0-1.5*0.5 ; c2 = 1.0-0.5 ; c3 = 2.0 ; slope = 0.5

      do j=jsq,jeq+1 ; do i=is-1,ie
        uhc = uh_center(i,j)
        uhm = uh(i,j,k)
        ! This sometimes matters with some types of open boundary conditions.
        if (g%dy_Cu(i,j) == 0.0) uhc = uhm

        if (abs(uhc) < 0.1*abs(uhm)) then
          uhm = 10.0*uhc
        elseif (abs(uhc) > c1*abs(uhm)) then
          if (abs(uhc) < c2*abs(uhm)) then ; uhc = (3.0*uhc+(1.0-c2*3.0)*uhm)
          elseif (abs(uhc) <= c3*abs(uhm)) then ; uhc = uhm
          else ; uhc = slope*uhc+(1.0-c3*slope)*uhm
          endif
        endif

        if (uhc > uhm) then
          uh_min(i,j) = uhm ; uh_max(i,j) = uhc
        else
          uh_max(i,j) = uhm ; uh_min(i,j) = uhc
        endif
      enddo ; enddo
      do j=js-1,je ; do i=isq,ieq+1
        vhc = vh_center(i,j)
        vhm = vh(i,j,k)
        ! This sometimes matters with some types of open boundary conditions.
        if (g%dx_Cv(i,j) == 0.0) vhc = vhm

        if (abs(vhc) < 0.1*abs(vhm)) then
          vhm = 10.0*vhc
        elseif (abs(vhc) > c1*abs(vhm)) then
          if (abs(vhc) < c2*abs(vhm)) then ; vhc = (3.0*vhc+(1.0-c2*3.0)*vhm)
          else if (abs(vhc) <= c3*abs(vhm)) then ; vhc = vhm
          else ; vhc = slope*vhc+(1.0-c3*slope)*vhm
          endif
        endif

        if (vhc > vhm) then
          vh_min(i,j) = vhm ; vh_max(i,j) = vhc
        else
          vh_max(i,j) = vhm ; vh_min(i,j) = vhc
        endif
      enddo ; enddo
    endif

    ! Calculate KE and the gradient of KE
    call gradke(u, v, h, ke, kex, key, k, obc, g, cs)

    ! Calculate the tendencies of zonal velocity due to the Coriolis
    ! force and momentum advection.  On a Cartesian grid, this is
    !     CAu =  q * vh - d(KE)/dx.
    if (cs%Coriolis_Scheme == sadourny75_energy) then
      if (cs%Coriolis_En_Dis) then
        ! Energy dissipating biased scheme, Hallberg 200x
        do j=js,je ; do i=isq,ieq
          if (q(i,j)*u(i,j,k) == 0.0) then
            temp1 = q(i,j) * ( (vh_max(i,j)+vh_max(i+1,j)) &
                             + (vh_min(i,j)+vh_min(i+1,j)) )*0.5
          elseif (q(i,j)*u(i,j,k) < 0.0) then
            temp1 = q(i,j) * (vh_max(i,j)+vh_max(i+1,j))
          else
            temp1 = q(i,j) * (vh_min(i,j)+vh_min(i+1,j))
          endif
          if (q(i,j-1)*u(i,j,k) == 0.0) then
            temp2 = q(i,j-1) * ( (vh_max(i,j-1)+vh_max(i+1,j-1)) &
                               + (vh_min(i,j-1)+vh_min(i+1,j-1)) )*0.5
          elseif (q(i,j-1)*u(i,j,k) < 0.0) then
            temp2 = q(i,j-1) * (vh_max(i,j-1)+vh_max(i+1,j-1))
          else
            temp2 = q(i,j-1) * (vh_min(i,j-1)+vh_min(i+1,j-1))
          endif
          cau(i,j,k) = 0.25 * g%IdxCu(i,j) * (temp1 + temp2)
        enddo ; enddo
      else
        ! Energy conserving scheme, Sadourny 1975
        do j=js,je ; do i=isq,ieq
          cau(i,j,k) = 0.25 * &
            (q(i,j) * (vh(i+1,j,k) + vh(i,j,k)) + &
             q(i,j-1) * (vh(i,j-1,k) + vh(i+1,j-1,k))) * g%IdxCu(i,j)
        enddo ; enddo
      endif
    elseif (cs%Coriolis_Scheme == sadourny75_enstro) then
      do j=js,je ; do i=isq,ieq
        cau(i,j,k) = 0.125 * (g%IdxCu(i,j) * (q(i,j) + q(i,j-1))) * &
                     ((vh(i+1,j,k) + vh(i,j,k)) + (vh(i,j-1,k) + vh(i+1,j-1,k)))
      enddo ; enddo
    elseif ((cs%Coriolis_Scheme == arakawa_hsu90) .or. &
            (cs%Coriolis_Scheme == arakawa_lamb81) .or. &
            (cs%Coriolis_Scheme == al_blend)) then
      ! (Global) Energy and (Local) Enstrophy conserving, Arakawa & Hsu 1990
      do j=js,je ; do i=isq,ieq
        cau(i,j,k) = ((a(i,j) * vh(i+1,j,k) + &
                       c(i,j) * vh(i,j-1,k))  &
                    + (b(i,j) * vh(i,j,k) +   &
                       d(i,j) * vh(i+1,j-1,k))) * g%IdxCu(i,j)
      enddo ; enddo
    elseif (cs%Coriolis_Scheme == robust_enstro) then
      ! An enstrophy conserving scheme robust to vanishing layers
      ! Note: Heffs are in lieu of h_at_v that should be returned by the
      !       continuity solver. AJA
      do j=js,je ; do i=isq,ieq
        heff1 = abs(vh(i,j,k)*g%IdxCv(i,j))/(eps_vel+abs(v(i,j,k)))
        heff1 = max(heff1,min(h(i,j,k),h(i,j+1,k)))
        heff1 = min(heff1,max(h(i,j,k),h(i,j+1,k)))
        heff2 = abs(vh(i,j-1,k)*g%IdxCv(i,j-1))/(eps_vel+abs(v(i,j-1,k)))
        heff2 = max(heff2,min(h(i,j-1,k),h(i,j,k)))
        heff2 = min(heff2,max(h(i,j-1,k),h(i,j,k)))
        heff3 = abs(vh(i+1,j,k)*g%IdxCv(i+1,j))/(eps_vel+abs(v(i+1,j,k)))
        heff3 = max(heff3,min(h(i+1,j,k),h(i+1,j+1,k)))
        heff3 = min(heff3,max(h(i+1,j,k),h(i+1,j+1,k)))
        heff4 = abs(vh(i+1,j-1,k)*g%IdxCv(i+1,j-1))/(eps_vel+abs(v(i+1,j-1,k)))
        heff4 = max(heff4,min(h(i+1,j-1,k),h(i+1,j,k)))
        heff4 = min(heff4,max(h(i+1,j-1,k),h(i+1,j,k)))
        if (cs%PV_Adv_Scheme == pv_adv_centered) then
          cau(i,j,k) = 0.5*(abs_vort(i,j)+abs_vort(i,j-1)) * &
                       ((vh(i  ,j  ,k)+vh(i+1,j-1,k)) +      &
                        (vh(i  ,j-1,k)+vh(i+1,j  ,k)) ) /    &
                       (h_tiny +((heff1+heff4) +(heff2+heff3)) ) * g%IdxCu(i,j)
        elseif (cs%PV_Adv_Scheme == pv_adv_upwind1) then
          vheff = ((vh(i  ,j  ,k)+vh(i+1,j-1,k)) +      &
                   (vh(i  ,j-1,k)+vh(i+1,j  ,k)) )
          qvheff = 0.5*( (abs_vort(i,j)+abs_vort(i,j-1))*vheff &
                        -(abs_vort(i,j)-abs_vort(i,j-1))*abs(vheff) )
          cau(i,j,k) = qvheff / &
                     (h_tiny +((heff1+heff4) +(heff2+heff3)) ) * g%IdxCu(i,j)
        endif
      enddo ; enddo
    endif
    ! Add in the additonal terms with Arakawa & Lamb.
    if ((cs%Coriolis_Scheme == arakawa_lamb81) .or. &
        (cs%Coriolis_Scheme == al_blend)) then ; do j=js,je ; do i=isq,ieq
      cau(i,j,k) = cau(i,j,k) + &
            (ep_u(i,j)*uh(i-1,j,k) - ep_u(i+1,j)*uh(i+1,j,k)) * g%IdxCu(i,j)
    enddo ; enddo ; endif


    if (cs%bound_Coriolis) then
      do j=js,je ; do i=isq,ieq
        max_fv = max(max_fvq(i,j),max_fvq(i,j-1))
        min_fv = min(min_fvq(i,j),min_fvq(i,j-1))
       ! CAu(I,j,k) = min( CAu(I,j,k), max_fv )
       ! CAu(I,j,k) = max( CAu(I,j,k), min_fv )
        if (cau(i,j,k) > max_fv) then
            cau(i,j,k) = max_fv
        else
          if (cau(i,j,k) < min_fv) cau(i,j,k) = min_fv
        endif
      enddo ; enddo
    endif

    ! Term - d(KE)/dx.
    do j=js,je ; do i=isq,ieq
      cau(i,j,k) = cau(i,j,k) - kex(i,j)
      if (ASSOCIATED(ad%gradKEu)) ad%gradKEu(i,j,k) = -kex(i,j)
    enddo ; enddo


    ! Calculate the tendencies of meridional velocity due to the Coriolis
    ! force and momentum advection.  On a Cartesian grid, this is
    !     CAv = - q * uh - d(KE)/dy.
    if (cs%Coriolis_Scheme == sadourny75_energy) then
      if (cs%Coriolis_En_Dis) then
        ! Energy dissipating biased scheme, Hallberg 200x
        do j=jsq,jeq ; do i=is,ie
          if (q(i-1,j)*v(i,j,k) == 0.0) then
            temp1 = q(i-1,j) * ( (uh_max(i-1,j)+uh_max(i-1,j+1)) &
                               + (uh_min(i-1,j)+uh_min(i-1,j+1)) )*0.5
          elseif (q(i-1,j)*v(i,j,k) > 0.0) then
            temp1 = q(i-1,j) * (uh_max(i-1,j)+uh_max(i-1,j+1))
          else
            temp1 = q(i-1,j) * (uh_min(i-1,j)+uh_min(i-1,j+1))
          endif
          if (q(i,j)*v(i,j,k) == 0.0) then
            temp2 = q(i,j) * ( (uh_max(i,j)+uh_max(i,j+1)) &
                             + (uh_min(i,j)+uh_min(i,j+1)) )*0.5
          elseif (q(i,j)*v(i,j,k) > 0.0) then
            temp2 = q(i,j) * (uh_max(i,j)+uh_max(i,j+1))
          else
            temp2 = q(i,j) * (uh_min(i,j)+uh_min(i,j+1))
          endif
          cav(i,j,k) = - 0.25 * g%IdyCv(i,j) * (temp1 + temp2)
        enddo ; enddo
      else
        ! Energy conserving scheme, Sadourny 1975
        do j=jsq,jeq ; do i=is,ie
          cav(i,j,k) = - 0.25* &
              (q(i-1,j)*(uh(i-1,j,k) + uh(i-1,j+1,k)) + &
               q(i,j)*(uh(i,j,k) + uh(i,j+1,k))) * g%IdyCv(i,j)
        enddo ; enddo
      endif
    elseif (cs%Coriolis_Scheme == sadourny75_enstro) then
      do j=jsq,jeq ; do i=is,ie
        cav(i,j,k) = -0.125 * (g%IdyCv(i,j) * (q(i-1,j) + q(i,j))) * &
                     ((uh(i-1,j,k) + uh(i-1,j+1,k)) + (uh(i,j,k) + uh(i,j+1,k)))
      enddo ; enddo
    elseif ((cs%Coriolis_Scheme == arakawa_hsu90) .or. &
            (cs%Coriolis_Scheme == arakawa_lamb81) .or. &
            (cs%Coriolis_Scheme == al_blend)) then
      ! (Global) Energy and (Local) Enstrophy conserving, Arakawa & Hsu 1990
      do j=jsq,jeq ; do i=is,ie
        cav(i,j,k) = - ((a(i-1,j)   * uh(i-1,j,k) + &
                         c(i,j+1)   * uh(i,j+1,k))  &
                      + (b(i,j)     * uh(i,j,k) +   &
                         d(i-1,j+1) * uh(i-1,j+1,k))) * g%IdyCv(i,j)
      enddo ; enddo
    elseif (cs%Coriolis_Scheme == robust_enstro) then
      ! An enstrophy conserving scheme robust to vanishing layers
      ! Note: Heffs are in lieu of h_at_u that should be returned by the
      !       continuity solver. AJA
      do j=jsq,jeq ; do i=is,ie
        heff1 = abs(uh(i,j,k)*g%IdyCu(i,j))/(eps_vel+abs(u(i,j,k)))
        heff1 = max(heff1,min(h(i,j,k),h(i+1,j,k)))
        heff1 = min(heff1,max(h(i,j,k),h(i+1,j,k)))
        heff2 = abs(uh(i-1,j,k)*g%IdyCu(i-1,j))/(eps_vel+abs(u(i-1,j,k)))
        heff2 = max(heff2,min(h(i-1,j,k),h(i,j,k)))
        heff2 = min(heff2,max(h(i-1,j,k),h(i,j,k)))
        heff3 = abs(uh(i,j+1,k)*g%IdyCu(i,j+1))/(eps_vel+abs(u(i,j+1,k)))
        heff3 = max(heff3,min(h(i,j+1,k),h(i+1,j+1,k)))
        heff3 = min(heff3,max(h(i,j+1,k),h(i+1,j+1,k)))
        heff4 = abs(uh(i-1,j+1,k)*g%IdyCu(i-1,j+1))/(eps_vel+abs(u(i-1,j+1,k)))
        heff4 = max(heff4,min(h(i-1,j+1,k),h(i,j+1,k)))
        heff4 = min(heff4,max(h(i-1,j+1,k),h(i,j+1,k)))
        if (cs%PV_Adv_Scheme == pv_adv_centered) then
          cav(i,j,k) = - 0.5*(abs_vort(i,j)+abs_vort(i-1,j)) * &
                         ((uh(i  ,j  ,k)+uh(i-1,j+1,k)) +      &
                          (uh(i-1,j  ,k)+uh(i  ,j+1,k)) ) /    &
                      (h_tiny + ((heff1+heff4) +(heff2+heff3)) ) * g%IdyCv(i,j)
        elseif (cs%PV_Adv_Scheme == pv_adv_upwind1) then
          uheff = ((uh(i  ,j  ,k)+uh(i-1,j+1,k)) +      &
                   (uh(i-1,j  ,k)+uh(i  ,j+1,k)) )
          quheff = 0.5*( (abs_vort(i,j)+abs_vort(i-1,j))*uheff &
                        -(abs_vort(i,j)-abs_vort(i-1,j))*abs(uheff) )
          cav(i,j,k) = - quheff / &
                       (h_tiny + ((heff1+heff4) +(heff2+heff3)) ) * g%IdyCv(i,j)
        endif
      enddo ; enddo
    endif
    ! Add in the additonal terms with Arakawa & Lamb.
    if ((cs%Coriolis_Scheme == arakawa_lamb81) .or. &
        (cs%Coriolis_Scheme == al_blend)) then ; do j=jsq,jeq ; do i=is,ie
      cav(i,j,k) = cav(i,j,k) + &
            (ep_v(i,j)*vh(i,j-1,k) - ep_v(i,j+1)*vh(i,j+1,k)) * g%IdyCv(i,j)
    enddo ; enddo ; endif

    if (cs%bound_Coriolis) then
      do j=jsq,jeq ; do i=is,ie
        max_fu = max(max_fuq(i,j),max_fuq(i-1,j))
        min_fu = min(min_fuq(i,j),min_fuq(i-1,j))
        if (cav(i,j,k) > max_fu) then
            cav(i,j,k) = max_fu
        else
          if (cav(i,j,k) < min_fu) cav(i,j,k) = min_fu
        endif
      enddo ; enddo
    endif

    ! Term - d(KE)/dy.
    do j=jsq,jeq ; do i=is,ie
      cav(i,j,k) = cav(i,j,k) - key(i,j)
      if (ASSOCIATED(ad%gradKEv)) ad%gradKEv(i,j,k) = -key(i,j)
    enddo ; enddo

    if (ASSOCIATED(ad%rv_x_u) .or. ASSOCIATED(ad%rv_x_v)) then
      ! Calculate the Coriolis-like acceleration due to relative vorticity.
      if (cs%Coriolis_Scheme == sadourny75_energy) then
        if (ASSOCIATED(ad%rv_x_u)) then
          do j=jsq,jeq ; do i=is,ie
            ad%rv_x_u(i,j,k) = - 0.25* &
              (q2(i-1,j)*(uh(i-1,j,k) + uh(i-1,j+1,k)) + &
               q2(i,j)*(uh(i,j,k) + uh(i,j+1,k))) * g%IdyCv(i,j)
          enddo ; enddo
        endif

        if (ASSOCIATED(ad%rv_x_v)) then
          do j=js,je ; do i=isq,ieq
            ad%rv_x_v(i,j,k) = 0.25 * &
              (q2(i,j) * (vh(i+1,j,k) + vh(i,j,k)) + &
               q2(i,j-1) * (vh(i,j-1,k) + vh(i+1,j-1,k))) * g%IdxCu(i,j)
          enddo ; enddo
        endif
      else
        if (ASSOCIATED(ad%rv_x_u)) then
          do j=jsq,jeq ; do i=is,ie
            ad%rv_x_u(i,j,k) = -g%IdyCv(i,j) * c1_12 * &
              ((q2(i,j) + q2(i-1,j) + q2(i-1,j-1)) * uh(i-1,j,k) + &
               (q2(i-1,j) + q2(i,j) + q2(i,j-1)) * uh(i,j,k) + &
               (q2(i-1,j) + q2(i,j+1) + q2(i,j)) * uh(i,j+1,k) + &
               (q2(i,j) + q2(i-1,j+1) + q2(i-1,j)) * uh(i-1,j+1,k))
          enddo ; enddo
        endif

        if (ASSOCIATED(ad%rv_x_v)) then
          do j=js,je ; do i=isq,ieq
            ad%rv_x_v(i,j,k) = g%IdxCu(i,j) * c1_12 * &
              ((q2(i+1,j) + q2(i,j) + q2(i,j-1)) * vh(i+1,j,k) + &
               (q2(i-1,j) + q2(i,j) + q2(i,j-1)) * vh(i,j,k) + &
               (q2(i-1,j-1) + q2(i,j) + q2(i,j-1)) * vh(i,j-1,k) + &
               (q2(i+1,j-1) + q2(i,j) + q2(i,j-1)) * vh(i+1,j-1,k))
          enddo ; enddo
        endif
      endif
    endif

  enddo ! k-loop.

  ! Here the various Coriolis-related derived quantities are offered for averaging.
  if (query_averaging_enabled(cs%diag)) then
    if (cs%id_rv > 0) call post_data(cs%id_rv, rv, cs%diag)
    if (cs%id_PV > 0) call post_data(cs%id_PV, pv, cs%diag)
    if (cs%id_gKEu>0) call post_data(cs%id_gKEu, ad%gradKEu, cs%diag)
    if (cs%id_gKEv>0) call post_data(cs%id_gKEv, ad%gradKEv, cs%diag)
    if (cs%id_rvxu > 0) call post_data(cs%id_rvxu, ad%rv_x_u, cs%diag)
    if (cs%id_rvxv > 0) call post_data(cs%id_rvxv, ad%rv_x_v, cs%diag)
  endif

Here is the call graph for this function:

Here is the caller graph for this function:

coriolisadv_end()

subroutine, public mom_coriolisadv::coriolisadv_end

(

type(coriolisadv_cs), pointer

CS

)

Destructor for coriolisadv_cs.

Parameters

cs	Control structure fro MOM_CoriolisAdv

Definition at line 1045 of file MOM_CoriolisAdv.F90.

  type(coriolisadv_cs), pointer :: cs !< Control structure fro MOM_CoriolisAdv
  deallocate(cs)

coriolisadv_init()

subroutine, public mom_coriolisadv::coriolisadv_init	(	type(time_type), intent(in), target	Time,
		type(ocean_grid_type), intent(in)	G,
		type(param_file_type), intent(in)	param_file,
		type(diag_ctrl), intent(inout), target	diag,
		type(accel_diag_ptrs), intent(inout), target	AD,
		type(coriolisadv_cs), pointer	CS
	)

Initializes the control structure for coriolisadv_cs.

Parameters

[in]	time	Current model time
[in]	g	Ocean grid structure
[in]	param_file	Runtime parameter handles
[in,out]	diag	Diagnostics control structure
[in,out]	ad	Strorage for acceleration diagnostics
	cs	Control structure fro MOM_CoriolisAdv

Definition at line 875 of file MOM_CoriolisAdv.F90.

References al_blend, al_blend_string, arakawa_hsu90, arakawa_hsu_string, arakawa_lamb81, arakawa_lamb_string, ke_arakawa, ke_arakawa_string, ke_gudonov, ke_gudonov_string, ke_simple_gudonov, ke_simple_gudonov_string, mom_error_handler::mom_error(), mom_error_handler::mom_mesg(), pv_adv_centered, pv_adv_centered_string, pv_adv_upwind1, pv_adv_upwind1_string, mom_diag_mediator::register_diag_field(), robust_enstro, robust_enstro_string, sadourny75_energy, sadourny75_energy_string, sadourny75_enstro, sadourny75_enstro_string, and mom_string_functions::uppercase().

  type(time_type), target, intent(in)    :: time !< Current model time
  type(ocean_grid_type),   intent(in)    :: g !< Ocean grid structure
  type(param_file_type),   intent(in)    :: param_file !< Runtime parameter handles
  type(diag_ctrl), target, intent(inout) :: diag !< Diagnostics control structure
  type(accel_diag_ptrs),   target, intent(inout) :: ad !< Strorage for acceleration diagnostics
  type(coriolisadv_cs),    pointer       :: cs !< Control structure fro MOM_CoriolisAdv
  ! Local variables
! This include declares and sets the variable "version".
#include "version_variable.h"
  character(len=40)  :: mdl = "MOM_CoriolisAdv" ! This module's name.
  character(len=20)  :: tmpstr
  character(len=400) :: mesg
  integer :: isd, ied, jsd, jed, isdb, iedb, jsdb, jedb, nz

  isd = g%isd ; ied = g%ied ; jsd = g%jsd ; jed = g%jed ; nz = g%ke
  isdb = g%IsdB ; iedb = g%IedB ; jsdb = g%JsdB ; jedb = g%JedB

  if (associated(cs)) then
    call mom_error(warning, "CoriolisAdv_init called with associated control structure.")
    return
  endif
  allocate(cs)

  cs%diag => diag ; cs%Time => time

  ! Read all relevant parameters and write them to the model log.
  call log_version(param_file, mdl, version, "")
  call get_param(param_file, mdl, "NOSLIP", cs%no_slip, &
                 "If true, no slip boundary conditions are used; otherwise \n"//&
                 "free slip boundary conditions are assumed. The \n"//&
                 "implementation of the free slip BCs on a C-grid is much \n"//&
                 "cleaner than the no slip BCs. The use of free slip BCs \n"//&
                 "is strongly encouraged, and no slip BCs are not used with \n"//&
                 "the biharmonic viscosity.", default=.false.)

  call get_param(param_file, mdl, "CORIOLIS_EN_DIS", cs%Coriolis_En_Dis, &
                 "If true, two estimates of the thickness fluxes are used \n"//&
                 "to estimate the Coriolis term, and the one that \n"//&
                 "dissipates energy relative to the other one is used.", &
                 default=.false.)

  ! Set %Coriolis_Scheme
  ! (Select the baseline discretization for the Coriolis term)
  call get_param(param_file, mdl, "CORIOLIS_SCHEME", tmpstr, &
                 "CORIOLIS_SCHEME selects the discretization for the \n"//&
                 "Coriolis terms. Valid values are: \n"//&
                 "\t SADOURNY75_ENERGY - Sadourny, 1975; energy cons. \n"//&
                 "\t ARAKAWA_HSU90     - Arakawa & Hsu, 1990 \n"//&
                 "\t SADOURNY75_ENSTRO - Sadourny, 1975; enstrophy cons. \n"//&
                 "\t ARAKAWA_LAMB81    - Arakawa & Lamb, 1981; En. + Enst.\n"//&
                 "\t ARAKAWA_LAMB_BLEND - A blend of Arakawa & Lamb with \n"//&
                 "\t                      Arakawa & Hsu and Sadourny energy", &
                 default=sadourny75_energy_string)
  tmpstr = uppercase(tmpstr)
  select case (tmpstr)
    case (sadourny75_energy_string)
      cs%Coriolis_Scheme = sadourny75_energy
    case (arakawa_hsu_string)
      cs%Coriolis_Scheme = arakawa_hsu90
    case (sadourny75_enstro_string)
      cs%Coriolis_Scheme = sadourny75_enstro
    case (arakawa_lamb_string)
      cs%Coriolis_Scheme = arakawa_lamb81
    case (al_blend_string)
      cs%Coriolis_Scheme = al_blend
    case (robust_enstro_string)
      cs%Coriolis_Scheme = robust_enstro
      cs%Coriolis_En_Dis = .false.
    case default
      call mom_mesg('CoriolisAdv_init: Coriolis_Scheme ="'//trim(tmpstr)//'"', 0)
      call mom_error(fatal, "CoriolisAdv_init: Unrecognized setting "// &
            "#define CORIOLIS_SCHEME "//trim(tmpstr)//" found in input file.")
  end select
  if (cs%Coriolis_Scheme == al_blend) then
    call get_param(param_file, mdl, "CORIOLIS_BLEND_WT_LIN", cs%wt_lin_blend, &
                 "A weighting value for the ratio of inverse thicknesses, \n"//&
                 "beyond which the blending between Sadourny Energy and \n"//&
                 "Arakawa & Hsu goes linearly to 0 when CORIOLIS_SCHEME \n"//&
                 "is ARAWAKA_LAMB_BLEND. This must be between 1 and 1e-16.", &
                 units="nondim", default=0.125)
    call get_param(param_file, mdl, "CORIOLIS_BLEND_F_EFF_MAX", cs%F_eff_max_blend, &
                 "The factor by which the maximum effective Coriolis \n"//&
                 "acceleration from any point can be increased when \n"//&
                 "blending different discretizations with the \n"//&
                 "ARAKAWA_LAMB_BLEND Coriolis scheme.  This must be \n"//&
                 "greater than 2.0 (the max value for Sadourny energy).", &
                 units="nondim", default=4.0)
    cs%wt_lin_blend = min(1.0, max(cs%wt_lin_blend,1e-16))
    if (cs%F_eff_max_blend < 2.0) call mom_error(warning, "CoriolisAdv_init: "//&
           "CORIOLIS_BLEND_F_EFF_MAX should be at least 2.")
  endif

  mesg = "If true, the Coriolis terms at u-points are bounded by \n"//&
         "the four estimates of (f+rv)v from the four neighboring \n"//&
         "v-points, and similarly at v-points."
  if (cs%Coriolis_En_Dis .and. (cs%Coriolis_Scheme == sadourny75_energy)) then
    mesg = trim(mesg)//"  This option is \n"//&
                 "always effectively false with CORIOLIS_EN_DIS defined and \n"//&
                 "CORIOLIS_SCHEME set to "//trim(sadourny75_energy_string)//"."
  else
    mesg = trim(mesg)//"  This option would \n"//&
                 "have no effect on the SADOURNY Coriolis scheme if it \n"//&
                 "were possible to use centered difference thickness fluxes."
  endif
  call get_param(param_file, mdl, "BOUND_CORIOLIS", cs%bound_Coriolis, mesg, &
                 default=.false.)
  if ((cs%Coriolis_En_Dis .and. (cs%Coriolis_Scheme == sadourny75_energy)) .or. &
      (cs%Coriolis_Scheme == robust_enstro)) cs%bound_Coriolis = .false.

  ! Set KE_Scheme (selects discretization of KE)
  call get_param(param_file, mdl, "KE_SCHEME", tmpstr, &
                 "KE_SCHEME selects the discretization for acceleration \n"//&
                 "due to the kinetic energy gradient. Valid values are: \n"//&
                 "\t KE_ARAKAWA, KE_SIMPLE_GUDONOV, KE_GUDONOV", &
                 default=ke_arakawa_string)
  tmpstr = uppercase(tmpstr)
  select case (tmpstr)
    case (ke_arakawa_string); cs%KE_Scheme = ke_arakawa
    case (ke_simple_gudonov_string); cs%KE_Scheme = ke_simple_gudonov
    case (ke_gudonov_string); cs%KE_Scheme = ke_gudonov
    case default
      call mom_mesg('CoriolisAdv_init: KE_Scheme ="'//trim(tmpstr)//'"', 0)
      call mom_error(fatal, "CoriolisAdv_init: "// &
               "#define KE_SCHEME "//trim(tmpstr)//" in input file is invalid.")
  end select

  ! Set PV_Adv_Scheme (selects discretization of PV advection)
  call get_param(param_file, mdl, "PV_ADV_SCHEME", tmpstr, &
                 "PV_ADV_SCHEME selects the discretization for PV \n"//&
                 "advection. Valid values are: \n"//&
                 "\t PV_ADV_CENTERED - centered (aka Sadourny, 75) \n"//&
                 "\t PV_ADV_UPWIND1  - upwind, first order", &
                 default=pv_adv_centered_string)
  select case (uppercase(tmpstr))
    case (pv_adv_centered_string)
      cs%PV_Adv_Scheme = pv_adv_centered
    case (pv_adv_upwind1_string)
      cs%PV_Adv_Scheme = pv_adv_upwind1
    case default
      call mom_mesg('CoriolisAdv_init: PV_Adv_Scheme ="'//trim(tmpstr)//'"', 0)
      call mom_error(fatal, "CoriolisAdv_init: "// &
                     "#DEFINE PV_ADV_SCHEME in input file is invalid.")
  end select

  cs%id_rv = register_diag_field('ocean_model', 'RV', diag%axesBL, time, &
     'Relative Vorticity', 'second-1')

  cs%id_PV = register_diag_field('ocean_model', 'PV', diag%axesBL, time, &
     'Potential Vorticity', 'meter-1 second-1')

  cs%id_gKEu = register_diag_field('ocean_model', 'gKEu', diag%axesCuL, time, &
     'Zonal Acceleration from Grad. Kinetic Energy', 'meter-1 second-2')
  if (cs%id_gKEu > 0) call safe_alloc_ptr(ad%gradKEu,isdb,iedb,jsd,jed,nz)

  cs%id_gKEv = register_diag_field('ocean_model', 'gKEv', diag%axesCvL, time, &
     'Meridional Acceleration from Grad. Kinetic Energy', 'meter-1 second-2')
  if (cs%id_gKEv > 0) call safe_alloc_ptr(ad%gradKEv,isd,ied,jsdb,jedb,nz)

  cs%id_rvxu = register_diag_field('ocean_model', 'rvxu', diag%axesCvL, time, &
     'Meridional Acceleration from Relative Vorticity', 'meter-1 second-2')
  if (cs%id_rvxu > 0) call safe_alloc_ptr(ad%rv_x_u,isd,ied,jsdb,jedb,nz)

  cs%id_rvxv = register_diag_field('ocean_model', 'rvxv', diag%axesCuL, time, &
     'Zonal Acceleration from Relative Vorticity', 'meter-1 second-2')
  if (cs%id_rvxv > 0) call safe_alloc_ptr(ad%rv_x_v,isdb,iedb,jsd,jed,nz)

Here is the call graph for this function:

gradke()

subroutine mom_coriolisadv::gradke ( real, dimension( g %isdb: g %iedb, g %jsd: g %jed, g %ke), intent(in)  u, real, dimension( g %isd: g %ied, g %jsdb: g %jedb, g %ke), intent(in)  v, real, dimension( g %isd: g %ied, g %jsd: g %jed, g %ke), intent(in)  h, real, dimension( g %isd: g %ied , g %jsd: g %jed ), intent(out)  KE, real, dimension( g %isdb: g %iedb, g %jsd: g %jed ), intent(out)  KEx, real, dimension( g %isd: g %ied , g %jsdb: g %jedb), intent(out)  KEy, integer, intent(in)  k, type(ocean_obc_type), pointer  OBC, type(ocean_grid_type), intent(in)  G, type(coriolisadv_cs), pointer  CS  )

private

Calculates the acceleration due to the gradient of kinetic energy.

Parameters

[in]	g	Ocen grid structure
[in]	u	Zonal velocity (m/s)
[in]	v	Meridional velocity (m/s)
[in]	h	Layer thickness (m or kg/m2)
[out]	ke	Kinetic energy (m2/s2)
[out]	kex	Zonal acceleration due to kinetic energy gradient (m/s2)
[out]	key	Meridional acceleration due to kinetic energy gradient (m/s2)
[in]	k	Layer number to calculate for
	obc	Open boundary control structure
	cs	Control structure for MOM_CoriolisAdv

Definition at line 791 of file MOM_CoriolisAdv.F90.

References ke_arakawa, ke_gudonov, and ke_simple_gudonov.

Referenced by coradcalc().

  type(ocean_grid_type),                      intent(in)  :: g !< Ocen grid structure
  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)),  intent(in)  :: u !< Zonal velocity (m/s)
  real, dimension(SZI_(G),SZJB_(G),SZK_(G)),  intent(in)  :: v !< Meridional velocity (m/s)
  real, dimension(SZI_(G),SZJ_(G),SZK_(G)),   intent(in)  :: h !< Layer thickness (m or kg/m2)
  real, dimension(SZI_(G) ,SZJ_(G) ),         intent(out) :: ke !< Kinetic energy (m2/s2)
  real, dimension(SZIB_(G),SZJ_(G) ),         intent(out) :: kex !< Zonal acceleration due to kinetic
                                                                 !! energy gradient (m/s2)
  real, dimension(SZI_(G) ,SZJB_(G)),         intent(out) :: key !< Meridional acceleration due to kinetic
                                                                 !! energy gradient (m/s2)
  integer,                                    intent(in)  :: k !< Layer number to calculate for
  type(ocean_obc_type),                       pointer     :: obc !< Open boundary control structure
  type(coriolisadv_cs),                       pointer     :: cs !< Control structure for MOM_CoriolisAdv
  ! Local variables
  real :: um, up, vm, vp         ! Temporary variables with units of m s-1.
  real :: um2, up2, vm2, vp2     ! Temporary variables with units of m2 s-2.
  real :: um2a, up2a, vm2a, vp2a ! Temporary variables with units of m4 s-2.
  integer :: i, j, is, ie, js, je, isq, ieq, jsq, jeq, nz, n

  is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
  isq = g%IscB ; ieq = g%IecB ; jsq = g%JscB ; jeq = g%JecB


  ! Calculate KE (Kinetic energy for use in the -grad(KE) acceleration term).
  if (cs%KE_Scheme.eq.ke_arakawa) then
    ! The following calculation of Kinetic energy includes the metric terms
    ! identified in Arakawa & Lamb 1982 as important for KE conservation.  It
    ! also includes the possibility of partially-blocked tracer cell faces.
    do j=jsq,jeq+1 ; do i=isq,ieq+1
      ke(i,j) = ( ( g%areaCu( i ,j)*(u( i ,j,k)*u( i ,j,k))   &
                   +g%areaCu(i-1,j)*(u(i-1,j,k)*u(i-1,j,k)) ) &
                 +( g%areaCv(i, j )*(v(i, j ,k)*v(i, j ,k))   &
                   +g%areaCv(i,j-1)*(v(i,j-1,k)*v(i,j-1,k)) ) &
                )*0.25*g%IareaT(i,j)
    enddo ; enddo
  elseif (cs%KE_Scheme.eq.ke_simple_gudonov) then
    ! The following discretization of KE is based on the one-dimensinal Gudonov
    ! scheme which does not take into account any geometric factors
    do j=jsq,jeq+1 ; do i=isq,ieq+1
      up = 0.5*( u(i-1,j,k) + abs( u(i-1,j,k) ) ) ; up2 = up*up
      um = 0.5*( u( i ,j,k) - abs( u( i ,j,k) ) ) ; um2 = um*um
      vp = 0.5*( v(i,j-1,k) + abs( v(i,j-1,k) ) ) ; vp2 = vp*vp
      vm = 0.5*( v(i, j ,k) - abs( v(i, j ,k) ) ) ; vm2 = vm*vm
      ke(i,j) = ( max(up2,um2) + max(vp2,vm2) ) *0.5
    enddo ; enddo
  elseif (cs%KE_Scheme.eq.ke_gudonov) then
    ! The following discretization of KE is based on the one-dimensinal Gudonov
    ! scheme but has been adapted to take horizontal grid factors into account
    do j=jsq,jeq+1 ; do i=isq,ieq+1
      up = 0.5*( u(i-1,j,k) + abs( u(i-1,j,k) ) ) ; up2a = up*up*g%areaCu(i-1,j)
      um = 0.5*( u( i ,j,k) - abs( u( i ,j,k) ) ) ; um2a = um*um*g%areaCu( i ,j)
      vp = 0.5*( v(i,j-1,k) + abs( v(i,j-1,k) ) ) ; vp2a = vp*vp*g%areaCv(i,j-1)
      vm = 0.5*( v(i, j ,k) - abs( v(i, j ,k) ) ) ; vm2a = vm*vm*g%areaCv(i, j )
      ke(i,j) = ( max(um2a,up2a) + max(vm2a,vp2a) )*0.5*g%IareaT(i,j)
    enddo ; enddo
  endif

  ! Term - d(KE)/dx.
  do j=js,je ; do i=isq,ieq
    kex(i,j) = (ke(i+1,j) - ke(i,j)) * g%IdxCu(i,j)
  enddo ; enddo

  ! Term - d(KE)/dy.
  do j=jsq,jeq ; do i=is,ie
    key(i,j) = (ke(i,j+1) - ke(i,j)) * g%IdyCv(i,j)
  enddo ; enddo

  if (associated(obc)) then
    do n=1,obc%number_of_segments
      if (obc%segment(n)%is_N_or_S) then
        do i=obc%segment(n)%HI%isd,obc%segment(n)%HI%ied
          key(i,obc%segment(n)%HI%JsdB) = 0.
        enddo
      elseif (obc%segment(n)%is_E_or_W) then
        do j=obc%segment(n)%HI%jsd,obc%segment(n)%HI%jed
          kex(obc%segment(n)%HI%IsdB,j) = 0.
        enddo
      endif
    enddo
  endif

Here is the caller graph for this function:

Variable Documentation

al_blend

integer, parameter mom_coriolisadv::al_blend = 6

private

Definition at line 83 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: al_blend          = 6

al_blend_string

character*(20), parameter mom_coriolisadv::al_blend_string = "ARAKAWA_LAMB_BLEND"

private

Definition at line 89 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: AL_BLEND_STRING = "ARAKAWA_LAMB_BLEND"

arakawa_hsu90

integer, parameter mom_coriolisadv::arakawa_hsu90 = 2

private

Definition at line 79 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: arakawa_hsu90     = 2

arakawa_hsu_string

character*(20), parameter mom_coriolisadv::arakawa_hsu_string = "ARAKAWA_HSU90"

private

Definition at line 85 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: ARAKAWA_HSU_STRING = "ARAKAWA_HSU90"

arakawa_lamb81

integer, parameter mom_coriolisadv::arakawa_lamb81 = 5

private

Definition at line 82 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: arakawa_lamb81    = 5

arakawa_lamb_string

character*(20), parameter mom_coriolisadv::arakawa_lamb_string = "ARAKAWA_LAMB81"

private

Definition at line 88 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: ARAKAWA_LAMB_STRING = "ARAKAWA_LAMB81"

ke_arakawa

integer, parameter mom_coriolisadv::ke_arakawa = 10

private

Definition at line 91 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init(), and gradke().

integer, parameter :: ke_arakawa        = 10

ke_arakawa_string

character*(20), parameter mom_coriolisadv::ke_arakawa_string = "KE_ARAKAWA"

private

Definition at line 94 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: KE_ARAKAWA_STRING = "KE_ARAKAWA"

ke_gudonov

integer, parameter mom_coriolisadv::ke_gudonov = 12

private

Definition at line 93 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init(), and gradke().

integer, parameter :: ke_gudonov        = 12

ke_gudonov_string

character*(20), parameter mom_coriolisadv::ke_gudonov_string = "KE_GUDONOV"

private

Definition at line 96 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: KE_GUDONOV_STRING = "KE_GUDONOV"

ke_simple_gudonov

integer, parameter mom_coriolisadv::ke_simple_gudonov = 11

private

Definition at line 92 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init(), and gradke().

integer, parameter :: ke_simple_gudonov = 11

ke_simple_gudonov_string

character*(20), parameter mom_coriolisadv::ke_simple_gudonov_string = "KE_SIMPLE_GUDONOV"

private

Definition at line 95 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: KE_SIMPLE_GUDONOV_STRING = "KE_SIMPLE_GUDONOV"

pv_adv_centered

integer, parameter mom_coriolisadv::pv_adv_centered = 21

private

Definition at line 98 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: pv_adv_centered   = 21

pv_adv_centered_string

character*(20), parameter mom_coriolisadv::pv_adv_centered_string = "PV_ADV_CENTERED"

private

Definition at line 100 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: PV_ADV_CENTERED_STRING = "PV_ADV_CENTERED"

pv_adv_upwind1

integer, parameter mom_coriolisadv::pv_adv_upwind1 = 22

private

Definition at line 99 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: pv_adv_upwind1    = 22

pv_adv_upwind1_string

character*(20), parameter mom_coriolisadv::pv_adv_upwind1_string = "PV_ADV_UPWIND1"

private

Definition at line 101 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: PV_ADV_UPWIND1_STRING = "PV_ADV_UPWIND1"

robust_enstro

integer, parameter mom_coriolisadv::robust_enstro = 3

private

Definition at line 80 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: robust_enstro     = 3

robust_enstro_string

character*(20), parameter mom_coriolisadv::robust_enstro_string = "ROBUST_ENSTRO"

private

Definition at line 86 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: ROBUST_ENSTRO_STRING = "ROBUST_ENSTRO"

sadourny75_energy

integer, parameter mom_coriolisadv::sadourny75_energy = 1

Definition at line 78 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: sadourny75_energy = 1

sadourny75_energy_string

character*(20), parameter mom_coriolisadv::sadourny75_energy_string = "SADOURNY75_ENERGY"

private

Definition at line 84 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: SADOURNY75_ENERGY_STRING = "SADOURNY75_ENERGY"

sadourny75_enstro

integer, parameter mom_coriolisadv::sadourny75_enstro = 4

private

Definition at line 81 of file MOM_CoriolisAdv.F90.

Referenced by coradcalc(), and coriolisadv_init().

integer, parameter :: sadourny75_enstro = 4

sadourny75_enstro_string

character*(20), parameter mom_coriolisadv::sadourny75_enstro_string = "SADOURNY75_ENSTRO"

private

Definition at line 87 of file MOM_CoriolisAdv.F90.

Referenced by coriolisadv_init().

character*(20), parameter :: SADOURNY75_ENSTRO_STRING = "SADOURNY75_ENSTRO"