Class | restart_surftemp_io |
In: |
io/restart_surftemp_io.f90
|
Note that Japanese and English are described in parallel.
地表面温度のリスタートデータもしくは初期値データの入出力を行います. 入力ファイル, 出力ファイル, データの出力の間隔は NAMELIST#restart_surftemp_io_nml で設定します.
リスタートデータの入力ファイルが指定されない場合, surface_data モジュールで生成された初期値データを取得します.
Restart data or initial data is input/output. Settings of input file, output file, and interval of data output is configured by "NAMELIST#restart_surftemp_io_nml".
If input file of restart data is not set, initial data is generated in "surface_data" module, and obtained data from the module.
RestartSurfTempOpen : | リスタートファイルのオープン |
RestartSurfTempOutput : | リスタートファイルへのデータ出力 |
RestartSurfTempClose : | リスタートファイルのクローズ |
RestartSurfTempGet : | リスタートファイルの入力 |
———— : | ———— |
RestartSurfTempOpen : | Open restart file |
RestartSurfTempOutput : | Data output to restart file |
RestartSurfTempClose : | Close restart file |
RestartSurfTempGet : | Input restart file |
Subroutine : |
リスタートデータファイル出力の終了処理を行います.
Terminate restart data files output.
subroutine RestartSurfTempClose ! ! リスタートデータファイル出力の終了処理を行います. ! ! Terminate restart data files output. ! モジュール引用 ; USE statements ! ! gtool4 データ出力 ! Gtool4 data output ! use gtool_history, only: HistoryClose ! 宣言文 ; Declaration statements ! ! 作業変数 ! Work variables ! ! 実行文 ; Executable statement ! if ( .not. restart_surftemp_opened ) return call HistoryClose( history = gthst_rst ) ! (inout) restart_surftemp_opened = .false. end subroutine RestartSurfTempClose
Subroutine : | |||
xy_SurfTemp(0:imax-1, 1:jmax) : | real(DP), intent(out)
| ||
xyz_SoilTemp(0:imax-1, 1:jmax, 1:kslmax) : | real(DP), intent(out) | ||
xyz_SOSeaIceTemp(0:imax-1, 1:jmax, 1:ksimax) : | real(DP), intent(out) | ||
xy_SOSeaIceMass(0:imax-1, 1:jmax) : | real(DP), intent(out) | ||
xy_SurfMajCompIceB(0:imax-1, 1:jmax) : | real(DP), intent(out) | ||
xy_SoilMoistB(0:imax-1, 1:jmax) : | real(DP), intent(out) | ||
xy_SurfSnowB(0:imax-1, 1:jmax) : | real(DP), intent(out) | ||
xy_SurfMajCompIceN(0:imax-1, 1:jmax) : | real(DP), intent(out) | ||
xy_SoilMoistN(0:imax-1, 1:jmax) : | real(DP), intent(out) | ||
xy_SurfSnowN(0:imax-1, 1:jmax) : | real(DP), intent(out) |
リスタートデータの入力を行います. リスタートデータが存在しない場合には, surface_data モジュールを用いて, 地表面温度データ生成を行います.
Input restart data. If restart data is not exist, surface temperature data is created by "surface_data".
subroutine RestartSurfTempGet( xy_SurfTemp, xyz_SoilTemp, xyz_SOSeaIceTemp, xy_SOSeaIceMass, xy_SurfMajCompIceB, xy_SoilMoistB, xy_SurfSnowB, xy_SurfMajCompIceN, xy_SoilMoistN, xy_SurfSnowN ) ! ! リスタートデータの入力を行います. ! リスタートデータが存在しない場合には, ! surface_data モジュールを用いて, 地表面温度データ生成を行います. ! ! Input restart data. ! If restart data is not exist, ! surface temperature data is created by "surface_data". ! モジュール引用 ; USE statements ! ! 地表面データ提供 ! Prepare surface data ! use surface_data, only: SetSurfData ! 時刻管理 ! Time control ! use timeset, only: RestartTime ! リスタート開始時刻. ! Retart time of calculation ! gtool4 データ入力 ! Gtool4 data input ! use gtool_history, only: HistoryGet ! 文字列操作 ! Character handling ! use dc_string, only: toChar ! NetCDF のラッパープログラム ! NetCDF wrapper ! use netcdf_wrapper, only: NWPresentAVarInFile ! 宣言文 ; Declaration statements ! real(DP), intent(out) :: xy_SurfTemp (0:imax-1, 1:jmax) ! 地表面温度. ! Surface temperature real(DP), intent(out) :: xyz_SoilTemp (0:imax-1, 1:jmax, 1:kslmax) real(DP), intent(out) :: xyz_SOSeaIceTemp (0:imax-1, 1:jmax, 1:ksimax) real(DP), intent(out) :: xy_SOSeaIceMass (0:imax-1, 1:jmax) real(DP), intent(out) :: xy_SurfMajCompIceB(0:imax-1, 1:jmax) real(DP), intent(out) :: xy_SoilMoistB (0:imax-1, 1:jmax) real(DP), intent(out) :: xy_SurfSnowB (0:imax-1, 1:jmax) real(DP), intent(out) :: xy_SurfMajCompIceN(0:imax-1, 1:jmax) real(DP), intent(out) :: xy_SoilMoistN (0:imax-1, 1:jmax) real(DP), intent(out) :: xy_SurfSnowN (0:imax-1, 1:jmax) ! 作業変数 ! Work variables ! character(TOKEN) :: time_range ! 時刻の指定. ! Specification of time logical :: flag_mpi_init integer :: k ! 実行文 ; Executable statement ! if ( .not. restart_surftemp_io_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if ! The variable name is fixed, 'SurfTemp' (yot, 2011/09/10). !!$ if ( trim(InputFile) == '' .or. trim(InputName) == '' ) then if ( trim(InputFile) == '' ) then ! データを surface_data モジュールから取得 ! Data is input from "surface_data" module ! call SetSurfData( xy_SurfTemp = xy_SurfTemp ) ! (out) optional do k = 1, kslmax xyz_SoilTemp(:,:,k) = xy_SurfTemp end do do k = 1, ksimax xyz_SOSeaIceTemp(:,:,k) = xy_SurfTemp end do xy_SOSeaIceMass = 0.0_DP xy_SurfMajCompIceB = 0.0_DP xy_SoilMoistB = 0.0_DP xy_SurfSnowB = 0.0_DP xy_SurfMajCompIceN = xy_SurfMajCompIceB xy_SoilMoistN = xy_SoilMoistB xy_SurfSnowN = xy_SurfSnowB else ! データを InputFile から取得 ! Data is input from InputFile ! ! 時刻情報の取得 ! Get time information ! time_range = time_name // '=' // toChar( RestartTime ) !!$ flag_mpi_init = .false. flag_mpi_init = .true. ! データ入力 ! Data input ! call HistoryGet( InputFile, 'SurfTemp', range = time_range, array = xy_SurfTemp, flag_mpi_split = flag_mpi_init ) if ( kslmax /= 0 ) then if ( NWPresentAVarInFile( InputFile, 'SoilTemp' ) ) then call HistoryGet( InputFile, 'SoilTemp', range = time_range, array = xyz_SoilTemp, flag_mpi_split = flag_mpi_init ) else do k = 1, kslmax xyz_SoilTemp(:,:,k) = xy_SurfTemp end do call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SoilTemp', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SoilTemp is assumed to be the same as SurfTemp.' ) end if end if if ( ksimax /= 0 ) then if ( NWPresentAVarInFile( InputFile, 'SOSeaIceTemp' ) ) then call HistoryGet( InputFile, 'SOSeaIceTemp', range = time_range, array = xyz_SOSeaIceTemp, flag_mpi_split = flag_mpi_init ) else do k = 1, ksimax xyz_SOSeaIceTemp(:,:,k) = xy_SurfTemp end do call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SOSeaIceTemp', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SOSeaIceTemp is assumed to be the same as SurfTemp.' ) end if end if if ( NWPresentAVarInFile( InputFile, 'SOSeaIceMass' ) ) then call HistoryGet( InputFile, 'SOSeaIceMass', range = time_range, array = xy_SOSeaIceMass, flag_mpi_split = flag_mpi_init ) else xy_SOSeaIceMass = 0.0_DP call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SOSeaIceMass', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SOSeaIceMass is assumed to be zero.' ) end if if ( NWPresentAVarInFile( InputFile, 'SurfMajCompIceB' ) ) then call HistoryGet( InputFile, 'SurfMajCompIceB', range = time_range, array = xy_SurfMajCompIceB, flag_mpi_split = flag_mpi_init ) else xy_SurfMajCompIceB = 0.0_DP call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SurfMajCompIceB', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SurfMajCompIceB is assumed to be zero.' ) end if if ( NWPresentAVarInFile( InputFile, 'SoilMoistB' ) ) then call HistoryGet( InputFile, 'SoilMoistB', range = time_range, array = xy_SoilMoistB, flag_mpi_split = flag_mpi_init ) else xy_SoilMoistB = 0.0_DP call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SoilMoistB', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SoilMoistB is assumed to be zero.' ) end if if ( NWPresentAVarInFile( InputFile, 'SurfSnowB' ) ) then call HistoryGet( InputFile, 'SurfSnowB', range = time_range, array = xy_SurfSnowB, flag_mpi_split = flag_mpi_init ) else xy_SurfSnowB = 0.0_DP call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SurfSnowB', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SurfSnowB is assumed to be zero.' ) end if if ( NWPresentAVarInFile( InputFile, 'SurfMajCompIceN' ) ) then call HistoryGet( InputFile, 'SurfMajCompIceN', range = time_range, array = xy_SurfMajCompIceN, flag_mpi_split = flag_mpi_init ) else xy_SurfMajCompIceN = xy_SurfMajCompIceB call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SurfMajCompIceN', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SurfMajCompIceN is assumed to be the same as SurfMajCompIceB.' ) end if if ( NWPresentAVarInFile( InputFile, 'SoilMoistN' ) ) then call HistoryGet( InputFile, 'SoilMoistN', range = time_range, array = xy_SoilMoistN, flag_mpi_split = flag_mpi_init ) else xy_SoilMoistN = xy_SoilMoistB call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SoilMoistN', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SoilMoistN is assumed to be the same as SoilMoistB.' ) end if if ( NWPresentAVarInFile( InputFile, 'SurfSnowN' ) ) then call HistoryGet( InputFile, 'SurfSnowN', range = time_range, array = xy_SurfSnowN, flag_mpi_split = flag_mpi_init ) else xy_SurfSnowN = xy_SurfSnowB call MessageNotify( 'M', module_name, 'Variable, %c, cannot be found in file, %c.', c1 = 'SurfSnowN', c2 = trim(InputFile) ) call MessageNotify( 'M', module_name, 'Values in SurfSnowN is assumed to be the same as SurfSnowB.' ) end if end if end subroutine RestartSurfTempGet
Subroutine : |
restart_surftemp_io モジュールの初期化を行います. NAMELIST#restart_surftemp_io_nml の読み込みはこの手続きで行われます.
"restart_surftemp_io" module is initialized. "NAMELIST#restart_surftemp_io_nml" is loaded in this procedure.
This procedure input/output NAMELIST#restart_surftemp_io_nml .
subroutine RestartSurfTempInit ! ! restart_surftemp_io モジュールの初期化を行います. ! NAMELIST#restart_surftemp_io_nml の読み込みはこの手続きで行われます. ! ! "restart_surftemp_io" module is initialized. ! "NAMELIST#restart_surftemp_io_nml" is loaded in this procedure. ! ! モジュール引用 ; USE statements ! ! リスタートデータ入出力 ! Restart data input/output ! use restart_file_io, only: restart_file_io_inited, RestartFileIntValue => IntValue, RestartFileIntUnit => IntUnit ! リスタートデータの出力間隔の単位. ! Unit for interval of restart data output ! 時刻管理 ! Time control ! use timeset, only: RestartTime ! リスタート開始時刻. ! Retart time of calculation ! NAMELIST ファイル入力に関するユーティリティ ! Utilities for NAMELIST file input ! use namelist_util, only: namelist_filename, NmlutilMsg ! 暦と日時の取り扱い ! Calendar and Date handler ! use dc_calendar, only: DCCalConvertByUnit ! ファイル入出力補助 ! File I/O support ! use dc_iounit, only: FileOpen ! 種別型パラメタ ! Kind type parameter ! use dc_types, only: STDOUT ! 標準出力の装置番号. Unit number of standard output ! 宣言文 ; Declaration statements ! ! 作業変数 ! 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 /restart_surftemp_io_nml/ InputFile, OutputFile, IntValue, IntUnit ! ! デフォルト値については初期化手続 "restart_surftemp_io#RestartSurfTempInit" ! のソースコードを参照のこと. ! ! Refer to source codes in the initialization procedure ! "restart_surftemp_io#RestartSurfTempInit" for the default values. ! ! 実行文 ; Executable statement ! if ( restart_surftemp_io_inited ) return ! デフォルト値の設定 ! Default values settings ! InputFile = '' ! The variable name is fixed, 'SurfTemp' (yot, 2011/09/10). !!$ InputName = 'SurfTemp' if ( .not. flag_init_data_save ) then OutputFile = 'rst_sst.nc' else OutputFile = 'sst.nc' end if if ( restart_file_io_inited ) then IntValue = RestartFileIntValue IntUnit = RestartFileIntUnit else IntValue = 365.0_DP IntUnit = 'day' end if ! NAMELIST の読み込み ! NAMELIST is input ! if ( trim(namelist_filename) /= '' ) then call FileOpen( unit_nml, namelist_filename, mode = 'r' ) ! (in) rewind( unit_nml ) read( unit_nml, nml = restart_surftemp_io_nml, iostat = iostat_nml ) ! (out) close( unit_nml ) call NmlutilMsg( iostat_nml, module_name ) ! (in) if ( iostat_nml == 0 ) write( STDOUT, nml = restart_surftemp_io_nml ) end if ! 出力時間間隔の設定 ! Configure time interval of output ! IntTime = DCCalConvertByUnit( IntValue, IntUnit, 'sec' ) ! (in) PrevOutputTime = RestartTime ! フラグの初期化 ! Initialize flag ! flag_output_end = .false. ! 印字 ; Print ! call MessageNotify( 'M', module_name, '----- Initialization Messages -----' ) call MessageNotify( 'M', module_name, 'Input:: ' ) call MessageNotify( 'M', module_name, ' InputFile = %c', c1 = trim(InputFile) ) ! The variable name is fixed, 'SurfTemp' (yot, 2011/09/10). !!$ call MessageNotify( 'M', module_name, ' InputName = %c', c1 = trim(InputName) ) call MessageNotify( 'M', module_name, 'Output:: ' ) call MessageNotify( 'M', module_name, ' OutputFile = %c', c1 = trim(OutputFile) ) call MessageNotify( 'M', module_name, ' IntTime = %f [%c]', d = (/ IntValue /), c1 = trim(IntUnit) ) call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) ) restart_surftemp_io_inited = .true. end subroutine RestartSurfTempInit
Subroutine : | |||
flag_init_data : | logical, intent(in), optional
|
リスタートファイルをオープンします.
A restart file is opened.
subroutine RestartSurfTempOpen( flag_init_data ) ! ! リスタートファイルをオープンします. ! ! A restart file is opened. ! ! モジュール引用 ; USE statements ! ! 出力ファイルの基本情報 ! Basic information for output files ! use fileset, only: FileTitle, FileSource, FileInstitution ! データファイルを最終的に変更した組織/個人. ! Institution or person that changes data files for the last time ! 物理・数学定数設定 ! Physical and mathematical constants settings ! use constants0, only: PI ! $ \pi $. ! 円周率. Circular constant ! 座標データ設定 ! Axes data settings ! use axesset, only: x_Lon, x_Lon_Weight, y_Lat, y_Lat_Weight, z_SSDepth, z_SIDepth ! 海氷の格子点の深さ ! sea ice grid at midpoint of layer ! 時刻管理 ! Time control ! use timeset, only: RestartTime ! リスタート開始時刻. ! Retart time of calculation ! gtool4 データ出力 ! Gtool4 data output ! use gtool_history, only: HistoryCreate, HistoryAddVariable, HistoryPut, HistoryAddAttr ! 文字列操作 ! Character handling ! use dc_string, only: StoA, CPrintf, LChar ! 組み込み関数 PRESENT の拡張版関数 ! Extended functions of intrinsic function "PRESENT" ! use dc_present, only: present_and_true ! 宣言文 ; Declaration statements ! logical, intent(in), optional:: flag_init_data ! 初期値データを作成する場合は, ! この引数に .true. を与えます. ! ! If initial data is created, ! give ".true." to this argument. ! 作業変数 ! Work variables ! integer :: NDims character(STRING), allocatable:: a_DimNames (:) integer , allocatable:: a_DimLens (:) character(STRING), allocatable:: a_DimLongNames(:) character(STRING), allocatable:: a_DimUnits (:) character(STRING):: title_msg ! 表題に付加するメッセージ. ! Message added to title real(DP):: origin_time ! 計算開始時刻. ! Start time of calculation character(STRING):: time_unit ! 日時の単位. Units of date and time logical:: flag_mpi_init integer:: n ! 実行文 ; Executable statement ! ! 初期化 ! Initialization ! flag_init_data_save = present_and_true( flag_init_data ) if ( .not. restart_surftemp_io_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if if ( restart_surftemp_opened ) return ! 表題に付加するメッセージの設定 ! Configure message added to title ! if ( .not. flag_init_data_save ) then title_msg = ' restart data of surface temperature' else title_msg = ' initial data of surface temperature' end if ! 時刻情報の取得 ! Get time information ! if ( .not. flag_init_data_save ) then origin_time = RestartTime + IntTime else origin_time = RestartTime end if time_unit = IntUnit flag_mpi_init = .true. ! リスタートファイルのオープン ! Open a restart file ! n = 3 ! lon, lat, time if ( kslmax /= 0 ) n = n + 1 ! dep if ( ksimax /= 0 ) n = n + 1 ! sidep NDims = n ! allocate( a_DimNames ( 1:NDims ) ) allocate( a_DimLens ( 1:NDims ) ) allocate( a_DimLongNames( 1:NDims ) ) allocate( a_DimUnits ( 1:NDims ) ) ! n = 1 ! lon a_DimNames (n) = 'lon' a_DimLens (n) = imax a_DimLongNames(n) = 'longitude' a_DimUnits (n) = 'degrees_east' n = n + 1 ! lat a_DimNames (n) = 'lat' a_DimLens (n) = jmax a_DimLongNames(n) = 'latitude' a_DimUnits (n) = 'degrees_north' if ( kslmax /= 0 ) then n = n + 1 ! dep a_DimNames (n) = 'dep' a_DimLens (n) = kslmax a_DimLongNames(n) = 'depth' a_DimUnits (n) = 'm' end if if ( ksimax /= 0 ) then n = n + 1 ! sidep a_DimNames (n) = 'sidep' a_DimLens (n) = ksimax a_DimLongNames(n) = 'depth' a_DimUnits (n) = 'm' end if n = n + 1 ! time a_DimNames (n) = time_name a_DimLens (n) = 0 a_DimLongNames(n) = time_name a_DimUnits (n) = time_unit ! !!$ call HistoryCreate( & !!$ & file = OutputFile, & ! (in) !!$ & title = trim(FileTitle) // trim(title_msg), & ! (in) !!$ & source = FileSource, institution = FileInstitution, & ! (in) !!$ & dims = StoA( 'lon ', 'lat ', 'dep ', 'sidep', time_name ), & ! (in) !!$ & dimsizes = (/ imax, jmax, kslmax, ksimax, 0 /), & ! (in) !!$ & longnames = StoA( 'longitude', 'latitude', 'depth ', 'depth ', time_name ), & ! (in) !!$ & units = StoA( 'degrees_east', 'degrees_north', 'm ', 'm ', & ! (in) !!$ & time_unit ), & ! (in) !!$ & origind = origin_time, & ! (in) optional !!$ & intervald = IntValue, & ! (in) optional !!$ & flag_mpi_split = flag_mpi_init, & ! (in) optional !!$ & history = gthst_rst & ! (out) optional !!$ & ) call HistoryCreate( file = OutputFile, title = trim(FileTitle) // trim(title_msg), source = FileSource, institution = FileInstitution, dims = a_DimNames, dimsizes = a_DimLens, longnames = a_DimLongNames, units = a_DimUnits, origind = origin_time, intervald = IntValue, flag_mpi_split = flag_mpi_init, history = gthst_rst ) ! deallocate( a_DimNames ) deallocate( a_DimLens ) deallocate( a_DimLongNames ) deallocate( a_DimUnits ) ! 座標データの設定 ! Axes data settings ! call HistoryAddAttr( varname = 'lon', attrname = 'standard_name', value = 'longitude', history = gthst_rst ) call HistoryAddAttr( varname = 'lat', attrname = 'standard_name', value = 'latitude', history = gthst_rst ) if ( kslmax /= 0 ) then call HistoryAddAttr( varname = 'dep', attrname = 'standard_name', value = 'depth', history = gthst_rst ) end if if ( ksimax /= 0 ) then call HistoryAddAttr( varname = 'sidep', attrname = 'standard_name', value = 'depth', history = gthst_rst ) end if call HistoryAddAttr( varname = time_name, attrname = 'standard_name', value = 'time', history = gthst_rst ) call HistoryPut( varname = 'lon', array = x_Lon / PI * 180.0_DP, history = gthst_rst ) call HistoryPut( varname = 'lat', array = y_Lat / PI * 180.0_DP, history = gthst_rst ) if ( kslmax /= 0 ) then call HistoryPut( varname = 'dep', array = z_SSDepth, history = gthst_rst ) end if if ( ksimax /= 0 ) then call HistoryPut( varname = 'sidep', array = z_SIDepth, history = gthst_rst ) end if ! 座標重みの設定 ! Axes weights settings ! call HistoryAddVariable( varname = 'lon_weight', dims = StoA('lon'), longname = 'weight for integration in longitude', units = 'radian', xtype = 'double', history = gthst_rst ) ! (inout) call HistoryAddAttr( varname = 'lon', attrname = 'gt_calc_weight', value = 'lon_weight', history = gthst_rst ) ! (inout) call HistoryPut( varname = 'lon_weight', array = x_Lon_Weight, history = gthst_rst ) ! (inout) call HistoryAddVariable( varname = 'lat_weight', dims = StoA('lat'), longname = 'weight for integration in latitude', units = 'radian', xtype = 'double', history = gthst_rst ) ! (inout) call HistoryAddAttr( varname = 'lat', attrname = 'gt_calc_weight', value = 'lat_weight', history = gthst_rst ) ! (inout) call HistoryPut( varname = 'lat_weight', array = y_Lat_Weight, history = gthst_rst ) ! (inout) ! 予報変数の設定 ! Predictional variables settings ! call HistoryAddVariable( varname = 'SurfTemp', dims = StoA('lon', 'lat', time_name), longname = 'surface temperature', units = 'K', xtype = 'double', history = gthst_rst ) if ( kslmax /= 0 ) then call HistoryAddVariable( varname = 'SoilTemp', dims = StoA('lon', 'lat', 'dep', time_name), longname = 'soil temperature', units = 'K', xtype = 'double', history = gthst_rst ) end if if ( ksimax /= 0 ) then call HistoryAddVariable( varname = 'SOSeaIceTemp', dims = StoA('lon', 'lat', 'sidep', time_name), longname = 'sea ice temperature', units = 'K', xtype = 'double', history = gthst_rst ) end if call HistoryAddVariable( varname = 'SOSeaIceMass', dims = StoA('lon', 'lat', time_name), longname = 'slab sea ice mass', units = 'kg m-2', xtype = 'double', history = gthst_rst ) call HistoryAddVariable( varname = 'SurfMajCompIceB', dims = StoA('lon', 'lat', time_name), longname = 'major component ice amount on the surface', units = 'kg m-2', xtype = 'double', history = gthst_rst ) call HistoryAddVariable( varname = 'SoilMoistB', dims = StoA('lon', 'lat', time_name), longname = 'soil moisture', units = 'kg m-2', xtype = 'double', history = gthst_rst ) call HistoryAddVariable( varname = 'SurfSnowB', dims = StoA('lon', 'lat', time_name), longname = 'snow amount on the surface', units = 'kg m-2', xtype = 'double', history = gthst_rst ) call HistoryAddVariable( varname = 'SurfMajCompIceN', dims = StoA('lon', 'lat', time_name), longname = 'major component ice amount on the surface', units = 'kg m-2', xtype = 'double', history = gthst_rst ) call HistoryAddVariable( varname = 'SoilMoistN', dims = StoA('lon', 'lat', time_name), longname = 'soil moisture', units = 'kg m-2', xtype = 'double', history = gthst_rst ) call HistoryAddVariable( varname = 'SurfSnowN', dims = StoA('lon', 'lat', time_name), longname = 'snow amount on the surface', units = 'kg m-2', xtype = 'double', history = gthst_rst ) restart_surftemp_opened = .true. end subroutine RestartSurfTempOpen
Subroutine : | |||
xy_SurfTemp(0:imax-1, 1:jmax) : | real(DP), intent(in)
| ||
xyz_SoilTemp(0:imax-1, 1:jmax, 1:kslmax) : | real(DP), intent(in), optional | ||
xyz_SOSeaIceTemp(0:imax-1, 1:jmax, 1:ksimax) : | real(DP), intent(in), optional | ||
xy_SOSeaIceMass(0:imax-1, 1:jmax) : | real(DP), intent(in), optional | ||
xy_SurfMajCompIceB(0:imax-1, 1:jmax) : | real(DP), intent(in), optional | ||
xy_SoilMoistB(0:imax-1, 1:jmax) : | real(DP), intent(in), optional | ||
xy_SurfSnowB(0:imax-1, 1:jmax) : | real(DP), intent(in), optional | ||
xy_SurfMajCompIceN(0:imax-1, 1:jmax) : | real(DP), intent(in), optional | ||
xy_SoilMoistN(0:imax-1, 1:jmax) : | real(DP), intent(in), optional | ||
xy_SurfSnowN(0:imax-1, 1:jmax) : | real(DP), intent(in), optional |
リスタートデータの出力を行います.
Output restart data
subroutine RestartSurfTempOutput( xy_SurfTemp, xyz_SoilTemp, xyz_SOSeaIceTemp, xy_SOSeaIceMass, xy_SurfMajCompIceB, xy_SoilMoistB, xy_SurfSnowB, xy_SurfMajCompIceN, xy_SoilMoistN, xy_SurfSnowN ) ! ! リスタートデータの出力を行います. ! ! Output restart data ! モジュール引用 ; USE statements ! ! gtool4 データ出力 ! Gtool4 data output ! use gtool_history, only: HistoryPut, HistorySetTime ! 時刻管理 ! Time control ! use timeset, only: TimeN, EndTime ! 計算終了時刻. ! End time of calculation ! 組み込み関数 PRESENT の拡張版関数 ! Extended functions of intrinsic function "PRESENT" ! use dc_present, only: present_and_true ! 宣言文 ; Declaration statements ! real(DP), intent(in) :: xy_SurfTemp (0:imax-1, 1:jmax) ! 地表面温度. ! Surface temperature real(DP), intent(in), optional :: xyz_SoilTemp (0:imax-1, 1:jmax, 1:kslmax) real(DP), intent(in), optional :: xyz_SOSeaIceTemp (0:imax-1, 1:jmax, 1:ksimax) real(DP), intent(in), optional :: xy_SOSeaIceMass (0:imax-1, 1:jmax) real(DP), intent(in), optional :: xy_SurfMajCompIceB(0:imax-1, 1:jmax) real(DP), intent(in), optional :: xy_SoilMoistB (0:imax-1, 1:jmax) real(DP), intent(in), optional :: xy_SurfSnowB (0:imax-1, 1:jmax) real(DP), intent(in), optional :: xy_SurfMajCompIceN(0:imax-1, 1:jmax) real(DP), intent(in), optional :: xy_SoilMoistN (0:imax-1, 1:jmax) real(DP), intent(in), optional :: xy_SurfSnowN (0:imax-1, 1:jmax) ! 作業変数 ! Work variables ! logical:: flag_output ! 出力のフラグ. ! Flag for output ! 実行文 ; Executable statement ! if ( .not. restart_surftemp_opened ) call RestartSurfTempOpen ! 出力タイミングのチェック ! Check output timing ! flag_output = TimeN - PrevOutputTime >= IntTime if ( TimeN >= EndTime .and. .not. flag_output_end ) then flag_output = .true. flag_output_end = .true. end if flag_output = ( .not. TimeN == PrevOutputTime ) .and. flag_output flag_output = flag_init_data_save .or. flag_output if ( .not. flag_output ) return ! 次回用に, 今回の出力 (希望) 時刻 を保存 ! Save output time (expected) in this time, for next time ! PrevOutputTime = PrevOutputTime + IntTime ! 時刻の設定 ! Set time ! call HistorySetTime( timed = TimeN, history = gthst_rst ) ! データ出力 ! Data output ! call HistoryPut( 'SurfTemp' , xy_SurfTemp , history = gthst_rst ) ! (in) if ( present( xyz_SoilTemp ) ) then if ( kslmax /= 0 ) then call HistoryPut( 'SoilTemp' , xyz_SoilTemp , history = gthst_rst ) ! (in) end if end if if ( present( xyz_SOSeaIceTemp ) ) then if ( ksimax /= 0 ) then call HistoryPut( 'SOSeaIceTemp' , xyz_SOSeaIceTemp , history = gthst_rst ) ! (in) end if end if if ( present( xy_SOSeaIceMass ) ) then call HistoryPut( 'SOSeaIceMass' , xy_SOSeaIceMass , history = gthst_rst ) ! (in) end if if ( present( xy_SurfMajCompIceB ) ) then call HistoryPut( 'SurfMajCompIceB', xy_SurfMajCompIceB, history = gthst_rst ) ! (in) end if if ( present( xy_SoilMoistB ) ) then call HistoryPut( 'SoilMoistB', xy_SoilMoistB, history = gthst_rst ) ! (in) end if if ( present( xy_SurfSnowB ) ) then call HistoryPut( 'SurfSnowB' , xy_SurfSnowB , history = gthst_rst ) ! (in) end if if ( present( xy_SurfMajCompIceN ) ) then call HistoryPut( 'SurfMajCompIceN', xy_SurfMajCompIceN, history = gthst_rst ) ! (in) end if if ( present( xy_SoilMoistN ) ) then call HistoryPut( 'SoilMoistN', xy_SoilMoistN, history = gthst_rst ) ! (in) end if if ( present( xy_SurfSnowN ) ) then call HistoryPut( 'SurfSnowN' , xy_SurfSnowN , history = gthst_rst ) ! (in) end if end subroutine RestartSurfTempOutput
Variable : | |||
restart_surftemp_io_inited = .false. : | logical, save, public
|
Variable : | |||
restart_surftemp_opened = .false. : | logical, save, public
|
Variable : | |||
InputFile : | character(STRING), save
|
Variable : | |||
IntUnit : | character(TOKEN)
|
Variable : | |||
OutputFile : | character(STRING), save
|
Constant : | |||
deltime_name = ‘deltime‘ : | character(*), parameter
|
Variable : | |||
flag_init_data_save : | logical, save
|
Variable : | |||
flag_output_end : | logical, save
|
Variable : | |||
gthst_rst : | type(GT_HISTORY), save
|
Constant : | |||
module_name = ‘restart_surftemp_io‘ : | character(*), parameter
|
Constant : | |||
time_name = ‘time‘ : | character(*), parameter
|
Constant : | |||
version = ’$Name: $’ // ’$Id: restart_surftemp_io.f90,v 1.11 2014/05/07 09:39:18 murashin Exp $’ : | character(*), parameter
|