From e11521a409c7ae280cb14694afec8bd854bd1eb1 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Tue, 6 Dec 2011 18:38:55 +0000
Subject: [PATCH] Modified pressure0 to be converted from positive to negative
from the input file, consistent with the DOUAR pressure sign convention
---
src/define_ov.f90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/define_ov.f90 b/src/define_ov.f90
index 452e4e73..b4705e60 100644
--- a/src/define_ov.f90
+++ b/src/define_ov.f90
@@ -120,10 +120,10 @@ if (params%irestart.eq.0) then
ov%wnode=0.d0
ov%wnodeiso=0.d0
!ov%pressure=0.d0
- ov%pressure=params%pressure0
+ ov%pressure=-params%pressure0
ov%spressure=0.d0
!ov%temporary_nodal_pressure=0.d0
- ov%temporary_nodal_pressure=params%pressure0
+ ov%temporary_nodal_pressure=-params%pressure0
call initialize_temperature (ov,params%ztemp)
if (iproc.eq.0) write(8,*) ov%nnode,' nodes in initial velocity octree'
--
GitLab