Detailed Description

Calculates energy input to the internal tides.

Data Types
type	int_tide_input_cs
	This control structure holds parameters that regulate internal tide energy inputs. More...

type	int_tide_input_type
	This type is used to exchange fields related to the internal tides. More...

Functions/Subroutines
subroutine, public	set_int_tide_input (u, v, h, tv, fluxes, itide, dt, G, GV, US, CS)
	Sets the model-state dependent internal tide energy sources. More...

subroutine	find_n2_bottom (h, tv, T_f, S_f, h2, fluxes, G, GV, US, N2_bot)
	Estimates the near-bottom buoyancy frequency (N^2). More...

subroutine, public	int_tide_input_init (Time, G, GV, US, param_file, diag, CS, itide)
	Initializes the data related to the internal tide input module. More...

subroutine, public	int_tide_input_end (CS)
	Deallocates any memory related to the internal tide input module. More...

Function/Subroutine Documentation

◆ find_n2_bottom()

subroutine mom_int_tide_input::find_n2_bottom	(	real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	h,
		type(thermo_var_ptrs), intent(in)	tv,
		real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	T_f,
		real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	S_f,
		real, dimension(szi_(g),szj_(g)), intent(in)	h2,
		type(forcing), intent(in)	fluxes,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US,
		real, dimension( g %isd: g %ied, g %jsd: g %jed), intent(out)	N2_bot
	)

private

Estimates the near-bottom buoyancy frequency (N^2).

Parameters

[in]	g	The ocean's grid structure
[in]	gv	The ocean's vertical grid structure
[in]	us	A dimensional unit scaling type
[in]	h	Layer thicknesses [H ~> m or kg m-2]
[in]	tv	A structure containing pointers to the thermodynamic fields
[in]	t_f	Temperature after vertical filtering to smooth out the values in thin layers [degC].
[in]	s_f	Salinity after vertical filtering to smooth out the values in thin layers [ppt].
[in]	h2	Bottom topographic roughness [Z2 ~> m2].
[in]	fluxes	A structure of thermodynamic surface fluxes
[out]	n2_bot	The squared buoyancy freqency at the ocean bottom [T-2 ~> s-2].

