anchor.diff:    deallocate(cap_eff_admit)
     deallocate(tube_eff_admit)
     deallocate(prop_const)
  +  if(cap_model == 2)then
  +    deallocate(sheet_admit_matrix)
  +  endif
     deallocate (Pressure, STAT = AllocateStatus)
     deallocate (l_a, STAT = AllocateStatus)
     deallocate (l_v, STAT = AllocateStatus)
  
positive.code: 
    sub_name = 'cap_flow_admit'
    call enter_exit(sub_name,1)
  
    !     Number of non-zero entries in solution matrix.
    NonZeros=3
  
negative.code:   complex(dp), allocatable :: tube_eff_admit(:,:)
    complex(dp), allocatable :: prop_const(:,:)
    complex(dp), allocatable :: prop_const_cap(:,:)
    complex(dp), allocatable :: sheet_admit_matrix(:,:,:)
  
    complex(dp) :: daughter_admit,sister_admit,reflect_coeff,sister_current
  
  