scale_land_phy.F90

Go to the documentation of this file.

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

  abstract interface
     subroutine lnd( &
           LAND_TEMP_t,       &
           LAND_WATER_t,      &
           LAND_TEMP,         &
           LAND_WATER,        &
           LAND_WaterLimit,   &
           LAND_ThermalCond,  &
           LAND_HeatCapacity, &
           LAND_WaterDiff,    &
           LAND_SFLX_GH,      &
           LAND_SFLX_prec,    &
           LAND_SFLX_evap,    &
           CDZ,               &
           dt                 )
       use scale_precision
       use scale_grid_index
       use scale_land_grid_index
       implicit none

       real(RP), intent(out) :: LAND_TEMP_t      (LKMAX,IA,JA)
       real(RP), intent(out) :: LAND_WATER_t     (LKMAX,IA,JA)

       real(RP), intent(in)  :: LAND_TEMP        (LKMAX,IA,JA)
       real(RP), intent(in)  :: LAND_WATER       (LKMAX,IA,JA)
       real(RP), intent(in)  :: LAND_WaterLimit  (IA,JA)
       real(RP), intent(in)  :: LAND_ThermalCond (IA,JA)
       real(RP), intent(in)  :: LAND_HeatCapacity(IA,JA)
       real(RP), intent(in)  :: LAND_WaterDiff   (IA,JA)
       real(RP), intent(in)  :: LAND_SFLX_GH     (IA,JA)
       real(RP), intent(in)  :: LAND_SFLX_prec   (IA,JA)
       real(RP), intent(in)  :: LAND_SFLX_evap   (IA,JA)
       real(RP), intent(in)  :: CDZ              (LKMAX)
       real(DP), intent(in)  :: dt
     end subroutine lnd
  end interface
  procedure(lnd), pointer :: land_phy => null()
  public :: land_phy

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine land_phy_setup( LAND_TYPE )
    use scale_process, only: &
       prc_mpistop
    use scale_land_phy_const, only: &
       land_phy_const_setup, &
       land_phy_const
    use scale_land_phy_slab, only: &
       land_phy_slab_setup, &
       land_phy_slab
    use scale_land_phy_matsiro, only: &
       land_phy_matsiro_setup, &
       land_phy_matsiro
    implicit none

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

    select case( land_type )
    case( 'CONST' )
       call land_phy_const_setup( land_type )
       land_phy => land_phy_const
    case( 'SLAB', 'THIN-SLAB', 'THICK-SLAB' )
       call land_phy_slab_setup( land_type )
       land_phy => land_phy_slab
    case( 'MATSIRO' )
       call land_phy_matsiro_setup( land_type )
       land_phy => land_phy_matsiro
    case default
       write(*,*) 'xxx invalid Land type(', trim(land_type), '). CHECK!'
       call prc_mpistop
    end select

    return
  end subroutine land_phy_setup

end module scale_land_phy