Definition at line 151 of file MOM_internal_tide_input.F90.

   type(ocean_grid_type),                    intent(in)  :: G    !< The ocean's grid structure
   type(verticalGrid_type),                  intent(in)  :: GV   !< The ocean's vertical grid structure
   type(unit_scale_type),                    intent(in)  :: US   !< A dimensional unit scaling type
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), intent(in)  :: h    !< Layer thicknesses [H ~> m or kg m-2]
   type(thermo_var_ptrs),                    intent(in)  :: tv   !< A structure containing pointers to the
                                                                 !! thermodynamic fields
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), intent(in)  :: T_f  !< Temperature after vertical filtering to
                                                                 !! smooth out the values in thin layers [degC].
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), intent(in)  :: S_f  !< Salinity after vertical filtering to
                                                                 !! smooth out the values in thin layers [ppt].
   real, dimension(SZI_(G),SZJ_(G)),         intent(in)  :: h2   !< Bottom topographic roughness [Z2 ~> m2].
   type(forcing),                            intent(in)  :: fluxes !< A structure of thermodynamic surface fluxes
   type(int_tide_input_CS),                  pointer     :: CS    !<  This module's control structure.
   real, dimension(SZI_(G),SZJ_(G)),         intent(out) :: N2_bot !< The squared buoyancy freqency at the
                                                                  !! ocean bottom [T-2 ~> s-2].
   ! Local variables
   real, dimension(SZI_(G),SZK_(G)+1) :: &
     dRho_int      ! The unfiltered density differences across interfaces [R ~> kg m-3].
   real, dimension(SZI_(G)) :: &
     pres, &       ! The pressure at each interface [R L2 T-2 ~> Pa].
     Temp_int, &   ! The temperature at each interface [degC].
     Salin_int, &  ! The salinity at each interface [ppt].
     drho_bot, &   ! The density difference at the bottom of a layer [R ~> kg m-3]
     h_amp, &      ! The amplitude of topographic roughness [Z ~> m].
     hb, &         ! The depth below a layer [Z ~> m].
     z_from_bot, & ! The height of a layer center above the bottom [Z ~> m].
     dRho_dT, &    ! The partial derivative of density with temperature [R degC-1 ~> kg m-3 degC-1]
     dRho_dS       ! The partial derivative of density with salinity [R ppt-1 ~> kg m-3 ppt-1].
  
   real :: dz_int  ! The thickness associated with an interface [Z ~> m].
   real :: G_Rho0  ! The gravitation acceleration divided by the Boussinesq
                   ! density [Z T-2 R-1 ~> m4 s-2 kg-1].
   logical :: do_i(SZI_(G)), do_any
   integer, dimension(2) :: EOSdom ! The i-computational domain for the equation of state
   integer :: i, j, k, is, ie, js, je, nz
  
   is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
   g_rho0 = (us%L_to_Z**2*gv%g_Earth) / gv%Rho0
   eosdom(:) = eos_domain(g%HI)
  
   ! Find the (limited) density jump across each interface.
   do i=is,ie
     drho_int(i,1) = 0.0 ; drho_int(i,nz+1) = 0.0
   enddo
 !$OMP parallel do default(none) shared(is,ie,js,je,nz,tv,fluxes,G,GV,US,h,T_f,S_f, &
 !$OMP                                  h2,N2_bot,G_Rho0,EOSdom) &
 !$OMP                          private(pres,Temp_Int,Salin_Int,dRho_dT,dRho_dS, &
 !$OMP                                  hb,dRho_bot,z_from_bot,do_i,h_amp,       &
 !$OMP                                  do_any,dz_int) &
 !$OMP                     firstprivate(dRho_int)
   do j=js,je
     if (associated(tv%eqn_of_state)) then
       if (associated(fluxes%p_surf)) then
         do i=is,ie ; pres(i) = fluxes%p_surf(i,j) ; enddo
       else
         do i=is,ie ; pres(i) = 0.0 ; enddo
       endif
       do k=2,nz
         do i=is,ie
           pres(i) = pres(i) + (gv%g_Earth*gv%H_to_RZ)*h(i,j,k-1)
           temp_int(i) = 0.5 * (t_f(i,j,k) + t_f(i,j,k-1))
           salin_int(i) = 0.5 * (s_f(i,j,k) + s_f(i,j,k-1))
         enddo
         call calculate_density_derivs(temp_int, salin_int, pres, drho_dt(:), drho_ds(:), &
                                       tv%eqn_of_state, eosdom)
         do i=is,ie
           drho_int(i,k) = max(drho_dt(i)*(t_f(i,j,k) - t_f(i,j,k-1)) + &
                               drho_ds(i)*(s_f(i,j,k) - s_f(i,j,k-1)), 0.0)
         enddo
       enddo
     else
       do k=2,nz ; do i=is,ie
         drho_int(i,k) = (gv%Rlay(k) - gv%Rlay(k-1))
       enddo ; enddo
     endif
  
     ! Find the bottom boundary layer stratification.
     do i=is,ie
       hb(i) = 0.0 ; drho_bot(i) = 0.0
       z_from_bot(i) = 0.5*gv%H_to_Z*h(i,j,nz)
       do_i(i) = (g%mask2dT(i,j) > 0.5)
       h_amp(i) = sqrt(h2(i,j))
     enddo
  
     do k=nz,2,-1
       do_any = .false.
       do i=is,ie ; if (do_i(i)) then
         dz_int = 0.5*gv%H_to_Z*(h(i,j,k) + h(i,j,k-1))
         z_from_bot(i) = z_from_bot(i) + dz_int ! middle of the layer above
  
         hb(i) = hb(i) + dz_int
         drho_bot(i) = drho_bot(i) + drho_int(i,k)
  
         if (z_from_bot(i) > h_amp(i)) then
           if (k>2) then
             ! Always include at least one full layer.
             hb(i) = hb(i) + 0.5*gv%H_to_Z*(h(i,j,k-1) + h(i,j,k-2))
             drho_bot(i) = drho_bot(i) + drho_int(i,k-1)
           endif
           do_i(i) = .false.
         else
           do_any = .true.
         endif
       endif ; enddo
       if (.not.do_any) exit
     enddo
  
     do i=is,ie
       if (hb(i) > 0.0) then
         n2_bot(i,j) = (g_rho0 * drho_bot(i)) / hb(i)
       else ;  n2_bot(i,j) = 0.0 ; endif
     enddo
   enddo
  

