Removed passing of viscosity0 and expon

src/make_pressure.f90

@@ -16,9 +16,8 @@
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
 
-subroutine make_pressure (params,icon,xg,yg,zg,viscosity0,penalty,expon,unode, &
-                         vnode,wnode,temp,pressure,nnode,r0,s0,t0,rst,         &
-                         eviscosity)
+subroutine make_pressure (params,icon,xg,yg,zg,penalty,unode,vnode,wnode,temp, &
+                         pressure,nnode,r0,s0,t0,rst,eviscosity)
 
 use definitions
 use invariants
@@ -36,9 +35,7 @@ include 'mpif.h'
 
 ! icon is connectivity array for the current element
 ! xg,yg,zg are the global coordinate arrays of length nnode
-! viscosity0 is viscosity
 ! penalty is incompressibility penalty factor
-! expon is nonlinear viscosity exponent
 ! unode,vnode,wnode is the velocity array (obtained from previous time step
 ! or at least containing the proper velocity at the fixed dofs)
 ! temp and pressure anre temperature and pressure
@@ -54,9 +51,7 @@ include 'mpif.h'
 type(parameters) params
 integer icon(params%mpe)
 double precision xg(nnode),yg(nnode),zg(nnode)
-double precision viscosity0
 double precision penalty
-double precision expon
 double precision unode(nnode),vnode(nnode),wnode(nnode)
 double precision temp(nnode)
 double precision pressure
@@ -89,6 +84,10 @@ double precision :: rr,ss,tt,ww
 call mpi_comm_size (mpi_comm_world,nproc,ierr)
 call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 
+
+!write (*,*) 'in make pressure'
+
+
 eps=tiny(eps)
 
 allocate (x(params%mpe),stat=err) ; if (err.ne.0) call stop_run ('Error alloc x in make_pressure$')