mom_shortwave_abs Module Reference

Data Types
type	optics_type

Functions/Subroutines
subroutine, public	absorbremainingsw (G, GV, h, opacity_band, nsw, j, dt, H_limit_fluxes, adjustAbsorptionProfile, absorbAllSW, T, Pen_SW_bnd, eps, ksort, htot, Ttot, TKE, dSV_dT)
	Apply shortwave heating below surface boundary layer. More...
subroutine, public	sumswoverbands (G, GV, h, opacity_band, nsw, j, dt, H_limit_fluxes, absorbAllSW, iPen_SW_bnd, netPen)

Function/Subroutine Documentation

absorbremainingsw()

subroutine, public mom_shortwave_abs::absorbremainingsw	(	type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		real, dimension( g %isd: g %ied, g %ke), intent(in)	h,
		real, dimension(:,:,:), intent(in)	opacity_band,
		integer, intent(in)	nsw,
		integer, intent(in)	j,
		real, intent(in)	dt,
		real, intent(in)	H_limit_fluxes,
		logical, intent(in)	adjustAbsorptionProfile,
		logical, intent(in)	absorbAllSW,
		real, dimension( g %isd: g %ied, g %ke), intent(inout)	T,
		real, dimension(:,:), intent(inout)	Pen_SW_bnd,
		real, dimension( g %isd: g %ied, g %ke), intent(in), optional	eps,
		integer, dimension( g %isd: g %ied, g %ke), intent(in), optional	ksort,
		real, dimension( g %isd: g %ied), intent(in), optional	htot,
		real, dimension( g %isd: g %ied), intent(inout), optional	Ttot,
		real, dimension( g %isd: g %ied, g %ke), intent(inout), optional	TKE,
		real, dimension( g %isd: g %ied, g %ke), intent(in), optional	dSV_dT
	)

Apply shortwave heating below surface boundary layer.

Parameters

[in]	g	The ocean's grid structure.
[in]	gv	The ocean's vertical grid structure.
[in]	h	Layer thicknesses, in H (usually m or kg m-2).
[in]	opacity_band	Opacity in each band of penetrating shortwave radiation (1/H). The indicies are band, i, k.
[in]	nsw	Number of bands of penetrating shortwave radiation.
[in]	j	j-index to work on.
[in]	dt	Time step (seconds).
[in]	h_limit_fluxes	If the total ocean depth is less than this, they are scaled away to avoid numerical instabilities. (H) This would not be necessary if a finite heat capacity mud-layer were added.
[in]	adjustabsorptionprofile	If true, apply heating above the layers in which it should have occurred to get the correct mean depth (and potential energy change) of the shortwave that should be absorbed by each layer.
[in]	absorballsw	If true, apply heating above the layers in which it should have occurred to get the correct mean depth (and potential energy change) of the shortwave that should be absorbed by each layer.
[in,out]	t	Layer potential/conservative temperatures (deg C)
[in,out]	pen_sw_bnd	Penetrating shortwave heating in each band that hits the bottom and will be redistributed through the water column (units of K*H), size nsw x G isd: G ied.
[in]	eps	Small thickness that must remain in
[in]	ksort	Density-sorted k-indicies.
[in]	htot	Total mixed layer thickness, in H .
[in,out]	ttot	Depth integrated mixed layer
[in]	dsv_dt	The partial derivative of specific
[in,out]	tke	The TKE sink from mixing the heating

Definition at line 60 of file MOM_shortwave_abs.F90.

Referenced by mom_diabatic_aux::applyboundaryfluxesinout(), and mom_bulk_mixed_layer::bulkmixedlayer().

