Skip to content
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.

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

reduce orbital parameters output to log #33

Merged
merged 2 commits into from
Sep 27, 2024
Merged
Changes from all commits
Commits
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
18 changes: 10 additions & 8 deletions cicecore/drivers/nuopc/cmeps/ice_comp_nuopc.F90
Original file line number Diff line number Diff line change
Expand Up @@ -404,6 +404,7 @@ subroutine InitializeAdvertise(gcomp, importState, exportState, clock, rc)
! Determine attributes - also needed in realize phase to get grid information
!----------------------------------------------------------------------------


! Get orbital values
! Note that these values are obtained in a call to init_orbit in ice_shortwave.F90
! if CESMCOUPLED is not defined
Expand Down Expand Up @@ -1426,8 +1427,9 @@ subroutine ice_orbital_init(gcomp, clock, logunit, mastertask, rc)
type(ESMF_Time) :: CurrTime ! current time
integer :: year ! model year at current time
integer :: orb_year ! orbital year for current orbital computation
integer, save :: prev_orb_year=0 ! orbital year for previous orbital computation
logical :: lprint
logical :: first_time = .true.
logical, save :: first_time = .true.
character(len=*) , parameter :: subname = "(cice_orbital_init)"
!-------------------------------------------------------------------------------

Expand Down Expand Up @@ -1509,24 +1511,24 @@ subroutine ice_orbital_init(gcomp, clock, logunit, mastertask, rc)
return ! bail out
endif
end if

lprint = .false.
if (trim(orb_mode) == trim(orb_variable_year)) then
call ESMF_ClockGet(clock, CurrTime=CurrTime, rc=rc)
if (chkerr(rc,__LINE__,u_FILE_u)) return
call ESMF_TimeGet(CurrTime, yy=year, rc=rc)
if (chkerr(rc,__LINE__,u_FILE_u)) return
orb_year = orb_iyear + (year - orb_iyear_align)
lprint = mastertask
else
orb_year = orb_iyear
if (first_time) then
lprint = mastertask
else
lprint = .false.
end if
end if

if (orb_year .ne. prev_orb_year) then
lprint = mastertask
! this prevents the orbital print happening before the log file is opened.
if (.not. first_time) prev_orb_year = orb_year
endif
eccen = orb_eccen

call shr_orb_params(orb_year, eccen, orb_obliq, orb_mvelp, obliqr, lambm0, mvelpp, lprint)

if ( eccen == SHR_ORB_UNDEF_REAL .or. obliqr == SHR_ORB_UNDEF_REAL .or. &
Expand Down
Loading