Class HistoryFileIO
In: io/historyfileio.f90

ファイル出力. 長い時間ステップの値を出力.

Methods

Included Modules

gt4_history gridset basicset fileset StoreSet average

Public Instance methods

Subroutine :

ヒストリファイルのクローズ

[Source]

  subroutine HistoryFile_Close
    !
    !ヒストリファイルのクローズ
    !
    
    !暗黙の型宣言禁止
    implicit none

    !変数定義
    integer :: s ! ループ添字

    !ファイルを閉じる
    do s = 1, FileNum
      call HistoryClose(gt_hist(s))
    end do

  end subroutine HistoryFile_Close
Subroutine :

ヒストリファイルの定義

[Source]

  subroutine HistoryFile_Open( )
    !
    !ヒストリファイルの定義
    !
    
    !暗黙の型宣言禁止
    implicit none
    
    !変数定義
    integer :: s ! ループ添字

    !-----------------------------------------------------------
    ! ヒストリー作成
    !-----------------------------------------------------------
    do s = 1, FileNum
      call HistoryCreate( file = HistoryFile(s), title = exptitle, source = expsrc, institution = expinst, dims=(/'x','z','t'/), dimsizes=(/FileNX, FileNZ, 0/), longnames=(/'X-coordinate', 'Z-coordinate', 'Time        '/), units=(/'m','m','s'/), origin=0.0, interval=0.0, history=gt_hist(s) )

      call HistoryPut('x', s_X( FileXMin: FileXMax ), gt_hist(s) )
      call HistoryPut('z', s_Z( FileZMin: FileZMax ), gt_hist(s) )
