anchor.diff:    if (AllocateStatus /= 0) STOP "*** Not enough memory for prop_const ***"
     allocate(prop_const_cap(ngen,no_freq), STAT = AllocateStatus)
     if (AllocateStatus /= 0) STOP "*** Not enough memory for prop_const_cap ***"
  -
  +  if(cap_model == 2)then !Allocating memory for sheet admittance matrix
  +    allocate(sheet_admit_matrix(ngen,no_freq,4), STAT = AllocateStatus)
  +    if (AllocateStatus /= 0) STOP "*** Not enough memory for prop_const_cap ***"
  +  endif
   
   
     cap_admit=cmplx(0.0_dp,0.0_dp,8) !initialise
  
positive.code:   deallocate(cap_eff_admit)
    deallocate(tube_eff_admit)
    deallocate(prop_const)
    deallocate (Pressure, STAT = AllocateStatus)
    deallocate (l_a, STAT = AllocateStatus)
    deallocate (l_v, STAT = AllocateStatus)
  
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
  
  