subroutine physics_ground( xy_SurfTemp , y_Lat ) ! (in) 緯度座標
!!!! AGCM5 で他に設定できるらしい地表面情報 !!!!!
! xy_SurfTemp ,
! xy_SurfHeatCapa ,
! xy_SurfAlbedo ,
! xy_SurfHumidCoeff ,
! xy_SurfRoughLength ,
! REAL GRZSD ( IDIM*JDIM ) !" 地表高度分散 = 0.0
! INTEGER IDSFC ( IDIM*JDIM ) !" 地表状態 = 0.0
! REAL GFLUX ( IDIM*JDIM ) !" 地中熱フラックス = 0.0
!==== Dependency
use type_mod, only: REKIND, DBKIND, INTKIND, TOKEN, STRING
use grid_3d_mod, only: im, jm, km
use dc_trace, only: SetDebug, BeginSub, EndSub, DbgMessage, DataDump
implicit none
!==== Input
!
real(DBKIND), intent(in) :: y_Lat(:) ! (in) 緯度座標
!==== Output
!
real(DBKIND), intent(out) :: xy_SurfTemp(im,jm) ! (out) 地表面温度
!----- 作業用内部変数 -----
character(STRING), parameter:: subname = "physics_ground"
real(DBKIND) :: xy_SeaSurfaceTemp(im,jm) ! 海表面温度
! Hosaka et al. (1998) SST 用パラメータ
real(DBKIND), parameter :: TEQ = 302.0d0
real(DBKIND), parameter :: ALAT0 = 0.0d0
real(DBKIND), parameter :: ALAT1 = 30.0d0
real(DBKIND), parameter :: ALPHA = 60.0d0
real(DBKIND), parameter :: BETA = 32.0d0
real(DBKIND), parameter :: GAMMA = 0.0d0
real(DBKIND), parameter :: ALACON = 7.0d0
continue
!----------------------------------------------------------------
! 開始処理
!----------------------------------------------------------------
call BeginSub(subname)
!----------------------------------------------------------------
! 地表面温度の設定
!----------------------------------------------------------------
!----- 海面温度の設定 -----
! 本当はファイルから読み出すべきだが暫定的処置.
! Hosaka et al. (1998) の SST 分布を与える.
call mksst( xy_SeaSurfaceTemp , y_Lat , TEQ , ALAT0 , ALAT1 , ALPHA , BETA , GAMMA , ALACON )
!----- 地表面温度を海面温度で置き換え-----
xy_SurfTemp = xy_SeaSurfaceTemp
!----------------------------------------------------------------
! 終了処理
!----------------------------------------------------------------
call EndSub(subname)
end subroutine physics_ground