scale_land_phy_matsiro.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_land_phy_matsiro
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_debug
  use scale_grid_index
  use scale_land_grid_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: land_phy_matsiro_setup
  public :: land_phy_matsiro

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine land_phy_matsiro_setup( LAND_TYPE )
    use scale_process, only: &
       prc_mpistop
    implicit none

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

    logical :: dummy

    namelist / param_land_matsiro / &
       dummy

    integer :: ierr
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[MATSIRO] / Categ[LAND PHY] / Origin[SCALElib]'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_land_matsiro,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_LAND_MATSIRO. Check!'
       call prc_mpistop
    endif
    if( io_lnml ) write(io_fid_log,nml=param_land_matsiro)

    return
  end subroutine land_phy_matsiro_setup

  !-----------------------------------------------------------------------------
  subroutine land_phy_matsiro( &
       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_grid_index
    implicit none

    ! arguments
    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
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Land step: Matsiro'

    land_temp_t(:,:,:) = 0.0_rp
    land_water_t(:,:,:) = 0.0_rp

    return
  end subroutine land_phy_matsiro

end module scale_land_phy_matsiro