!< This subroutine applies shortwave heating below the boundary layer (when running
!! with the bulk mixed layer from GOLD) or throughout the water column.  In
!! addition, it causes all of the remaining SW radiation to be absorbed,
!! provided that the total water column thickness is greater than
!! H_limit_fluxes.  For thinner water columns, the heating is scaled down
!! proportionately, the assumption being that the remaining heating (which is
!! left in Pen_SW) should go into an (absent for now) ocean bottom sediment layer.

  type(ocean_grid_type),            intent(in)    :: g     !< The ocean's grid structure.
  type(verticalgrid_type),          intent(in)    :: gv    !< The ocean's vertical grid structure.
  real, dimension(SZI_(G),SZK_(G)), intent(in)    :: h     !< Layer thicknesses, in H (usually m or
                                                           !! kg m-2).
  real, dimension(:,:,:),           intent(in)    :: opacity_band !< Opacity in each band of
                                                           !! penetrating shortwave radiation (1/H).
                                                           !! The indicies are band, i, k.
  integer,                          intent(in)    :: nsw   !< Number of bands of penetrating
                                                           !! shortwave radiation.
  integer,                          intent(in)    :: j     !< j-index to work on.
  real,                             intent(in)    :: dt    !< Time step (seconds).
  real,                             intent(in)    :: h_limit_fluxes !< If the total ocean depth is
                                                           !! less than this, they are scaled away
                                                           !! to avoid numerical instabilities. (H)
                                                           !! This would not be necessary if a
                                                           !! finite heat capacity mud-layer
                                                           !! were added.
  logical,                          intent(in)    :: adjustabsorptionprofile !< If true, apply
                                                           !! heating above the layers in which it
                                                           !! should have occurred to get the
                                                           !! correct mean depth (and potential
                                                           !! energy change) of the shortwave that
                                                           !! should be absorbed by each layer.
  logical,                          intent(in)    :: absorballsw !< If true, apply heating above the
                                                           !! layers in which it should have occurred
                                                           !! to get the correct mean depth (and
                                                           !! potential energy change) of the
                                                           !! shortwave that should be absorbed by
                                                           !! each layer.
  real, dimension(SZI_(G),SZK_(G)), intent(inout) :: t     !< Layer potential/conservative
                                                           !! temperatures (deg C)
  real, dimension(:,:),             intent(inout) :: pen_sw_bnd !< Penetrating shortwave heating in
                                                           !! each band that hits the bottom and
                                                           !! will be redistributed through the
                                                           !! water column (units of K*H), size
                                                           !! nsw x SZI_(G).
  real, dimension(SZI_(G),SZK_(G)),    &
                          optional, intent(in)    :: eps   !< Small thickness that must remain in
                                                           !! each layer, and which will not be
                                                           !! subject to heating (units of H)
  integer, dimension(SZI_(G),SZK_(G)), &
                          optional, intent(in)    :: ksort !< Density-sorted k-indicies.
  real, dimension(SZI_(G)),            &
                          optional, intent(in)    :: htot  !< Total mixed layer thickness, in H .
  real, dimension(SZI_(G)),            &
                          optional, intent(inout) :: ttot  !< Depth integrated mixed layer
                                                           !! temperature (units of K H).
  real, dimension(SZI_(G),SZK_(G)),    &
                          optional, intent(in)    :: dsv_dt !< The partial derivative of specific
                                                           !! volume with temperature, in m3 kg-1
                                                           !! K-1.
  real, dimension(SZI_(G),SZK_(G)),    &
                          optional, intent(inout) :: tke   !< The TKE sink from mixing the heating
                                                           !! throughout a layer, in J m-2.

