mom_opacity Module Reference

Data Types
type	opacity_cs

Functions/Subroutines
subroutine, public	set_opacity (optics, fluxes, G, GV, CS)
subroutine	opacity_from_chl (optics, fluxes, G, CS, chl_in)
real function, public	opacity_morel (chl_data)
real function	sw_pen_frac_morel (chl_data)
real function, public	opacity_manizza (chl_data)
subroutine, public	opacity_init (Time, G, param_file, diag, tracer_flow, CS, optics)
subroutine, public	opacity_end (CS, optics)

Variables
integer, parameter	no_scheme = 0
integer, parameter	manizza_05 = 1
integer, parameter	morel_88 = 2
integer, parameter	single_exp = 3
integer, parameter	double_exp = 4
character *(10), parameter	manizza_05_string = "MANIZZA_05"
character *(10), parameter	morel_88_string = "MOREL_88"
character *(10), parameter	single_exp_string = "SINGLE_EXP"
character *(10), parameter	double_exp_string = "DOUBLE_EXP"

Function/Subroutine Documentation

opacity_end()

subroutine, public mom_opacity::opacity_end	(	type(opacity_cs), pointer	CS,
		type(optics_type), optional, pointer	optics
	)

Definition at line 692 of file MOM_opacity.F90.

Referenced by mom_diabatic_driver::diabatic_driver_end().

  type(opacity_cs),  pointer           :: cs
  type(optics_type), pointer, optional :: optics

  if (associated(cs%id_opacity)) deallocate(cs%id_opacity)
  if (associated(cs)) deallocate(cs)

  if (present(optics)) then ; if (associated(optics)) then
    if (ASSOCIATED(optics%opacity_band)) deallocate(optics%opacity_band)
    if (ASSOCIATED(optics%sw_pen_band)) deallocate(optics%sw_pen_band)
  endif ; endif

Here is the caller graph for this function:

opacity_from_chl()

subroutine mom_opacity::opacity_from_chl ( type(optics_type), intent(inout)  optics, type(forcing), intent(in)  fluxes, type(ocean_grid_type), intent(in)  G, type(opacity_cs), pointer  CS, real, dimension( g %isd: g %ied, g %jsd: g %jed, g %ke), intent(in), optional  chl_in  )

private

Parameters

[in]	fluxes	A structure containing pointers to any possible forcing fields. Unused fields have NULL ptrs.
[in]	g	The ocean's grid structure.
	cs	The control structure.
[in]	chl_in	A 3-d field of chlorophyll A,

Definition at line 251 of file MOM_opacity.F90.

References manizza_05, mom_error_handler::mom_error(), morel_88, opacity_morel(), and sw_pen_frac_morel().

Referenced by set_opacity().

  type(optics_type),              intent(inout)  :: optics
  type(forcing),                  intent(in)     :: fluxes !< A structure containing pointers to any
                                                           !! possible forcing fields. Unused fields
                                                           !! have NULL ptrs.
  type(ocean_grid_type),          intent(in)     :: g      !< The ocean's grid structure.
  type(opacity_cs),               pointer        :: cs     !< The control structure.
  real, dimension(SZI_(G),SZJ_(G),SZK_(G)), &
                      intent(in), optional       :: chl_in !< A 3-d field of chlorophyll A,
                                                           !! in mg m-3.
! Arguments: fluxes - A structure containing pointers to any possible
!                     forcing fields.  Unused fields have NULL ptrs.
!  (out)     opacity - The inverse of the vertical absorption decay
!                           scale for penetrating shortwave radiation, in m-1.
!  (in)      G - The ocean's grid structure.
!  (in)      chl_in - A 3-d field of chlorophyll A, in mg m-3.

  real :: chl_data(szi_(g),szj_(g)) ! The chlorophyll A concentrations in
                                    ! a layer, in mg/m^3.
  real :: inv_nbands        ! The inverse of the number of bands of penetrating
                            ! shortwave radiation.
  real :: inv_nbands_nir    ! The inverse of the number of bands of penetrating
                            ! near-infrafed radiation.
  real :: sw_pen_tot        ! The sum across the bands of the penetrating
                            ! shortwave radiation, in W m-2.
  real :: sw_vis_tot        ! The sum across the visible bands of shortwave
                            ! radiation, in W m-2.
  real :: sw_nir_tot        ! The sum across the near infrared bands of shortwave
                            ! radiation, in W m-2.
  type(time_type) :: day
  character(len=128) :: mesg
  integer :: days, seconds
  integer :: i, j, k, n, is, ie, js, je, nz, nbands
  logical :: multiband_vis_input, multiband_nir_input

  is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke

