From d5c94c68e13b8d2b75765a7b40f9c391ce174301 Mon Sep 17 00:00:00 2001
From: holger <holger.grosshans@ptb.de>
Date: Mon, 5 Oct 2020 11:19:52 +0200
Subject: [PATCH] modified: src/pressure.f90
---
src/pressure.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/pressure.f90 b/src/pressure.f90
index a873402..3224269 100644
--- a/src/pressure.f90
+++ b/src/pressure.f90
@@ -32,7 +32,7 @@
call mass2(ra,i,j,l)
de= (defect_c(i,j,l)+ra)/h2sum*urfp
- dpa= de*(1._pr/dt+2._pr*nuf/((xf(i)-xf(i-1))**2+(yf(j)-yf(j-1))**2+(zf(l)-zf(l-1))**2))
+ dpa= de*rhof*(1._pr/dt+2._pr*nuf/((xf(i)-xf(i-1))**2+(yf(j)-yf(j-1))**2+(zf(l)-zf(l-1))**2))
dex= de/(xf(i)-xf(i-1))
dey= de/(yf(j)-yf(j-1))
dez= de/(zf(l)-zf(l-1))
--
GitLab