! Arguments:
!  (in)    G            = the ocean grid structure.
!  (in)    GV           = The ocean's vertical grid structure.
!  (in)    h            = the layer thicknesses, in m or kg m-2.
!                         units of h are referred to as "H" below.
!  (in)    opacity_band = opacity in each band of penetrating shortwave
!                         radiation (1/H). The indicies are band, i, k.
!  (in)    nsw          = number of bands of penetrating shortwave radiation
!  (in)    j            = j-index to work on
!  (in)    dt           = time step (seconds)
!  (in)    H_limit_fluxes = if the total ocean depth is less than this, they
!                         are scaled away to avoid numerical instabilities. (H)
!                         This would not be necessary if a finite heat
!                         capacity mud-layer were added.
!  (in)    adjustAbsorptionProfile = if true, apply heating above the layers
!                         in which it should have occurred to get the correct
!                         mean depth (and potential energy change) of the
!                         shortwave that should be absorbed by each layer.
!  (in)    absorbAllSW  = if true, any shortwave radiation that hits the
!                         bottom is absorbed uniformly over the water column.
!  (inout) T            = layer potential/conservative temperatures (deg C)
!  (inout) Pen_SW_bnd   = penetrating shortwave heating in each band that
!                         hits the bottom and will be redistributed through
!                         the water column (units of K*H), size nsw x SZI_(G).

