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, 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.

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().

     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

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,
		real(rp), dimension (oia,oja), intent(out)	TC_dz
	)

Definition at line 84 of file scale_ocean_phy_tc.F90.

References scale_const::const_eps, scale_ocean_phy_ice_simple::ocean_phy_ice_density, ocean_phy_thermalcond_max, and ocean_phy_thermalcond_seaice.

Referenced by mod_ocean_driver::ocean_driver_calc_tendency().

     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]
     real(RP), intent(out) :: tc_dz   (oia,oja) ! thermal conductivity / depth [J/m2/s/K]
 
     real(RP) :: ice_depth
     integer  :: i, j
     !---------------------------------------------------------------------------
 
     do j = ojs, oje
     do i = ois, oie
        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 )
     enddo
     enddo
 
     return