module OCEAN / Physics Slab model More...

Functions/Subroutines
subroutine, public	ocean_phy_slab_setup (OCEAN_TYPE)
	Setup. More...

subroutine, public	ocean_phy_slab (OCEAN_TEMP_t, OCEAN_TEMP, OCEAN_SFLX_WH, OCEAN_SFLX_prec, OCEAN_SFLX_evap, dt)
	Slab ocean model. More...

Detailed Description

module OCEAN / Physics Slab model

Description: ocean physics module, slab model

Author: Team SCALE

NAMELIST

PARAM_OCEAN_PHY_SLAB

name type default value comment

OCEAN_PHY_SLAB_DEPTH real(RP) 10.0_RP water depth of slab ocean [m]

History Output: No history output

Function/Subroutine Documentation

◆ ocean_phy_slab_setup()

subroutine, public scale_ocean_phy_slab::ocean_phy_slab_setup ( character(len=*), intent(in) OCEAN_TYPE )

Setup.

Definition at line 54 of file scale_ocean_phy_slab.F90.

References scale_const::const_cl, scale_const::const_dwatr, scale_grid_index::ia, scale_grid_index::ie, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, scale_stdio::io_lnml, scale_grid_index::is, scale_grid_index::ja, scale_grid_index::je, scale_grid_index::js, scale_landuse::landuse_fact_ocean, and scale_process::prc_mpistop().

Referenced by scale_ocean_phy::ocean_phy_setup().

     use scale_process, only: &
        prc_mpistop
     use scale_const, only: &
        dwatr => const_dwatr, &
        cl    => const_cl
     use scale_landuse, only: &
        landuse_fact_ocean
     implicit none
 
     character(len=*), intent(in) :: ocean_type
 
     namelist / param_ocean_phy_slab / &
        ocean_phy_slab_depth
 
     integer :: i, j
     integer :: ierr
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[SLAB] / Categ[OCEAN PHY] / Origin[SCALElib]'
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_ocean_phy_slab,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        write(*,*) 'xxx Not appropriate names in namelist PARAM_OCEAN_PHY_SLAB. Check!'
        call prc_mpistop
     endif
     if( io_lnml ) write(io_fid_log,nml=param_ocean_phy_slab)
 
     if( ocean_type == 'CONST' ) then
        ocean_phy_slab_fixedsst = .true.
     else if( ocean_type == 'SLAB' ) then
        ocean_phy_slab_fixedsst = .false.
     else
        write(*,*) 'xxx wrong OCEAN_TYPE. Check!'
        call prc_mpistop
     end if
 
     ocean_phy_slab_heatcapacity = dwatr * cl * ocean_phy_slab_depth
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** Prevent SST change?          : ', ocean_phy_slab_fixedsst
     if( io_l ) write(io_fid_log,*) '*** Slab ocean depth [m]         : ', ocean_phy_slab_depth
     if( io_l ) write(io_fid_log,*) '*** Ocean heat capacity [J/K/m2] : ', ocean_phy_slab_heatcapacity
 
     ! judge to run slab ocean model
     allocate( is_ocn(ia,ja) )
 
     do j = js, je
     do i = is, ie
       if( landuse_fact_ocean(i,j) > 0.0_rp ) then
         is_ocn(i,j) = .true.
       else
         is_ocn(i,j) = .false.
       end if
     end do
     end do
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ocean_phy_slab()

subroutine, public scale_ocean_phy_slab::ocean_phy_slab	(	real(rp), dimension (ia,ja), intent(out)	OCEAN_TEMP_t,
		real(rp), dimension (ia,ja), intent(in)	OCEAN_TEMP,
		real(rp), dimension (ia,ja), intent(in)	OCEAN_SFLX_WH,
		real(rp), dimension(ia,ja), intent(in)	OCEAN_SFLX_prec,
		real(rp), dimension(ia,ja), intent(in)	OCEAN_SFLX_evap,
		real(dp), intent(in)	dt
	)

Slab ocean model.

Definition at line 127 of file scale_ocean_phy_slab.F90.

References scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, and scale_grid_index::js.