! These optional arguments apply when the bulk mixed layer is used
! but are unnecessary with other schemes.
!  (in,opt)    eps      = small thickness that must remain in each layer, and
!                         which will not be subject to heating (units of H)
!  (inout,opt) ksort    = density-sorted k-indicies
!  (in,opt)    htot     = total mixed layer thickness, in H
!  (inout,opt) Ttot     = depth integrated mixed layer temperature (units of K H)
!  (in,opt)    dSV_dT   = the partial derivative of specific volume with temperature, in m3 kg-1 K-1.
!  (inout,opt) TKE      = the TKE sink from mixing the heating throughout a layer, in J m-2.

  real, dimension(SZI_(G),SZK_(G)) :: &
    t_chg_above    ! A temperature change that will be applied to all the thick
                   ! layers above a given layer, in K.  This is only nonzero if
                   ! adjustAbsorptionProfile is true, in which case the net
                   ! change in the temperature of a layer is the sum of the
                   ! direct heating of that layer plus T_chg_above from all of
                   ! the layers below, plus any contribution from absorbing
                   ! radiation that hits the bottom.
  real, dimension(SZI_(G)) :: &
    h_heat, &      ! The thickness of the water column that will be heated by
                   ! any remaining shortwave radiation (H units).
    t_chg, &       ! The temperature change of thick layers due to the remaining
                   ! shortwave radiation and contributions from T_chg_above, in K.
    pen_sw_rem     ! The sum across all wavelength bands of the penetrating shortwave
                   ! heating that hits the bottom and will be redistributed through
                   ! the water column (in units of K H)
  real :: sw_trans          ! fraction of shortwave radiation that is not
                            ! absorbed in a layer (nondimensional)
  real :: unabsorbed        ! fraction of the shortwave radiation that
                            ! is not absorbed because the layers are too thin
  real :: ih_limit          ! inverse of the total depth at which the
                            ! surface fluxes start to be limited (1/H)
  real :: h_min_heat        ! minimum thickness layer that should get heated (H)
  real :: opt_depth         ! optical depth of a layer (non-dim)
  real :: exp_od            ! exp(-opt_depth) (non-dim)
  real :: heat_bnd          ! heating due to absorption in the current
                            ! layer by the current band, including any piece that
                            ! is moved upward (K H units)
  real :: swa               ! fraction of the absorbed shortwave that is
                            ! moved to layers above with adjustAbsorptionProfile (non-dim)
  real :: coswa_frac        ! The fraction of SWa that is actually moved upward.
  real :: min_sw_heating    ! A minimum remaining shortwave heating rate that will be
                            ! simply absorbed in the next layer for computational
                            ! efficiency, instead of continuing to penetrate, in units
                            ! of K H s-1.  The default, 2.5e-11, is about 0.08 K m / century.
  real :: epsilon           ! A small thickness that must remain in each
                            ! layer, and which will not be subject to heating (units of H)
  real :: i_g_earth
  real :: g_hconv2
  logical :: sw_remains     ! If true, some column has shortwave radiation that
                            ! was not entirely absorbed.
  logical :: tke_calc       ! If true, calculate the implications to the
                            ! TKE budget of the shortwave heating.
  real :: c1_6, c1_60
  integer :: is, ie, nz, i, k, ks, n
  sw_remains = .false.

  min_sw_heating = 2.5e-11

  h_min_heat = 2.0*gv%Angstrom + gv%H_subroundoff
  is = g%isc ; ie = g%iec ; nz = g%ke
  c1_6 = 1.0 / 6.0 ; c1_60 = 1.0 / 60.0

  tke_calc = (present(tke) .and. present(dsv_dt))
  g_hconv2 = gv%g_Earth * gv%H_to_kg_m2**2

  h_heat(:) = 0.0
  if (present(htot)) then ; do i=is,ie ; h_heat(i) = htot(i) ; enddo ; endif

  ! Apply penetrating SW radiation to remaining parts of layers.
  ! Excessively thin layers are not heated to avoid runaway temps.
  do ks=1,nz ; do i=is,ie
    k = ks
    if (present(ksort)) then
      if (ksort(i,ks) <= 0) cycle
      k = ksort(i,ks)
    endif
    epsilon = 0.0 ; if (present(eps)) epsilon = eps(i,k)

    t_chg_above(i,k) = 0.0

    if (h(i,k) > 1.5*epsilon) then
      do n=1,nsw ; if (pen_sw_bnd(n,i) > 0.0) then
        ! SW_trans is the SW that is transmitted THROUGH the layer
        opt_depth = h(i,k) * opacity_band(n,i,k)
        exp_od = exp(-opt_depth)
        sw_trans = exp_od

        ! Heating at a rate of less than 10-4 W m-2 = 10-3 K m / Century,
        ! and of the layer in question less than 1 K / Century, can be
        ! absorbed without further penetration.
        ! ###Make these numbers into parameters!
        if (nsw*pen_sw_bnd(n,i)*sw_trans < &
            dt*min_sw_heating*min(gv%m_to_H,1e3*h(i,k)) ) sw_trans = 0.0

        heat_bnd = pen_sw_bnd(n,i) * (1.0 - sw_trans)
        if (adjustabsorptionprofile .and. (h_heat(i) > 0.0)) then
          !   In this case, a fraction of the heating is applied to the
          ! overlying water so that the mean pressure at which the shortwave
          ! heating occurs is exactly what it would have been with a careful
          ! pressure-weighted averaging of the exponential heating profile,
          ! hence there should be no TKE budget requirements due to this
          ! layer.  Very clever, but this is also limited so that the
          ! water above is not heated at a faster rate than the layer
          ! actually being heated, i.e., SWA <= h_heat / (h_heat + h(i,k))
          ! and takes the energetics of the rest of the heating into account.
          ! (-RWH, ~7 years later.)
          if (opt_depth > 1e-5) then
            swa = ((opt_depth + (opt_depth + 2.0)*exp_od) - 2.0) / &
              ((opt_depth + opacity_band(n,i,k) * h_heat(i)) * &
               (1.0 - exp_od))
          else
            ! Use Taylor series expansion of the expression above for a
            ! more accurate form with very small layer optical depths.
            swa = h(i,k) * (opt_depth * (1.0 - opt_depth)) / &
              ((h_heat(i) + h(i,k)) * (6.0 - 3.0*opt_depth))
          endif
          coswa_frac = 0.0
          if (swa*(h_heat(i) + h(i,k)) > h_heat(i)) then
            coswa_frac = (swa*(h_heat(i) + h(i,k)) - h_heat(i) ) / &
                         (swa*(h_heat(i) + h(i,k)))
            swa = h_heat(i) / (h_heat(i) + h(i,k))
          endif

          t_chg_above(i,k) = t_chg_above(i,k) + (swa * heat_bnd) / h_heat(i)
          t(i,k) = t(i,k) + ((1.0 - swa) * heat_bnd) / h(i,k)
        else
          coswa_frac = 1.0
          t(i,k) = t(i,k) + pen_sw_bnd(n,i) * (1.0 - sw_trans) / h(i,k)
        endif

        if (tke_calc) then
          if (opt_depth > 1e-2) then
            tke(i,k) = tke(i,k) - coswa_frac*heat_bnd*dsv_dt(i,k)* &
               (0.5*h(i,k)*g_hconv2) * &
               (opt_depth*(1.0+exp_od) - 2.0*(1.0-exp_od)) / (opt_depth*(1.0-exp_od))
          else
            ! Use Taylor series-derived approximation to the above expression
            ! that is well behaved and more accurate when opt_depth is small.
            tke(i,k) = tke(i,k) - coswa_frac*heat_bnd*dsv_dt(i,k)* &
               (0.5*h(i,k)*g_hconv2) * &
               (c1_6*opt_depth * (1.0 - c1_60*opt_depth**2))
          endif
        endif

        pen_sw_bnd(n,i) = pen_sw_bnd(n,i) * sw_trans
      endif ; enddo
    endif

    ! Add to the accumulated thickness above that could be heated.
    ! Only layers greater than h_min_heat thick should get heated.
    if (h(i,k) >= 2.0*h_min_heat) then
      h_heat(i) = h_heat(i) + h(i,k)
    elseif (h(i,k) > h_min_heat) then
      h_heat(i) = h_heat(i) + (2.0*h(i,k) - 2.0*h_min_heat)
    endif
  enddo ; enddo ! i & k loops


