From efd1b513cf556de9cd80e7badce67156e0aaf0b2 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 27 Aug 2012 15:43:17 +0000 Subject: [PATCH] fixed bug in divergenceFDM --- code/math.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/code/math.f90 b/code/math.f90 index 18f3dae0f..4b320a552 100644 --- a/code/math.f90 +++ b/code/math.f90 @@ -3212,7 +3212,7 @@ function math_divergenceFDM(geomdim,order,field) + (/1_pInt,1_pInt,1_pInt/) do l = 1_pInt, vec_tens math_divergenceFDM(i+1_pInt,j+1_pInt,k+1_pInt,l) = math_divergenceFDM(i+1_pInt,j+1_pInt,k+1_pInt,l) & - + FDcoefficient(m,order) * & + + FDcoefficient(m,order+1_pInt) * & ((field(coordinates(1,1),coordinates(1,2),coordinates(1,3),l,1)- & field(coordinates(2,1),coordinates(2,2),coordinates(2,3),l,1))*real(res(1),pReal)/geomdim(1) +& (field(coordinates(3,1),coordinates(3,2),coordinates(3,3),l,2)- &