SCALE-RM
mod_admin_time.f90
Go to the documentation of this file.
1 !-------------------------------------------------------------------------------
10 module mod_admin_time
11  !-----------------------------------------------------------------------------
12  !
13  !++ used modules
14  !
15  use scale_precision
16  use scale_stdio
17  use scale_prof
18  !-----------------------------------------------------------------------------
19  implicit none
20  private
21  !-----------------------------------------------------------------------------
22  !
23  !++ Public procedure
24  !
25  public :: admin_time_setup
26  public :: admin_time_checkstate
27  public :: admin_time_advance
28 
29  !-----------------------------------------------------------------------------
30  !
31  !++ Public parameters & variables
32  !
33  real(DP), public :: time_dtsec_atmos_restart
34  real(DP), public :: time_dtsec_ocean_restart
35  real(DP), public :: time_dtsec_land_restart
36  real(DP), public :: time_dtsec_urban_restart
37  real(DP), public :: time_dtsec_resume
38 
39  integer, public :: time_dstep_atmos_restart
40  integer, public :: time_dstep_ocean_restart
41  integer, public :: time_dstep_land_restart
42  integer, public :: time_dstep_urban_restart
43  integer, public :: time_dstep_resume
44 
45  logical, public :: time_doatmos_step
46  logical, public :: time_doatmos_dyn
47  logical, public :: time_doatmos_phy_cp
48  logical, public :: time_doatmos_phy_mp
49  logical, public :: time_doatmos_phy_rd
50  logical, public :: time_doatmos_phy_sf
51  logical, public :: time_doatmos_phy_tb
52  logical, public :: time_doatmos_phy_ch
53  logical, public :: time_doatmos_phy_ae
54  logical, public :: time_doatmos_restart
55  logical, public :: time_doocean_step
56  logical, public :: time_doocean_restart
57  logical, public :: time_doland_step
58  logical, public :: time_doland_restart
59  logical, public :: time_dourban_step
60  logical, public :: time_dourban_restart
61  logical, public :: time_doresume
62  logical, public :: time_doend
63 
64  !-----------------------------------------------------------------------------
65  !
66  !++ Private procedure
67  !
68  !-----------------------------------------------------------------------------
69  !
70  !++ Private parameters & variables
71  !
72  integer, private :: time_startdate(6) = (/ -999, 1, 1, 0, 0, 0 /)
73  real(DP), private :: time_startms = 0.0_dp
74  integer, private :: time_startday
75  real(DP), private :: time_startsec
76 
77  integer, private :: time_enddate(6)
78  real(DP), private :: time_endms
79  integer, private :: time_endday
80  real(DP), private :: time_endsec
81 
82  integer, private :: time_res_atmos_dyn = 0
83  integer, private :: time_res_atmos_phy_cp = 0
84  integer, private :: time_res_atmos_phy_mp = 0
85  integer, private :: time_res_atmos_phy_rd = 0
86  integer, private :: time_res_atmos_phy_sf = 0
87  integer, private :: time_res_atmos_phy_tb = 0
88  integer, private :: time_res_atmos_phy_ch = 0
89  integer, private :: time_res_atmos_phy_ae = 0
90  integer, private :: time_res_atmos_restart = 0
91  integer, private :: time_res_ocean = 0
92  integer, private :: time_res_ocean_restart = 0
93  integer, private :: time_res_land = 0
94  integer, private :: time_res_land_restart = 0
95  integer, private :: time_res_urban = 0
96  integer, private :: time_res_urban_restart = 0
97  integer, private :: time_res_resume
98 
99  real(DP), private :: time_wallclock_start ! Start time of wall clock [sec]
100  real(DP), private :: time_wallclock_limit = -1.0_dp ! Elapse time limit of wall clock time [sec]
101  real(DP), private :: time_wallclock_safe = 0.9_dp ! Safety coefficient for elapse time limit
102 
103  real(DP), private, parameter :: eps = 1.e-6_dp
104 
105  !-----------------------------------------------------------------------------
106 contains
107  !-----------------------------------------------------------------------------
109  subroutine admin_time_setup( &
110  setup_TimeIntegration )
111  use gtool_file, only: &
112  filegetdatainfo
113  use scale_process, only: &
114  prc_myrank, &
115  prc_mpistop, &
116  prc_mpitime
117  use scale_const, only: &
118  undef8 => const_undef8
119  use scale_calendar, only: &
120  calendar_date2daysec, &
121  calendar_daysec2date, &
122  calendar_adjust_daysec, &
123  calendar_combine_daysec, &
124  calendar_unit2sec, &
125  calendar_cfunits2sec, &
126  calendar_date2char
127  use scale_time, only: &
128  time_dtsec, &
129  time_nowdate, &
130  time_nowms, &
131  time_nowday, &
132  time_nowsec, &
133  time_nowdaysec, &
134  time_nowstep, &
135  time_nstep, &
136  time_dtsec_atmos_dyn, &
137  nstep_dyn => time_nstep_atmos_dyn, &
138  time_dtsec_atmos_phy_cp, &
139  time_dtsec_atmos_phy_mp, &
140  time_dtsec_atmos_phy_rd, &
141  time_dtsec_atmos_phy_sf, &
142  time_dtsec_atmos_phy_tb, &
143  time_dtsec_atmos_phy_ch, &
144  time_dtsec_atmos_phy_ae, &
145  time_dtsec_ocean, &
146  time_dtsec_land, &
147  time_dtsec_urban, &
148  time_dtsec_wallclock_check, &
149  time_dstep_atmos_dyn, &
150  time_dstep_atmos_phy_cp, &
151  time_dstep_atmos_phy_mp, &
152  time_dstep_atmos_phy_rd, &
153  time_dstep_atmos_phy_sf, &
154  time_dstep_atmos_phy_tb, &
155  time_dstep_atmos_phy_ch, &
156  time_dstep_atmos_phy_ae, &
157  time_dstep_ocean, &
158  time_dstep_land, &
159  time_dstep_urban, &
160  time_dstep_wallclock_check, &
161  time_offset_year, &
162  time_startdaysec
163  use mod_atmos_vars, only: &
164  restart_in_basename => atmos_restart_in_basename
165  implicit none
166 
167  logical, intent(in) :: setup_TimeIntegration
168 
169  real(DP) :: TIME_DURATION = undef8
170  character(len=H_SHORT) :: TIME_DURATION_UNIT = "SEC"
171  real(DP) :: TIME_DT = undef8
172  character(len=H_SHORT) :: TIME_DT_UNIT = "SEC"
173 
174  real(DP) :: TIME_DT_ATMOS_DYN = undef8
175  character(len=H_SHORT) :: TIME_DT_ATMOS_DYN_UNIT = "SEC"
176  integer :: TIME_NSTEP_ATMOS_DYN = -1
177  real(DP) :: TIME_DT_ATMOS_PHY_CP = undef8
178  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_CP_UNIT = ""
179  real(DP) :: TIME_DT_ATMOS_PHY_MP = undef8
180  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_MP_UNIT = ""
181  real(DP) :: TIME_DT_ATMOS_PHY_RD = undef8
182  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_RD_UNIT = ""
183  real(DP) :: TIME_DT_ATMOS_PHY_SF = undef8
184  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_SF_UNIT = ""
185  real(DP) :: TIME_DT_ATMOS_PHY_TB = undef8
186  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_TB_UNIT = ""
187  real(DP) :: TIME_DT_ATMOS_PHY_CH = undef8
188  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_CH_UNIT = ""
189  real(DP) :: TIME_DT_ATMOS_PHY_AE = undef8
190  character(len=H_SHORT) :: TIME_DT_ATMOS_PHY_AE_UNIT = ""
191  real(DP) :: TIME_DT_ATMOS_RESTART = undef8
192  character(len=H_SHORT) :: TIME_DT_ATMOS_RESTART_UNIT = ""
193 
194  real(DP) :: TIME_DT_OCEAN = undef8
195  character(len=H_SHORT) :: TIME_DT_OCEAN_UNIT = ""
196  real(DP) :: TIME_DT_OCEAN_RESTART = undef8
197  character(len=H_SHORT) :: TIME_DT_OCEAN_RESTART_UNIT = ""
198 
199  real(DP) :: TIME_DT_LAND = undef8
200  character(len=H_SHORT) :: TIME_DT_LAND_UNIT = ""
201  real(DP) :: TIME_DT_LAND_RESTART = undef8
202  character(len=H_SHORT) :: TIME_DT_LAND_RESTART_UNIT = ""
203 
204  real(DP) :: TIME_DT_URBAN = undef8
205  character(len=H_SHORT) :: TIME_DT_URBAN_UNIT = ""
206  real(DP) :: TIME_DT_URBAN_RESTART = undef8
207  character(len=H_SHORT) :: TIME_DT_URBAN_RESTART_UNIT = ""
208 
209  real(DP) :: TIME_DT_RESUME = undef8
210  character(len=H_SHORT) :: TIME_DT_RESUME_UNIT = ""
211 
212  real(DP) :: TIME_DT_WALLCLOCK_CHECK = undef8
213  character(len=H_SHORT) :: TIME_DT_WALLCLOCK_CHECK_UNIT = ""
214 
215  namelist / param_time / &
216  time_startdate, &
217  time_startms, &
218  time_duration, &
219  time_duration_unit, &
220  time_dt, &
221  time_dt_unit, &
222  time_dt_atmos_dyn, &
223  time_dt_atmos_dyn_unit, &
224  time_nstep_atmos_dyn, &
225  time_dt_atmos_phy_cp, &
226  time_dt_atmos_phy_cp_unit, &
227  time_dt_atmos_phy_mp, &
228  time_dt_atmos_phy_mp_unit, &
229  time_dt_atmos_phy_rd, &
230  time_dt_atmos_phy_rd_unit, &
231  time_dt_atmos_phy_sf, &
232  time_dt_atmos_phy_sf_unit, &
233  time_dt_atmos_phy_tb, &
234  time_dt_atmos_phy_tb_unit, &
235  time_dt_atmos_phy_ch, &
236  time_dt_atmos_phy_ch_unit, &
237  time_dt_atmos_phy_ae, &
238  time_dt_atmos_phy_ae_unit, &
239  time_dt_atmos_restart, &
240  time_dt_atmos_restart_unit, &
241  time_dt_ocean, &
242  time_dt_ocean_unit, &
243  time_dt_ocean_restart, &
244  time_dt_ocean_restart_unit, &
245  time_dt_land, &
246  time_dt_land_unit, &
247  time_dt_land_restart, &
248  time_dt_land_restart_unit, &
249  time_dt_urban, &
250  time_dt_urban_unit, &
251  time_dt_urban_restart, &
252  time_dt_urban_restart_unit, &
253  time_dt_wallclock_check, &
254  time_dt_wallclock_check_unit, &
255  time_dt_resume, &
256  time_dt_resume_unit, &
257  time_wallclock_limit, &
258  time_wallclock_safe
259 
260  integer :: dateday
261  real(DP) :: datesec
262  real(DP) :: cftime
263  character(len=H_MID) :: cfunits
264 
265  real(DP) :: TIME_DURATIONSEC
266  character(len=27) :: startchardate
267  character(len=27) :: endchardate
268 
269  integer :: ierr
270  !---------------------------------------------------------------------------
271 
272  if( io_l ) write(io_fid_log,*)
273  if( io_l ) write(io_fid_log,*) '++++++ Module[TIME] / Categ[COMMON] / Origin[SCALElib]'
274 
275  time_nstep_atmos_dyn = -1
276 
277  !--- read namelist
278  rewind(io_fid_conf)
279  read(io_fid_conf,nml=param_time,iostat=ierr)
280  if( ierr < 0 ) then !--- missing
281  if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
282  elseif( ierr > 0 ) then !--- fatal error
283  write(*,*) 'xxx Not appropriate names in namelist PARAM_TIME. Check!'
284  call prc_mpistop
285  endif
286  if( io_lnml ) write(io_fid_log,nml=param_time)
287 
288  ! check time setting
289  if ( setup_timeintegration ) then
290  if ( time_dt == undef8 ) then
291  write(*,*) 'xxx Not found TIME_DT.'
292  call prc_mpistop
293  endif
294  if ( time_duration == undef8 ) then
295  write(*,*) 'xxx Not found TIME_DURATION.'
296  call prc_mpistop
297  endif
298 
299  ! DYN
300  if ( time_dt_atmos_dyn == undef8 ) then
301  if ( time_nstep_atmos_dyn < 0 ) then
302  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_DYN. TIME_DT is used.'
303  time_dt_atmos_dyn = time_dt
304  endif
305  endif
306  if ( time_dt_atmos_dyn_unit == '' ) then
307  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_DYN_UNIT. TIME_DT_UNIT is used.'
308  time_dt_atmos_dyn_unit = time_dt_unit
309  endif
310  ! PHY_CP
311  if ( time_dt_atmos_phy_cp == undef8 ) then
312  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_CP. TIME_DT is used.'
313  time_dt_atmos_phy_cp = time_dt
314  endif
315  if ( time_dt_atmos_phy_cp_unit == '' ) then
316  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_CP_UNIT. TIME_DT_UNIT is used.'
317  time_dt_atmos_phy_cp_unit = time_dt_unit
318  endif
319  ! PHY_MP
320  if ( time_dt_atmos_phy_mp == undef8 ) then
321  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_MP. TIME_DT is used.'
322  time_dt_atmos_phy_mp = time_dt
323  endif
324  if ( time_dt_atmos_phy_mp_unit == '' ) then
325  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_MP_UNIT. TIME_DT_UNIT is used.'
326  time_dt_atmos_phy_mp_unit = time_dt_unit
327  endif
328  ! PHY_RD
329  if ( time_dt_atmos_phy_rd == undef8 ) then
330  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_RD. TIME_DT is used.'
331  time_dt_atmos_phy_rd = time_dt
332  endif
333  if ( time_dt_atmos_phy_rd_unit == '' ) then
334  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_RD_UNIT. TIME_DT_UNIT is used.'
335  time_dt_atmos_phy_rd_unit = time_dt_unit
336  endif
337  ! PHY_SF
338  if ( time_dt_atmos_phy_sf == undef8 ) then
339  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_SF. TIME_DT is used.'
340  time_dt_atmos_phy_sf = time_dt
341  endif
342  if ( time_dt_atmos_phy_sf_unit == '' ) then
343  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_SF_UNIT. TIME_DT_UNIT is used.'
344  time_dt_atmos_phy_sf_unit = time_dt_unit
345  endif
346  ! PHY_TB
347  if ( time_dt_atmos_phy_tb == undef8 ) then
348  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_TB. TIME_DT is used.'
349  time_dt_atmos_phy_tb = time_dt
350  endif
351  if ( time_dt_atmos_phy_tb_unit == '' ) then
352  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_TB_UNIT. TIME_DT_UNIT is used.'
353  time_dt_atmos_phy_tb_unit = time_dt_unit
354  endif
355  ! PHY_CH
356  if ( time_dt_atmos_phy_ch == undef8 ) then
357  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_CH. TIME_DT is used.'
358  time_dt_atmos_phy_ch = time_dt
359  endif
360  if ( time_dt_atmos_phy_ch_unit == '' ) then
361  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_CH_UNIT. TIME_DT_UNIT is used.'
362  time_dt_atmos_phy_ch_unit = time_dt_unit
363  endif
364  ! PHY_AE
365  if ( time_dt_atmos_phy_ae == undef8 ) then
366  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_AE. TIME_DT is used.'
367  time_dt_atmos_phy_ae = time_dt
368  endif
369  if ( time_dt_atmos_phy_ae_unit == '' ) then
370  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_PHY_AE_UNIT. TIME_DT_UNIT is used.'
371  time_dt_atmos_phy_ae_unit = time_dt_unit
372  endif
373  ! ATMOS RESTART
374  if ( time_dt_atmos_restart == undef8 ) then
375  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_RESTART. TIME_DURATION is used.'
376  time_dt_atmos_restart = time_duration
377  endif
378  if ( time_dt_atmos_restart_unit == '' ) then
379  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_ATMOS_RESTART_UNIT. TIME_DURATION_UNIT is used.'
380  time_dt_atmos_restart_unit = time_duration_unit
381  endif
382  ! OCEAN
383  if ( time_dt_ocean == undef8 ) then
384  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_OCEAN. TIME_DT is used.'
385  time_dt_ocean = time_dt
386  endif
387  if ( time_dt_ocean_unit == '' ) then
388  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_OCEAN_UNIT. TIME_DT_UNIT is used.'
389  time_dt_ocean_unit = time_dt_unit
390  endif
391  ! OCEAN RESTART
392  if ( time_dt_ocean_restart == undef8 ) then
393  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_OCEAN_RESTART. TIME_DURATION is used.'
394  time_dt_ocean_restart = time_duration
395  endif
396  if ( time_dt_ocean_restart_unit == '' ) then
397  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_OCEAN_RESTART_UNIT. TIME_DURATION_UNIT is used.'
398  time_dt_ocean_restart_unit = time_duration_unit
399  endif
400  ! LAND
401  if ( time_dt_land == undef8 ) then
402  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_LAND. TIME_DT is used.'
403  time_dt_land = time_dt
404  endif
405  if ( time_dt_land_unit == '' ) then
406  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_LAND_UNIT. TIME_DT_UNIT is used.'
407  time_dt_land_unit = time_dt_unit
408  endif
409  ! LAND RESTART
410  if ( time_dt_land_restart == undef8 ) then
411  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_LAND_RESTART. TIME_DURATION is used.'
412  time_dt_land_restart = time_duration
413  endif
414  if ( time_dt_land_restart_unit == '' ) then
415  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_LAND_RESTART_UNIT. TIME_DURATION_UNIT is used.'
416  time_dt_land_restart_unit = time_duration_unit
417  endif
418  ! URBAN
419  if ( time_dt_urban == undef8 ) then
420  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_URBAN. TIME_DT is used.'
421  time_dt_urban = time_dt
422  endif
423  if ( time_dt_urban_unit == '' ) then
424  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_URBAN_UNIT. TIME_DT_UNIT is used.'
425  time_dt_urban_unit = time_dt_unit
426  endif
427  ! URBAN RESTART
428  if ( time_dt_urban_restart == undef8 ) then
429  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_URBAN_RESTART. TIME_DURATION is used.'
430  time_dt_urban_restart = time_duration
431  endif
432  if ( time_dt_urban_restart_unit == '' ) then
433  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_URBAN_RESTART_UNIT. TIME_DURATION_UNIT is used.'
434  time_dt_urban_restart_unit = time_duration_unit
435  endif
436  ! Resume
437  if ( time_dt_resume == undef8 ) then
438  time_dt_resume = time_duration
439  endif
440  if ( time_dt_resume_unit == '' ) then
441  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_RESUME_UNIT. TIME_DURATION_UNIT is used.'
442  time_dt_resume_unit = time_duration_unit
443  endif
444  endif
445 
446  !--- calculate time
447  if ( time_startdate(1) == -999 ) then
448  if ( restart_in_basename /= '' ) then ! read start time from the restart data
449  call filegetdatainfo( restart_in_basename, & ! [IN]
450  'DENS', & ! [IN]
451  prc_myrank, & ! [IN]
452  0, & ! [IN] step
453  time_start = cftime, & ! [OUT]
454  time_units = cfunits ) ! [OUT]
455 
456  dateday = 0
457  datesec = calendar_cfunits2sec( cftime, cfunits, 0 )
458 
459  call calendar_adjust_daysec( dateday, datesec )
460 
461  call calendar_daysec2date( time_startdate, & ! [OUT]
462  time_startms, & ! [OUT]
463  dateday, & ! [IN]
464  datesec, & ! [IN]
465  0 ) ! [IN]
466  else
467  time_startdate = (/ 0, 1, 1, 0, 0, 0 /)
468  time_startms = 0.0_dp
469  endif
470  else
471  time_startms = time_startms * 1.e-3_dp
472  endif
473 
474  time_offset_year = time_startdate(1)
475 
476  call calendar_date2daysec( time_startday, & ! [OUT]
477  time_startsec, & ! [OUT]
478  time_startdate(:), & ! [IN]
479  time_startms, & ! [IN]
480  time_offset_year ) ! [IN]
481 
482  call calendar_date2char( startchardate, & ! [OUT]
483  time_startdate(:), & ! [IN]
484  time_startms ) ! [IN]
485 
486  time_startdaysec = calendar_combine_daysec( time_startday, time_startsec )
487 
488  time_nowdate(:) = time_startdate(:)
489  time_nowms = time_startms
490  time_nowday = time_startday
491  time_nowsec = time_startsec
492  time_nowdaysec = calendar_combine_daysec( time_nowday, time_nowsec )
493 
494  time_endday = time_startday
495 
496  if ( setup_timeintegration ) then
497  call calendar_unit2sec( time_durationsec, time_duration, time_duration_unit )
498  time_endsec = time_startsec + time_durationsec
499  else
500  time_endsec = time_startsec
501  endif
502 
503  call calendar_adjust_daysec( time_endday, time_endsec ) ! [INOUT]
504 
505  call calendar_daysec2date( time_enddate(:), & ! [OUT]
506  time_endms, & ! [OUT]
507  time_endday, & ! [IN]
508  time_endsec, & ! [IN]
509  time_offset_year ) ! [IN]
510 
511  call calendar_date2char( endchardate, & ! [OUT]
512  time_enddate(:), & ! [IN]
513  time_endms ) ! [IN]
514 
515  if( io_l ) write(io_fid_log,*)
516  if( io_l ) write(io_fid_log,*) '*** Date/time setting ***'
517  if( io_l ) write(io_fid_log,'(1x,A,A)') '*** START Date : ', startchardate
518  if( io_l ) write(io_fid_log,'(1x,A,A)') '*** END Date : ', endchardate
519 
520  if ( setup_timeintegration ) then
521 
522  call calendar_unit2sec( time_dtsec, time_dt, time_dt_unit )
523 
524  time_nstep = int( time_durationsec / time_dtsec )
525  time_nowstep = 1
526 
527  if( io_l ) write(io_fid_log,'(1x,A,F10.3)') '*** delta t (sec.) :', time_dtsec
528  if( io_l ) write(io_fid_log,'(1x,A,I10)' ) '*** No. of steps :', time_nstep
529 
530  !--- calculate intervals for atmosphere
531  if ( time_dt_atmos_dyn /= undef8 ) then
532  call calendar_unit2sec( time_dtsec_atmos_dyn, time_dt_atmos_dyn, time_dt_atmos_dyn_unit )
533  else
534  time_dtsec_atmos_dyn = time_dtsec / time_nstep_atmos_dyn
535  endif
536  call calendar_unit2sec( time_dtsec_atmos_phy_cp, time_dt_atmos_phy_cp, time_dt_atmos_phy_cp_unit )
537  call calendar_unit2sec( time_dtsec_atmos_phy_mp, time_dt_atmos_phy_mp, time_dt_atmos_phy_mp_unit )
538  call calendar_unit2sec( time_dtsec_atmos_phy_rd, time_dt_atmos_phy_rd, time_dt_atmos_phy_rd_unit )
539  call calendar_unit2sec( time_dtsec_atmos_phy_sf, time_dt_atmos_phy_sf, time_dt_atmos_phy_sf_unit )
540  call calendar_unit2sec( time_dtsec_atmos_phy_tb, time_dt_atmos_phy_tb, time_dt_atmos_phy_tb_unit )
541  call calendar_unit2sec( time_dtsec_atmos_phy_ch, time_dt_atmos_phy_ch, time_dt_atmos_phy_ch_unit )
542  call calendar_unit2sec( time_dtsec_atmos_phy_ae, time_dt_atmos_phy_ae, time_dt_atmos_phy_ae_unit )
543  call calendar_unit2sec( time_dtsec_atmos_restart, time_dt_atmos_restart, time_dt_atmos_restart_unit )
544  call calendar_unit2sec( time_dtsec_ocean, time_dt_ocean, time_dt_ocean_unit )
545  call calendar_unit2sec( time_dtsec_ocean_restart, time_dt_ocean_restart, time_dt_ocean_restart_unit )
546  call calendar_unit2sec( time_dtsec_land, time_dt_land, time_dt_land_unit )
547  call calendar_unit2sec( time_dtsec_land_restart, time_dt_land_restart, time_dt_land_restart_unit )
548  call calendar_unit2sec( time_dtsec_urban, time_dt_urban, time_dt_urban_unit )
549  call calendar_unit2sec( time_dtsec_urban_restart, time_dt_urban_restart, time_dt_urban_restart_unit )
550  call calendar_unit2sec( time_dtsec_resume, time_dt_resume, time_dt_resume_unit )
551 
552  time_nstep_atmos_dyn = max( nint( time_dtsec / time_dtsec_atmos_dyn ), 1 )
553  nstep_dyn = time_nstep_atmos_dyn
554 
555  time_dtsec_atmos_dyn = max( time_dtsec_atmos_dyn, time_dtsec /time_nstep_atmos_dyn )
556  time_dtsec_atmos_phy_cp = max( time_dtsec_atmos_phy_cp, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
557  time_dtsec_atmos_phy_mp = max( time_dtsec_atmos_phy_mp, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
558  time_dtsec_atmos_phy_rd = max( time_dtsec_atmos_phy_rd, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
559  time_dtsec_atmos_phy_sf = max( time_dtsec_atmos_phy_sf, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
560  time_dtsec_atmos_phy_tb = max( time_dtsec_atmos_phy_tb, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
561  time_dtsec_atmos_phy_ch = max( time_dtsec_atmos_phy_ch, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
562  time_dtsec_atmos_phy_ae = max( time_dtsec_atmos_phy_ae, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
563  time_dtsec_atmos_restart = max( time_dtsec_atmos_restart, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
564  time_dtsec_ocean = max( time_dtsec_ocean, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
565  time_dtsec_ocean_restart = max( time_dtsec_ocean_restart, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
566  time_dtsec_land = max( time_dtsec_land, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
567  time_dtsec_land_restart = max( time_dtsec_land_restart, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
568  time_dtsec_urban = max( time_dtsec_urban, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
569  time_dtsec_urban_restart = max( time_dtsec_urban_restart, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
570  time_dtsec_resume = max( time_dtsec_resume, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
571 
572  time_dstep_atmos_dyn = nint( time_dtsec_atmos_dyn / time_dtsec * time_nstep_atmos_dyn )
573  time_dstep_atmos_phy_cp = nint( time_dtsec_atmos_phy_cp / time_dtsec )
574  time_dstep_atmos_phy_mp = nint( time_dtsec_atmos_phy_mp / time_dtsec )
575  time_dstep_atmos_phy_rd = nint( time_dtsec_atmos_phy_rd / time_dtsec )
576  time_dstep_atmos_phy_sf = nint( time_dtsec_atmos_phy_sf / time_dtsec )
577  time_dstep_atmos_phy_tb = nint( time_dtsec_atmos_phy_tb / time_dtsec )
578  time_dstep_atmos_phy_ch = nint( time_dtsec_atmos_phy_ch / time_dtsec )
579  time_dstep_atmos_phy_ae = nint( time_dtsec_atmos_phy_ae / time_dtsec )
580  time_dstep_ocean = nint( time_dtsec_ocean / time_dtsec )
581  time_dstep_land = nint( time_dtsec_land / time_dtsec )
582  time_dstep_urban = nint( time_dtsec_urban / time_dtsec )
583  time_dstep_atmos_restart = nint( time_dtsec_atmos_restart / time_dtsec )
584  time_dstep_ocean_restart = nint( time_dtsec_ocean_restart / time_dtsec )
585  time_dstep_land_restart = nint( time_dtsec_land_restart / time_dtsec )
586  time_dstep_urban_restart = nint( time_dtsec_urban_restart / time_dtsec )
587  time_dstep_resume = nint( time_dtsec_resume / time_dtsec )
588 
589  time_res_resume = time_dstep_resume - 1
590 
591  if ( abs( real(time_nstep_atmos_dyn,kind=dp)*time_dtsec_atmos_dyn &
592  - real(time_dstep_atmos_dyn,kind=dp)*time_dtsec ) > eps ) then
593  write(*,*) 'xxx delta t(ATMOS_DYN) must be a multiple of delta t ', &
594  time_dtsec_atmos_dyn, real(time_dstep_atmos_dyn,kind=dp)*time_dtsec
595  call prc_mpistop
596  endif
597  if ( abs(time_dtsec_atmos_phy_cp-real(time_dstep_atmos_phy_cp,kind=dp)*time_dtsec) > eps ) then
598  write(*,*) 'xxx delta t(ATMOS_PHY_CP) must be a multiple of delta t ', &
599  time_dtsec_atmos_phy_cp, real(time_dstep_atmos_phy_cp,kind=dp)*time_dtsec
600  call prc_mpistop
601  endif
602  if ( abs(time_dtsec_atmos_phy_mp-real(time_dstep_atmos_phy_mp,kind=dp)*time_dtsec) > eps ) then
603  write(*,*) 'xxx delta t(ATMOS_PHY_MP) must be a multiple of delta t ', &
604  time_dtsec_atmos_phy_mp, real(time_dstep_atmos_phy_mp,kind=dp)*time_dtsec
605  call prc_mpistop
606  endif
607  if ( abs(time_dtsec_atmos_phy_rd-real(time_dstep_atmos_phy_rd,kind=dp)*time_dtsec) > eps ) then
608  write(*,*) 'xxx delta t(ATMOS_PHY_RD) must be a multiple of delta t ', &
609  time_dtsec_atmos_phy_rd, real(time_dstep_atmos_phy_rd,kind=dp)*time_dtsec
610  call prc_mpistop
611  endif
612  if ( abs(time_dtsec_atmos_phy_sf-real(time_dstep_atmos_phy_sf,kind=dp)*time_dtsec) > eps ) then
613  write(*,*) 'xxx delta t(ATMOS_PHY_SF) must be a multiple of delta t ', &
614  time_dtsec_atmos_phy_sf, real(time_dstep_atmos_phy_sf,kind=dp)*time_dtsec
615  call prc_mpistop
616  endif
617  if ( abs(time_dtsec_atmos_phy_tb-real(time_dstep_atmos_phy_tb,kind=dp)*time_dtsec) > eps ) then
618  write(*,*) 'xxx delta t(ATMOS_PHY_TB) must be a multiple of delta t ', &
619  time_dtsec_atmos_phy_tb, real(time_dstep_atmos_phy_tb,kind=dp)*time_dtsec
620  call prc_mpistop
621  endif
622  if ( abs(time_dtsec_atmos_phy_ch-real(time_dstep_atmos_phy_ch,kind=dp)*time_dtsec) > eps ) then
623  write(*,*) 'xxx delta t(ATMOS_PHY_CH) must be a multiple of delta t ', &
624  time_dtsec_atmos_phy_ch, real(time_dstep_atmos_phy_ch,kind=dp)*time_dtsec
625  call prc_mpistop
626  endif
627  if ( abs(time_dtsec_atmos_phy_ae-real(time_dstep_atmos_phy_ae,kind=dp)*time_dtsec) > eps ) then
628  write(*,*) 'xxx delta t(ATMOS_PHY_AE) must be a multiple of delta t ', &
629  time_dtsec_atmos_phy_ae, real(time_dstep_atmos_phy_ae,kind=dp)*time_dtsec
630  call prc_mpistop
631  endif
632  if ( abs(time_dtsec_ocean-real(time_dstep_ocean,kind=dp)*time_dtsec) > eps ) then
633  write(*,*) 'xxx delta t(OCEAN) must be a multiple of delta t ', &
634  time_dtsec_ocean, real(time_dstep_ocean,kind=dp)*time_dtsec
635  call prc_mpistop
636  endif
637  if ( abs(time_dtsec_land-real(time_dstep_land,kind=dp)*time_dtsec) > eps ) then
638  write(*,*) 'xxx delta t(LAND) must be a multiple of delta t ', &
639  time_dtsec_land, real(time_dstep_land,kind=dp)*time_dtsec
640  call prc_mpistop
641  endif
642  if ( abs(time_dtsec_urban-real(time_dstep_urban,kind=dp)*time_dtsec) > eps ) then
643  write(*,*) 'xxx delta t(URBAN) must be a multiple of delta t ', &
644  time_dtsec_urban, real(time_dstep_urban,kind=dp)*time_dtsec
645  call prc_mpistop
646  endif
647  if ( abs(time_dtsec_atmos_restart-real(time_dstep_atmos_restart,kind=dp)*time_dtsec) > eps ) then
648  write(*,*) 'xxx delta t(ATMOS_RESTART) must be a multiple of delta t ', &
649  time_dtsec_atmos_restart, real(time_dstep_atmos_restart,kind=dp)*time_dtsec
650  call prc_mpistop
651  endif
652  if ( abs(time_dtsec_ocean_restart-real(time_dstep_ocean_restart,kind=dp)*time_dtsec) > eps ) then
653  write(*,*) 'xxx delta t(OCEAN_RESTART) must be a multiple of delta t ', &
654  time_dtsec_ocean_restart, real(time_dstep_ocean_restart,kind=dp)*time_dtsec
655  call prc_mpistop
656  endif
657  if ( abs(time_dtsec_land_restart-real(time_dstep_land_restart,kind=dp)*time_dtsec) > eps ) then
658  write(*,*) 'xxx delta t(LAND_RESTART) must be a multiple of delta t ', &
659  time_dtsec_land_restart, real(time_dstep_land_restart,kind=dp)*time_dtsec
660  call prc_mpistop
661  endif
662  if ( abs(time_dtsec_urban_restart-real(time_dstep_urban_restart,kind=dp)*time_dtsec) > eps ) then
663  write(*,*) 'xxx delta t(URBAN_RESTART) must be a multiple of delta t ', &
664  time_dtsec_urban_restart, real(time_dstep_urban_restart,kind=dp)*time_dtsec
665  call prc_mpistop
666  endif
667  if ( abs(time_dtsec_urban_restart-real(time_dstep_urban_restart,kind=dp)*time_dtsec) > eps ) then
668  write(*,*) 'xxx delta t(URBAN_RESTART) must be a multiple of delta t ', &
669  time_dtsec_urban_restart, real(time_dstep_urban_restart,kind=dp)*time_dtsec
670  call prc_mpistop
671  endif
672  if ( abs(time_dtsec_resume-real(time_dstep_resume,kind=dp)*time_dtsec) > eps ) then
673  write(*,*) 'xxx delta t(RESUME) must be a multiple of delta t ', &
674  time_dtsec_resume, real(time_dstep_resume,kind=dp)*time_dtsec
675  call prc_mpistop
676  endif
677 
678  if( io_l ) write(io_fid_log,*)
679  if( io_l ) write(io_fid_log,*) '*** Time interval for each processes (sec.)'
680  if( io_l ) write(io_fid_log,*) '*** Atmosphere'
681  if( io_l ) write(io_fid_log,'(1x,A,F10.3)') '*** Dynamics (time) : ', time_dtsec_atmos_dyn
682  if( io_l ) write(io_fid_log,'(1x,A,I10,A,I8,A)') '*** (step) : ', time_nstep_atmos_dyn, &
683  ' (step interval=', time_dstep_atmos_dyn, ')'
684  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Cumulus : ', time_dtsec_atmos_phy_cp, &
685  ' (step interval=', time_dstep_atmos_phy_cp, ')'
686  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Cloud Microphysics : ', time_dtsec_atmos_phy_mp, &
687  ' (step interval=', time_dstep_atmos_phy_mp, ')'
688  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Radiation : ', time_dtsec_atmos_phy_rd, &
689  ' (step interval=', time_dstep_atmos_phy_rd, ')'
690  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Surface Flux : ', time_dtsec_atmos_phy_sf, &
691  ' (step interval=', time_dstep_atmos_phy_sf, ')'
692  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Turbulence : ', time_dtsec_atmos_phy_tb, &
693  ' (step interval=', time_dstep_atmos_phy_tb, ')'
694  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Chemistry : ', time_dtsec_atmos_phy_ch, &
695  ' (step interval=', time_dstep_atmos_phy_ch, ')'
696  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Physics, Aerosol : ', time_dtsec_atmos_phy_ae, &
697  ' (step interval=', time_dstep_atmos_phy_ae, ')'
698  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Ocean : ', time_dtsec_ocean, &
699  ' (step interval=', time_dstep_ocean, ')'
700  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Land : ', time_dtsec_land, &
701  ' (step interval=', time_dstep_land, ')'
702  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Urban : ', time_dtsec_urban, &
703  ' (step interval=', time_dstep_urban, ')'
704  if( io_l ) write(io_fid_log,*)
705  if( io_l ) write(io_fid_log,*) '*** Time interval for restart (sec.)'
706  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Atmospheric Variables : ', time_dtsec_atmos_restart, &
707  ' (step interval=', time_dstep_atmos_restart, ')'
708  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Ocean Variables : ', time_dtsec_ocean_restart, &
709  ' (step interval=', time_dstep_ocean_restart, ')'
710  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Land Variables : ', time_dtsec_land_restart, &
711  ' (step interval=', time_dstep_land_restart, ')'
712  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Urban Variables : ', time_dtsec_urban_restart, &
713  ' (step interval=', time_dstep_urban_restart, ')'
714  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Resume : ', time_dtsec_resume, &
715  ' (step interval=', time_dstep_resume, ')'
716  else
717  time_dtsec = 1.0_rp
718  endif
719 
720  ! WALLCLOCK TERMINATOR SETUP
721  if ( time_wallclock_limit > 0.0_dp ) then
722  if( io_l ) write(io_fid_log,*)
723  if( io_l ) write(io_fid_log,*) '*** Wall clock time limit of execution is specified.'
724 
725  if ( time_dt_wallclock_check == undef8 ) then
726  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_WALLCLOCK_CHECK. largest time step interval is used.'
727  time_dtsec_wallclock_check = max( time_dtsec_atmos_dyn*time_nstep_atmos_dyn, &
728  time_dtsec_atmos_phy_cp, &
729  time_dtsec_atmos_phy_mp, &
730  time_dtsec_atmos_phy_rd, &
731  time_dtsec_atmos_phy_sf, &
732  time_dtsec_atmos_phy_tb, &
733  time_dtsec_atmos_phy_ch, &
734  time_dtsec_atmos_phy_ae, &
735  time_dtsec_ocean, &
736  time_dtsec_land, &
737  time_dtsec_urban )
738  else
739  if ( time_dt_wallclock_check_unit == '' ) then
740  if( io_l ) write(io_fid_log,*) '*** Not found TIME_DT_WALLCLOCK_CHECK_UNIT. TIME_DURATION_UNIT is used.'
741  time_dt_wallclock_check_unit = time_duration_unit
742  endif
743  call calendar_unit2sec( time_dtsec_wallclock_check, time_dt_wallclock_check, time_dt_wallclock_check_unit )
744  time_dtsec_wallclock_check = max( time_dtsec_wallclock_check, time_dtsec_atmos_dyn*time_nstep_atmos_dyn )
745  endif
746 
747  time_dstep_wallclock_check = int( time_dtsec_wallclock_check / time_dtsec )
748 
749  time_wallclock_safe = max( min( time_wallclock_safe, 1.0_dp ), 0.0_dp )
750  time_wallclock_start = prc_mpitime()
751 
752  if( io_l ) write(io_fid_log,'(1x,A,F10.1,A)') '*** This job stops after ', &
753  time_wallclock_limit * time_wallclock_safe, ' seconds.'
754  if( io_l ) write(io_fid_log,'(1x,A,F10.3,A,I8,A)') '*** Step interval for check : ', time_dtsec_wallclock_check, &
755  ' (step interval=', time_dstep_wallclock_check, ')'
756  endif
757 
758  return
759  end subroutine admin_time_setup
760 
761  !-----------------------------------------------------------------------------
763  subroutine admin_time_checkstate
764  use scale_process, only: &
765  prc_mpitime
766  use scale_calendar, only: &
767  calendar_date2char
768  use scale_time, only: &
769  time_nowdate, &
770  time_nowms, &
771  time_nowstep, &
772  time_nstep, &
773  time_dstep_atmos_dyn, &
774  time_dstep_atmos_phy_cp, &
775  time_dstep_atmos_phy_mp, &
776  time_dstep_atmos_phy_rd, &
777  time_dstep_atmos_phy_sf, &
778  time_dstep_atmos_phy_tb, &
779  time_dstep_atmos_phy_ch, &
780  time_dstep_atmos_phy_ae, &
781  time_dstep_ocean, &
782  time_dstep_land, &
783  time_dstep_urban
784  implicit none
785 
786  real(DP) :: wallclock_elapse
787  character(len=27) :: nowchardate
788  !---------------------------------------------------------------------------
789 
790  time_doatmos_step = .false.
791  time_doatmos_dyn = .false.
792  time_doatmos_phy_cp = .false.
793  time_doatmos_phy_mp = .false.
794  time_doatmos_phy_rd = .false.
795  time_doatmos_phy_sf = .false.
796  time_doatmos_phy_tb = .false.
797  time_doatmos_phy_ch = .false.
798  time_doatmos_phy_ae = .false.
799  time_doocean_step = .false.
800  time_doland_step = .false.
801  time_dourban_step = .false.
802  time_doresume = .false.
803 
804  time_res_atmos_dyn = time_res_atmos_dyn + 1
805  time_res_atmos_phy_cp = time_res_atmos_phy_cp + 1
806  time_res_atmos_phy_mp = time_res_atmos_phy_mp + 1
807  time_res_atmos_phy_rd = time_res_atmos_phy_rd + 1
808  time_res_atmos_phy_sf = time_res_atmos_phy_sf + 1
809  time_res_atmos_phy_tb = time_res_atmos_phy_tb + 1
810  time_res_atmos_phy_ch = time_res_atmos_phy_ch + 1
811  time_res_atmos_phy_ae = time_res_atmos_phy_ae + 1
812  time_res_ocean = time_res_ocean + 1
813  time_res_land = time_res_land + 1
814  time_res_urban = time_res_urban + 1
815  time_res_resume = time_res_resume + 1
816 
817  if ( time_res_atmos_dyn == time_dstep_atmos_dyn ) then
818  time_doatmos_step = .true.
819  time_doatmos_dyn = .true.
820  time_res_atmos_dyn = 0
821  endif
822  if ( time_res_atmos_phy_cp == time_dstep_atmos_phy_cp ) then
823  time_doatmos_step = .true.
824  time_doatmos_phy_cp = .true.
825  time_res_atmos_phy_cp = 0
826  endif
827  if ( time_res_atmos_phy_mp == time_dstep_atmos_phy_mp ) then
828  time_doatmos_step = .true.
829  time_doatmos_phy_mp = .true.
830  time_res_atmos_phy_mp = 0
831  endif
832  if ( time_res_atmos_phy_rd == time_dstep_atmos_phy_rd ) then
833  time_doatmos_step = .true.
834  time_doatmos_phy_rd = .true.
835  time_res_atmos_phy_rd = 0
836  endif
837  if ( time_res_atmos_phy_sf == time_dstep_atmos_phy_sf ) then
838  time_doatmos_step = .true.
839  time_doatmos_phy_sf = .true.
840  time_res_atmos_phy_sf = 0
841  endif
842  if ( time_res_atmos_phy_tb == time_dstep_atmos_phy_tb ) then
843  time_doatmos_step = .true.
844  time_doatmos_phy_tb = .true.
845  time_res_atmos_phy_tb = 0
846  endif
847  if ( time_res_atmos_phy_ch == time_dstep_atmos_phy_ch ) then
848  time_doatmos_step = .true.
849  time_doatmos_phy_ch = .true.
850  time_res_atmos_phy_ch = 0
851  endif
852  if ( time_res_atmos_phy_ae == time_dstep_atmos_phy_ae ) then
853  time_doatmos_step = .true.
854  time_doatmos_phy_ae = .true.
855  time_res_atmos_phy_ae = 0
856  endif
857 
858  if ( time_res_ocean == time_dstep_ocean ) then
859  time_doocean_step = .true.
860  time_res_ocean = 0
861  endif
862  if ( time_res_land == time_dstep_land ) then
863  time_doland_step = .true.
864  time_res_land = 0
865  endif
866  if ( time_res_urban == time_dstep_urban ) then
867  time_dourban_step = .true.
868  time_res_urban = 0
869  endif
870  if ( time_res_resume == time_dstep_resume ) then
871  time_doresume = .true.
872  time_res_resume = 0
873  endif
874 
875  call calendar_date2char( nowchardate, & ! [OUT]
876  time_nowdate(:), & ! [IN]
877  time_nowms ) ! [IN]
878 
879  if ( time_wallclock_limit > 0.0_dp ) then
880  wallclock_elapse = prc_mpitime() - time_wallclock_start
881  if( io_l ) write(io_fid_log,'(1x,3A,I7,A,I7,A,F10.1)') '*** TIME: ', nowchardate,' STEP:',time_nowstep, '/', time_nstep, &
882  ' WCLOCK:', wallclock_elapse
883  else
884  if( io_l ) write(io_fid_log,'(1x,3A,I7,A,I7)') '*** TIME: ', nowchardate,' STEP:',time_nowstep, '/', time_nstep
885  endif
886 
887  return
888  end subroutine admin_time_checkstate
889 
890  !-----------------------------------------------------------------------------
892  subroutine admin_time_advance
893  use scale_process, only: &
894  prc_ismaster, &
895  prc_mpitime
896  use scale_calendar, only: &
897  calendar_daysec2date, &
898  calendar_adjust_daysec, &
899  calendar_combine_daysec
900  use scale_time, only: &
901  time_dtsec, &
902  time_nowdate, &
903  time_nowms, &
904  time_nowday, &
905  time_nowsec, &
906  time_nowdaysec, &
907  time_nowstep, &
908  time_nstep, &
909  time_offset_year, &
910  time_dstep_wallclock_check
911  implicit none
912 
913  real(DP) :: WALLCLOCK_elapse
914  logical :: exists
915  !---------------------------------------------------------------------------
916 
917  time_doend = .false.
918 
919  time_nowstep = time_nowstep + 1
920  time_nowday = time_startday
921  time_nowsec = time_startsec + real(TIME_NOWSTEP-1,kind=DP) * TIME_DTSEC
922 
923  ! reallocate day & sub-day
924  call calendar_adjust_daysec( time_nowday, time_nowsec ) ! [INOUT]
925 
926  call calendar_daysec2date( time_nowdate(:), & ! [OUT]
927  time_nowms, & ! [OUT]
928  time_nowday, & ! [IN]
929  time_nowsec, & ! [IN]
930  time_offset_year ) ! [IN]
931 
932  time_nowdaysec = calendar_combine_daysec( time_nowday, time_nowsec )
933 
934  if ( time_nowstep > time_nstep ) then
935  time_doend = .true.
936  endif
937 
938  if ( time_wallclock_limit > 0.0_dp &
939  .AND. mod(time_nowstep-1,time_dstep_wallclock_check) == 0 ) then
940  wallclock_elapse = prc_mpitime() - time_wallclock_start
941 
942  if ( wallclock_elapse > time_wallclock_limit * time_wallclock_safe ) then
943  if( io_l ) write(io_fid_log,*) '*** Elapse time limit is detected. Termination operation starts.'
944  time_doend = .true.
945  endif
946  endif
947 
948  ! QUIT file control
949  if ( prc_ismaster ) then ! master node
950  inquire(file='QUIT', exist=exists)
951 
952  if( exists ) then
953  if( io_l ) write(io_fid_log,*) '*** QUIT file is found. Termination operation starts.'
954  time_doend = .true.
955  endif
956  endif
957 
958  time_doatmos_restart = .false.
959  time_doocean_restart = .false.
960  time_doland_restart = .false.
961  time_dourban_restart = .false.
962 
963  time_res_atmos_restart = time_res_atmos_restart + 1
964  time_res_ocean_restart = time_res_ocean_restart + 1
965  time_res_land_restart = time_res_land_restart + 1
966  time_res_urban_restart = time_res_urban_restart + 1
967 
968  if ( time_res_atmos_restart == time_dstep_atmos_restart ) then
969  time_doatmos_restart = .true.
970  time_res_atmos_restart = 0
971  elseif( time_doend ) then
972  time_doatmos_restart = .true.
973  endif
974 
975  if ( time_res_ocean_restart == time_dstep_ocean_restart ) then
976  time_doocean_restart = .true.
977  time_res_ocean_restart = 0
978  elseif( time_doend ) then
979  time_doocean_restart = .true.
980  endif
981 
982  if ( time_res_land_restart == time_dstep_land_restart ) then
983  time_doland_restart = .true.
984  time_res_land_restart = 0
985  elseif( time_doend ) then
986  time_doland_restart = .true.
987  endif
988 
989  if ( time_res_urban_restart == time_dstep_urban_restart ) then
990  time_dourban_restart = .true.
991  time_res_urban_restart = 0
992  elseif( time_doend ) then
993  time_dourban_restart = .true.
994  endif
995 
996  return
997  end subroutine admin_time_advance
998 
999 end module mod_admin_time
scale_time::time_nowstep
integer, public time_nowstep
current step [number]
Definition: scale_time.F90:70
scale_time::time_dtsec_ocean
real(dp), public time_dtsec_ocean
time interval of ocean step [sec]
Definition: scale_time.F90:47
scale_time::time_dstep_atmos_phy_mp
integer, public time_dstep_atmos_phy_mp
step interval of physics(microphysics)
Definition: scale_time.F90:54
gtool_file
module GTOOL_FILE
Definition: gtool_file.f90:17
scale_calendar::calendar_unit2sec
subroutine, public calendar_unit2sec(second, value, unit)
Convert several units to second.
Definition: scale_calendar.F90:366
scale_process::prc_ismaster
logical, public prc_ismaster
master process in local communicator?
Definition: scale_process.F90:79
scale_process::prc_mpistop
subroutine, public prc_mpistop
Abort MPI.
Definition: scale_process.F90:234
mod_admin_time::time_doatmos_phy_rd
logical, public time_doatmos_phy_rd
execute physics in this step? (radiation )
Definition: mod_admin_time.f90:49
mod_admin_time::time_dourban_step
logical, public time_dourban_step
execute urban component in this step?
Definition: mod_admin_time.f90:59
scale_calendar::calendar_combine_daysec
real(dp) function, public calendar_combine_daysec(absday, abssec)
Combine day and second.
Definition: scale_calendar.F90:347
scale_time::time_dtsec_atmos_phy_rd
real(dp), public time_dtsec_atmos_phy_rd
time interval of physics(radiation ) [sec]
Definition: scale_time.F90:42
mod_admin_time::time_dstep_atmos_restart
integer, public time_dstep_atmos_restart
interval of atmosphere restart [step]
Definition: mod_admin_time.f90:39
scale_time::time_nstep_atmos_dyn
integer, public time_nstep_atmos_dyn
small step of dynamics
Definition: scale_time.F90:39
scale_time::time_nowms
real(dp), public time_nowms
subsecond part of current time [millisec]
Definition: scale_time.F90:66
mod_admin_time::time_doend
logical, public time_doend
finish program in this step?
Definition: mod_admin_time.f90:62
scale_time::time_nowsec
real(dp), public time_nowsec
subday part of current time [sec]
Definition: scale_time.F90:68
scale_calendar::calendar_date2char
subroutine, public calendar_date2char(chardate, ymdhms, subsec)
Convert from gregorian date to absolute day/second.
Definition: scale_calendar.F90:497
scale_time::time_nstep
integer, public time_nstep
total steps [number]
Definition: scale_time.F90:71
mod_admin_time::time_dstep_urban_restart
integer, public time_dstep_urban_restart
interval of urban restart [step]
Definition: mod_admin_time.f90:42
scale_time::time_dtsec_atmos_phy_mp
real(dp), public time_dtsec_atmos_phy_mp
time interval of physics(microphysics) [sec]
Definition: scale_time.F90:41
scale_stdio::io_l
logical, public io_l
output log or not? (this process)
Definition: scale_stdio.F90:59
scale_time::time_dtsec_atmos_phy_ch
real(dp), public time_dtsec_atmos_phy_ch
time interval of physics(chemistry ) [sec]
Definition: scale_time.F90:45
mod_atmos_vars
module ATMOSPHERIC Variables
Definition: mod_atmos_vars.f90:17
scale_time::time_dstep_atmos_dyn
integer, public time_dstep_atmos_dyn
step interval of dynamics
Definition: scale_time.F90:52
scale_time::time_nowdaysec
real(dp), public time_nowdaysec
second of current time [sec]
Definition: scale_time.F90:69
scale_stdio
module STDIO
Definition: scale_stdio.F90:12
scale_time::time_startdaysec
real(dp), public time_startdaysec
second of start time [sec]
Definition: scale_time.F90:74
mod_admin_time::time_doatmos_restart
logical, public time_doatmos_restart
execute atmosphere restart output in this step?
Definition: mod_admin_time.f90:54
mod_admin_time::time_doocean_restart
logical, public time_doocean_restart
execute ocean restart output in this step?
Definition: mod_admin_time.f90:56
scale_time::time_nowday
integer, public time_nowday
absolute day of current time [day]
Definition: scale_time.F90:67
scale_time::time_dstep_urban
integer, public time_dstep_urban
step interval of urban step
Definition: scale_time.F90:62
mod_admin_time::time_dtsec_atmos_restart
real(dp), public time_dtsec_atmos_restart
time interval of atmosphere restart [sec]
Definition: mod_admin_time.f90:33
mod_admin_time::time_doatmos_phy_mp
logical, public time_doatmos_phy_mp
execute physics in this step? (microphysics)
Definition: mod_admin_time.f90:48
scale_time::time_dtsec_land
real(dp), public time_dtsec_land
time interval of land step [sec]
Definition: scale_time.F90:48
mod_admin_time::time_dstep_land_restart
integer, public time_dstep_land_restart
interval of land restart [step]
Definition: mod_admin_time.f90:41
scale_time::time_dstep_wallclock_check
integer, public time_dstep_wallclock_check
step interval of wallclock terminator
Definition: scale_time.F90:63
mod_admin_time::time_dtsec_resume
real(dp), public time_dtsec_resume
time interval for resume [sec]
Definition: mod_admin_time.f90:37
mod_admin_time::admin_time_setup
subroutine, public admin_time_setup(setup_TimeIntegration)
Setup.
Definition: mod_admin_time.f90:111
scale_calendar::calendar_adjust_daysec
subroutine, public calendar_adjust_daysec(absday, abssec)
Adjust day and second.
Definition: scale_calendar.F90:327
scale_time::time_dstep_ocean
integer, public time_dstep_ocean
step interval of ocean step
Definition: scale_time.F90:60
mod_admin_time::time_dstep_ocean_restart
integer, public time_dstep_ocean_restart
interval of ocean restart [step]
Definition: mod_admin_time.f90:40
mod_admin_time::time_doocean_step
logical, public time_doocean_step
execute ocean component in this step?
Definition: mod_admin_time.f90:55
scale_time::time_dtsec
real(dp), public time_dtsec
time interval of model [sec]
Definition: scale_time.F90:36
mod_admin_time::time_doresume
logical, public time_doresume
resume in this step?
Definition: mod_admin_time.f90:61
scale_calendar::calendar_cfunits2sec
real(dp) function, public calendar_cfunits2sec(cftime, cfunits, offset_year, startdaysec)
Convert time in units of the CF convention to second.
Definition: scale_calendar.F90:397
scale_time::time_dstep_atmos_phy_sf
integer, public time_dstep_atmos_phy_sf
step interval of physics(surface flux)
Definition: scale_time.F90:56
dc_types::dp
integer, parameter, public dp
Definition: dc_types.f90:27
scale_time::time_offset_year
integer, public time_offset_year
time offset [year]
Definition: scale_time.F90:73
scale_time::time_dstep_atmos_phy_tb
integer, public time_dstep_atmos_phy_tb
step interval of physics(turbulence )
Definition: scale_time.F90:57
scale_time::time_dtsec_wallclock_check
real(dp), public time_dtsec_wallclock_check
time interval of wallclock terminator [sec]
Definition: scale_time.F90:50
scale_time::time_dstep_atmos_phy_rd
integer, public time_dstep_atmos_phy_rd
step interval of physics(radiation )
Definition: scale_time.F90:55
scale_time
module TIME
Definition: scale_time.F90:15
scale_time::time_dtsec_atmos_dyn
real(dp), public time_dtsec_atmos_dyn
time interval of dynamics [sec]
Definition: scale_time.F90:38
scale_process::prc_mpitime
real(dp) function, public prc_mpitime()
Get MPI time.
Definition: scale_process.F90:715
scale_process
module PROCESS
Definition: scale_process.F90:14
scale_time::time_dtsec_atmos_phy_sf
real(dp), public time_dtsec_atmos_phy_sf
time interval of physics(surface flux) [sec]
Definition: scale_time.F90:43
scale_time::time_dtsec_atmos_phy_tb
real(dp), public time_dtsec_atmos_phy_tb
time interval of physics(turbulence ) [sec]
Definition: scale_time.F90:44
scale_time::time_dstep_land
integer, public time_dstep_land
step interval of land step
Definition: scale_time.F90:61
mod_admin_time::time_doatmos_phy_ch
logical, public time_doatmos_phy_ch
execute physics in this step? (chemistry )
Definition: mod_admin_time.f90:52
mod_admin_time::time_dtsec_urban_restart
real(dp), public time_dtsec_urban_restart
time interval of urban restart [sec]
Definition: mod_admin_time.f90:36
mod_admin_time::time_doatmos_phy_cp
logical, public time_doatmos_phy_cp
execute physics in this step? (cumulus )
Definition: mod_admin_time.f90:47
scale_time::time_dtsec_atmos_phy_ae
real(dp), public time_dtsec_atmos_phy_ae
time interval of physics(aerosol ) [sec]
Definition: scale_time.F90:46
gtool_file::filegetdatainfo
subroutine, public filegetdatainfo(basename, varname, myrank, istep, single, description, units, datatype, dim_rank, dim_name, dim_size, time_start, time_end, time_units)
Definition: gtool_file.f90:1397
scale_const
module CONSTANT
Definition: scale_const.F90:14
scale_time::time_dstep_atmos_phy_cp
integer, public time_dstep_atmos_phy_cp
step interval of physics(cumulus )
Definition: scale_time.F90:53
mod_admin_time::time_dtsec_land_restart
real(dp), public time_dtsec_land_restart
time interval of land restart [sec]
Definition: mod_admin_time.f90:35
scale_time::time_dtsec_atmos_phy_cp
real(dp), public time_dtsec_atmos_phy_cp
time interval of physics(cumulus ) [sec]
Definition: scale_time.F90:40
scale_process::prc_myrank
integer, public prc_myrank
process num in local communicator
Definition: scale_process.F90:78
mod_admin_time::time_dstep_resume
integer, public time_dstep_resume
interval for resume [step]
Definition: mod_admin_time.f90:43
mod_admin_time::time_dourban_restart
logical, public time_dourban_restart
execute urban restart output in this step?
Definition: mod_admin_time.f90:60
scale_const::const_undef8
real(dp), parameter, public const_undef8
undefined value (REAL8)
Definition: scale_const.F90:42
scale_time::time_dstep_atmos_phy_ae
integer, public time_dstep_atmos_phy_ae
step interval of physics(aerosol )
Definition: scale_time.F90:59
scale_calendar::calendar_daysec2date
subroutine, public calendar_daysec2date(ymdhms, subsec, absday, abssec, offset_year)
Convert from gregorian date to absolute day/second.
Definition: scale_calendar.F90:149
scale_time::time_dstep_atmos_phy_ch
integer, public time_dstep_atmos_phy_ch
step interval of physics(chemistry )
Definition: scale_time.F90:58
scale_prof
module profiler
Definition: scale_prof.F90:10
mod_admin_time::time_dtsec_ocean_restart
real(dp), public time_dtsec_ocean_restart
time interval of ocean restart [sec]
Definition: mod_admin_time.f90:34
mod_admin_time::time_doatmos_phy_ae
logical, public time_doatmos_phy_ae
execute physics in this step? (aerosol )
Definition: mod_admin_time.f90:53
mod_admin_time::admin_time_advance
subroutine, public admin_time_advance
Advance the time & evaluate restart & stop.
Definition: mod_admin_time.f90:893
scale_stdio::io_lnml
logical, public io_lnml
output log or not? (for namelist, this process)
Definition: scale_stdio.F90:60
mod_atmos_vars::atmos_restart_in_basename
character(len=h_long), public atmos_restart_in_basename
basename of the restart file
Definition: mod_atmos_vars.f90:59
mod_admin_time::time_doatmos_step
logical, public time_doatmos_step
execute atmosphere component in this step?
Definition: mod_admin_time.f90:45
scale_precision
module PRECISION
Definition: scale_precision.F90:13
mod_admin_time::time_doatmos_phy_tb
logical, public time_doatmos_phy_tb
execute physics in this step? (turbulence )
Definition: mod_admin_time.f90:51
mod_admin_time::time_doland_restart
logical, public time_doland_restart
execute land restart output in this step?
Definition: mod_admin_time.f90:58
mod_admin_time
module ADMIN TIME
Definition: mod_admin_time.f90:10
scale_calendar
module CALENDAR
Definition: scale_calendar.F90:15
scale_time::time_nowdate
integer, dimension(6), public time_nowdate
current time [YYYY MM DD HH MM SS]
Definition: scale_time.F90:65
scale_stdio::io_fid_conf
integer, public io_fid_conf
Config file ID.
Definition: scale_stdio.F90:55
mod_admin_time::time_doatmos_dyn
logical, public time_doatmos_dyn
execute dynamics in this step?
Definition: mod_admin_time.f90:46
scale_stdio::io_fid_log
integer, public io_fid_log
Log file ID.
Definition: scale_stdio.F90:56
mod_admin_time::time_doatmos_phy_sf
logical, public time_doatmos_phy_sf
execute physics in this step? (surface flux)
Definition: mod_admin_time.f90:50
scale_time::time_dtsec_urban
real(dp), public time_dtsec_urban
time interval of urban step [sec]
Definition: scale_time.F90:49
mod_admin_time::time_doland_step
logical, public time_doland_step
execute land component in this step?
Definition: mod_admin_time.f90:57
scale_calendar::calendar_date2daysec
subroutine, public calendar_date2daysec(absday, abssec, ymdhms, subsec, offset_year)
Convert from gregorian date to absolute day/second.
Definition: scale_calendar.F90:117
mod_admin_time::admin_time_checkstate
subroutine, public admin_time_checkstate
Evaluate component execution.
Definition: mod_admin_time.f90:764
Generated by   doxygen 1.8.13