anchor.diff:      use arrays,only: dp,num_elems,num_nodes,elem_field
       use diagnostics, only: enter_exit
       !local variables
  -    integer :: mesh_dof,variable_types
  -    integer, allocatable :: mesh_from_variables(:,:,:)
  -    integer, allocatable :: variables_at_node(:,:,:)
  -    integer, allocatable :: variables_at_elem(:,:,:)
  -    integer, dimension(0:2,2) :: variable_totals
  +    integer :: mesh_dof,depvar_types
  +    integer, allocatable :: mesh_from_depvar(:,:,:)
  +    integer, allocatable :: depvar_at_node(:,:,:)
  +    integer, allocatable :: depvar_at_elem(:,:,:)
  +    integer, dimension(0:2,2) :: depvar_totals
       integer, allocatable :: SparseCol(:)
       integer, allocatable :: SparseRow(:)
       integer, allocatable :: update_resistance_entries(:)
  
positive.code: end subroutine calc_press_area
  
  end module pressure_resistance_flow
  
negative.code:     real(dp) :: viscosity,density,inletbc,outletbc,grav_vect(3),gamma,total_resistance,ERR
      logical, allocatable :: FIX(:)
      logical :: ADD=.FALSE.,CONVERGED=.FALSE.
      character(len=60) :: sub_name,mesh_type,vessel_type,mechanics_type,bc_type,grav_type
      integer :: grav_dirn,no,depvar,KOUNT,nz,ne,SOLVER_FLAG
      real(dp) :: MIN_ERR,N_MIN_ERR,ptrans,beta,Go_artery,Ptm_max,grav_factor,pleural_density
  
      sub_name = 'evaluate_prq'
      call enter_exit(sub_name,1)
  