◆ int_tide_input_end()

subroutine, public mom_int_tide_input::int_tide_input_end ( type(int_tide_input_cs), pointer CS )

Deallocates any memory related to the internal tide input module.

Parameters

cs	This module's control structure, which is deallocated here.

Definition at line 426 of file MOM_internal_tide_input.F90.

   type(int_tide_input_CS), pointer :: CS !< This module's control structure, which is deallocated here.
  
   if (associated(cs)) deallocate(cs)

◆ int_tide_input_init()

subroutine, public mom_int_tide_input::int_tide_input_init	(	type(time_type), intent(in)	Time,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US,
		type(param_file_type), intent(in)	param_file,
		type(diag_ctrl), intent(inout), target	diag,
		type(int_tide_input_cs), pointer	CS,
		type(int_tide_input_type), pointer	itide
	)

Initializes the data related to the internal tide input module.

Parameters

[in]	time	The current model time
[in]	g	The ocean's grid structure
[in]	gv	The ocean's vertical grid structure
[in]	us	A dimensional unit scaling type
[in]	param_file	A structure to parse for run-time parameters
[in,out]	diag	structure used to regulate diagnostic output.
	cs	This module's control structure, which is initialized here.
	itide	A structure containing fields related to the internal tide sources.

Definition at line 269 of file MOM_internal_tide_input.F90.

   type(time_type),           intent(in)    :: Time !< The current model time
   type(ocean_grid_type),     intent(in)    :: G    !< The ocean's grid structure
   type(verticalGrid_type),   intent(in)    :: GV   !< The ocean's vertical grid structure
   type(unit_scale_type),     intent(in)    :: US   !< A dimensional unit scaling type
   type(param_file_type),     intent(in)    :: param_file !< A structure to parse for run-time parameters
   type(diag_ctrl),   target, intent(inout) :: diag !< structure used to regulate diagnostic output.
   type(int_tide_input_CS),   pointer       :: CS   !< This module's control structure, which is initialized here.
   type(int_tide_input_type), pointer       :: itide !< A structure containing fields related
                                                    !! to the internal tide sources.
   ! Local variables
   type(vardesc) :: vd
   logical :: read_tideamp
   ! This include declares and sets the variable "version".
 # include "version_variable.h"
   character(len=40)  :: mdl = "MOM_int_tide_input"  ! This module's name.
   character(len=20)  :: tmpstr
   character(len=200) :: filename, tideamp_file, h2_file
  
   real :: mask_itidal        ! A multiplicative land mask, 0 or 1 [nondim]
   real :: max_frac_rough     ! The fraction relating the maximum topographic roughness
                              ! to the mean depth [nondim]
   real :: utide              ! constant tidal amplitude [L T-1 ~> m s-1] to be used if
                              ! tidal amplitude file is not present.
   real :: kappa_h2_factor    ! factor for the product of wavenumber * rms sgs height [nondim].
   real :: kappa_itides       ! topographic wavenumber and non-dimensional scaling [L-1 ~> m-1]
   real :: min_zbot_itides    ! Minimum ocean depth for internal tide conversion [Z ~> m].
   integer :: tlen_days       !< Time interval from start for adding wave source
                              !! for testing internal tides (BDM)
   integer :: i, j, is, ie, js, je, isd, ied, jsd, jed
  
   if (associated(cs)) then
     call mom_error(warning, "int_tide_input_init called with an associated "// &
                             "control structure.")
     return
   endif
   if (associated(itide)) then
     call mom_error(warning, "int_tide_input_init called with an associated "// &
                             "internal tide input type.")
     return
   endif
   allocate(cs)
   allocate(itide)
  
   is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec
   isd = g%isd ; ied = g%ied ; jsd = g%jsd ; jed = g%jed
  
   cs%diag => diag
  
   ! Read all relevant parameters and write them to the model log.
   call log_version(param_file, mdl, version, "")
  
   call get_param(param_file, mdl, "INPUTDIR", cs%inputdir, default=".")
   cs%inputdir = slasher(cs%inputdir)
  
   call get_param(param_file, mdl, "DEBUG", cs%debug, default=.false., do_not_log=.true.)
  
   call get_param(param_file, mdl, "MIN_ZBOT_ITIDES", min_zbot_itides, &
                "Turn off internal tidal dissipation when the total "//&
                "ocean depth is less than this value.", units="m", default=0.0, scale=us%m_to_Z)
   call get_param(param_file, mdl, "KD_SMOOTH", cs%kappa_fill, &
                  "A diapycnal diffusivity that is used to interpolate "//&
                  "more sensible values of T & S into thin layers.", &
                  units="m2 s-1", default=1.0e-6, scale=us%m2_s_to_Z2_T)
  
   call get_param(param_file, mdl, "UTIDE", utide, &
                "The constant tidal amplitude used with INT_TIDE_DISSIPATION.", &
                units="m s-1", default=0.0, scale=us%m_s_to_L_T)
  
   allocate(itide%Nb(isd:ied,jsd:jed))  ; itide%Nb(:,:) = 0.0
   allocate(itide%h2(isd:ied,jsd:jed))  ; itide%h2(:,:) = 0.0
   allocate(itide%TKE_itidal_input(isd:ied,jsd:jed)) ; itide%TKE_itidal_input(:,:) = 0.0
   allocate(itide%tideamp(isd:ied,jsd:jed)) ; itide%tideamp(:,:) = utide
   allocate(cs%TKE_itidal_coef(isd:ied,jsd:jed)) ; cs%TKE_itidal_coef(:,:) = 0.0
  
   call get_param(param_file, mdl, "KAPPA_ITIDES", kappa_itides, &
                "A topographic wavenumber used with INT_TIDE_DISSIPATION. "//&
                "The default is 2pi/10 km, as in St.Laurent et al. 2002.", &
                units="m-1", default=8.e-4*atan(1.0), scale=us%L_to_m)
  
   call get_param(param_file, mdl, "KAPPA_H2_FACTOR", kappa_h2_factor, &
                "A scaling factor for the roughness amplitude with n"//&
                "INT_TIDE_DISSIPATION.",  units="nondim", default=1.0)
   call get_param(param_file, mdl, "TKE_ITIDE_MAX", cs%TKE_itide_max, &
                "The maximum internal tide energy source available to mix "//&
                "above the bottom boundary layer with INT_TIDE_DISSIPATION.", &
                units="W m-2", default=1.0e3, scale=us%W_m2_to_RZ3_T3)
  
   call get_param(param_file, mdl, "READ_TIDEAMP", read_tideamp, &
                "If true, read a file (given by TIDEAMP_FILE) containing "//&
                "the tidal amplitude with INT_TIDE_DISSIPATION.", default=.false.)
   if (read_tideamp) then
     call get_param(param_file, mdl, "TIDEAMP_FILE", tideamp_file, &
                "The path to the file containing the spatially varying "//&
                "tidal amplitudes with INT_TIDE_DISSIPATION.", default="tideamp.nc")
     filename = trim(cs%inputdir) // trim(tideamp_file)
     call log_param(param_file, mdl, "INPUTDIR/TIDEAMP_FILE", filename)
     call mom_read_data(filename, 'tideamp', itide%tideamp, g%domain, timelevel=1, scale=us%m_s_to_L_T)
   endif
  
   call get_param(param_file, mdl, "H2_FILE", h2_file, &
                "The path to the file containing the sub-grid-scale "//&
                "topographic roughness amplitude with INT_TIDE_DISSIPATION.", &
                fail_if_missing=.true.)
   filename = trim(cs%inputdir) // trim(h2_file)
   call log_param(param_file, mdl, "INPUTDIR/H2_FILE", filename)
   call mom_read_data(filename, 'h2', itide%h2, g%domain, timelevel=1, scale=us%m_to_Z**2)
  
   call get_param(param_file, mdl, "FRACTIONAL_ROUGHNESS_MAX", max_frac_rough, &
                  "The maximum topographic roughness amplitude as a fraction of the mean depth, "//&
                  "or a negative value for no limitations on roughness.", &
                  units="nondim", default=0.1)
  
   ! The following parameters are used in testing the internal tide code.
   call get_param(param_file, mdl, "INTERNAL_TIDE_SOURCE_TEST", cs%int_tide_source_test, &
                  "If true, apply an arbitrary generation site for internal tide testing", &
                  default=.false.)
   if (cs%int_tide_source_test)then
     call get_param(param_file, mdl, "INTERNAL_TIDE_SOURCE_X", cs%int_tide_source_x, &
                  "X Location of generation site for internal tide", default=1.)
     call get_param(param_file, mdl, "INTERNAL_TIDE_SOURCE_Y", cs%int_tide_source_y, &
                  "Y Location of generation site for internal tide", default=1.)
     call get_param(param_file, mdl, "INTERNAL_TIDE_SOURCE_TLEN_DAYS", tlen_days, &
                  "Time interval from start of experiment for adding wave source", &
                  units="days", default=0)
     cs%time_max_source = time + set_time(0, days=tlen_days)
   endif
  
   do j=js,je ; do i=is,ie
     mask_itidal = 1.0
     if (g%bathyT(i,j) < min_zbot_itides) mask_itidal = 0.0
  
     itide%tideamp(i,j) = itide%tideamp(i,j) * mask_itidal * g%mask2dT(i,j)
  
     ! Restrict rms topo to a fraction (often 10 percent) of the column depth.
     if (max_frac_rough >= 0.0) &
       itide%h2(i,j) = min((max_frac_rough*g%bathyT(i,j))**2, itide%h2(i,j))
  
     ! Compute the fixed part of internal tidal forcing; units are [R Z3 T-2 ~> J m-2] here.
     cs%TKE_itidal_coef(i,j) = 0.5*us%L_to_Z*kappa_h2_factor*gv%Rho0*&
          kappa_itides * itide%h2(i,j) * itide%tideamp(i,j)**2
   enddo ; enddo
  
  
   cs%id_TKE_itidal = register_diag_field('ocean_model','TKE_itidal_itide',diag%axesT1,time, &
       'Internal Tide Driven Turbulent Kinetic Energy', &
       'W m-2', conversion=us%RZ3_T3_to_W_m2)
  
   cs%id_Nb = register_diag_field('ocean_model','Nb_itide',diag%axesT1,time, &
        'Bottom Buoyancy Frequency', 's-1', conversion=us%s_to_T)
  
   cs%id_N2_bot = register_diag_field('ocean_model','N2_b_itide',diag%axesT1,time, &
        'Bottom Buoyancy frequency squared', 's-2', conversion=us%s_to_T**2)
  

