anchor.diff:       call exit(0)
   endif
   
  -!!---------DESCRIPTION OF IMPORTANT PARAMETERS-----------
  +!!---------PHYSICAL PARAMETERS-----------
   !viscosity: fluid viscosity
   !density:fluid density
  -!
  -
  -!!!set the default values for the parameters that control the prq simulation, these should be controlled by user input
  -    call read_params_evaluate_prq(viscosity,density,gamma)
  -
  +!gamma:Pedley correction factor
  +density=0.10500e-02_dp !kg/cm3
  +viscosity=0.33600e-02_dp !Pa.s
  +gamma = 0.327_dp !=1.85/(4*sqrt(2))
   
   !! Allocate memory to depvar arrays
       mesh_dof=num_elems+num_nodes
  
positive.code: !need to write solution to element/nodal fields for export
      call map_solution_to_mesh(prq_solution,depvar_at_elem,depvar_at_node,mesh_dof)
  
      write(*,*) prq_solution
  
      deallocate (mesh_from_depvar, STAT = AllocateStatus)
      deallocate (depvar_at_elem, STAT = AllocateStatus)
      deallocate (depvar_at_node, STAT = AllocateStatus)
  
negative.code:     deallocate (update_resistance_entries, STAT=AllocateStatus)
      call enter_exit(sub_name,2)
    end subroutine evaluate_prq
  !
  !###################################################################################
  !
  