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