| Class | HistoryFileIO |
| In: |
io/historyfileio.f90
|
ファイル出力. 長い時間ステップの値を出力.
| Subroutine : |
ヒストリファイルのクローズ
subroutine HistoryFile_Close
!
!ヒストリファイルのクローズ
!
!暗黙の型宣言禁止
implicit none
!ファイルを閉じる
call HistoryClose
end subroutine HistoryFile_Close
| Subroutine : |
ヒストリファイルの定義
subroutine HistoryFile_Open( )
!
!ヒストリファイルの定義
!
!暗黙の型宣言禁止
implicit none
!変数定義
real(4) :: SpcID(SpcNum)
integer :: s
SpcID = 0.0d0
do s = 1, SpcNum
SpcID(s) = real(s, 4)
end do
!-----------------------------------------------------------
! ヒストリー作成
!-----------------------------------------------------------
call HistoryCreate( file = HistoryFile, title = exptitle, source = expsrc, institution = expinst, dims=(/'x','z','s','t'/), dimsizes=(/FileNX, FileNZ, SpcNum, 0/), longnames=(/'X-coordinate', 'Z-coordinate', 'Species Num ', 'Time '/), units=(/'m','m','1','s'/), origin=0.0, interval=0.0 )
!-----------------------------------------------------------
! 軸の出力
!-----------------------------------------------------------
call HistoryPut('x', s_X( FileXMin: FileXMax ) )
call HistoryPut('z', s_Z( FileZMin: FileZMax ) )
call HistoryPut('s', real(SpcID, 4))
!-----------------------------------------------------------
! 予報変数の出力
!-----------------------------------------------------------
!無次元圧力の擾乱
call HistoryAddVariable( varname='Exner', dims=(/'x','z','t'/), longname='disturbunce of nondimensional pressure', units=' ', xtype='double' )
!温位の擾乱
call HistoryAddVariable( varname='PotTemp', dims=(/'x','z','t'/), longname='disturbunce of potential temperature', units='K', xtype='double' )
!水平速度
call HistoryAddVariable( varname='VelX', dims=(/'x','z','t'/), longname='zonal velocity', units='m s|-1"', xtype='double' )
!鉛直速度
call HistoryAddVariable( varname='VelZ', dims=(/'x','z','t'/), longname='vertical velocity', units='m s|-1"', xtype='double' )
!渦粘性係数
call HistoryAddVariable( varname='Km', dims=(/'x','z','t'/), longname='Km', units='1', xtype='double' )
!渦粘性係数
call HistoryAddVariable( varname='Kh', dims=(/'x','z','t'/), longname='Kh', units='1', xtype='double' )
!混合比
call HistoryAddVariable( varname='MixRt', dims=(/'x','z','s','t'/), longname='Mixing Ratio', units='kg kg|-1"', xtype='double' )
!----------------------------------------------------------------
! 静的安定度
!----------------------------------------------------------------
call HistoryAddVariable( varname='Stab', dims=(/'z','t'/), longname='Static Stability', units='s|-2"', xtype='double' )
call HistoryAddVariable( varname='StabTemp', dims=(/'z','t'/), longname='Static Stability contributed by temperature', units='s|-2"', xtype='double' )
call HistoryAddVariable( varname='StabMolWt', dims=(/'z','t'/), longname='Static Stability contributed by Molecular Weight', units='s|-2"', xtype='double' )
!----------------------------------------------------------------
! 温位の時間変化
!----------------------------------------------------------------
call HistoryAddVariable( varname='PotTempAdv', dims=(/'z','t'/), longname='Advection term of potential temperature', units='K day|-1"', xtype='double' )
call HistoryAddVariable( varname='PotTempTurb', dims=(/'z','t'/), longname='Turbulence term of potential temperature', units='K day|-1"', xtype='double' )
call HistoryAddVariable( varname='PotTempDisp', dims=(/'z','t'/), longname='Dissipation term of potential temperature', units='K day|-1"', xtype='double' )
call HistoryAddVariable( varname='PotTempRad', dims=(/'z','t'/), longname='Radiation term of potential temperature', units='K day|-1"', xtype='double' )
call HistoryAddVariable( varname='PotTempDiff', dims=(/'z','t'/), longname='Numerical diffusion term of potential temperature', units='K day|-1"', xtype='double' )
call HistoryAddVariable( varname='PotTempCond', dims=(/'z','t'/), longname='Latent heat term of potential temperature', units='K day|-1"', xtype='double' )
call HistoryAddVariable( varname='PotTempDamp', dims=(/'z','t'/), longname='Newtonian Cooling term of potential temperature', units='K day|-1"', xtype='double' )
!-----------------------------------------------------------
! 基本場の出力
!-----------------------------------------------------------
!無次元圧力の基本場
call HistoryAddVariable( varname='ExnerBasicZ', dims=(/'x','z'/), longname='nondimensional pressure', units='1', xtype='double' )
!温位の基本場
call HistoryAddVariable( varname='PotTempBasicZ', dims=(/'x','z'/), longname='potential temperature', units='K', xtype='double' )
!密度の基本場
call HistoryAddVariable( varname='DensBasicZ', dims=(/'x','z'/), longname='density', units='Kg/m^3', xtype='double' )
!音波速度の基本場
call HistoryAddVariable( varname='VelSoundBasicZ', dims=(/'x','z'/), longname='sound velocity', units='m/s|2', xtype='double' )
!温度の基本場
call HistoryAddVariable( varname='TempBasicZ', dims=(/'x','z'/), longname='Temperature of basic state', units='K', xtype='double' )
!圧力の基本場
call HistoryAddVariable( varname='PressBasicZ', dims=(/'x','z'/), longname='Pressure of basic state', units='Pa', xtype='double' )
!水蒸気混合比の基本場
call HistoryAddVariable( varname='MixRtBasicZ', dims=(/'x','z','s'/), longname='Mixing ratio of Condensible volatiles', units='kg/kg', xtype='double' )
!分子量効果
call HistoryAddVariable( varname='EffMolWtBasicZ', dims=(/'x','z'/), longname='Effect of Mole Weight', units='1', xtype='double' )
!-------------------------------------------------------------
! 基本場のファイル出力
!-------------------------------------------------------------
call HistoryPut( 'DensBasicZ', xz_DensBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
call HistoryPut( 'ExnerBasicZ', xz_ExnerBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
call HistoryPut( 'PotTempBasicZ', xz_PotTempBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
call HistoryPut( 'VelSoundBasicZ', xz_VelSoundBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
call HistoryPut( 'TempBasicZ', xz_TempBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
call HistoryPut( 'PressBasicZ', xz_PressBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
call HistoryPut( 'MixRtBasicZ', xza_MixRtBasicZ(FileXMin:FileXMax, FileZMin:FileZMax, 1:SpcNum))
call HistoryPut( 'EffMolWtBasicZ', xz_EffMolWtBasicZ(FileXMin:FileXMax, FileZMin:FileZMax))
end subroutine HistoryFile_Open
| Subroutine : | |
| Time : | real(8), intent(in) |
| xz_PotTemp(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in) |
| xz_Exner(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in) |
| pz_VelX(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in) |
| xr_VelZ(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in) |
| xza_MixRt(DimXMin:DimXMax, DimZMin:DimZMax, SpcNum) : | real(8), intent(in) |
| xz_Km(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in) |
| xz_Kh(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in) |
予報変数のヒストリファイルへの出力. 出力時には半格子点の位置でプロット.
subroutine HistoryFile_OutPut( Time, xz_PotTemp, xz_Exner, pz_VelX, xr_VelZ, xza_MixRt, xz_Km, xz_Kh )
!
!予報変数のヒストリファイルへの出力. 出力時には半格子点の位置でプロット.
!
!モジュール読み込み
use average, only: xz_avr_pz, xz_avr_xr
!暗黙の型宣言禁止
implicit none
!変数定義
real(8), intent(in) :: Time
real(8), intent(in) :: pz_VelX(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(in) :: xr_VelZ(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(in) :: xz_Exner(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(in) :: xz_PotTemp(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(in) :: xz_Km(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(in) :: xz_Kh(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(in) :: xza_MixRt(DimXMin:DimXMax, DimZMin:DimZMax, SpcNum)
real(8) :: xz_VelX(DimXMin:DimXMax, DimZMin:DimZMax)
real(8) :: xz_VelZ(DimXMin:DimXMax, DimZMin:DimZMax)
!----------------------------------------------------------------
! 格子点位置を変換
!----------------------------------------------------------------
xz_VelX = xz_avr_pz( pz_VelX )
xz_VelZ = xz_avr_xr( xr_VelZ )
!----------------------------------------------------------------
! 値を出力
!----------------------------------------------------------------
call HistoryPut( 't', Time )
call HistoryPut('Exner', xz_Exner(FileXMin:FileXMax, FileZMin:FileZMax) )
call HistoryPut('PotTemp', xz_PotTemp(FileXMin:FileXMax, FileZMin:FileZMax) )
call HistoryPut( 'VelX', xz_VelX(FileXMin:FileXMax, FileZMin:FileZMax) )
call HistoryPut( 'VelZ', xz_VelZ(FileXMin:FileXMax, FileZMin:FileZMax) )
call HistoryPut( 'Km', xz_Km(FileXMin:FileXMax, FileZMin:FileZMax) )
call HistoryPut( 'Kh', xz_Kh(FileXMin:FileXMax, FileZMin:FileZMax) )
call HistoryPut( 'MixRt', xza_MixRt(FileXMin:FileXMax, FileZMin:FileZMax, 1:SpcNum) )
!----------------------------------------------------------------
! 解析値を出力
!----------------------------------------------------------------
call StoreMeanX()
call HistoryPut( 'PotTempAdv', z_Adv(FileZMin:FileZMax) )
call HistoryPut( 'PotTempTurb', z_Turb(FileZMin:FileZMax) )
call HistoryPut( 'PotTempDisp', z_Disp(FileZMin:FileZMax) )
call HistoryPut( 'PotTempDiff', z_Diff(FileZMin:FileZMax) )
call HistoryPut( 'PotTempRad', z_Rad(FileZMin:FileZMax) )
call HistoryPut( 'PotTempCond', z_Cond(FileZMin:FileZMax) )
call HistoryPut( 'PotTempDamp', z_Damp(FileZMin:FileZMax) )
call HistoryPut( 'Stab', z_Stab(FileZMin:FileZMax))
call HistoryPut( 'StabTemp', z_StabTemp(FileZMin:FileZMax))
call HistoryPut( 'StabMolWt', z_StabMolWt(FileZMin:FileZMax))
end subroutine HistoryFile_OutPut