Path: | ../src/main/arare_init-data.f90 |
Last Update: | Tue Jun 21 15:14:52 +0900 2011 |
Authors: | SUGIYAMA Ko-ichiro (杉山耕一朗), ODAKA Masatsugu (小高正嗣) |
Version: | $Id: arare_init-data.f90,v 1.15 2011-06-21 06:14:52 sugiyama Exp $ |
Tag Name: | $Name: arare5-20110623-2 $ |
Copyright: | Copyright (C) GFD Dennou Club, 2006. All rights reserved. |
License: | See COPYRIGHT |
非静力学モデル deepconv/arare.
* 方程式系は準圧縮系.
Main Program : |
非静力学モデル deepconv/arare.
This procedure input/output NAMELIST#initialdata_nml, NAMELIST#initialdata_basic_nml, NAMELIST#initialdata_disturb_gauss_nml, NAMELIST#initialdata_disturb_random_nml, NAMELIST#initialdata_disturb_dryreg_nml, NAMELIST#initialdata_disturb_moist_nml .
program ArareInitData ! !非静力学モデル deepconv/arare. ! !----- モジュール読み込み ------ !----- 型宣言, 文字列処理 ---- use dc_types, only : STRING, DP use dc_string, only : StoA use gtool_history, only : HistoryPut !----- メッセージ出力 ----- use dc_message, only: MessageNotify ! コマンドライン引数解釈 use argset, only : argset_init ! MPI の初期化 ! Initialize MPI wrapper ! use mpi_wrapper, only : MPIWrapperInit, MPIWrapperFinalize, myrank !----- 管理モジュール ----- ! 入出力ファイル名管理モジュール use fileset, only : fileset_init ! 格子点管理モジュール use gridset, only : gridset_init, imin, imax, jmin, jmax, kmin, kmax, nx, ny, nz, ncmax ! 基本場設定モジュール use basicset, only : basicset_init use axesset, only : axesset_init use constants, only : constants_init, molwtdry, Grav, TempSfc, PressSfc, pressbasis, gasrdry, cpdry, cvdry ! 湿潤ルーチン設定モジュール use composition, only: composition_init, molwtwet, SpcWetMolFr, GasNum, CondNum, IdxCC, IdxCG, SpcWetID ! 化学計算 use ChemCalc, only: ChemCalc_init, xyz_SvapPress ! 境界条件適用モジュール use setmargin use axesset, only: z_Z, z_dz !Z 方向の格子点間隔 use namelist_util, only: NmlutilInit, namelist_filename use eccm, only: ECCM_Dry, ECCM_Wet use initialdata_disturb, only: initialdata_disturb_random, initialdata_disturb_gaussXZ, initialdata_disturb_gaussXY, initialdata_disturb_gaussXYZ, initialdata_disturb_dryreg, initialdata_disturb_moist use initialdata_yamasaki1983, only: initialdata_yamasaki1983_basic use initialdata_takemi2007, only: initialdata_takemi2007_init, initialdata_takemi2007_basic, initialdata_takemi2007_wind use initialdata_toon2002, only: initialdata_toon2002_basic !----- 入出力モジュール ----- ! リスタートファイル入出力モジュール use RestartFileIO, only : ReStartFileIO_Init, ReStartFileIO_Finalize, rstat use Arare4FileIO, only : Arare4FileIO_Init, Arare4FileIO_Var_Get, Arare4FileIO_BZ_Get, Arare4FileIO_MMC_Var_Get, Arare4FileIO_MMC_BZ_Get !暗黙の型宣言禁止 implicit none !内部変数 real(DP), allocatable :: pyz_VelX(:,:,:) real(DP), allocatable :: xqz_VelY(:,:,:) real(DP), allocatable :: xyr_VelZ(:,:,:) real(DP), allocatable :: xyz_Exner(:,:,:) real(DP), allocatable :: xyz_PTemp(:,:,:) real(DP), allocatable :: xyz_Km(:,:,:) real(DP), allocatable :: xyz_Kh(:,:,:) real(DP), allocatable :: xyz_CDens(:,:,:) real(DP), allocatable :: xyzf_QMix(:,:,:,:) real(DP), allocatable :: xyz_DensBZ(:,:,:) real(DP), allocatable :: xyz_PressBZ(:,:,:) real(DP), allocatable :: xyz_ExnerBZ(:,:,:) real(DP), allocatable :: xyz_TempBZ(:,:,:) real(DP), allocatable :: xyz_PTempBZ(:,:,:) real(DP), allocatable :: xyz_VelSoundBZ(:,:,:) real(DP), allocatable :: xyzf_QMixBZ(:,:,:,:) real(DP), allocatable :: xyz_EffMolWtBZ(:,:,:) real(DP), allocatable :: z_TempBZ(:) real(DP), allocatable :: z_PressBZ(:) real(DP), allocatable :: xyzf_MolFr(:,:,:,:) real(DP), allocatable :: zf_MolFr(:,:) real(DP), allocatable :: xyzf_QMixDivMolWt(:,:,:,:) real(DP), allocatable :: xyzf_HumBZ(:,:,:,:) real(DP) :: Time integer :: i, j, k, s !変数定義 ! Moist 用 real(DP) :: Humidity = 0.0d0 !相対湿度 ! Gauss 用 real(DP) :: PTempMax = 0.0d0 !最大値 real(DP) :: ExnerMax = 0.0d0 !最大値 real(DP) :: QMixMax = 0.0d0 !最大値 real(DP) :: Xc = 0.0d0 !擾乱の中心位置(X方向) real(DP) :: Yc = 0.0d0 !擾乱の中心位置(Y方向) real(DP) :: Zc = 0.0d0 !擾乱の中心位置(鉛直方向) real(DP) :: Xr = 0.0d0 !擾乱の半径(X方向) real(DP) :: Yr = 0.0d0 !擾乱の半径(Y方向) real(DP) :: Zr = 0.0d0 !擾乱の半径(鉛直方向) ! Therma-Random 用 real(DP) :: Zpos = 0.0d0 !擾乱の Z 座標 [m] ! Dryreg 用 real(DP) :: XposMin = 0.0d0 !乾燥域の X 座標 [m] real(DP) :: YposMin = 0.0d0 !乾燥域の Y 座標 [m] real(DP) :: ZposMin = 0.0d0 !乾燥域の Z 座標 [m] real(DP) :: XposMax = 0.0d0 !乾燥域の X 座標 [m] real(DP) :: YposMax = 0.0d0 !乾燥域の Y 座標 [m] real(DP) :: ZposMax = 0.0d0 !乾燥域の Z 座標 [m] real(DP) :: TempTr = 100.0d0 !対流圏界面の温度 [k] real(DP) :: Dhight = 10.0d3 !重み関数のパラメータ [m] integer :: IDBasic = 0 integer, parameter :: IDBasicArare4 = 1 integer, parameter :: IDBasicArare4mmc = 2 integer, parameter :: IDBasicYamasaki1983 = 3 integer, parameter :: IDBasicTakemi2007 = 4 integer, parameter :: IDBasicIsoThermal = 5 integer, parameter :: IDBasicDry = 6 integer, parameter :: IDBasicMoist = 7 integer, parameter :: IDBasicToon2002 = 8 integer :: IDDisturbPTemp = 0 integer, parameter :: IDDisturbPTempGaussXY = 1 integer, parameter :: IDDisturbPTempGaussXZ = 2 integer, parameter :: IDDisturbPTempGaussXYZ = 3 integer, parameter :: IDDisturbPTempRandom = 4 integer :: IDDisturbExner = 0 integer, parameter :: IDDisturbExnerGaussXY = 1 integer, parameter :: IDDisturbExnerGaussXZ = 2 integer, parameter :: IDDisturbExnerGaussXYZ = 3 integer :: IDDisturbQMix = 0 integer, parameter :: IDDisturbQMixGaussXY = 1 integer, parameter :: IDDisturbQMixGaussXZ = 2 integer, parameter :: IDDisturbQMixGaussXYZ = 3 integer, parameter :: IDDisturbQMixDryreg = 4 integer, parameter :: IDDisturbQMixMoist = 5 integer :: IDDisturbWind = 0 integer, parameter :: IDDisturbWindTakemi2007= 1 integer :: IDDisturb = 0 integer, parameter :: IDDisturbArare4 = 1 integer, parameter :: IDDisturbArare4mmc = 2 !------------------------------------------ ! 初期化手続き ; Initialize procedure ! call MainInit Time = 0.0d0 call MessageNotify( "M", "main", "Making Initial data...." ) !--------------------------------------------------------------- ! 基本場を作成する. ! select case(IDBasic) case (IDBasicArare4) ! deepconv/arare4 のヒストリーファイルから初期値を読み込む ! call MessageNotify( "M", "main", "Making Initial data using arare4 files ...." ) call Arare4fileio_BZ_Get( xyz_PressBZ, xyz_ExnerBZ, xyz_TempBZ, xyz_PTempBZ, xyz_DensBZ, xyz_VelSoundBZ, xyzf_QMixBZ, xyz_EffMolWtBZ ) case (IDBasicArare4mmc) ! deepconv/arare4 火星版のヒストリーファイルから初期値を読み込む ! call MessageNotify( "M", "main", "Making Initial data using arare4 files ...." ) call Arare4fileio_MMC_BZ_Get( xyz_PressBZ, xyz_ExnerBZ, xyz_TempBZ, xyz_PTempBZ, xyz_DensBZ, xyz_VelSoundBZ, xyzf_QMixBZ, xyz_EffMolWtBZ ) case(IDBasicIsoThermal) ! 等温/重力考慮しない場合 ! call MessageNotify( "M", "main", "Making Initial data (basic) named IsoThermal..." ) z_TempBZ = TempSfc z_PressBZ = PressSfc zf_MolFr = 0.0d0 ! 残りの基本場の値を決める call DetermineBZ case(IDBasicDry) ! 乾燥断熱的な初期場 ! call MessageNotify( "M", "main", "Making Initial data (basic) named DRY..." ) call eccm_dry( SpcWetMolFr(1:ncmax), Humidity, z_TempBZ, z_PressBZ, zf_MolFr ) if (minval(z_TempBZ(1:nz)) < TempTr) then call initialdata_basic_strat( z_TempBZ, z_PressBZ ) !(inout) end if ! 残りの基本場の値を決める call DetermineBZ ! デバッグ済んでいない ! ! case(IDBasicMoist) ! ! 湿潤断熱的な初期場 ! ! ! call MessageNotify( "M", "main", "Making Initial data (basic) named MOIST..." ) ! call eccm_wet( SpcWetMolFr(1:ncmax), Humidity, z_TempBZ, z_PressBZ, zf_MolFr ) ! call initialdata_basic_strat( z_TempBZ, z_PressBZ ) !(inout) ! ! ! 残りの基本場の値を決める ! call DetermineBZ case(IDBasicToon2002) ! 火星用 ! call MessageNotify( "M", "main", "Making Initial data (basic) named Toon et al. 2002..." ) call initialdata_toon2002_basic( z_TempBZ, z_PressBZ ) !(out) ! call initialdata_basic_strat( z_TempBZ, z_PressBZ ) !(inout) ! 残りの基本場の値を決める call DetermineBZ case(IDBasicTakemi2007) ! Takemi(2007)の基本場を使用する場合 ! call MessageNotify( "M", "main", "Making Initial data named Takemi2007..." ) call initialdata_takemi2007_basic( z_TempBZ, z_PressBZ, zf_MolFr ) ! 残りの基本場の値を決める call DetermineBZ case(IDBasicYamasaki1983) ! Yamasaki(1983)の温度と相対湿度の観測値を使用する場合 ! call MessageNotify( "M", "main", "Making Initial data named Yamasaki1983..." ) call initialdata_yamasaki1983_basic( z_TempBZ, z_PressBZ, zf_MolFr ) ! 残りの基本場の値を決める call DetermineBZ end select !--------------------------------------------------------- ! 擾乱場の値を決める ! ! 温位: 擾乱場の値を決める ! select case(IDDisturbPTemp) case(IDDisturbPTempGaussXY) ! 温位: ガウシアンな分布を与える (XY) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXY..." ) call initialdata_disturb_gaussXY(PTempMax, Xc, Xr, Yc, Yr, xyz_PTemp) call SetMargin_xyz( xyz_PTemp ) case(IDDisturbPTempGaussXZ) ! 温位: ガウシアンな分布を与える (XZ) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXZ..." ) call initialdata_disturb_gaussXZ(PTempMax, Xc, Xr, Zc, Zr, xyz_PTemp) call SetMargin_xyz( xyz_PTemp ) case(IDDisturbPTempGaussXYZ) ! 温位: ガウシアンな分布を与える (XYZ) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXYZ..." ) call initialdata_disturb_gaussXYZ(PTempMax, Xc, Xr, Yc, Yr, Zc, Zr, xyz_PTemp) call SetMargin_xyz( xyz_PTemp ) case(IDDisturbPTempRandom) ! 温位: 指定された高度にランダムな分布を与える. ! call MessageNotify( "M", "main", "Making Initial data (disturb) named random..." ) call initialdata_disturb_random(PTempMax, Zpos, xyz_PTemp) call SetMargin_xyz( xyz_PTemp ) end select ! エクスナー関数: 擾乱場の値を決める ! select case(IDDisturbExner) case(IDDisturbExnerGaussXY) ! エクスナー関数: ガウシアンな分布を与える (XY) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXY..." ) call initialdata_disturb_gaussXY(ExnerMax, Xc, Xr, Yc, Yr, xyz_Exner) call SetMargin_xyz( xyz_Exner ) case(IDDisturbExnerGaussXZ) ! エクスナー関数: ガウシアンな分布を与える (XZ) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXZ..." ) call initialdata_disturb_gaussXZ(ExnerMax, Xc, Xr, Zc, Zr, xyz_Exner) call SetMargin_xyz( xyz_Exner ) case(IDDisturbExnerGaussXYZ) ! エクスナー関数: ガウシアンな分布を与える (XYZ) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXYZ..." ) call initialdata_disturb_gaussXYZ(ExnerMax, Xc, Xr, Yc, Yr, Zc, Zr, xyz_Exner) call SetMargin_xyz( xyz_Exner ) end select ! 混合比: 擾乱場の値を決める ! select case(IDDisturbQMix) case(IDDisturbQMixGaussXY) ! 混合比: ガウシアンな分布を与える (XY) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXY..." ) call initialdata_disturb_gaussXY(QMixMax, Xc, Xr, Yc, Yr, xyzf_QMix(:,:,:,1)) call SetMargin_xyzf( xyzf_QMix ) case(IDDisturbQMixGaussXZ) ! 混合比: ガウシアンな分布を与える (XZ) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXZ..." ) call initialdata_disturb_gaussXZ(QMixMax, Xc, Xr, Zc, Zr, xyzf_QMix(:,:,:,1)) call SetMargin_xyzf( xyzf_QMix ) case(IDDisturbQMixGaussXYZ) ! 混合比: ガウシアンな分布を与える (XYZ) ! call MessageNotify( "M", "main", "Making Initial data (disturb) named GaussXYZ..." ) call initialdata_disturb_gaussXYZ(QMixMax, Xc, Xr, Yc, Yr, Zc, Zr, xyzf_QMix(:,:,:,1)) call SetMargin_xyzf( xyzf_QMix ) case(IDDisturbQMixDryreg) ! 混合比: 乾燥領域を作る ! call MessageNotify( "M", "main", "Making Initial data (disturb) named Dryreg..." ) call initialdata_disturb_dryreg( XposMin, XposMax, YposMin, YposMax, ZposMin, ZposMax, xyzf_QMix) call SetMargin_xyzf( xyzf_QMix ) case(IDDisturbQMixMoist) ! 混合比: 湿潤断熱的な分布を与える ! call MessageNotify( "M", "main", "Making Initial data (disturb) named MOIST..." ) call initialdata_disturb_moist(Humidity, xyzf_QMix) call SetMargin_xyzf( xyzf_QMix ) end select ! 混合比: 速度場の値を決める ! select case(IDDisturbWind) case (IDDisturbWindTakemi2007) ! Takemi (2007) の風速分布 ! call MessageNotify( "M", "main", "Making Initial data (disturb) named Wind takemi2007..." ) call initialdata_takemi2007_wind( pyz_VelX ) call SetMargin_pyz( pyz_VelX ) end select ! 一括して擾乱場を与える場合 ! select case(IDDisturb) case (IDDisturbArare4) call MessageNotify( "M", "main", "Making Initial data using arare4 files ...." ) call Arare4fileio_Var_Get( pyz_VelX, xqz_VelY, xyr_VelZ, xyz_PTemp, xyz_Exner, xyzf_QMix, xyz_Km, xyz_Kh, xyz_CDens ) ! (out) case (IDDisturbArare4mmc) call MessageNotify( "M", "main", "Making Initial data using arare4 files ...." ) call Arare4fileio_MMC_Var_Get( pyz_VelX, xqz_VelY, xyr_VelZ, xyz_PTemp, xyz_Exner, xyzf_QMix, xyz_Km, xyz_Kh, xyz_CDens ) ! (out) end select !------------------------------------------------ ! ファイル出力 ! call MessageNotify( "M", "main", "Output variables into netCDF file..." ) ! リスタートファイル作成. 基本場と擾乱場を出力. ! call ReStartFileIO_Init call HistoryPut( 't', Time, rstat) call HistoryPut( 'VelX', pyz_VelX, rstat) call HistoryPut( 'VelY', xqz_VelY, rstat) call HistoryPut( 'VelZ', xyr_VelZ, rstat) call HistoryPut( 'Exner', xyz_Exner, rstat) call HistoryPut( 'PTemp', xyz_PTemp, rstat) call HistoryPut( 'Km', xyz_Km, rstat) call HistoryPut( 'Kh', xyz_Kh, rstat) call HistoryPut( 'QMix', xyzf_QMix, rstat) call HistoryPut( 'CDens', xyz_CDens, rstat) ! 基本場のファイル出力 ! call HistoryPut( 'DensBZ', xyz_DensBZ , rstat) call HistoryPut( 'ExnerBZ', xyz_ExnerBZ , rstat) call HistoryPut( 'PTempBZ', xyz_PTempBZ , rstat) call HistoryPut( 'VelSoundBZ', xyz_VelSoundBZ, rstat) call HistoryPut( 'TempBZ', xyz_TempBZ , rstat) call HistoryPut( 'PressBZ', xyz_PressBZ , rstat) call HistoryPut( 'QMixBZ', xyzf_QMixBZ , rstat) call HistoryPut( 'EffMolWtBZ', xyz_EffMolWtBZ, rstat) call HistoryPut( 'HumBZ', xyzf_HumBZ , rstat) call ReStartFileIO_Finalize !---------------------------------------------------- ! MPI END ! call MPIWrapperFinalize contains subroutine ArareAlloc ! !初期化として, 配列を定義し, 値としてゼロを代入する. ! !暗黙の型宣言禁止 implicit none !配列割り当て allocate( pyz_VelX(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xqz_VelY(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyr_VelZ(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyz_Exner(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyz_PTemp(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyz_Km(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyz_Kh(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyz_CDens(imin:imax, jmin:jmax, kmin:kmax) ) allocate( xyzf_QMix(imin:imax, jmin:jmax, kmin:kmax, ncmax) ) allocate( xyz_DensBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( xyz_PressBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( xyz_ExnerBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( xyz_TempBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( xyz_PTempBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( xyz_VelSoundBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( xyzf_QMixBZ(imin:imax,jmin:jmax,kmin:kmax,ncmax)) allocate( xyz_EffMolWtBZ(imin:imax,jmin:jmax,kmin:kmax)) allocate( z_TempBZ(kmin:kmax)) allocate( z_PressBZ(kmin:kmax)) allocate( xyzf_MolFr(imin:imax, jmin:jmax, kmin:kmax, ncmax)) allocate( zf_MolFr(kmin:kmax, ncmax)) allocate( xyzf_QMixDivMolWt(imin:imax,jmin:jmax,kmin:kmax, ncmax)) allocate( xyzf_HumBZ(imin:imax,jmin:jmax,kmin:kmax, ncmax)) pyz_VelX = 0.0d0 xqz_VelY = 0.0d0 xyr_VelZ = 0.0d0 xyz_PTemp = 0.0d0 xyz_Exner = 0.0d0 xyz_Km = 0.0d0 xyz_Kh = 0.0d0 xyz_CDens = 0.0d0 xyzf_QMix = 0.0d0 xyz_DensBZ = 0.0d0 xyz_PressBZ = 0.0d0 xyz_ExnerBZ = 0.0d0 xyz_TempBZ = 0.0d0 xyz_PTempBZ = 0.0d0 xyz_VelSoundBZ = 0.0d0 xyzf_QMixBZ = 0.0d0 xyz_EffMolWtBZ = 0.0d0 z_TempBZ = 0.0d0 z_PressBZ = 0.0d0 xyzf_MolFr = 0.0d0 zf_MolFr = 0.0d0 xyzf_QMixDivMolWt = 0.0d0 xyzf_HumBZ = 0.0d0 end subroutine ArareAlloc subroutine DetermineBZ !--------------------------------------------------------------- ! 温度・圧力 ! ! 3 次元配列に格納 do k = 1, nz do j = 1, ny do i = 1, nx xyz_TempBZ(i,j,k) = z_TempBZ(k) xyz_PressBZ(i,j,k) = z_PressBZ(k) end do end do end do !のり代の値を決める ! call SetMargin_xyz( xyz_TempBZ ) call SetMargin_xyz( xyz_PressBZ) !--------------------------------------------------------------- ! 混合比 ! !水平方向には一様 do k = 1, nz do j = 1, ny do i = 1, nx xyzf_MolFr(i,j,k,:) = zf_MolFr(k,:) end do end do end do !のり代の値を決める ! call SetMargin_xyzf( xyzf_MolFr ) !気相のモル比を混合比に変換 do s = 1, ncmax xyzf_QMixBZ(:,:,:,s) = xyzf_MolFr(:,:,:,s) * MolWtWet(s) / MolWtDry end do ! !値が小さくなりすぎないように最低値を与える ! where (xyzf_QMixBZ <= 1.0d-20 ) ! xyzf_QMixBZ = 1.0d-20 ! end where !のり代の値を決める ! call SetMargin_xyzf( xyzf_QMixBZ ) !--------------------------------------------------------------- ! 分子量の効果 ! do s = 1, ncmax xyzf_QMixDivMolWt(:,:,:,s) = xyzf_QMixBZ(:,:,:,s) / MolWtWet(s) end do xyz_EffMolWtBZ = (1.0d0 + sum(xyzf_QMixBZ,4) ) / ( MolWtDry * ((1.0d0 / MolWtDry) + sum(xyzf_QMixDivMolWt,4)) ) !のり代の値を決める ! call SetMargin_xyz( xyz_EffMolWtBZ ) !--------------------------------------------------------------- ! 温位 ! xyz_PTempBZ = xyz_TempBZ * (PressBasis / xyz_PressBZ) ** (GasRDry / CpDry) !のり代の値を決める ! call SetMargin_xyz( xyz_PTempBZ ) !--------------------------------------------------------------- ! エクスナー関数 ! xyz_ExnerBZ = xyz_TempBZ / xyz_PTempBZ !のり代の値を決める ! call SetMargin_xyz( xyz_ExnerBZ ) !--------------------------------------------------------------- ! 密度 ! xyz_DensBZ = PressBasis * (xyz_ExnerBZ ** (CvDry / GasRDry)) / (GasRDry * xyz_PTempBZ / xyz_EffMolWtBZ) !のり代の値を決める ! call SetMargin_xyz( xyz_DensBZ ) !--------------------------------------------------------------- ! 音速 ! xyz_VelSoundBZ = sqrt ( CpDry * GasRDry * xyz_ExnerBZ * xyz_PTempBZ / (CvDry * xyz_EffMolWtBZ) ) !のり代の値を決める ! call SetMargin_xyz( xyz_VelSoundBZ ) !--------------------------------------------------------------- ! 湿度 ! if (CondNum >= 1) then do s = 1, CondNum xyzf_HumBZ(:,:,:,IdxCG(s)) = xyz_PressBZ * xyzf_MolFr(:,:,:,IdxCG(s)) / xyz_SvapPress(SpcWetID(IdxCC(s)), xyz_TempBZ) * 100.0d0 end do end if !---------------------------------------------------------- ! BasicSet モジュールに値を設定 ! call BasicSet_Init( xyz_PressBZ, xyz_ExnerBZ, xyz_TempBZ, xyz_PTempBZ, xyz_DensBZ, xyz_VelSoundBZ, xyzf_QMixBZ, xyz_EffMolWtBZ ) end subroutine DetermineBZ subroutine Initialdata_init use dc_iounit, only : FileOpen implicit none integer :: unit !装置番号 character(STRING) :: FlagBasic = "" character(STRING) :: FlagDisturbPTemp = "" character(STRING) :: FlagDisturbExner = "" character(STRING) :: FlagDisturbQMix = "" character(STRING) :: FlagDisturbWind = "" character(STRING) :: FlagDisturb = "" NAMELIST /initialdata_nml/ FlagBasic, FlagDisturb, FlagDisturbPTemp, FlagDisturbExner, FlagDisturbQMix, FlagDisturbWind NAMELIST /initialdata_basic_nml/ Humidity, TempTr, DHight NAMELIST /initialdata_disturb_gauss_nml/ PTempMax, ExnerMax, QMixMax, Xc, Xr, Yc, Yr, Zc, Zr NAMELIST /initialdata_disturb_random_nml/ PTempMax, ExnerMax, QMixMax, Zpos NAMELIST /initialdata_disturb_dryreg_nml/ XposMin, YposMin, ZposMin, XposMax, YposMax, ZposMax NAMELIST /initialdata_disturb_moist_nml/ Humidity call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_nml) close(unit) if (FlagBasic == "Arare4") then IDBasic = IDBasicArare4 elseif(FlagBasic == "Arare4mmc") then IDBasic = IDBasicArare4mmc elseif(FlagBasic == "Yamasaki1983") then IDBasic = IDBasicYamasaki1983 elseif(FlagBasic == "Takemi2007") then IDBasic = IDBasicTakemi2007 elseif(FlagBasic == "IsoThermal") then IDBasic = IDBasicIsoThermal elseif(FlagBasic == "Dry" ) then IDBasic = IDBasicDry elseif(FlagBasic == "Moist" ) then IDBasic = IDBasicMoist elseif(FlagBasic == "Toon2002" ) then IDBasic = IDBasicToon2002 end if if(FlagDisturbPTemp == "GaussXY" ) then IDDisturbPTemp = IDDisturbPTempGaussXY elseif(FlagDisturbPTemp == "GaussXZ" ) then IDDisturbPTemp = IDDisturbPTempGaussXZ elseif(FlagDisturbPTemp == "GaussXYZ") then IDDisturbPTemp = IDDisturbPTempGaussXYZ elseif(FlagDisturbPTemp == "Random" ) then IDDisturbPTemp = IDDisturbPTempRandom end if if(FlagDisturbExner == "GaussXY" ) then IDDisturbExner = IDDisturbExnerGaussXY elseif(FlagDisturbExner == "GaussXZ" ) then IDDisturbExner = IDDisturbExnerGaussXZ elseif(FlagDisturbExner == "GaussXYZ") then IDDisturbExner = IDDisturbExnerGaussXYZ end if if(FlagDisturbQMix == "GaussXY" ) then IDDisturbQMix = IDDisturbQMixGaussXY elseif(FlagDisturbQMix == "GaussXZ" ) then IDDisturbQMix = IDDisturbQMixGaussXZ elseif(FlagDisturbQMix == "GaussXYZ" ) then IDDisturbQMix = IDDisturbQMixGaussXYZ elseif(FlagDisturbQMix == "Dryreg" ) then IDDisturbQMix = IDDisturbQMixDryreg elseif(FlagDisturbQMix == "Moist" ) then IDDisturbQMix = IDDisturbQMixMoist end if if (FlagDisturbWind == "Takemi2007") then IDDisturbWind = IDDisturbWindTakemi2007 end if if(FlagDisturb == "Arare4" ) then IDDisturb = IDDisturbArare4 elseif(FlagDisturb == "Arare4mmc" ) then IDDisturb = IDDisturbArare4mmc end if !---------------------------------------------------- ! モジュールの初期化 ! ! Yamasaki, シンプル擾乱場は, 初期化する必要ない. ! select case (IDDisturbPTemp) case (IDDisturbPTempGaussXY:IDDisturbPTempGaussXYZ) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_disturb_gauss_nml) close(unit) case (IDDisturbPTempRandom) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_disturb_random_nml) close(unit) end select select case (IDDisturbExner) case (IDDisturbExnerGaussXY:IDDisturbExnerGaussXYZ) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_disturb_gauss_nml) close(unit) end select select case (IDDisturbQMix) case (IDDisturbQMixGaussXY:IDDisturbQMixGaussXYZ) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_disturb_gauss_nml) close(unit) case(IDDisturbQMixDryreg) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_disturb_dryreg_nml) close(unit) case(IDDisturbQMixMoist) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_disturb_moist_nml) close(unit) end select select case (IDBasic) case (IDBasicDry, IDBasicMoist) call FileOpen(unit, file=namelist_filename, mode='r') read(unit, NML=initialdata_basic_nml) close(unit) end select if (FlagBasic == "Takemi2007" .OR. FlagDisturbWind == "Takemi2007") then call initialdata_takemi2007_init end if if ( FlagBasic == "Arare4" .OR. FlagDisturb == "Arare4" .OR. FlagBasic == "Arare4mmc" .OR. FlagDisturb == "Arare4mmc") then call arare4fileio_init end if if (myrank == 0) then call MessageNotify( "M", "main", "FlagBasic = %c", c1=trim(FlagBasic)) call MessageNotify( "M", "main", "IDBasic = %d", i=(/IDBasic/)) call MessageNotify( "M", "main", "FlagDisturbPTemp = %c", c1=trim(FlagDisturbPTemp)) call MessageNotify( "M", "main", "IDDisturbPTemp = %d", i=(/IDDisturbPTemp/)) call MessageNotify( "M", "main", "FlagDisturbExner = %c", c1=trim(FlagDisturbExner)) call MessageNotify( "M", "main", "IDDisturbExner = %d", i=(/IDDisturbExner/)) call MessageNotify( "M", "main", "FlagDisturbQMix = %c", c1=trim(FlagDisturbQMix)) call MessageNotify( "M", "main", "IDDisturbQMix = %d", i=(/IDDisturbQMix/)) call MessageNotify( "M", "main", "FlagDisturb = %c", c1=trim(FlagDisturb)) call MessageNotify( "M", "main", "IDDisturb = %d", i=(/IDDisturb/)) end if end subroutine Initialdata_init !!!-------------------------------------------------------------- subroutine initialdata_basic_strat(z_TempBZ, z_PressBZ) ! implicit none real(DP), intent(inout) :: z_TempBZ(kmin:kmax) real(DP), intent(inout) :: z_PressBZ(kmin:kmax) real(DP) :: HeightTr real(DP) :: z_DTempDZ(kmin:kmax) real(DP) :: DTempDZ real(DP) :: Weight integer :: k do k = 1, kmax-1 z_DTempDZ(k) = (z_TempBZ(k) - z_TempBZ(k-1)) / z_dz(k-1) end do ! 対流圏界面 ! HeightTr = minval(z_z(1:nz), 1, z_TempBZ(1:nz) < TempTr) ! 対流圏界面より上の扱い ! 圏界面より上は等温大気とする. 等温位大気から等温大気への遷移は ! tanh を用いてなめらかにつなぐ do k = 2, kmax !重みつけの関数を用意. tanh を用いる Weight = ( tanh( (z_Z(k) - HeightTr ) / Dhight ) + 1.0d0 ) * 5.0d-1 !仮値として温度を計算する. 圏界面より上では TempTr の等温大気に近づける z_TempBZ(k) = z_TempBZ(k) * ( 1.0d0 - Weight ) + TempTr * Weight !温度減率が断熱温度減率より小さくならないように DTempDZ = max( z_DTempDZ(k), (z_TempBZ(k) - z_TempBZ(k-1)) / z_dz(k-1) ) !基本場の温度を決める z_TempBZ(k) = z_TempBZ(k-1) + DTempDZ * z_dz(k-1) !圧力を静水圧平衡から計算 z_PressBZ(k) = z_PressBZ(k-1) * ( ( z_TempBZ(k-1) / z_TempBZ(k) ) ** (Grav / ( DTempDZ * GasRDry ) ) ) end do end subroutine initialdata_basic_strat subroutine MainInit implicit none character(STRING) :: cfgfile = "" ! MPI ! call MPIWrapperInit !NAMELIST ファイル名の読み込み ! call argset_init(cfgfile) !(out) ! NAMELIST ファイル名のモジュールへの登録 ! Loading NAMELIST file. ! call NmlutilInit( cfgfile ) !(in) !格子点情報の初期化 ! NAMELIST から情報を得て, 格子点を計算する ! call gridset_init ! 化学計算ルーチンの初期化 ! Initialization of chemical routines. ! call chemcalc_init ! 軸の情報の初期化 ! Initialization of axis variables. ! call axesset_init ! 定数の情報の初期化 ! Initialization of constant variables. ! call constants_init ! 湿潤過程共有変数の初期化 ! Initialization of common variables for moist process. ! call composition_init ! I/O ファイル名の初期化 ! Initialization of output file name. ! call fileset_init ! nml から情報を取り出す (内部サブルーチン) ! call InitialData_init !内部変数の初期化. とりあえずゼロを入れて値を確定させておく. ! call ArareAlloc end subroutine MainInit end program ArareInitData