!   In this model, the Morel (modified) and Manizza (modified) schemes
! use the "blue" band in the paramterizations to determine the e-folding
! depth of the incoming shortwave attenuation. The red portion is lumped
! into the net heating at the surface.
!
! Morel, A., Optical modeling of the upper ocean in relation to its biogenous
!   matter content (case-i waters).,J. Geo. Res., {93}, 10,749--10,768, 1988.
!

! Manizza, M., C.~L. Quere, A.~Watson, and E.~T. Buitenhuis, Bio-optical
!   feedbacks amoung phytoplankton, upper ocean physics and sea-ice in a
!   global model, Geophys. Res. Let., , L05,603, 2005.

  nbands = optics%nbands

  if (nbands <= 1) then ; inv_nbands = 1.0
  else ; inv_nbands = 1.0 / real(nbands) ; endif

  if (nbands <= 2) then ; inv_nbands_nir = 0.0
  else ; inv_nbands_nir = 1.0 / real(nbands - 2.0) ; endif

  multiband_vis_input = (associated(fluxes%sw_vis_dir) .and. &
                         associated(fluxes%sw_vis_dif))
  multiband_nir_input = (associated(fluxes%sw_nir_dir) .and. &
                         associated(fluxes%sw_nir_dif))

  chl_data(:,:) = 0.0
  if(present(chl_in)) then
    do j=js,je ; do i=is,ie ; chl_data(i,j) = chl_in(i,j,1) ; enddo ; enddo
    do k=1,nz; do j=js,je ; do i=is,ie
      if ((g%mask2dT(i,j) > 0.5) .and. (chl_in(i,j,k) < 0.0)) then
          write(mesg,'(" Negative chl_in of ",(1pe12.4)," found at i,j,k = ", &
                    & 3(1x,i3), " lon/lat = ",(1pe12.4)," E ", (1pe12.4), " N.")') &
                     chl_in(i,j,k), i, j, k, g%geoLonT(i,j), g%geoLatT(i,j)
          call mom_error(fatal,"MOM_opacity opacity_from_chl: "//trim(mesg))
      endif
    enddo; enddo; enddo
  else
    ! Only the 2-d surface chlorophyll can be read in from a file.  The
    ! same value is assumed for all layers.
    call get_time(cs%Time,seconds,days)
    call time_interp_external(cs%sbc_chl, cs%Time, chl_data)
    do j=js,je ; do i=is,ie
      if ((g%mask2dT(i,j) > 0.5) .and. (chl_data(i,j) < 0.0)) then
        write(mesg,'(" Time_interp negative chl of ",(1pe12.4)," at i,j = ",&
                  & 2(i3), "lon/lat = ",(1pe12.4)," E ", (1pe12.4), " N.")') &
                   chl_data(i,j), i, j, g%geoLonT(i,j), g%geoLatT(i,j)
        call mom_error(fatal,"MOM_opacity opacity_from_chl: "//trim(mesg))
      endif
    enddo ; enddo
  endif

  if (cs%id_chl > 0) then
    if(present(chl_in)) then
      call post_data(cs%id_chl, chl_in(:,:,1), cs%diag)
    else
      call post_data(cs%id_chl, chl_data, cs%diag)
    endif
  endif

  select case (cs%opacity_scheme)
    case (manizza_05)
!$OMP parallel do default(none) shared(is,ie,js,je,fluxes,optics,CS,G,multiband_nir_input, &
!$OMP                                  nbands,Inv_nbands_nir,multiband_vis_input )         &
!$OMP                          private(SW_vis_tot,SW_nir_tot)
      do j=js,je ; do i=is,ie
        sw_vis_tot = 0.0 ; sw_nir_tot = 0.0
        if (g%mask2dT(i,j) > 0.5) then
          if (multiband_vis_input) then
            sw_vis_tot = fluxes%sw_vis_dir(i,j) + fluxes%sw_vis_dif(i,j)
          else  ! Follow Manizza 05 in assuming that 42% of SW is visible.
            sw_vis_tot = 0.42 * fluxes%sw(i,j)
          endif
          if (multiband_nir_input) then
            sw_nir_tot = fluxes%sw_nir_dir(i,j) + fluxes%sw_nir_dif(i,j)
          else
            sw_nir_tot = fluxes%sw(i,j) - sw_vis_tot
          endif
        endif

        ! Band 1 is Manizza blue.
        optics%sw_pen_band(1,i,j) = cs%blue_frac*sw_vis_tot
        ! Band 2 (if used) is Manizza red.
        if (nbands > 1) &
          optics%sw_pen_band(2,i,j) = (1.0-cs%blue_frac)*sw_vis_tot
        ! All remaining bands are NIR, for lack of something better to do.
        do n=3,nbands
          optics%sw_pen_band(n,i,j) = inv_nbands_nir * sw_nir_tot
        enddo
      enddo ; enddo
    case (morel_88)
!$OMP parallel do default(none) shared(is,ie,js,je,G,multiband_vis_input,chl_data, &
!$OMP                                  fluxes,nbands,optics,Inv_nbands) &
!$OMP                          private(SW_pen_tot)
      do j=js,je ; do i=is,ie
        sw_pen_tot = 0.0
        if (g%mask2dT(i,j) > 0.5) then ; if (multiband_vis_input) then
            sw_pen_tot = sw_pen_frac_morel(chl_data(i,j)) * &
                (fluxes%sw_vis_dir(i,j) + fluxes%sw_vis_dif(i,j))
          else
            sw_pen_tot = sw_pen_frac_morel(chl_data(i,j)) * &
                0.5*fluxes%sw(i,j)
        endif ; endif

        do n=1,nbands
          optics%sw_pen_band(n,i,j) = inv_nbands*sw_pen_tot
        enddo
      enddo ; enddo
    case default
        call mom_error(fatal, "opacity_from_chl: CS%opacity_scheme is not valid.")
    end select

!$OMP parallel do default(none) shared(nz,is,ie,js,je,CS,G,chl_in,optics,nbands) &
!$OMP                     firstprivate(chl_data)
  do k=1,nz
    if (present(chl_in)) then
      do j=js,je ; do i=is,ie ; chl_data(i,j) = chl_in(i,j,k) ; enddo ; enddo
    endif

    select case (cs%opacity_scheme)
      case (manizza_05)
        do j=js,je ; do i=is,ie
          if (g%mask2dT(i,j) <= 0.5) then
            do n=1,optics%nbands
              optics%opacity_band(n,i,j,k) = cs%opacity_land_value
            enddo
          else
            ! Band 1 is Manizza blue.
            optics%opacity_band(1,i,j,k) = 0.0232 + 0.074*chl_data(i,j)**0.674
            if (nbands >= 2) &  !  Band 2 is Manizza red.
              optics%opacity_band(2,i,j,k) = 0.225 + 0.037*chl_data(i,j)**0.629
            ! All remaining bands are NIR, for lack of something better to do.
            do n=3,nbands ; optics%opacity_band(n,i,j,k) = 2.86 ; enddo
          endif
        enddo ; enddo
      case (morel_88)
        do j=js,je ; do i=is,ie
          optics%opacity_band(1,i,j,k) = cs%opacity_land_value
          if (g%mask2dT(i,j) > 0.5) &
            optics%opacity_band(1,i,j,k) = opacity_morel(chl_data(i,j))

          do n=2,optics%nbands
            optics%opacity_band(n,i,j,k) = optics%opacity_band(1,i,j,k)
          enddo
        enddo ; enddo

      case default
        call mom_error(fatal, "opacity_from_chl: CS%opacity_scheme is not valid.")
    end select
  enddo

Here is the call graph for this function:

Here is the caller graph for this function:

opacity_init()

subroutine, public mom_opacity::opacity_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(tracer_flow_control_cs), intent(in), target	tracer_flow,
		type(opacity_cs), pointer	CS,
		type(optics_type), pointer	optics
	)

Parameters

[in]	time	The current model time.
[in]	g	The ocean's grid structure.
[in]	param_file	A structure to parse for run-time parameters.
[in,out]	diag	A structure that is used to regulate diagnostic output.
	cs	A pointer that is set to point to the control structure for this module.

Definition at line 488 of file MOM_opacity.F90.

References double_exp, double_exp_string, manizza_05, manizza_05_string, morel_88, morel_88_string, no_scheme, single_exp, and single_exp_string.

Referenced by mom_diabatic_driver::diabatic_driver_init().

  type(time_type), target, intent(in)    :: time !< The current model time.
  type(ocean_grid_type),   intent(in)    :: g    !< The ocean's grid structure.
  type(param_file_type),   intent(in)    :: param_file !< A structure to parse for run-time
                                                 !! parameters.
  type(diag_ctrl), target, intent(inout) :: diag !< A structure that is used to regulate diagnostic
                                                 !! output.
  type(tracer_flow_control_cs), &
                  target, intent(in)     :: tracer_flow
  type(opacity_cs),        pointer       :: cs   !< A pointer that is set to point to the control
                                                 !! structure for this module.
  type(optics_type),       pointer       :: optics
! Arguments: Time - The current model time.
!  (in)      G - The ocean's grid structure.
!  (in)      param_file - A structure indicating the open file to parse for
!                         model parameter values.
!  (in)      diag - A structure that is used to regulate diagnostic output.
!
!  (in/out)  CS - A pointer that is set to point to the control structure
!                  for this module
! This include declares and sets the variable "version".
#include "version_variable.h"
  character(len=200) :: inputdir   ! The directory where NetCDF input files
  character(len=240) :: filename
  character(len=200) :: tmpstr
  character(len=40)  :: mdl = "MOM_opacity"
  character(len=40)  :: bandnum, shortname
  character(len=200) :: longname
  character(len=40)  :: scheme_string
  logical :: use_scheme
  integer :: isd, ied, jsd, jed, nz, n
  isd = g%isd ; ied = g%ied ; jsd = g%jsd ; jed = g%jed ; nz = g%ke

  if (associated(cs)) then
    call mom_error(warning, "opacity_init called with an associated"// &
                             "associated control structure.")
    return
  else ; allocate(cs) ; endif

  cs%diag => diag
  cs%Time => time
  cs%tracer_flow_CSp => tracer_flow

  ! Read all relevant parameters and write them to the model log.
  call log_version(param_file, mdl, version, "")

! parameters for CHL_A routines
  call get_param(param_file, mdl, "VAR_PEN_SW", cs%var_pen_sw, &
                 "If true, use one of the CHL_A schemes specified by \n"//&
                 "OPACITY_SCHEME to determine the e-folding depth of \n"//&
                 "incoming short wave radiation.", default=.false.)

  cs%opacity_scheme = no_scheme ; scheme_string = ""
  if (cs%var_pen_sw) then
    call get_param(param_file, mdl, "OPACITY_SCHEME", tmpstr, &
                 "This character string specifies how chlorophyll \n"//&
                 "concentrations are translated into opacities. Currently \n"//&
                 "valid options include:\n"//&
                 " \t\t  MANIZZA_05 - Use Manizza et al., GRL, 2005. \n"//&
                 " \t\t  MOREL_88 - Use Morel, JGR, 1988.", &
                 default=manizza_05_string)
    if (len_trim(tmpstr) > 0) then
      tmpstr = uppercase(tmpstr)
      select case (tmpstr)
        case (manizza_05_string)
          cs%opacity_scheme = manizza_05 ; scheme_string = manizza_05_string
        case (morel_88_string)
          cs%opacity_scheme = morel_88 ; scheme_string = morel_88_string
        case default
          call mom_error(fatal, "opacity_init: #DEFINE OPACITY_SCHEME "//&
                                  trim(tmpstr) // "in input file is invalid.")
      end select
      call mom_mesg('opacity_init: opacity scheme set to "'//trim(tmpstr)//'".', 5)
    endif
    if (cs%opacity_scheme == no_scheme) then
      call mom_error(warning, "opacity_init: No scheme has successfully "//&
               "been specified for the opacity.  Using the default MANIZZA_05.")
      cs%opacity_scheme = manizza_05 ; scheme_string = manizza_05_string
    endif

    call get_param(param_file, mdl, "CHL_FROM_FILE", cs%chl_from_file, &
                 "If true, chl_a is read from a file.", default=.true.)
    if (cs%chl_from_file) then
      call time_interp_external_init()

      call get_param(param_file, mdl, "INPUTDIR", inputdir, default=".")
      call get_param(param_file, mdl, "CHL_FILE", cs%chl_file, &
                 "CHL_FILE is the file containing chl_a concentrations in \n"//&
                 "the variable CHL_A. It is used when VAR_PEN_SW and \n"//&
                 "CHL_FROM_FILE are true.", fail_if_missing=.true.)

      filename = trim(slasher(inputdir))//trim(cs%chl_file)
      call log_param(param_file, mdl, "INPUTDIR/CHL_FILE", filename)
      cs%sbc_chl = init_external_field(filename,'CHL_A',domain=g%Domain%mpp_domain)
    endif

    call get_param(param_file, mdl, "BLUE_FRAC_SW", cs%blue_frac, &
                 "The fraction of the penetrating shortwave radiation \n"//&
                 "that is in the blue band.", default=0.5, units="nondim")
  else
    call get_param(param_file, mdl, "EXP_OPACITY_SCHEME", tmpstr, &
                 "This character string specifies which exponential \n"//&
                 "opacity scheme to utilize. Currently \n"//&
                 "valid options include:\n"//&
                 " \t\t  SINGLE_EXP - Single Exponent decay. \n"//&
                 " \t\t  DOUBLE_EXP - Double Exponent decay.", &
                 default=single_exp_string)!New default for "else" above (non-Chl scheme)
    if (len_trim(tmpstr) > 0) then
      tmpstr = uppercase(tmpstr)
      select case (tmpstr)
        case (single_exp_string)
          cs%opacity_scheme = single_exp ; scheme_string = single_exp_string
        case (double_exp_string)
          cs%opacity_scheme = double_exp ; scheme_string = double_exp_string
      end select
      call mom_mesg('opacity_init: opacity scheme set to "'//trim(tmpstr)//'".', 5)
    endif
    call get_param(param_file, mdl, "PEN_SW_SCALE", cs%pen_sw_scale, &
                 "The vertical absorption e-folding depth of the \n"//&
                 "penetrating shortwave radiation.", units="m", default=0.0)
    !BGR/ Added for opacity_scheme==double_exp read in 2nd exp-decay and fraction
    if (cs%Opacity_scheme == double_exp ) then
      call get_param(param_file, mdl, "PEN_SW_SCALE_2ND", cs%pen_sw_scale_2nd, &
                 "The (2nd) vertical absorption e-folding depth of the \n"//&
                 "penetrating shortwave radiation \n"//&
                 "(use if SW_EXP_MODE==double.)",&
                 units="m", default=0.0)
      call get_param(param_file, mdl, "SW_1ST_EXP_RATIO", cs%sw_1st_exp_ratio, &
                 "The fraction of 1st vertical absorption e-folding depth \n"//&
                 "penetrating shortwave radiation if SW_EXP_MODE==double.",&
                  units="m", default=0.0)
    elseif (cs%OPACITY_SCHEME == single_exp) then
      !/Else disable 2nd_exp scheme
      cs%pen_sw_scale_2nd = 0.0
      cs%sw_1st_exp_ratio = 1.0
    endif
    call get_param(param_file, mdl, "PEN_SW_FRAC", cs%pen_sw_frac, &
                 "The fraction of the shortwave radiation that penetrates \n"//&
                 "below the surface.", units="nondim", default=0.0)

  endif
  call get_param(param_file, mdl, "PEN_SW_NBANDS", optics%nbands, &
                 "The number of bands of penetrating shortwave radiation.", &
                 default=1)
  if (cs%Opacity_scheme == double_exp ) then
    if (optics%nbands.ne.2) then
      call mom_error(fatal, "set_opacity: "// &
         "Cannot use a double_exp opacity scheme with nbands!=2.")
    endif
  elseif (cs%Opacity_scheme == single_exp ) then
    if (optics%nbands.ne.1) then
      call mom_error(fatal, "set_opacity: "// &
         "Cannot use a single_exp opacity scheme with nbands!=1.")
    endif
  endif
  if (.not.ASSOCIATED(optics%min_wavelength_band)) &
    allocate(optics%min_wavelength_band(optics%nbands))
  if (.not.ASSOCIATED(optics%max_wavelength_band)) &
    allocate(optics%max_wavelength_band(optics%nbands))

  if (cs%opacity_scheme == manizza_05) then
    optics%min_wavelength_band(1) =0
    optics%max_wavelength_band(1) =550
    if (optics%nbands >= 2) then
      optics%min_wavelength_band(2)=550
      optics%max_wavelength_band(2)=700
    endif
    if (optics%nbands > 2) then
      do n=3,optics%nbands
        optics%min_wavelength_band(n) =700
        optics%max_wavelength_band(n) =2800
      enddo
    endif
  endif

  call get_param(param_file, mdl, "OPACITY_LAND_VALUE", cs%opacity_land_value, &
                 "The value to use for opacity over land. The default is \n"//&
                 "10 m-1 - a value for muddy water.", units="m-1", default=10.0)

  if (.not.ASSOCIATED(optics%opacity_band)) &
    allocate(optics%opacity_band(optics%nbands,isd:ied,jsd:jed,nz))
  if (.not.ASSOCIATED(optics%sw_pen_band)) &
    allocate(optics%sw_pen_band(optics%nbands,isd:ied,jsd:jed))
  allocate(cs%id_opacity(optics%nbands)) ; cs%id_opacity(:) = -1

  cs%id_sw_pen = register_diag_field('ocean_model', 'SW_pen', diag%axesT1, time, &
      'Penetrating shortwave radiation flux into ocean', 'Watt meter-2')
  cs%id_sw_vis_pen = register_diag_field('ocean_model', 'SW_vis_pen', diag%axesT1, time, &
      'Visible penetrating shortwave radiation flux into ocean', 'Watt meter-2')
  do n=1,optics%nbands
    write(bandnum,'(i3)') n
    shortname = 'opac_'//trim(adjustl(bandnum))
    longname = 'Opacity for shortwave radiation in band '//trim(adjustl(bandnum))
    cs%id_opacity(n) = register_diag_field('ocean_model', shortname, diag%axesTL, time, &
      longname, 'meter-1')
  enddo
  if (cs%var_pen_sw) &
    cs%id_chl = register_diag_field('ocean_model', 'Chl_opac', diag%axesT1, time, &
        'Surface chlorophyll A concentration used to find opacity', 'mg meter-3')

Here is the caller graph for this function:

opacity_manizza()

real function, public mom_opacity::opacity_manizza

(

real, intent(in)

chl_data

)

Definition at line 478 of file MOM_opacity.F90.

  real, intent(in)  :: chl_data
  real :: opacity_manizza
! Argument : chl_data - The chlorophyll-A concentration in mg m-3.
!   This sets the blue-wavelength opacity according to the scheme proposed by
! Manizza, M. et al, 2005.

  opacity_manizza = 0.0232 + 0.074*chl_data**0.674

opacity_morel()

real function, public mom_opacity::opacity_morel

(

real, intent(in)

chl_data

)

Definition at line 442 of file MOM_opacity.F90.

Referenced by opacity_from_chl().

  real, intent(in)  :: chl_data
  real :: opacity_morel
! Argument : chl_data - The chlorophyll-A concentration in mg m-3.
!   The following are coefficients for the optical model taken from Morel and
! Antoine (1994). These coeficients represent a non uniform distribution of
! chlorophyll-a through the water column.  Other approaches may be more
! appropriate when using an interactive ecosystem model that predicts
! three-dimensional chl-a values.
  real, dimension(6), parameter :: &
       z2_coef=(/7.925, -6.644, 3.662, -1.815, -0.218,  0.502/)
  real :: chl, chl2 ! The log10 of chl_data (in mg m-3), and Chl^2.

  chl = log10(min(max(chl_data,0.02),60.0)) ; chl2 = chl*chl
  opacity_morel = 1.0 / ( (z2_coef(1) + z2_coef(2)*chl) + chl2 * &
      ((z2_coef(3) + chl*z2_coef(4)) + chl2*(z2_coef(5) + chl*z2_coef(6))) )

Here is the caller graph for this function:

set_opacity()

subroutine, public mom_opacity::set_opacity	(	type(optics_type), intent(inout)	optics,
		type(forcing), intent(in)	fluxes,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(opacity_cs), pointer	CS
	)

Parameters

[in]	fluxes	A structure containing pointers to any possible forcing fields. Unused fields have NULL ptrs.
[in]	g	The ocean's grid structure.
[in]	gv	The ocean's vertical grid structure.
	cs	The control structure earlier set up by opacity_init.

Definition at line 122 of file MOM_opacity.F90.

References double_exp, mom_tracer_flow_control::get_chl_from_model(), manizza_05, mom_error_handler::mom_error(), opacity_from_chl(), and mom_diag_mediator::query_averaging_enabled().

Referenced by mom_diabatic_driver::diabatic().

  type(optics_type),         intent(inout) :: optics
  type(forcing),             intent(in)    :: fluxes !< A structure containing pointers to any
                                                     !! possible forcing fields. Unused fields
                                                     !! have NULL ptrs.
  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(opacity_cs),          pointer       :: cs     !< The control structure earlier set up by
                                                     !! opacity_init.

! Arguments: (inout) opacity - The inverse of the vertical absorption decay
!                     scale for penetrating shortwave radiation, in m-1.
!            (inout) fluxes - A structure containing pointers to any possible
!                     forcing fields.  Unused fields have NULL ptrs.
!            (in)    G - The ocean's grid structure.
!  (in)      GV - The ocean's vertical grid structure.
!            (in)    CS - The control structure earlier set up by opacity_init.

! local variables
  integer :: i, j, k, n, is, ie, js, je, nz
  real :: inv_sw_pen_scale  ! The inverse of the e-folding scale, in m-1.
  real :: inv_nbands        ! The inverse of the number of bands of penetrating
                            ! shortwave radiation.
  logical :: call_for_surface  ! if horizontal slice is the surface layer
  real :: tmp(szi_(g),szj_(g),szk_(g))  ! A 3-d temporary array.
  real :: chl(szi_(g),szj_(g),szk_(g))  ! The concentration of chlorophyll-A
                                        ! in mg m-3.
  real :: pen_sw_tot(szi_(g),szj_(g))   ! The penetrating shortwave radiation
                                        ! summed across all bands, in W m-2.
  is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke

  if (.not. associated(cs)) call mom_error(fatal, "set_opacity: "// &
         "Module must be initialized via opacity_init before it is used.")

  if (cs%var_pen_sw) then
    if (cs%chl_from_file) then
      call opacity_from_chl(optics, fluxes, g, cs)
    else
      call get_chl_from_model(chl, g, cs%tracer_flow_CSp)
      call opacity_from_chl(optics, fluxes, g, cs, chl)
    endif
  else ! Use sw e-folding scale set by MOM_input
    if (optics%nbands <= 1) then ; inv_nbands = 1.0
    else ; inv_nbands = 1.0 / real(optics%nbands) ; endif

    ! Make sure there is no division by 0.
    inv_sw_pen_scale = 1.0 / max(cs%pen_sw_scale, 0.1*gv%Angstrom_z, &
                                 gv%H_to_m*gv%H_subroundoff)
!$OMP parallel default(none) shared(is,ie,js,je,nz,optics,inv_sw_pen_scale,fluxes,CS,Inv_nbands,GV)
    if ( cs%Opacity_scheme == double_exp ) then
!$OMP do
      do k=1,nz ; do j=js,je ; do i=is,ie
        optics%opacity_band(1,i,j,k) = inv_sw_pen_scale
        optics%opacity_band(2,i,j,k) = 1.0 / max(cs%pen_sw_scale_2nd, &
             0.1*gv%Angstrom_z,gv%H_to_m*gv%H_subroundoff)
      enddo ; enddo ; enddo
      if (.not.associated(fluxes%sw) .or. (cs%pen_SW_scale <= 0.0)) then
!$OMP do
        do j=js,je ; do i=is,ie ; do n=1,optics%nbands
          optics%sw_pen_band(n,i,j) = 0.0
        enddo ; enddo ; enddo
      else
!$OMP do
        do j=js,je ; do i=is,ie ;
          optics%sw_pen_band(1,i,j) = (cs%SW_1st_EXP_RATIO) * fluxes%sw(i,j)
          optics%sw_pen_band(2,i,j) = (1.-cs%SW_1st_EXP_RATIO) * fluxes%sw(i,j)
        enddo ; enddo ;
      endif
    else
      do k=1,nz ; do j=js,je ; do i=is,ie  ; do n=1,optics%nbands
        optics%opacity_band(n,i,j,k) = inv_sw_pen_scale
      enddo ; enddo ; enddo ; enddo
      if (.not.associated(fluxes%sw) .or. (cs%pen_SW_scale <= 0.0)) then
!$OMP do
        do j=js,je ; do i=is,ie ; do n=1,optics%nbands
          optics%sw_pen_band(n,i,j) = 0.0
        enddo ; enddo ; enddo
      else
!$OMP do
        do j=js,je ; do i=is,ie ; do n=1,optics%nbands
          optics%sw_pen_band(n,i,j) = cs%pen_SW_frac * inv_nbands * fluxes%sw(i,j)
        enddo ; enddo ; enddo
      endif
    endif
!$OMP end parallel
  endif
  if (query_averaging_enabled(cs%diag)) then
    if (cs%id_sw_pen > 0) then
!$OMP parallel do default(none) shared(is,ie,js,je,Pen_SW_tot,optics)
      do j=js,je ; do i=is,ie
        pen_sw_tot(i,j) = 0.0
        do n=1,optics%nbands
          pen_sw_tot(i,j) = pen_sw_tot(i,j) + optics%sw_pen_band(n,i,j)
        enddo
      enddo ; enddo
      call post_data(cs%id_sw_pen, pen_sw_tot, cs%diag)
    endif
    if (cs%id_sw_vis_pen > 0) then
      if (cs%opacity_scheme == manizza_05) then
!$OMP parallel do default(none) shared(is,ie,js,je,Pen_SW_tot,optics)
        do j=js,je ; do i=is,ie
          pen_sw_tot(i,j) = 0.0
          do n=1,min(optics%nbands,2)
            pen_sw_tot(i,j) = pen_sw_tot(i,j) + optics%sw_pen_band(n,i,j)
          enddo
        enddo ; enddo
      else
!$OMP parallel do default(none) shared(is,ie,js,je,Pen_SW_tot,optics)
        do j=js,je ; do i=is,ie
          pen_sw_tot(i,j) = 0.0
          do n=1,optics%nbands
            pen_sw_tot(i,j) = pen_sw_tot(i,j) + optics%sw_pen_band(n,i,j)
          enddo
        enddo ; enddo
      endif
      call post_data(cs%id_sw_vis_pen, pen_sw_tot, cs%diag)
    endif
    do n=1,optics%nbands ; if (cs%id_opacity(n) > 0) then
!$OMP parallel do default(none) shared(nz,is,ie,js,je,tmp,optics,n)
      do k=1,nz ; do j=js,je ; do i=is,ie
        tmp(i,j,k) = optics%opacity_band(n,i,j,k)
      enddo ; enddo ; enddo
      call post_data(cs%id_opacity(n), tmp, cs%diag)
    endif ; enddo
  endif

Here is the call graph for this function:

Here is the caller graph for this function:

sw_pen_frac_morel()

real function mom_opacity::sw_pen_frac_morel ( real, intent(in)  chl_data )

private

Definition at line 460 of file MOM_opacity.F90.

Referenced by opacity_from_chl().

  real, intent(in)  :: chl_data
  real :: sw_pen_frac_morel
! Argument : chl_data - The chlorophyll-A concentration in mg m-3.
!   The following are coefficients for the optical model taken from Morel and
! Antoine (1994). These coeficients represent a non uniform distribution of
! chlorophyll-a through the water column.  Other approaches may be more
! appropriate when using an interactive ecosystem model that predicts
! three-dimensional chl-a values.
  real :: chl, chl2         ! The log10 of chl_data in mg m-3, and Chl^2.
  real, dimension(6), parameter :: &
       v1_coef=(/0.321,  0.008, 0.132,  0.038, -0.017, -0.007/)

  chl = log10(min(max(chl_data,0.02),60.0)) ; chl2 = chl*chl
  sw_pen_frac_morel = 1.0 - ( (v1_coef(1) + v1_coef(2)*chl) + chl2 * &
       ((v1_coef(3) + chl*v1_coef(4)) + chl2*(v1_coef(5) + chl*v1_coef(6))) )

Here is the caller graph for this function:

Variable Documentation

double_exp

integer, parameter mom_opacity::double_exp = 4

private

Definition at line 111 of file MOM_opacity.F90.

Referenced by opacity_init(), and set_opacity().

double_exp_string

character*(10), parameter mom_opacity::double_exp_string = "DOUBLE_EXP"

private

Definition at line 117 of file MOM_opacity.F90.

Referenced by opacity_init().

character*(10), parameter :: DOUBLE_EXP_STRING = "DOUBLE_EXP"

manizza_05

integer, parameter mom_opacity::manizza_05 = 1

private

Definition at line 111 of file MOM_opacity.F90.

Referenced by opacity_from_chl(), opacity_init(), and set_opacity().

manizza_05_string

character*(10), parameter mom_opacity::manizza_05_string = "MANIZZA_05"

private

Definition at line 114 of file MOM_opacity.F90.

Referenced by opacity_init().

character*(10), parameter :: MANIZZA_05_STRING = "MANIZZA_05"

morel_88

integer, parameter mom_opacity::morel_88 = 2

private

Definition at line 111 of file MOM_opacity.F90.

Referenced by opacity_from_chl(), and opacity_init().

morel_88_string

character*(10), parameter mom_opacity::morel_88_string = "MOREL_88"

private

Definition at line 115 of file MOM_opacity.F90.

Referenced by opacity_init().

character*(10), parameter :: MOREL_88_STRING = "MOREL_88"

no_scheme

integer, parameter mom_opacity::no_scheme = 0

Definition at line 111 of file MOM_opacity.F90.

Referenced by opacity_init().

integer, parameter :: no_scheme = 0, manizza_05 = 1, morel_88 = 2, &
                      single_exp = 3, double_exp = 4

single_exp

integer, parameter mom_opacity::single_exp = 3

private

Definition at line 111 of file MOM_opacity.F90.

Referenced by opacity_init().

single_exp_string

character*(10), parameter mom_opacity::single_exp_string = "SINGLE_EXP"

private

Definition at line 116 of file MOM_opacity.F90.

Referenced by opacity_init().

character*(10), parameter :: SINGLE_EXP_STRING = "SINGLE_EXP"