!      call HistoryPut('z', xz_PressBasicZ( 1, FileZMin: FileZMax ), gt_hist(s) )

    end do

    !-----------------------------------------------------------  
    ! 予報変数の出力
    !-----------------------------------------------------------  
    !無次元圧力の擾乱
    call HistoryAddVariable( varname='Exner', dims=(/'x','z','t'/), longname='disturbunce of nondimensional pressure', units=' ', xtype='double', history=gt_hist(1) )
    
    !温位の擾乱
    call HistoryAddVariable( varname='PotTemp', dims=(/'x','z','t'/), longname='disturbunce of potential temperature', units='K', xtype='double', history=gt_hist(2) )

    !水平速度
    call HistoryAddVariable( varname='VelX', dims=(/'x','z','t'/), longname='zonal velocity', units='m s|-1"', xtype='double', history=gt_hist(3) )

    !鉛直速度
    call HistoryAddVariable( varname='VelZ', dims=(/'x','z','t'/), longname='vertical velocity', units='m s|-1"', xtype='double', history=gt_hist(4) )

    !渦粘性係数(運動量)
    call HistoryAddVariable( varname='Km', dims=(/'x','z','t'/), longname='turbulet diffusion coefficient', units='1', xtype='double', history=gt_hist(5) )
  
    !渦粘性係数(熱)
    call HistoryAddVariable( varname='Kh', dims=(/'x','z','t'/), longname='turbulet diffusion coefficient for heat', units='1', xtype='double', history=gt_hist(6) )
  
    !混合比
    do s = 1, SpcNum
      call HistoryAddVariable( varname=trim(SpcWetSymbol(s)), dims=(/'x','z','t'/), longname=trim(SpcWetSymbol(s))//' Mixing Ratio', units='kg kg|-1"', xtype='double', history=gt_hist(8+s))

      call HistoryAddVariable( varname=trim(SpcWetSymbol(s))//'_Mean', dims=(/'z','t'/), longname='Horiizontal Mean of ' //trim(SpcWetSymbol(s))//' Mixing Ratio', units='kg kg|-1"', xtype='double', history=gt_hist(8+s) ) 
    end do
    
    !----------------------------------------------------------------
    ! Stored Mixing Ratio
    !----------------------------------------------------------------

    !----------------------------------------------------------------
    ! 静的安定度
    !----------------------------------------------------------------
    call HistoryAddVariable( varname='Stab', dims=(/'z','t'/), longname='Static Stability', units='s|-2"', xtype='double', history=gt_hist(8) ) 

    call HistoryAddVariable( varname='StabTemp', dims=(/'z','t'/), longname='Static Stability contributed by temperature', units='s|-2"', xtype='double', history=gt_hist(8) ) 

    call HistoryAddVariable( varname='StabMolWt', dims=(/'z','t'/), longname='Static Stability contributed by Molecular Weight', units='s|-2"', xtype='double', history=gt_hist(8) ) 

    !----------------------------------------------------------------
    ! 温位の時間変化
    !----------------------------------------------------------------
    call HistoryAddVariable( varname='PotTempAdv', dims=(/'z','t'/), longname='Advection term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 

    call HistoryAddVariable( varname='PotTempTurb', dims=(/'z','t'/), longname='Turbulence term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 

    call HistoryAddVariable( varname='PotTempDisp', dims=(/'z','t'/), longname='Dissipation term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 

    call HistoryAddVariable( varname='PotTempRad', dims=(/'z','t'/), longname='Radiation term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 
    
    call HistoryAddVariable( varname='PotTempDiff', dims=(/'z','t'/), longname='Numerical diffusion term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 
    
    call HistoryAddVariable( varname='PotTempCond', dims=(/'z','t'/), longname='Latent heat term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 

    call HistoryAddVariable( varname='PotTempDamp', dims=(/'z','t'/), longname='Newtonian Cooling term of potential temperature', units='K day|-1"', xtype='double', history=gt_hist(8) ) 

    !-----------------------------------------------------------  
    ! 基本場の出力
    !-----------------------------------------------------------  
    !無次元圧力の基本場
    call HistoryAddVariable( varname='ExnerBasicZ', dims=(/'x','z'/), longname='nondimensional pressure', units='1', xtype='double', history=gt_hist(7) ) 
    
    !温位の基本場
    call HistoryAddVariable( varname='PotTempBasicZ', dims=(/'x','z'/), longname='potential temperature', units='K', xtype='double', history=gt_hist(7) ) 
    
    !密度の基本場
    call HistoryAddVariable( varname='DensBasicZ', dims=(/'x','z'/), longname='density', units='Kg/m^3', xtype='double', history=gt_hist(7) ) 
    
    !音波速度の基本場
    call HistoryAddVariable( varname='VelSoundBasicZ', dims=(/'x','z'/), longname='sound velocity', units='m/s|2', xtype='double', history=gt_hist(7) ) 
    
    !温度の基本場
    call HistoryAddVariable( varname='TempBasicZ', dims=(/'x','z'/), longname='Temperature of basic state', units='K', xtype='double', history=gt_hist(7) ) 
    
    !圧力の基本場
    call HistoryAddVariable( varname='PressBasicZ', dims=(/'x','z'/), longname='Pressure of basic state', units='Pa', xtype='double', history=gt_hist(7) ) 
    
    !水蒸気混合比の基本場
    do s = 1, SpcNum
      call HistoryAddVariable( varname=trim(SpcWetSymbol(s))//'BasicZ', dims=(/'x','z'/), longname=trim(SpcWetSymbol(s))//' Mixing Ratio of basic state', units='kg/kg', xtype='double', history=gt_hist(7) ) 
    end do

    !分子量効果
    call HistoryAddVariable( varname='EffMolWtBasicZ', dims=(/'x','z'/), longname='Effect of Mole Weight', units='1', xtype='double', history=gt_hist(7) ) 

    !-------------------------------------------------------------
    ! 基本場のファイル出力
    !-------------------------------------------------------------
    call HistoryPut( 'DensBasicZ', xz_DensBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )

    call HistoryPut( 'ExnerBasicZ', xz_ExnerBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )

    call HistoryPut( 'PotTempBasicZ', xz_PotTempBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )

    call HistoryPut( 'VelSoundBasicZ', xz_VelSoundBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )

    call HistoryPut( 'TempBasicZ', xz_TempBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )

    call HistoryPut( 'PressBasicZ', xz_PressBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )

    do s = 1, SpcNum
      call HistoryPut( trim(SpcWetSymbol(s))//'BasicZ', xza_MixRtBasicZ(FileXMin:FileXMax, FileZMin:FileZMax, s), gt_hist(7) )
    end do

    call HistoryPut( 'EffMolWtBasicZ', xz_EffMolWtBasicZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(7) )
    
  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)

予報変数のヒストリファイルへの出力. 出力時には半格子点の位置でプロット.

[Source]

  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)
    integer             :: s

    !----------------------------------------------------------------
    ! 格子点位置を変換
    !----------------------------------------------------------------
    xz_VelX = xz_avr_pz( pz_VelX )
    xz_VelZ = xz_avr_xr( xr_VelZ )

    !----------------------------------------------------------------
    ! 値を出力
    !----------------------------------------------------------------

    do s = 1, 8 + SpcNum
      call HistoryPut( 't', Time, gt_hist(s) )
    end do
    
    call HistoryPut( 'Exner', xz_Exner(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(1) )

    call HistoryPut( 'PotTemp', xz_PotTemp(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(2) )

    call HistoryPut( 'VelX', xz_VelX(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(3) )
    
    call HistoryPut( 'VelZ', xz_VelZ(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(4) )

    call HistoryPut( 'Km', xz_Km(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(5) )
    
    call HistoryPut( 'Kh', xz_Kh(FileXMin:FileXMax, FileZMin:FileZMax), gt_hist(6) )
    
    do s = 1, SpcNum
      call HistoryPut( trim(SpcWetSymbol(s)), xza_MixRt(FileXMin:FileXMax, FileZMin:FileZMax, s), gt_hist(8+s) )
    end do

    !----------------------------------------------------------------
    ! 解析値を出力
    !----------------------------------------------------------------
    call StoreMeanX()

    do s = 1, SpcNum
      call HistoryPut( trim(SpcWetSymbol(s))//'_Mean', za_MixRt(FileZMin:FileZMax, s), gt_hist(8+s) )
    end do

    call HistoryPut( 'PotTempAdv', z_Adv(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'PotTempTurb', z_Turb(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'PotTempDisp', z_Disp(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'PotTempDiff', z_Diff(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'PotTempRad', z_Rad(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'PotTempCond', z_Cond(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'PotTempDamp', z_Damp(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'Stab', z_Stab(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'StabTemp', z_StabTemp(FileZMin:FileZMax), gt_hist(8) )

    call HistoryPut( 'StabMolWt', z_StabMolWt(FileZMin:FileZMax), gt_hist(8) )

  end subroutine HistoryFile_OutPut

[Validate]