Skip to content

Fix zero out accumulated total precip #979

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

Closed
junwang-noaa wants to merge 6 commits into from
Closed

Fix zero out accumulated total precip #979

Changes from 3 commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
dc7b73e
fix zero out accumulated total precip
junwang-noaa Jun 9, 2025
6a180d0
add a print line
junwang-noaa Jun 9, 2025
4741e65
set iauwindow_center default value
junwang-noaa Jun 9, 2025
532a2af
accumulation time change
junwang-noaa Jun 18, 2025
95b3d95
prints
junwang-noaa Jun 18, 2025
6e6ea30
remove if statement
junwang-noaa Jun 18, 2025
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
13 changes: 12 additions & 1 deletion ccpp/CCPP_driver.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -59,6 +59,8 @@ subroutine CCPP_step (step, nblks, ierr)
! Local variables
integer :: nb, nt, ntX
integer :: ierr2
integer :: kdt_iau
logical :: iauwindow_center
! DH* 20210104 - remove kdt_rad when code to clear diagnostic buckets is removed
integer :: kdt_rad

Expand Down Expand Up @@ -167,8 +169,17 @@ subroutine CCPP_step (step, nblks, ierr)
endif

!--- determine if physics diagnostics buckets need to be cleared
iauwindow_center = .false.
if (GFS_control%iau_offset > 0) then
kdt_iau = nint(GFS_control%iau_offset*3600./GFS_control%dtp)
if (GFS_control%kdt-1 == kdt_iau) then
iauwindow_center = .true.
if( GFS_control%me == 0)print *,'in ccpp step vary, iauwindow_center=',iauwindow_center,&
'kdt=',GFS_control%kdt,'dtp=',GFS_control%dtp,'iau_offset=',GFS_control%iau_offset
endif
endif
if ((mod(GFS_control%kdt-1,GFS_control%nszero)) == 0) then
call GFS_Intdiag%phys_zero(GFS_control)
call GFS_Intdiag%phys_zero(GFS_control, iauwindow_center=iauwindow_center)
Comment thread
DusanJovic-NOAA marked this conversation as resolved.
endif

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Expand Down