! if (.not.absorbAllSW .and. .not.adjustAbsorptionProfile) return

  ! Unless modified, there is no temperature change due to fluxes from the bottom.
  do i=is,ie ; t_chg(i) = 0.0 ; enddo

  if (absorballsw) then
    ! If there is still shortwave radiation at this point, it could go into
    ! the bottom (with a bottom mud model), or it could be redistributed back
    ! through the water column.
    do i=is,ie
      pen_sw_rem(i) = pen_sw_bnd(1,i)
      do n=2,nsw ; pen_sw_rem(i) = pen_sw_rem(i) + pen_sw_bnd(n,i) ; enddo
    enddo
    do i=is,ie ; if (pen_sw_rem(i) > 0.0) sw_remains = .true. ; enddo

    ih_limit = 1.0 / h_limit_fluxes
    do i=is,ie ; if ((pen_sw_rem(i) > 0.0) .and. (h_heat(i) > 0.0)) then
      if (h_heat(i)*ih_limit >= 1.0) then
        t_chg(i) = pen_sw_rem(i) / h_heat(i) ; unabsorbed = 0.0
      else
        t_chg(i) = pen_sw_rem(i) * ih_limit
        unabsorbed = 1.0 - h_heat(i)*ih_limit
      endif
      do n=1,nsw ; pen_sw_bnd(n,i) = unabsorbed * pen_sw_bnd(n,i) ; enddo
    endif ; enddo
  endif ! absorbAllSW

  if (absorballsw .or. adjustabsorptionprofile) then
    do ks=nz,1,-1 ; do i=is,ie
      k = ks
      if (present(ksort)) then
        if (ksort(i,ks) <= 0) cycle
        k = ksort(i,ks)
      endif

      if (t_chg(i) > 0.0) then
        ! Only layers greater than h_min_heat thick should get heated.
        if (h(i,k) >= 2.0*h_min_heat) then ; t(i,k) = t(i,k) + t_chg(i)
        elseif (h(i,k) > h_min_heat) then
          t(i,k) = t(i,k) + t_chg(i) * (2.0 - 2.0*h_min_heat/h(i,k))
        endif
      endif
      ! Increase the heating for layers above.
      t_chg(i) = t_chg(i) + t_chg_above(i,k)
    enddo ; enddo
    if (present(htot) .and. present(ttot)) then
      do i=is,ie ; ttot(i) = ttot(i) + t_chg(i) * htot(i) ; enddo
    endif
  endif ! absorbAllSW .or. adjustAbsorptionProfile

