anchor.diff:    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
   
  
positive.code:   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 ***"
  
  
  
    cap_admit=cmplx(0.0_dp,0.0_dp,8) !initialise
  
negative.code:   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)
  