Skip to content
Snippets Groups Projects
Commit e11521a4 authored by Dave Whipp's avatar Dave Whipp
Browse files

Modified pressure0 to be converted from positive to negative from the input...

Modified pressure0 to be converted from positive to negative from the input file, consistent with the DOUAR pressure sign convention
parent 52ce93d5
No related branches found
No related tags found
No related merge requests found
...@@ -120,10 +120,10 @@ if (params%irestart.eq.0) then ...@@ -120,10 +120,10 @@ if (params%irestart.eq.0) then
ov%wnode=0.d0 ov%wnode=0.d0
ov%wnodeiso=0.d0 ov%wnodeiso=0.d0
!ov%pressure=0.d0 !ov%pressure=0.d0
ov%pressure=params%pressure0 ov%pressure=-params%pressure0
ov%spressure=0.d0 ov%spressure=0.d0
!ov%temporary_nodal_pressure=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) call initialize_temperature (ov,params%ztemp)
if (iproc.eq.0) write(8,*) ov%nnode,' nodes in initial velocity octree' if (iproc.eq.0) write(8,*) ov%nnode,' nodes in initial velocity octree'
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment