module ocean / physics / surface thermal conductivity More...

Functions/Subroutines
subroutine, public	ocean_phy_tc_seaice_setup

subroutine, public	ocean_phy_tc_seaice (OIA, OIS, OIE, OJA, OJS, OJE, ICE_MASS, ICE_FRAC, mask, TC_dz)

Variables
real(rp), public	ocean_phy_thermalcond_max = 10.0_RP

real(rp), public	ocean_phy_thermalcond_seaice = 2.0_RP

Detailed Description

module ocean / physics / surface thermal conductivity

Description: surface thermal conductivity common module

Author: Team SCALE

NAMELIST

PARAM_OCEAN_PHY_TC_seaice

name	type	default value	comment
OCEAN_PHY_THERMALCOND_MAX	real(RP)	10.0_RP	maximum thermal conductivity / depth [J/m2/s/K]
OCEAN_PHY_THERMALCOND_SEAICE	real(RP)	2.0_RP	thermal conductivity of sea ice [J/m/s/K]

History Output: No history output

Function/Subroutine Documentation

◆ ocean_phy_tc_seaice_setup()

subroutine, public scale_ocean_phy_tc::ocean_phy_tc_seaice_setup

Definition at line 49 of file scale_ocean_phy_tc.F90.

     use scale_prc, only: &
        prc_abort
     implicit none
  
     namelist / param_ocean_phy_tc_seaice / &
        ocean_phy_thermalcond_max,   &
        ocean_phy_thermalcond_seaice
  
     integer :: ierr
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("OCEAN_PHY_TC_seaice_setup",*) 'Setup'
  
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_ocean_phy_tc_seaice,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("OCEAN_PHY_TC_seaice_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("OCEAN_PHY_TC_seaice_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_PHY_TC_seaice. Check!'
        call prc_abort
     endif
     log_nml(param_ocean_phy_tc_seaice)
  
     return

References scale_io::io_fid_conf, ocean_phy_thermalcond_max, ocean_phy_thermalcond_seaice, and scale_prc::prc_abort().

Referenced by mod_ocean_driver::ocean_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ocean_phy_tc_seaice()

subroutine, public scale_ocean_phy_tc::ocean_phy_tc_seaice	(	integer, intent(in)	OIA,
		integer, intent(in)	OIS,
		integer, intent(in)	OIE,
		integer, intent(in)	OJA,
		integer, intent(in)	OJS,
		integer, intent(in)	OJE,
		real(rp), dimension(oia,oja), intent(in)	ICE_MASS,
		real(rp), dimension(oia,oja), intent(in)	ICE_FRAC,
		logical, dimension (oia,oja), intent(in)	mask,
		real(rp), dimension (oia,oja), intent(out)	TC_dz
	)

Definition at line 85 of file scale_ocean_phy_tc.F90.

     use scale_const, only: &
        eps => const_eps
     use scale_ocean_phy_ice_simple, only: &
        ocean_phy_ice_density
     implicit none
  
     integer,  intent(in)  :: OIA, OIS, OIE
     integer,  intent(in)  :: OJA, OJS, OJE
     real(RP), intent(in)  :: ICE_MASS(OIA,OJA) ! sea ice amount               [kg/m2]
     real(RP), intent(in)  :: ICE_FRAC(OIA,OJA) ! sea ice area fraction        [1]
     logical,  intent(in)  :: mask    (OIA,OJA)
     real(RP), intent(out) :: TC_dz   (OIA,OJA) ! thermal conductivity / depth [J/m2/s/K]
  
     real(RP) :: ice_depth
     integer  :: i, j
     !---------------------------------------------------------------------------
     !$acc data copyin(ICE_MASS,ICE_FRAC,mask) copyout(TC_dz)
  
     !$acc kernels
     !$omp parallel do private(ice_depth)
     do j = ojs, oje
     do i = ois, oie
        if ( mask(i,j) ) then
           ice_depth  = ice_mass(i,j) / ocean_phy_ice_density / max(ice_frac(i,j),eps)
  
           tc_dz(i,j) = ocean_phy_thermalcond_seaice / max(ice_depth*0.5_rp,eps) ! at the middle point of the layer
  
           tc_dz(i,j) = min( tc_dz(i,j), ocean_phy_thermalcond_max )
        end if
     enddo
     enddo
     !$acc end kernels
  
     !$acc end data
     return