Detailed Description

Initialization of the boundary-forced-basing configuration.

Functions/Subroutines
subroutine, public	bfb_set_coord (Rlay, g_prime, GV, US, param_file, eqn_of_state)
	This subroutine specifies the vertical coordinate in terms of temperature at the surface and at the bottom. This case is set up in such a way that the temperature of the topmost layer is equal to the SST at the southern edge of the domain. The temperatures are then converted to densities of the top and bottom layers and linearly interpolated for the intermediate layers. More...

subroutine, public	bfb_initialize_sponges_southonly (G, GV, US, use_temperature, tv, param_file, CSp, h)
	This subroutine sets up the sponges for the southern bouundary of the domain. Maximum damping occurs within 2 degrees lat of the boundary. The damping linearly decreases northward over the next 2 degrees. More...

subroutine	write_bfb_log (param_file)
	Write output about the parameter values being used. More...

Variables
logical	first_call = .true.
	Unsafe model variable. More...

Function/Subroutine Documentation

◆ bfb_initialize_sponges_southonly()

subroutine, public bfb_initialization::bfb_initialize_sponges_southonly	(	type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US,
		logical, intent(in)	use_temperature,
		type(thermo_var_ptrs), intent(in)	tv,
		type(param_file_type), intent(in)	param_file,
		type(sponge_cs), pointer	CSp,
		real, dimension( g %isd: g %ied, g %jsd: g %jed, gv %ke), intent(in)	h
	)

This subroutine sets up the sponges for the southern bouundary of the domain. Maximum damping occurs within 2 degrees lat of the boundary. The damping linearly decreases northward over the next 2 degrees.

Parameters

[in]	g	The ocean's grid structure
[in]	gv	The ocean's vertical grid structure.
[in]	us	A dimensional unit scaling type
[in]	use_temperature	If true, temperature and salinity are used as state variables.
[in]	tv	A structure pointing to various thermodynamic variables
[in]	param_file	A structure to parse for run-time parameters
	csp	A pointer to the sponge control structure
[in]	h	Layer thicknesses [H ~> m or kg m-2]

Definition at line 79 of file BFB_initialization.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
   logical,                 intent(in) :: use_temperature !< If true, temperature and salinity are used as
                                             !! state variables.
   type(thermo_var_ptrs),   intent(in) :: tv   !< A structure pointing to various thermodynamic variables
   type(param_file_type),   intent(in) :: param_file !< A structure to parse for run-time parameters
   type(sponge_CS),         pointer    :: CSp  !< A pointer to the sponge control structure
   real, dimension(SZI_(G),SZJ_(G),SZK_(GV)), &
                            intent(in) :: h    !< Layer thicknesses [H ~> m or kg m-2]
  
   ! Local variables
   real :: eta(SZI_(G),SZJ_(G),SZK_(GV)+1) ! A temporary array for eta, in depth units [Z ~> m].
   real :: Idamp(SZI_(G),SZJ_(G))    ! The inverse damping rate [T-1 ~> s-1].
   real :: H0(SZK_(GV))              ! Resting layer thicknesses in depth units [Z ~> m].
   real :: min_depth                 ! The minimum ocean depth in depth units [Z ~> m].
   real :: slat, wlon, lenlat, lenlon, nlat
   real :: max_damping               ! The maximum damping rate [T-1 ~> s-1]
   character(len=40)  :: mdl = "BFB_initialize_sponges_southonly" ! This subroutine's name.
   integer :: i, j, k, is, ie, js, je, isd, ied, jsd, jed, nz
  
   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
  
   eta(:,:,:) = 0.0 ; idamp(:,:) = 0.0
  
 !  Here the inverse damping time [T-1 ~> s-1], is set. Set Idamp to 0
 !  wherever there is no sponge, and the subroutines that are called
 !  will automatically set up the sponges only where Idamp is positive
 !  and mask2dT is 1.
  
 !   Set up sponges for DOME configuration
   call get_param(param_file, mdl, "MINIMUM_DEPTH", min_depth, &
                  "The minimum depth of the ocean.", units="m", default=0.0, scale=us%m_to_Z)
  
   call get_param(param_file, mdl, "SOUTHLAT", slat, &
                  "The southern latitude of the domain.", units="degrees")
   call get_param(param_file, mdl, "LENLAT", lenlat, &
                  "The latitudinal length of the domain.", units="degrees")
   call get_param(param_file, mdl, "WESTLON", wlon, &
                  "The western longitude of the domain.", units="degrees", default=0.0)
   call get_param(param_file, mdl, "LENLON", lenlon, &
                  "The longitudinal length of the domain.", units="degrees")
   nlat = slat + lenlat
   do k=1,nz ; h0(k) = -g%max_depth * real(k-1) / real(nz) ; enddo
  
   ! Use for meridional thickness profile initialization
 !  do k=1,nz ; H0(k) = -G%max_depth * real(k-1) / real(nz-1) ; enddo
  
   max_damping = 1.0  / (86400.0*us%s_to_T)
  
   do i=is,ie; do j=js,je
     if (g%bathyT(i,j) <= min_depth) then ; idamp(i,j) = 0.0
     elseif (g%geoLatT(i,j) < slat+2.0) then ; idamp(i,j) = max_damping
     elseif (g%geoLatT(i,j) < slat+4.0) then
       idamp(i,j) = max_damping * (slat+4.0-g%geoLatT(i,j))/2.0
     else ; idamp(i,j) = 0.0
     endif
  
     ! These will be streched inside of apply_sponge, so they can be in
     ! depth space for Boussinesq or non-Boussinesq models.
  
     ! This section is used for uniform thickness initialization
     do k = 1,nz; eta(i,j,k) = h0(k); enddo
  
     ! The below section is used for meridional temperature profile thickness initiation
     ! do k = 1,nz; eta(i,j,k) = H0(k); enddo
     ! if (G%geoLatT(i,j) > 40.0) then
     !   do k = 1,nz
     !     eta(i,j,k) = -G%Angstrom_Z*(k-1)
     !   enddo
     ! elseif (G%geoLatT(i,j) > 20.0) then
     !   do k = 1,nz
     !     eta(i,j,k) = min(H0(k) + (G%geoLatT(i,j) - 20.0)*(G%max_depth - nz*G%Angstrom_Z)/20.0, &
     !                      -(k-1)*G%Angstrom_Z)
     !   enddo
     ! endif
     eta(i,j,nz+1) = -g%max_depth
  
   enddo ; enddo
  
 !  This call sets up the damping rates and interface heights.
 !  This sets the inverse damping timescale fields in the sponges.    !
   call initialize_sponge(idamp, eta, g, param_file, csp, gv)
  
 !   Now register all of the fields which are damped in the sponge.   !
 ! By default, momentum is advected vertically within the sponge, but !
 ! momentum is typically not damped within the sponge.                !
  
   if (first_call) call write_bfb_log(param_file)
  

