Modified to use either an input penalty value or bulk viscosity

src/make_pressure.f90

@@ -16,9 +16,9 @@
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
 
-subroutine make_pressure (params,icon,xg,yg,zg,viscosity0,penalty,expon, &
-                          unode,vnode,wnode,temp,pressure,strain,nnode,        &
-                          r0,s0,t0,rst,plasticity_type,plasticity_parameters,eviscosity)
+subroutine make_pressure (params,icon,xg,yg,zg,viscosity0,penalty,expon,unode, &
+                         vnode,wnode,temp,pressure,strain,nnode,r0,s0,t0,rst,  &
+                         plasticity_type,plasticity_parameters,eviscosity)
 use definitions 
 use invariants
 
@@ -93,10 +93,8 @@ INCLUDE 'mpif.h'
 call mpi_comm_size (mpi_comm_world,nproc,ierr)
 call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 
-
 eps=tiny(eps)
 
-
 allocate (x(params%mpe),stat=err) ; if (err.ne.0) call stop_run ('Error alloc x in make_pressure$')
 allocate (y(params%mpe),stat=err) ; if (err.ne.0) call stop_run ('Error alloc y in make_pressure$')
 allocate (z(params%mpe),stat=err) ; if (err.ne.0) call stop_run ('Error alloc z in make_pressure$')
@@ -155,7 +153,11 @@ enddo
    div=exx+eyy+ezz
 ! end of change
 
-pressure=div*penalty*eviscosity
+if (params%bulkvisc) then
+  pressure=div*penalty
+else
+  pressure=div*penalty*eviscosity
+endif
 
 deallocate (x)
 deallocate (y)