scale_ocean_phy_slab.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_ocean_phy_slab
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_debug
  use scale_grid_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_phy_slab_setup
  public :: ocean_phy_slab

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  real(RP), private :: OCEAN_PHY_SLAB_DEPTH = 10.0_rp 
  real(RP), private :: OCEAN_PHY_SLAB_HeatCapacity

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine ocean_phy_slab_setup( OCEAN_TYPE )
    use scale_process, only: &
       prc_mpistop
    use scale_const, only: &
       dwatr => const_dwatr, &
       cl    => const_cl
    implicit none

    character(len=*), intent(in) :: ocean_type

    namelist / param_ocean_phy_slab / &
       ocean_phy_slab_depth

    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_nml ) write(io_fid_nml,nml=param_ocean_phy_slab)

    ocean_phy_slab_heatcapacity = dwatr * cl * ocean_phy_slab_depth

    if( io_l ) write(io_fid_log,*)
    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

    return
  end subroutine ocean_phy_slab_setup

  !-----------------------------------------------------------------------------
  subroutine ocean_phy_slab( &
       OCEAN_TEMP_t,    &
       OCEAN_TEMP,      &
       OCEAN_SFLX_WH,   &
       OCEAN_SFLX_prec, &
       OCEAN_SFLX_evap, &
       dt               )
    use scale_landuse, only: &
       landuse_fact_ocean
    implicit none

    real(RP), intent(out) :: ocean_temp_t   (ia,ja)
    real(RP), intent(in)  :: ocean_temp     (ia,ja)
    real(RP), intent(in)  :: ocean_sflx_wh  (ia,ja)
    real(RP), intent(in)  :: ocean_sflx_prec(ia,ja)
    real(RP), intent(in)  :: ocean_sflx_evap(ia,ja)
    real(DP), intent(in)  :: dt

    integer :: i, j
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Ocean physics step: Slab'

    do j = js, je
    do i = is, ie
      if( landuse_fact_ocean(i,j) > 0.0_rp ) then
        ocean_temp_t(i,j) = - ocean_sflx_wh(i,j) / ocean_phy_slab_heatcapacity
      else
        ocean_temp_t(i,j) = 0.0_rp
      endif
    enddo
    enddo

    return
  end subroutine ocean_phy_slab

end module scale_ocean_phy_slab