From 58e3aade1dc422b254405ac79279ec0ac3f86db9 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Wed, 20 Jul 2022 10:53:48 +0200
Subject: [PATCH] typoo... Latest version works with shear, EM not benchmarked
yet
---
src/processing_mod.F90 | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/src/processing_mod.F90 b/src/processing_mod.F90
index 79d80f27..5b64b26e 100644
--- a/src/processing_mod.F90
+++ b/src/processing_mod.F90
@@ -317,8 +317,10 @@ SUBROUTINE compute_nadiab_moments_z_gradients_and_interp
+ qe_taue*kernel_e(ij,:,:,:,0)*phi(:,:,:)
ENDDO
ELSEIF(parray_e(ip) .EQ. 1) THEN
+ DO ij=ijgs_e,ijge_e
nadiab_moments_e(ip,ij,:,:,:) = moments_e(ip,ij,:,:,:,updatetlevel) &
- q_o_sqrt_tau_sigma_e*kernel_e(ij,:,:,:,0)*psi(:,:,:)
+ ENDDO
ELSE
DO ij=ijgs_e,ijge_e
nadiab_moments_e(ip,ij,:,:,:) = moments_e(ip,ij,:,:,:,updatetlevel)
@@ -334,8 +336,10 @@ SUBROUTINE compute_nadiab_moments_z_gradients_and_interp
+ qi_taui*kernel_i(ij,:,:,:,0)*phi(:,:,:)
ENDDO
ELSEIF(parray_i(ip) .EQ. 1) THEN
- nadiab_moments_e(ip,ij,:,:,:) = moments_i(ip,ij,:,:,:,updatetlevel) &
+ DO ij=ijgs_i,ijge_i
+ nadiab_moments_i(ip,ij,:,:,:) = moments_i(ip,ij,:,:,:,updatetlevel) &
- q_o_sqrt_tau_sigma_i*kernel_i(ij,:,:,:,0)*psi(:,:,:)
+ ENDDO
ELSE
DO ij=ijgs_i,ijge_i
nadiab_moments_i(ip,ij,:,:,:) = moments_i(ip,ij,:,:,:,updatetlevel)
--
GitLab