scale_ocean_phy.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_ocean_phy
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_phy_setup

  abstract interface
     subroutine ocn( &
           OCEAN_TEMP_t,    &
           OCEAN_TEMP,      &
           OCEAN_SFLX_WH,   &
           OCEAN_SFLX_prec, &
           OCEAN_SFLX_evap, &
           dt               )
       use scale_precision
       use scale_grid_index
       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
     end subroutine ocn
  end interface
  procedure(ocn), pointer :: ocean_phy => null()
  public :: ocean_phy

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine ocean_phy_setup( OCEAN_TYPE )
    use scale_process, only: &
       prc_mpistop
    use scale_ocean_phy_slab, only: &
       ocean_phy_slab_setup, &
       ocean_phy_slab
    use scale_ocean_phy_file, only: &
       ocean_phy_file_setup, &
       ocean_phy_file
    implicit none

    character(len=*), intent(in) :: OCEAN_TYPE
    !---------------------------------------------------------------------------

    select case ( ocean_type )
    case ( 'CONST' )
       call ocean_phy_slab_setup( ocean_type )
       ocean_phy => ocean_phy_slab
    case ( 'SLAB' )
       call ocean_phy_slab_setup( ocean_type )
       ocean_phy => ocean_phy_slab
    case ( 'FILE' )
       call ocean_phy_file_setup( ocean_type )
       ocean_phy => ocean_phy_file
    case default
       write(*,*) 'xxx invalid Ocean type(', trim(ocean_type), '). CHECK!'
       call prc_mpistop
    end select

    return
  end subroutine ocean_phy_setup

end module scale_ocean_phy