Here is the caller graph for this function:

sumswoverbands()

subroutine, public mom_shortwave_abs::sumswoverbands	(	type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		real, dimension(szi_(g),szk_(g)), intent(in)	h,
		real, dimension(:,:,:), intent(in)	opacity_band,
		integer, intent(in)	nsw,
		integer, intent(in)	j,
		real, intent(in)	dt,
		real, intent(in)	H_limit_fluxes,
		logical, intent(in)	absorbAllSW,
		real, dimension(:,:), intent(in)	iPen_SW_bnd,
		real, dimension(szi_(g),szk_(g)+1), intent(inout)	netPen
	)

Parameters

[in]	g	The ocean's grid structure.
[in]	gv	The ocean's vertical grid structure.
[in]	h	Layer thicknesses, in H (usually m or kg m-2).
[in]	opacity_band	opacity in each band of penetrating shortwave radiation, in m-1. The indicies are band, i, k.
[in]	nsw	number of bands of penetrating shortwave radiation.
[in]	j	j-index to work on.
[in]	dt	Time step (seconds).

Definition at line 363 of file MOM_shortwave_abs.F90.

Referenced by mom_diabatic_aux::applyboundaryfluxesinout(), and mom_forcing_type::calculatebuoyancyflux1d().

!< This subroutine calculates the total shortwave heat flux integrated over
!! bands as a function of depth.  This routine is only called for computing
!! buoyancy fluxes for use in KPP. This routine does not updat e the state.
  type(ocean_grid_type),              intent(in)    :: g   !< The ocean's grid structure.
  type(verticalgrid_type),            intent(in)    :: gv  !< The ocean's vertical grid structure.
  real, dimension(SZI_(G),SZK_(G)),   intent(in)    :: h   !< Layer thicknesses, in H (usually m
                                                           !! or kg m-2).
  real, dimension(:,:,:),             intent(in)    :: opacity_band !< opacity in each band of
                                                           !! penetrating shortwave radiation,
                                                           !! in m-1. The indicies are band, i, k.
  integer,                            intent(in)    :: nsw !< number of bands of penetrating
                                                           !! shortwave radiation.
  integer,                            intent(in)    :: j   !< j-index to work on.
  real,                               intent(in)    :: dt  !< Time step (seconds).
  real,                               intent(in)    :: h_limit_fluxes
  logical,                            intent(in)    :: absorballsw
  real, dimension(:,:),               intent(in)    :: ipen_sw_bnd
  real, dimension(SZI_(G),SZK_(G)+1), intent(inout) :: netpen ! Units of K H.

