Class | subsurface_diffusion_heat |
In: |
surface_flux/subsurface_diffusion_heat.f90
|
Note that Japanese and English are described in parallel.
PhyImplTendency : | 時間変化率の計算 |
———— : | ———— |
PhyImplTendency : | Calculate tendency |
Subroutine : | |||
xy_GroundTempFlux(0:imax-1, 1:jmax) : | real(DP), intent(in )
| ||
xy_SurfTemp(0:imax-1, 1:jmax) : | real(DP), intent(in )
| ||
xyz_SoilTemp(0:imax-1, 1:jmax, 1:kslmax) : | real(DP), intent(in )
| ||
xyr_SoilTempTransCoef(0:imax-1, 1:jmax, 0:kslmax) : | real(DP), intent(out)
| ||
xyr_SoilTempFlux(0:imax-1, 1:jmax, 0:kslmax) : | real(DP), intent(out)
|
時間変化率の計算を行います.
Calculate tendencies.
subroutine SubsurfaceDiffusion( xy_GroundTempFlux, xy_SurfTemp, xyz_SoilTemp, xyr_SoilTempTransCoef, xyr_SoilTempFlux ) ! ! 時間変化率の計算を行います. ! ! Calculate tendencies. ! ! モジュール引用 ; USE statements ! ! 時刻管理 ! Time control ! use timeset, only: TimesetClockStart, TimesetClockStop ! 座標データ設定 ! Axes data settings ! use axesset, only: r_SSDepth, z_SSDepth ! subsurface grid at midpoint of layer ! 時刻管理 ! Time control ! ! 宣言文 ; Declaration statements ! implicit none real(DP), intent(in ):: xy_GroundTempFlux (0:imax-1, 1:jmax) ! 地中熱フラックス. ! Ground temperature flux real(DP), intent(in ):: xy_SurfTemp (0:imax-1, 1:jmax) ! 地表面温度. ! Surface temperature real(DP), intent(in ):: xyz_SoilTemp (0:imax-1, 1:jmax, 1:kslmax) ! 土壌温度 (K) ! Soil temperature (K) real(DP), intent(out):: xyr_SoilTempTransCoef (0:imax-1, 1:jmax, 0:kslmax) ! 輸送係数:土壌温度. ! Transfer coefficient: soil temperature real(DP), intent(out):: xyr_SoilTempFlux (0:imax-1, 1:jmax, 0:kslmax) ! 土壌の熱フラックス. ! Heat flux in sub-surface soil ! 作業変数 ! Work variables ! integer:: k ! 鉛直方向に回る DO ループ用作業変数 ! Work variables for DO loop in vertical direction ! 実行文 ; Executable statement ! ! 計算時間計測開始 ! Start measurement of computation time ! call TimesetClockStart( module_name ) ! 初期化 ! Initialization ! if ( .not. subsurface_diffusion_inited ) call SubsurfaceDiffusionInit ! 土壌温度計算用の輸送係数の計算 ! Calculate transfer coefficient for heat diffusion in the soil ! k = 0 if ( kslmax == 0 ) then ! This line is used when kslmax == 0, because z_SSDepth(k+1) does not exist. xyr_SoilTempTransCoef(:,:,k) = 0.0d0 else xyr_SoilTempTransCoef(:,:,k) = xy_SoilHeatDiffCoef(:,:) / ( z_SSDepth(k+1) - 0.0d0 ) end if do k = 1, kslmax-1 xyr_SoilTempTransCoef(:,:,k) = xy_SoilHeatDiffCoef(:,:) / ( z_SSDepth(k+1) - z_SSDepth(k) ) end do k = kslmax xyr_SoilTempTransCoef(:,:,k) = 0.0d0 ! 土壌中の熱フラックスの計算 ! Calculate heat flux in sub-surface soil ! k = 0 if ( kslmax == 0 ) then ! This line is used when kslmax == 0, because xyz_SoilTemp(:,:,k+1) does not exist. xyr_SoilTempFlux(:,:,k) = 0.0d0 else xyr_SoilTempFlux(:,:,k) = - xyr_SoilTempTransCoef(:,:,k) * ( xyz_SoilTemp(:,:,1) - xy_SurfTemp(:,:) ) end if do k = 1, kslmax-1 xyr_SoilTempFlux(:,:,k) = - xyr_SoilTempTransCoef(:,:,k) * ( xyz_SoilTemp(:,:,k+1) - xyz_SoilTemp(:,:,k) ) end do k = kslmax xyr_SoilTempFlux(:,:,k) = xy_GroundTempFlux ! 計算時間計測一時停止 ! Pause measurement of computation time ! call TimesetClockStop( module_name ) end subroutine SubsurfaceDiffusion
Variable : | |||
subsurface_diffusion_inited = .false. : | logical, save, public
|
Variable : | |||
xy_SoilHeatCap(:,:) : | real(DP), allocatable, save, public
|
Variable : | |||
xy_SoilHeatDiffCoef(:,:) : | real(DP), allocatable, save, public
|
Variable : | |||
xy_SoilRho(:,:) : | real(DP), allocatable, save, public
|
Variable : | |||
xy_SoilSpecHeat(:,:) : | real(DP), allocatable, save, public
|
Subroutine : |
subsurface_diffusion_heat モジュールの初期化を行います. NAMELIST#subsurface_diffusion_heat_nml の読み込みはこの手続きで行われます.
"subsurface_diffusion_heat" module is initialized. "NAMELIST#subsurface_diffusion_heat_nml" is loaded in this procedure.
subroutine SubsurfaceDiffusionInit ! ! subsurface_diffusion_heat モジュールの初期化を行います. ! NAMELIST#subsurface_diffusion_heat_nml の読み込みはこの手続きで行われます. ! ! "subsurface_diffusion_heat" module is initialized. ! "NAMELIST#subsurface_diffusion_heat_nml" is loaded in this procedure. ! ! モジュール引用 ; USE statements ! ! NAMELIST ファイル入力に関するユーティリティ ! Utilities for NAMELIST file input ! use namelist_util, only: namelist_filename, NmlutilMsg, NmlutilAryValid ! ファイル入出力補助 ! File I/O support ! use dc_iounit, only: FileOpen ! 種別型パラメタ ! Kind type parameter ! use dc_types, only: STDOUT ! 標準出力の装置番号. Unit number of standard output ! 文字列操作 ! Character handling ! use dc_string, only: StoA ! ヒストリデータ出力 ! History data output ! use gtool_historyauto, only: HistoryAutoAddVariable ! 宣言文 ; Declaration statements ! implicit none ! 作業変数 ! Work variables ! !!$ integer:: unit_nml ! NAMELIST ファイルオープン用装置番号. !!$ ! Unit number for NAMELIST file open !!$ integer:: iostat_nml ! NAMELIST 読み込み時の IOSTAT. !!$ ! IOSTAT of NAMELIST read ! NAMELIST 変数群 ! NAMELIST group name ! !!$ namelist /phy_implicit_nml/ !!$ ! !!$ ! デフォルト値については初期化手続 "phy_implicit#PhyImplInit" !!$ ! のソースコードを参照のこと. !!$ ! !!$ ! Refer to source codes in the initialization procedure !!$ ! "phy_implicit#PhyImplInit" for the default values. !!$ ! ! 実行文 ; Executable statement ! if ( subsurface_diffusion_inited ) return !!$ call InitCheck ! デフォルト値の設定 ! Default values settings ! ! NAMELIST の読み込み ! NAMELIST is input ! !!$ if ( trim(namelist_filename) /= '' ) then !!$ call FileOpen( unit_nml, & ! (out) !!$ & namelist_filename, mode = 'r' ) ! (in) !!$ !!$ rewind( unit_nml ) !!$ read( unit_nml, & ! (in) !!$ & nml = phy_implicit_nml, & ! (out) !!$ & iostat = iostat_nml ) ! (out) !!$ close( unit_nml ) !!$ !!$ call NmlutilMsg( iostat_nml, module_name ) ! (in) !!$ end if allocate( xy_SoilSpecHeat ( 0:imax-1, 1:jmax ) ) allocate( xy_SoilRho ( 0:imax-1, 1:jmax ) ) allocate( xy_SoilHeatCap ( 0:imax-1, 1:jmax ) ) allocate( xy_SoilHeatDiffCoef( 0:imax-1, 1:jmax ) ) xy_SoilSpecHeat = 1.0d100 ! This variable is not used. xy_SoilRho = 1.0d100 ! This variable is not used. xy_SoilHeatCap = 2.1d6 ! J m-3 K-1 ! Value of Clay for porosity f=0.4, volumetric wetness theta=0.2 in Table 12.3 by ! Hillel (2004). ! Note that the unit of Table 12.3 of Hillel (2004) would be wrong. Although the ! unit in the Table is wrong, the volumetric heat capacity of 2.1d6 J m-3 K-1 is ! within the range of typical value of it. xy_SoilHeatDiffCoef = 1.2d0 ! W m-1 K-1 ! Value of Clay for porosity f=0.4, volumetric wetness theta=0.2 in Table 12.3 by ! Hillel (2004). ! Reference ! ! Hillet, D., ! Introduction to Environmental Soil Physics, ! Elsevier Academic Press, pp494, 2004. ! 印字 ; Print ! call MessageNotify( 'M', module_name, '----- Initialization Messages -----' ) call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) ) subsurface_diffusion_inited = .true. end subroutine SubsurfaceDiffusionInit
Constant : | |||
module_name = ‘subsurface_diffusion_heat‘ : | character(*), parameter
|
Constant : | |||
version = ’$Name: dcpam5-20100224 $’ // ’$Id: subsurface_diffusion_heat.f90,v 1.3 2009-08-04 09:42:36 yot Exp $’ : | character(*), parameter
|