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.

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

Bug fix in MOSART_physics_mod.F90 - numSubSteps defined as real instead of integer #6110

Merged
merged 2 commits into from
Jan 5, 2024
Merged

Bug fix in MOSART_physics_mod.F90 - numSubSteps defined as real instead of integer #6110

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
961d458
Defining numSubSteps as an integer
liho745 Dec 6, 2023
4416b05
Fixing indentation
liho745 Dec 7, 2023
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
11 changes: 6 additions & 5 deletions components/mosart/src/riverroute/MOSART_physics_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,7 +64,8 @@ subroutine Euler
integer :: iunit, idam, m, k, unitUp, cnt, ier, dd, nSubStep !local index
real(r8) :: temp_erout, localDeltaT, temp_haout, temp_Tt, temp_Tr, temp_T, temp_ha
real(r8) :: mud_erout, san_erout, temp_ehexch, temp_etexch, temp_erexch
real(r8) :: negchan, numSubSteps
real(r8) :: negchan
integer :: numSubSteps
integer :: yr,mon,day,tod
real(r8) :: myTINYVALUE
character(len=*),parameter :: subname = '(Euler)'
Expand Down Expand Up @@ -993,7 +994,7 @@ subroutine Routing_DW(iunit, nt, theDeltaT)
TRunoff%erout(iunit,nt) = -TRunoff%erin(iunit,nt)-TRunoff%erlateral(iunit,nt)
else
!TODO. If this channel is at basin outlet (downstream is ocean), use the KW method
if(rtmCTL%mask(iunit) .eq. 3) then
if(rtmCTL%mask(iunit) .eq. 3) then
call Routing_KW(iunit, nt, theDeltaT)
else
if(nt == nt_nliq) then
Expand Down Expand Up @@ -1229,9 +1230,9 @@ subroutine updateState_hillslope(iunit,nt)

if(nt==nt_nliq) then
TRunoff%yh(iunit,nt) = TRunoff%wh(iunit,nt) !/ TUnit%area(iunit) / TUnit%frac(iunit)
else
TRunoff%yh(iunit,nt) = 0._r8
end if
else
TRunoff%yh(iunit,nt) = 0._r8
end if

end subroutine updateState_hillslope

Expand Down