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_const, only: &
       ocean_phy_const_setup, &
       ocean_phy_const
    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_const_setup( ocean_type )
       ocean_phy => ocean_phy_const
    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