From c4cfeddb422fccc569dc2f908ca75f113191c516 Mon Sep 17 00:00:00 2001 From: Philip Pegion Date: Fri, 8 Nov 2019 16:53:17 +0000 Subject: [PATCH 1/2] zero out dissipation estimate in halo region, needed for SAR --- model/dyn_core.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/model/dyn_core.F90 b/model/dyn_core.F90 index 763654714..de9d7b7da 100644 --- a/model/dyn_core.F90 +++ b/model/dyn_core.F90 @@ -366,8 +366,8 @@ subroutine dyn_core(npx, npy, npz, ng, sphum, nq, bdt, n_split, zvir, cp, akap, endif !$OMP parallel do default(none) shared(is,ie,js,je,npz,diss_est) do k=1,npz - do j=js,je - do i=is,ie + do j=jsd,jed + do i=isd,ied diss_est(i,j,k) = 0. enddo enddo From b400a5c1b20abf0ffe83eea3460ec120f5ac3a11 Mon Sep 17 00:00:00 2001 From: Philip Pegion Date: Fri, 8 Nov 2019 17:58:47 +0000 Subject: [PATCH 2/2] fix of openmp directives for dissipation estimate --- model/dyn_core.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/model/dyn_core.F90 b/model/dyn_core.F90 index de9d7b7da..32e5ca97e 100644 --- a/model/dyn_core.F90 +++ b/model/dyn_core.F90 @@ -364,7 +364,7 @@ subroutine dyn_core(npx, npy, npz, ng, sphum, nq, bdt, n_split, zvir, cp, akap, allocate( dv(isd:ied+1,jsd:jed, npz) ) call init_ijk_mem(isd,ied+1, jsd,jed , npz, dv, 0.) endif -!$OMP parallel do default(none) shared(is,ie,js,je,npz,diss_est) +!$OMP parallel do default(none) shared(isd,ied,jsd,jed,npz,diss_est) do k=1,npz do j=jsd,jed do i=isd,ied