Add scrape of outer edge of domain so water doesn't pile up. (NCAR#658)

Add scrape of outer edge of domain so water doesn't pile up

src/Routing/Noah_distr_routing_overland.F

@@ -256,9 +256,6 @@ subroutine ov_rtng( &
  REAL :: DT_FRAC,SUM_INFXS,sum_head
  !INTEGER SO8RT_D(IXRT,JXRT,3), rt_option
 
-
-
-
  !DJG ----------------------------------------------------------------------
  ! DJG BEGIN 1-D or 2-D OVERLAND FLOW ROUTING LOOP
  !DJG ---------------------------------------------------------------------
@@ -473,6 +470,7 @@ SUBROUTINE ROUTE_OVERLAND1(dt, &
  REAL, INTENT(IN), DIMENSION(XX,YY,8) :: SO8RT
  REAL*8, DIMENSION(XX,YY) :: QBDRY_tmp, DH
  REAL*8, DIMENSION(XX,YY) :: DH_tmp
+ REAL, DIMENSION(XX,YY) :: edge_adjust ! mm
 
  !!! Declare Local Variables
 
@@ -554,6 +552,7 @@ SUBROUTINE ROUTE_OVERLAND1(dt, &
 
  !yw changed as following:
  tmp_adjust=qqsfc*1000
+
  if((h(i,j) - tmp_adjust) <0 ) then
 #ifdef HYDRO_D
  print*, "Error Warning: surface head is negative: ",i,j,ixx0,jyy0, &
@@ -581,10 +580,13 @@ SUBROUTINE ROUTE_OVERLAND1(dt, &
  QBDRY_tmp(IXX0,JYY0)=QBDRY_tmp(IXX0,JYY0) - qqsfc*1000.
  QBDRYT=QBDRYT - qqsfc
  DH_tmp(IXX0,JYY0)= DH_tmp(IXX0,JYY0)-tmp_adjust
+
  end if
  end if
+
  !! End loop to route sfc water
  end if
+
  end do
  end do
 
@@ -603,6 +605,54 @@ SUBROUTINE ROUTE_OVERLAND1(dt, &
 
  H = H + DH
 
+!!! Scrape the outermost edges
+edge_adjust = 0.0
+do j=1,YY,YY-1
+ do i=1,XX
+#ifdef MPP_LAND
+ if( ((i.eq.XX).and.(right_id .lt. 0)) .or. &
+ ((i.eq.1) .and.(left_id .lt. 0)) .or. &
+ ((j.eq.1) .and.(down_id .lt. 0)) .or. &
+ ((j.eq.YY).and.(up_id .lt. 0)) ) then
+#else
+ if ((i.eq.XX).or.(i.eq.1).or.(j.eq.1) &
+ .or.(j.eq.YY )) then
+#endif
+ if (h(i,j) .GT. retent_dep(i,j)) then
+ edge_adjust(i,j) = h(i,j) - retent_dep(i,j) ! positive mm
+ end if
+
+ end if
+ end do
+end do
+
+do i=1,XX,XX-1
+ do j=1,YY
+#ifdef MPP_LAND
+ if( ((i.eq.XX).and.(right_id .lt. 0)) .or. &
+ ((i.eq.1) .and.(left_id .lt. 0)) .or. &
+ ((j.eq.1) .and.(down_id .lt. 0)) .or. &
+ ((j.eq.YY).and.(up_id .lt. 0)) ) then
+#else
+ if ((i.eq.XX).or.(i.eq.1).or.(j.eq.1) &
+ .or.(j.eq.YY )) then
+#endif
+ if (h(i,j) .GT. retent_dep(i,j)) then
+ edge_adjust(i,j) = h(i,j) - retent_dep(i,j) ! positive mm
+ end if
+
+ end if
+ end do
+end do
+
+
+#ifdef MPP_LAND
+ call MPP_LAND_COM_REAL(edge_adjust,XX,YY,99)
+#endif
+ QBDRY = QBDRY - edge_adjust ! making this negative term more negative
+ H = H - edge_adjust ! making this positive term less positive
+!!! End outermost edge scrape
+
  return
 
  !DJG ----------------------------------------------------------------------