anchor.diff:      character(len=60) :: sub_name
       sub_name = 'calc_press_area'
       call enter_exit(sub_name,1)
  -    if(KOUNT.EQ.1)then !store initial, unstressed radius values
  +    if(KOUNT.eq.1)then !store initial, unstressed radius values
         do  ne=1,num_elems
           elem_field(ne_radius_in0,ne)=elem_field(ne_radius_in,ne)
           elem_field(ne_radius_out0,ne)=elem_field(ne_radius_out,ne)
  
positive.code:           endif
          endif
        elseif(vessel_type.eq.'elastic_alpha')then
           if(Ptm.LT.elasticity_parameters(2))THEN
            if(nn.eq.1) elem_field(ne_radius_in,ne)=R0*((Ptm*elasticity_parameters(1))+1.d0)
            if(nn.eq.2) elem_field(ne_radius_out,ne)=R0*((Ptm*elasticity_parameters(1))+1.d0)
          elseif(Ptm.lt.0.0_dp)THEN
            if(Ptm.lt.0)write(*,*) 'Transmural pressure < zero',ne,Ptm,Pblood,Ppl
            if(nn.eq.1) elem_field(ne_radius_in,ne)=R0
            if(nn.eq.2) elem_field(ne_radius_out,ne)=R0
  
negative.code:     NonZeros = num_elems*3 - fixed_pressures - fixed_flows
      !count of conservation of flow equations = sum of elements connected to nodes which have at least 2 connected elements - fixed flows
      do np=1, num_nodes
          if(elems_at_node(np,0).gt.1)then
              NonZeros = NonZeros + elems_at_node(np,0)
          endif
      enddo
  