From 961d458419757f05e967f8e390f91f018987df37 Mon Sep 17 00:00:00 2001 From: Hong-Yi Li Date: Wed, 6 Dec 2023 15:55:17 -0800 Subject: [PATCH 1/2] Defining numSubSteps as an integer --- components/mosart/src/riverroute/MOSART_physics_mod.F90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/components/mosart/src/riverroute/MOSART_physics_mod.F90 b/components/mosart/src/riverroute/MOSART_physics_mod.F90 index 6f19ac18c9a7..c78fab58d495 100644 --- a/components/mosart/src/riverroute/MOSART_physics_mod.F90 +++ b/components/mosart/src/riverroute/MOSART_physics_mod.F90 @@ -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)' From 4416b05dd42c0ec92167f7548452438f49276372 Mon Sep 17 00:00:00 2001 From: Hong-Yi Li Date: Wed, 6 Dec 2023 21:31:34 -0800 Subject: [PATCH 2/2] Fixing indentation --- .../mosart/src/riverroute/MOSART_physics_mod.F90 | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/components/mosart/src/riverroute/MOSART_physics_mod.F90 b/components/mosart/src/riverroute/MOSART_physics_mod.F90 index c78fab58d495..355a1041b129 100644 --- a/components/mosart/src/riverroute/MOSART_physics_mod.F90 +++ b/components/mosart/src/riverroute/MOSART_physics_mod.F90 @@ -65,7 +65,7 @@ subroutine Euler 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 - integer :: numSubSteps + integer :: numSubSteps integer :: yr,mon,day,tod real(r8) :: myTINYVALUE character(len=*),parameter :: subname = '(Euler)' @@ -994,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 @@ -1230,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