ocean_munk_delta_m Function

public pure function ocean_munk_delta_m(nu_h, beta) result(delta_m)

Munk boundary-layer width delta_M = (nu_h/beta)^(1/3). beta<=0 (f-plane — no meridional PV gradient, no Munk boundary layer) or nu_h<=0 returns huge(1.0_wp) (no constraint — never trips the >= 2-cell criterion).

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: nu_h
real(kind=wp), intent(in) :: beta

Return Value real(kind=wp)


Called by

proc~~ocean_munk_delta_m~~CalledByGraph proc~ocean_munk_delta_m ocean_munk_delta_m proc~ocean_munk_worst_case ocean_munk_worst_case proc~ocean_munk_worst_case->proc~ocean_munk_delta_m proc~ocean_stability_audit ocean_stability_audit proc~ocean_stability_audit->proc~ocean_munk_delta_m proc~ocean_stability_audit->proc~ocean_munk_worst_case proc~engine_setup engine_setup proc~engine_setup->proc~ocean_stability_audit proc~complete_ocean_create complete_ocean_create proc~complete_ocean_create->proc~engine_setup proc~driver_run_ocean driver_run_ocean proc~driver_run_ocean->proc~engine_setup proc~driver_validate driver_validate proc~driver_validate->proc~engine_setup proc~driver_run driver_run proc~driver_run->proc~driver_run_ocean proc~rdb_ocean_create_finalize rdb_ocean_create_finalize proc~rdb_ocean_create_finalize->proc~complete_ocean_create proc~rdb_ocean_create_from_string rdb_ocean_create_from_string proc~rdb_ocean_create_from_string->proc~complete_ocean_create

Source Code

   pure function ocean_munk_delta_m(nu_h, beta) result(delta_m)
      !! Munk boundary-layer width `delta_M = (nu_h/beta)^(1/3)`.
      !! `beta<=0` (f-plane — no meridional PV gradient, no Munk
      !! boundary layer) or `nu_h<=0` returns `huge(1.0_wp)` (no
      !! constraint — never trips the >= 2-cell criterion).
      real(wp), intent(in) :: nu_h, beta
      real(wp) :: delta_m
      if (beta > 0.0_wp .and. nu_h > 0.0_wp) then
         delta_m = (nu_h/beta)**(1.0_wp/3.0_wp)
      else
         delta_m = huge(1.0_wp)
      end if
   end function ocean_munk_delta_m