◆ set_int_tide_input()

subroutine, public mom_int_tide_input::set_int_tide_input	(	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,
		type(thermo_var_ptrs), intent(in)	tv,
		type(forcing), intent(in)	fluxes,
		type(int_tide_input_type), intent(inout)	itide,
		real, intent(in)	dt,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US,
		type(int_tide_input_cs), pointer	CS
	)

Sets the model-state dependent internal tide energy sources.

Parameters

[in]	g	The ocean's grid structure
[in]	gv	The ocean's vertical grid structure
[in]	us	A dimensional unit scaling type
[in]	u	The zonal velocity [L T-1 ~> m s-1]
[in]	v	The meridional velocity [L T-1 ~> m s-1]
[in]	h	Layer thicknesses [H ~> m or kg m-2]
[in]	tv	A structure containing pointers to the thermodynamic fields
[in]	fluxes	A structure of thermodynamic surface fluxes
[in,out]	itide	A structure containing fields related to the internal tide sources.
[in]	dt	The time increment [T ~> s].
	cs	This module's control structure.

Definition at line 75 of file MOM_internal_tide_input.F90.

   type(ocean_grid_type),                     intent(in)    :: G  !< The ocean's grid structure
   type(verticalGrid_type),                   intent(in)    :: GV !< The ocean's vertical grid structure
   type(unit_scale_type),                     intent(in)    :: US !< A dimensional unit scaling type
   real, dimension(SZIB_(G),SZJ_(G),SZK_(G)), intent(in)    :: u  !< The zonal velocity [L T-1 ~> m s-1]
   real, dimension(SZI_(G),SZJB_(G),SZK_(G)), intent(in)    :: v  !< The meridional velocity [L T-1 ~> m s-1]
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)),  intent(in)    :: h  !< Layer thicknesses [H ~> m or kg m-2]
   type(thermo_var_ptrs),                     intent(in)    :: tv !< A structure containing pointers to the
                                                                  !! thermodynamic fields
   type(forcing),                             intent(in)    :: fluxes !< A structure of thermodynamic surface fluxes
   type(int_tide_input_type),                 intent(inout) :: itide !< A structure containing fields related
                                                                  !! to the internal tide sources.
   real,                                      intent(in)    :: dt !< The time increment [T ~> s].
   type(int_tide_input_CS),                   pointer       :: CS !< This module's control structure.
   ! Local variables
   real, dimension(SZI_(G),SZJ_(G)) :: &
     N2_bot        ! The bottom squared buoyancy frequency [T-2 ~> s-2].
  
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)) :: &
     T_f, S_f      ! The temperature and salinity in [degC] and [ppt] with the values in
                   ! the massless layers filled vertically by diffusion.
   logical :: use_EOS    ! If true, density is calculated from T & S using an
                         ! equation of state.
   logical :: avg_enabled  ! for testing internal tides (BDM)
   type(time_type) :: time_end        !< For use in testing internal tides (BDM)
  
   integer :: i, j, k, is, ie, js, je, nz, isd, ied, jsd, jed
  
   is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
   isd = g%isd ; ied = g%ied ; jsd = g%jsd ; jed = g%jed
  
   if (.not.associated(cs)) call mom_error(fatal,"set_diffusivity: "//&
          "Module must be initialized before it is used.")
  
   use_eos = associated(tv%eqn_of_state)
  
   ! Smooth the properties through massless layers.
   if (use_eos) then
     call vert_fill_ts(h, tv%T, tv%S, cs%kappa_fill*dt, t_f, s_f, g, gv, larger_h_denom=.true.)
   endif
  
   call find_n2_bottom(h, tv, t_f, s_f, itide%h2, fluxes, g, gv, us, n2_bot)
  
   !$OMP parallel do default(shared)
   do j=js,je ; do i=is,ie
     itide%Nb(i,j) = g%mask2dT(i,j) * sqrt(n2_bot(i,j))
     itide%TKE_itidal_input(i,j) = min(cs%TKE_itidal_coef(i,j)*itide%Nb(i,j), cs%TKE_itide_max)
   enddo ; enddo
  
   if (cs%int_tide_source_test) then
     itide%TKE_itidal_input(:,:) = 0.0
     avg_enabled = query_averaging_enabled(cs%diag, time_end=time_end)
     if (time_end <= cs%time_max_source) then
       do j=js,je ; do i=is,ie
         ! Input  an arbitrary energy point source.id_
         if (((g%geoLonCu(i-1,j)-cs%int_tide_source_x) * (g%geoLonBu(i,j)-cs%int_tide_source_x) <= 0.0) .and. &
             ((g%geoLatCv(i,j-1)-cs%int_tide_source_y) * (g%geoLatCv(i,j)-cs%int_tide_source_y) <= 0.0)) then
           itide%TKE_itidal_input(i,j) = 1.0*us%kg_m3_to_R*us%m_to_Z**3*us%T_to_s**3
         endif
       enddo ; enddo
     endif
   endif
  
   if (cs%debug) then
     call hchksum(n2_bot,"N2_bot",g%HI,haloshift=0, scale=us%s_to_T**2)
     call hchksum(itide%TKE_itidal_input,"TKE_itidal_input",g%HI,haloshift=0, &
                  scale=us%RZ3_T3_to_W_m2)
   endif
  
   if (cs%id_TKE_itidal > 0) call post_data(cs%id_TKE_itidal, itide%TKE_itidal_input, cs%diag)
   if (cs%id_Nb > 0) call post_data(cs%id_Nb, itide%Nb, cs%diag)
   if (cs%id_N2_bot > 0 ) call post_data(cs%id_N2_bot, n2_bot, cs%diag)
  

Detailed Description

Data Types

Functions/Subroutines

Function/Subroutine Documentation

◆ find_n2_bottom()

◆ int_tide_input_end()

◆ int_tide_input_init()

◆ set_int_tide_input()