◆ bfb_set_coord()

subroutine, public bfb_initialization::bfb_set_coord	(	real, dimension(gv%ke), intent(out)	Rlay,
		real, dimension(gv%ke+1), intent(out)	g_prime,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US,
		type(param_file_type), intent(in)	param_file,
		type(eos_type), pointer	eqn_of_state
	)

This subroutine specifies the vertical coordinate in terms of temperature at the surface and at the bottom. This case is set up in such a way that the temperature of the topmost layer is equal to the SST at the southern edge of the domain. The temperatures are then converted to densities of the top and bottom layers and linearly interpolated for the intermediate layers.

Parameters

[in]	gv	The ocean's vertical grid structure
[out]	rlay	Layer potential density [R ~> kg m-3].
[out]	g_prime	The reduced gravity at each interface [L2 Z-1 T-2 ~> m s-2].
[in]	us	A dimensional unit scaling type
[in]	param_file	A structure to parse for run-time parameters
	eqn_of_state	Equation of state structure

Definition at line 38 of file BFB_initialization.F90.

   type(verticalGrid_type),  intent(in)  :: GV      !< The ocean's vertical grid structure
   real, dimension(GV%ke),   intent(out) :: Rlay    !< Layer potential density [R ~> kg m-3].
   real, dimension(GV%ke+1), intent(out) :: g_prime !< The reduced gravity at each
                                                    !! interface [L2 Z-1 T-2 ~> m s-2].
   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(EOS_type),           pointer     :: eqn_of_state !< Equation of state structure
   ! Local variables
   real                                 :: drho_dt, SST_s, T_bot, rho_top, rho_bot
   integer                              :: k, nz
   character(len=40)  :: mdl = "BFB_set_coord" ! This subroutine's name.
  
   call get_param(param_file, mdl, "DRHO_DT", drho_dt, &
           "Rate of change of density with temperature.", &
            units="kg m-3 K-1", default=-0.2, scale=us%kg_m3_to_R)
   call get_param(param_file, mdl, "SST_S", sst_s, &
           "SST at the suothern edge of the domain.", units="C", default=20.0)
   call get_param(param_file, mdl, "T_BOT", t_bot, &
                  "Bottom Temp", units="C", default=5.0)
   rho_top = gv%Rho0 + drho_dt*sst_s
   rho_bot = gv%Rho0 + drho_dt*t_bot
   nz = gv%ke
  
   do k = 1,nz
     rlay(k) = (rho_bot - rho_top)/(nz-1)*real(k-1) + rho_top
     if (k >1) then
       g_prime(k) = (rlay(k) - rlay(k-1)) * gv%g_Earth / (gv%Rho0)
     else
       g_prime(k) = gv%g_Earth
     endif
     !Rlay(:) = 0.0
     !g_prime(:) = 0.0
   enddo
  
   if (first_call) call write_bfb_log(param_file)
  

◆ write_bfb_log()

subroutine bfb_initialization::write_bfb_log ( type(param_file_type), intent(in) param_file )

private

Write output about the parameter values being used.

Parameters

[in] param_file A structure indicating the open file to parse for model parameter values.

Definition at line 174 of file BFB_initialization.F90.

   type(param_file_type), intent(in) :: param_file !< A structure indicating the
                                                   !! open file to parse for model
                                                   !! parameter values.
  
 ! This include declares and sets the variable "version".
 #include "version_variable.h"
   character(len=40)  :: mdl = "BFB_initialization" ! This module's name.
  
   call log_version(param_file, mdl, version)
   first_call = .false.
  

Variable Documentation

◆ first_call

logical bfb_initialization::first_call = .true.

private

Unsafe model variable.

Todo:: Remove this module variable

Definition at line 30 of file BFB_initialization.F90.

30 logical :: first_call = .true.

Detailed Description

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ bfb_initialize_sponges_southonly()

◆ bfb_set_coord()

◆ write_bfb_log()

Variable Documentation

◆ first_call