Changeset 595 for codes/icosagcm/devel/src/diagnostics/diagflux.F90

Timestamp:

10/19/17 10:32:52 (7 years ago)

Author:

dubos

Message:

devel : fix halo issues with computation of energy fluxes

File:

• codes/icosagcm/devel/src/diagnostics/diagflux.F90 (modified) (1 diff)

codes/icosagcm/devel/src/diagnostics/diagflux.F90

@@ -150 +150 @@
     INTEGER :: ij, l, ij_omp_begin_ext, ij_omp_end_ext
     Rd = kappa*cpp
-    CALL distrib_level(ij_end_ext-ij_begin_ext+1,ij_omp_begin_ext,ij_omp_end_ext)
-    ij_omp_begin_ext = ij_omp_begin_ext+ij_begin_ext-1
-    ij_omp_end_ext = ij_omp_end_ext+ij_begin_ext-1
+    ! even if loops are of the _ext variant, we still need halo exchanges before reconstructing fluxes at cell centers
+    ! => loop over interior region
+    CALL distrib_level(ij_end-ij_begin+1,ij_omp_begin_ext,ij_omp_end_ext)
+    ij_omp_begin_ext = ij_omp_begin_ext+ij_begin-1
+    ij_omp_end_ext = ij_omp_end_ext+ij_begin-1
 #include "../kernels/energy_fluxes.k90"
   END SUBROUTINE compute_diagflux_energy