! Arguments:
!  (in)      G             = ocean grid structure
!  (in)      GV            = The ocean's vertical grid structure.
!  (in)      h             = layer thickness (units of m or kg/m^2);
!                            units of h are referred to as H below.
!  (in)      opacity_band  = opacity in each band of penetrating shortwave
!                            radiation, in m-1. The indicies are band, i, k.
!  (in)      nsw           = number of bands of penetrating shortwave radiation
!  (in)      j             = j-index to work on
!  (in)      dt            = time step (seconds)
!  (inout)   Pen_SW_bnd    = penetrating shortwave heating in each band that
!                            hits the bottom and will be redistributed through
!                            the water column (K H units); size nsw x SZI_(G).
!  (out)     netPen        = attenuated flux at interfaces, summed over bands (K H units)

  real :: h_heat(szi_(g))     ! thickness of the water column that receives
                              ! remaining shortwave radiation, in H.
  real :: pen_sw_rem(szi_(g)) ! sum across all wavelength bands of the
                              ! penetrating shortwave heating that hits the bottom
                              ! and will be redistributed through the water column
                              ! (K H units)

  real, dimension(size(iPen_SW_bnd,1),size(iPen_SW_bnd,2)) :: pen_sw_bnd
  real :: sw_trans        ! fraction of shortwave radiation not
                          ! absorbed in a layer (nondimensional)
  real :: unabsorbed      ! fraction of the shortwave radiation
                          ! not absorbed because the layers are too thin.
  real :: ih_limit        ! inverse of the total depth at which the
                          ! surface fluxes start to be limited (1/H units)
  real :: h_min_heat      ! minimum thickness layer that should get heated (H units)
  real :: opt_depth       ! optical depth of a layer (non-dim)
  real :: exp_od          ! exp(-opt_depth) (non-dim)
  logical :: sw_remains   ! If true, some column has shortwave radiation that
                          ! was not entirely absorbed.

  integer :: is, ie, nz, i, k, ks, n
  sw_remains = .false.

  h_min_heat = 2.0*gv%Angstrom + gv%H_subroundoff
  is = g%isc ; ie = g%iec ; nz = g%ke

  pen_sw_bnd(:,:) = ipen_sw_bnd(:,:)
  do i=is,ie ; h_heat(i) = 0.0 ; enddo
  netpen(:,1) = sum( pen_sw_bnd(:,:), dim=1 ) ! Surface interface

  ! Apply penetrating SW radiation to remaining parts of layers.
  ! Excessively thin layers are not heated to avoid runaway temps.
  do k=1,nz

    do i=is,ie
      netpen(i,k+1) = 0.

      if (h(i,k) > 0.0) then
        do n=1,nsw ; if (pen_sw_bnd(n,i) > 0.0) then
          ! SW_trans is the SW that is transmitted THROUGH the layer
          opt_depth = h(i,k)*gv%H_to_m * opacity_band(n,i,k)
          exp_od = exp(-opt_depth)
          sw_trans = exp_od

          ! Heating at a rate of less than 10-4 W m-2 = 10-3 K m / Century,
          ! and of the layer in question less than 1 K / Century, can be
          ! absorbed without further penetration.
          if ((nsw*pen_sw_bnd(n,i)*sw_trans < gv%m_to_H*2.5e-11*dt) .and. &
              (nsw*pen_sw_bnd(n,i)*sw_trans < h(i,k)*dt*2.5e-8)) &
            sw_trans = 0.0

          pen_sw_bnd(n,i) = pen_sw_bnd(n,i) * sw_trans
          netpen(i,k+1)   = netpen(i,k+1) + pen_sw_bnd(n,i)
        endif ; enddo
      endif ! h(i,k) > 0.0

      ! Add to the accumulated thickness above that could be heated.
      ! Only layers greater than h_min_heat thick should get heated.
      if (h(i,k) >= 2.0*h_min_heat) then
        h_heat(i) = h_heat(i) + h(i,k)
      elseif (h(i,k) > h_min_heat) then
        h_heat(i) = h_heat(i) + (2.0*h(i,k) - 2.0*h_min_heat)
      endif
    enddo ! i loop
  enddo ! k loop

  if (absorballsw) then

    ! If there is still shortwave radiation at this point, it could go into
    ! the bottom (with a bottom mud model), or it could be redistributed back
    ! through the water column.
    do i=is,ie
      pen_sw_rem(i) = pen_sw_bnd(1,i)
      do n=2,nsw ; pen_sw_rem(i) = pen_sw_rem(i) + pen_sw_bnd(n,i) ; enddo
    enddo
    do i=is,ie ; if (pen_sw_rem(i) > 0.0) sw_remains = .true. ; enddo

    ih_limit = 1.0 / h_limit_fluxes
    do i=is,ie ; if ((pen_sw_rem(i) > 0.0) .and. (h_heat(i) > 0.0)) then
      if (h_heat(i)*ih_limit < 1.0) then
        unabsorbed = 1.0 - h_heat(i)*ih_limit
      else
        unabsorbed = 0.0
      endif
      do n=1,nsw ; pen_sw_bnd(n,i) = unabsorbed * pen_sw_bnd(n,i) ; enddo
    endif ; enddo

  endif ! absorbAllSW

Here is the caller graph for this function: