From f934010f85bf2c96f7c51c1f91e5063e33ec3317 Mon Sep 17 00:00:00 2001 From: Nicholas Hannah Date: Mon, 3 Aug 2015 07:52:29 -0700 Subject: [PATCH] Initialize domore_u, domore_v variables. Closes #22 --- SIS_tracer_advect.F90 | 18 +++++++++++++++--- 1 file changed, 15 insertions(+), 3 deletions(-) diff --git a/SIS_tracer_advect.F90 b/SIS_tracer_advect.F90 index f01f7ea876..25b453dc11 100644 --- a/SIS_tracer_advect.F90 +++ b/SIS_tracer_advect.F90 @@ -485,9 +485,21 @@ subroutine advect_scalar(scalar, h_prev, h_end, uhtr, vhtr, dt, G, CS) ! (, OBC) !$OMP do do k=1,ncat domore_k(k)=1 - ! Put the remaining (total) thickness fluxes into uhr and vhr. - do j=js,je ; do I=is-1,ie ; uhr(I,j,k) = dt*uhtr(I,j,k) ; enddo ; enddo - do J=js-1,je ; do i=is,ie ; vhr(i,J,k) = dt*vhtr(i,J,k) ; enddo ; enddo + + ! Put the remaining (total) thickness fluxes into uhr and vhr. + ! Initialise domore_u, domore_v + do j=js,je + domore_u(j, k) = .false. + do I=is-1,ie + uhr(I,j,k) = dt*uhtr(I,j,k) + enddo + enddo + do J=js-1,je + domore_v(j, k) = .false. + do i=is,ie + vhr(i,J,k) = dt*vhtr(i,J,k) + enddo + enddo ! Find the previous total mass (or volume) of ice, but in the case that this ! category is now dramatically thinner than it was previously, add a tiny ! bit of extra mass to avoid nonsensical tracer concentrations. This will