Skip to content

Correct the use iau_offset and add new functions to output restart files #176

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
MingjingTong-NOAA wants to merge 2 commits into
base: main
Choose a base branch
from
Open

Correct the use iau_offset and add new functions to output restart files #176

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
219e1a9
Correct the use iau_offset; Added the function to output the first time
MingjingTong-NOAA Mar 5, 2026
336ed4f
Remove iau_offset from full_coupler_mod.F90
MingjingTong-NOAA Mar 12, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions full/coupler_main.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -660,9 +660,9 @@ program coupler_main
endif

!> write out intermediate restart file when needead.
if (Time >= Time_restart) &
if ((Time >= Time_restart) .or. (nc == 1 .and. restart_first_time_step)) &
call coupler_intermediate_restart(Atm, Ice, Ocean, Ocean_state, Ocn_bc_restart, Ice_bc_restart, &
Time, Time_restart, Time_restart_current, Time_start)
Time, Time_restart, Time_restart_current, Time_start, nc)

call coupler_summarize_timestep(nc, num_cpld_calls, coupler_chksum_obj, Atm%pe, omp_sec, imb_sec)

Expand Down
23 changes: 17 additions & 6 deletions full/full_coupler_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -142,6 +142,8 @@ module full_coupler_mod
!! The format is (yr,mo,day,hr,min,sec). When restart_interval
!! is all zero, no intermediate restart file will be written out
integer, dimension(6), public :: restart_interval = (/ 0, 0, 0, 0, 0, 0/)
integer, dimension(6), public :: restart_start = (/ 0, 0, 0, 0, 0, 0/)
logical, public :: restart_first_time_step = .false.

!> The date that the current integration starts with. (See
!! force_date_from_namelist.)
Expand Down Expand Up @@ -230,7 +232,8 @@ module full_coupler_mod
concurrent, do_concurrent_radiation, use_lag_fluxes, &
check_stocks, restart_interval, do_debug, do_chksum, &
use_hyper_thread, concurrent_ice, slow_ice_with_ocean, &
do_endpoint_chksum, combined_ice_and_ocean
do_endpoint_chksum, combined_ice_and_ocean, &
restart_start, restart_first_time_step

!> coupler_clock_type derived type consist of all clock ids that will be set and used
!! in full coupler_main.
Expand Down Expand Up @@ -810,8 +813,13 @@ subroutine coupler_init(Atm, Ocean, Land, Ice, Ocean_state, Atmos_land_boundary,
else
Time_restart = fms_time_manager_set_date(date_restart(1), date_restart(2), date_restart(3), &
date_restart(4), date_restart(5), date_restart(6) )
Time_restart = fms_time_manager_increment_date(Time_restart, restart_interval(1), restart_interval(2), &
restart_interval(3), restart_interval(4), restart_interval(5), restart_interval(6) )
if (ALL(restart_start ==0)) then
Time_restart = fms_time_manager_increment_date(Time_restart, restart_interval(1), restart_interval(2), &
restart_interval(3), restart_interval(4), restart_interval(5), restart_interval(6) )
else
Time_restart = fms_time_manager_increment_date(Time_restart, restart_start(1), restart_start(2), &
restart_start(3), restart_start(4), restart_start(5), restart_start(6) )
endif
if (Time_restart <= Time) call fms_mpp_error(FATAL, &
'==>Error from program coupler: The first intermediate restart time is no larger than the start time')
endif
Expand Down Expand Up @@ -2350,7 +2358,7 @@ end subroutine coupler_update_ocean_model
!! is produced in the latter calls. Time_restart is the next timestep where the intermediate restart
!! file will be written out. Time_restart_current records the current restart time.
subroutine coupler_intermediate_restart(Atm, Ice, Ocean, Ocean_state, Ocn_bc_restart, Ice_bc_restart,&
Time_current, Time_restart, Time_restart_current, Time_start)
Time_current, Time_restart, Time_restart_current, Time_start, nc)

implicit none
type(atmos_data_type), intent(inout) :: Atm !< Atm
Expand All @@ -2364,6 +2372,7 @@ subroutine coupler_intermediate_restart(Atm, Ice, Ocean, Ocean_state, Ocn_bc_res
!! Time_restart_current records the current timestep the restart file is being written.
!! Time_restart_current does not necessary = Time_restart.
type(FmsTime_type), intent(inout) :: Time_restart, Time_restart_current
integer, intent(in) :: nc !< coupled timestep loop index
character(len=32) :: timestamp !< Time in string
integer :: outunit !< stdout

Expand All @@ -2383,8 +2392,10 @@ subroutine coupler_intermediate_restart(Atm, Ice, Ocean, Ocean_state, Ocn_bc_res
call coupler_restart(Atm, Ice, Ocean, Ocn_bc_restart, Ice_bc_restart, &
Time_current, Time_restart_current, Time_start, timestamp)

Time_restart = fms_time_manager_increment_date(Time_current, restart_interval(1), restart_interval(2), &
restart_interval(3), restart_interval(4), restart_interval(5), restart_interval(6) )
if (nc > 1) then
Time_restart = fms_time_manager_increment_date(Time_current, restart_interval(1), restart_interval(2), &
restart_interval(3), restart_interval(4), restart_interval(5), restart_interval(6) )
endif

end subroutine coupler_intermediate_restart

Expand Down
Loading