diff --git a/.gitignore b/.gitignore
index 2957d22..816b589 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,6 +1,10 @@
 dist-newstyle
-*.nec
-*.gao
-*.csv
+./*.nec
+./*.gao
+./*.csv
 .debpkg
 .debpkg/*
+tools/*
+report.html
+examples/*/*gao_*.nec
+examples/*/*.nec.csv
diff --git a/README.md b/README.md
index c944845..a32f3d8 100644
--- a/README.md
+++ b/README.md
@@ -26,6 +26,17 @@
 Xnec2c-gao is an external optimizer for the antenna modeling software [Xnec2c](https://www.xnec2c.org/) written by [Neoklis Kyriazis,
 5B4AZ](http://www.5b4az.org/). It uses an genetic algorithm to find solutions in a multidimensional search space.
 
+### Features
+
+  * automatically runs Xnec2c
+  * symbols in your model
+  * mathematic expression evaluation
+  * genetic algorithm control
+  * unlimited number of genes for genetic optimization
+  * ability to specify multiple "bands" as optimization targets
+  * one output nec file per survivor
+
+
 ## Installing xnec2c-gao
 
 ### Install .deb package
@@ -158,6 +169,37 @@ To vary the parameters defined in terms of SYM cards you can use the GSYM symbol
     GSYM  scale := [0,5...1,5]
     SYM   lambda := scale * vf * c / freq
 
+## BND Cards
+
+In xnec2c frequency ranges are specified as one ore more FR Cards. You can use these and
+xnec2c-gao will automatically sweep over the bands they define and does calculation of fittnes
+for all obtained steps. However if for example you specify two FR Cards for antenna model which isn't
+resonant on both bands, the AVSWR calculation show a distorted reading of the average over data points from both bands,
+resonant and non resonant. This makes it somewhat harder for the algorithm to rate the variations of the model and hence to
+produce good results fast.
+
+To prevent this and help the algorithm optimizing your model we have BND Cards in the form :
+
+    BND <label>   <lowfreq>   <highfreq>    <steps>
+
+for 2m ham band it would be:
+
+    BND 2m  144,000  146,000  10
+
+Each BND card instructs xnec2c-gao, to run for each BND Band in isolation for each variant of your model.
+In presence of BND Cards, FR Cards are not considered in optimization. On the resulting output .nec files, only the original
+FR Cards are included , BND Cards disapear.
+
+So you could have a FR card that sweeps the entire HF band from 3.5 MHzto 30 MHz and several BND Cards for only the HAM Bands.
+That would instruct xnec2c-gao to optimize for the HAM Bands and if you run the winner .nec files with xnec2c, it will show
+the performance for the entire 3.5 Mhz to 30Mhz range.
+
+The fewer and narrower the BND Cards are the faster xnec2c-gao will run. Also the with \<steps\> specified count of calculation points,
+will speed up the calculation if \<step\> is lower. So you could also define one FR Card that sweeps the 2m band very detailed, and
+a BND Card for 2m Band that has only a few steps, to make optimization fast and final display detailed.
+
+Example : examples/awx/
+
 # From .gao file to .gao.nec files
 
 After you prepared your .gao file for your antenna model, you can invoke the optimizer.
@@ -169,12 +211,16 @@ This will yield:
      
     xnec2c-gao - a genetic algorithm optimizer for your antenna model
 
-    Usage: xnec2c-gao (-f|--gaofile FILENAME) [-v] [-d|--select-distinct] 
-                      [-s|--population-size INT] [-c|--generation-count INT]
+    Usage: xnec2c-gao [--version] (-f|--gaofile FILENAME) [-v] 
+                      [-d|--select-distinct] [-s|--population-size INT] 
+                      [-c|--generation-count INT] [-o|--optimization-mode omode] 
+                      [-y|--directional-mode dmode]
 
       Run an optimizer for GAOModel FILENAME
 
     Available options:
+      -h,--help                Show this help text
+      --version                Show version
       -f,--gaofile FILENAME    GAO Model to optimize
       -v                       How verbose to optimize (can be specified multiple
                                times)
@@ -184,13 +230,17 @@ This will yield:
       -c,--generation-count INT
                                How many generations to run the optimizer
                                (default: 10)
-      -h,--help                Show this help text
+      -o,--optimization-mode omode
+                               What should we optimize for: vswr, gain, vswr+gain
+                               (default: vswr+gain)
+      -y,--directional-mode dmode
+                               Are we optimizing a symmetrical or directive antenna:
+                               symmetrical, directive (default: symetrical)
 
     Copyright 2022 Maurizio Di Pietro DC1MDP. Program is provided "as is". Author is
     not responsible for any havoc caused by the usage of this software. Use at own
     risk.
 
-
 ## Running your optimization
 
 So to execute your optimization invoke:
@@ -230,3 +280,11 @@ This error doesn't indicate that you used the wrong number delimeter , but at le
 xnec2c-gao will run for n generations and present you the evolved survivors, you can choose to continue optimizing 
 or if you are satisfied quit the optimization in which case xnec2c-gao will write one .nec file for each survivor.
 The filename of the necfile will also contain the AVSWR calculated by xnec2c.
+
+# TODO
+
+  * examples are not clean
+  * genetic parameters command line options
+  * sanity checks for radius, segment length, ratios
+  * verbosity logging
+  * ...
diff --git a/app/Display.hs b/app/Display.hs
new file mode 100644
index 0000000..1f815b5
--- /dev/null
+++ b/app/Display.hs
@@ -0,0 +1,127 @@
+{-# LANGUAGE OverloadedStrings #-}
+
+module Display where
+
+import Control.Concurrent
+import Control.Monad.State
+import qualified Data.Map as M
+import Data.Maybe
+import Data.Text (Text)
+import qualified Data.Text as T
+import qualified Data.Text.IO as T
+import Data.These
+import GHC.Float
+import Text.Builder
+import Types
+import Utils
+
+tab :: Builder
+tab = char '\t'
+
+nl :: Builder
+nl = char '\n'
+
+padr :: Builder -> Builder
+padr = padFromRight 3 ' '
+
+padl :: Builder -> Builder
+padl = padFromLeft 9 ' '
+
+toText :: Show a => a -> Builder
+toText a = padl $ text $ T.replace "." "," $ T.pack $ show a
+
+renderFitness :: Text -> Fitness -> Text
+renderFitness _ None = ""
+renderFitness lab (Fitness swr gain fbr) =
+  run $
+    (if T.null lab then mempty else tab <> padl (text lab))
+      <> tab
+      <> string "VSWR "
+      <> padl (fixedDouble 4 $ float2Double swr)
+      <> tab
+      <> string "Raw Gain"
+      <> padl (fixedDouble 2 $ float2Double gain)
+      <> string " dBi  "
+      <> tab
+      <> string "F/B Ratio "
+      <> padl (fixedDouble 2 $ float2Double fbr)
+      <> string " dB"
+
+renderScore :: OptFun -> Phenotype -> Text
+renderScore o pt =
+  run $
+    tab <> string "Score used for optimization (higher is better) : "
+      <> tab
+      <> padl (fixedDouble 2 $ float2Double $ score o pt)
+
+renderOptModes :: GAO Text
+renderOptModes = do
+  s <- get
+  pure $
+    run $
+      tab <> string "Optimizing for " <> string (omodeShow (omode $ opts s)) <> nl
+        <> tab
+        <> string "Optimizing a "
+        <> string (dmodeShow (dmode $ opts s))
+        <> " antenna."
+
+printGenotype :: Individual -> IO ()
+printGenotype i = do
+  T.putStr $
+    run $
+      tab <> string "Genotype is " <> text (renderGenotype i) <> nl
+
+renderGenotype :: Individual -> Text
+renderGenotype i = do
+  run $
+    tab <> string "<|"
+      <> foldl (<>) mempty (fmap (\(k, v) -> text k <> ": " <> fixedDouble 4 (float2Double v) <> "|") $ M.assocs $ let Genotype g = genotype i in g)
+      <> string ">"
+
+printGeneration :: GAO ()
+printGeneration = do
+  s <- get
+  liftIO $ T.putStrLn $ run $ string "Generation " <> decimal (genNum s) <> " of " <> decimal (genCount s)
+
+printGenerationSummary :: [Individual] -> GAO ()
+printGenerationSummary is = do
+  s <- get
+  let summary =
+        run $
+          tab <> text "Generation " <> decimal (genNum s - 1) <> " selected survivors" <> nl
+            <> tab
+            <> text "========================================"
+            <> nl
+            <> nl
+            <> mconcat
+              ( flip map (zip [(1 :: Int) ..] is) $ \(n, i) ->
+                  let p = fromJust . phenotype $ i
+                      scr = score (optfun s) p
+                      linetail None = mempty
+                      linetail f@(Fitness _ gain fbr) =
+                        tab
+                          <> "Raw Gain:  "
+                          <> padl (fixedDouble 2 (float2Double gain))
+                          <> " dBi  "
+                          <> tab
+                          <> "F/B Ratio: "
+                          <> padl (fixedDouble 2 (float2Double fbr))
+                          <> " dB"
+                          <> tab
+                          <> "Score:     "
+                          <> padl (fixedDouble 2 (float2Double (let (OF h) = optfun s in h f)))
+                          <> nl
+                      entryhead = tab <> decimal n <> tab <> text "Genotype is " <> padl (text (renderGenotype i)) <> nl <> nl
+                      entryfoot = nl <> tab <> tab <> "with Score " <> padl (fixedDouble 2 (float2Double scr)) <> nl <> nl
+                      entrybody = case getPhenotype p of
+                        This (PhenotypeData _ f@(Fitness swr _ _)) ->
+                          tab <> tab <> "AVSWR:  " <> padl (fixedDouble 2 (float2Double swr)) <> linetail f
+                        These _ bs ->
+                          let ls = (\(Band lab _ _, PhenotypeData _ f@(Fitness swr _ _)) -> tab <> tab <> padl (text lab) <> tab <> "AVSWR: " <> padl (fixedDouble 2 (float2Double swr)) <> linetail f) <$> M.assocs bs
+                           in foldl (<>) mempty ls
+                        _ -> mempty
+                   in entryhead <> entrybody <> entryfoot
+              )
+  liftIO $ do
+    T.putStrLn summary
+    threadDelay 6000000
diff --git a/app/GAOParser.hs b/app/GAOParser.hs
index 12998b4..f75d5ba 100644
--- a/app/GAOParser.hs
+++ b/app/GAOParser.hs
@@ -12,7 +12,6 @@ import Control.Monad.IO.Class
 import Control.Monad.State
 import qualified Data.Char as Char
 import Data.List (foldl')
-import Data.Maybe
 import Data.Proxy
 import Data.Scientific (Scientific, toRealFloat)
 import qualified Data.Scientific as Sci
@@ -33,45 +32,48 @@ type Parser = ParsecT Void Text GAO
 sc :: Parser ()
 sc = L.space hspace1 empty empty
 
+lexeme :: Parser a -> Parser a
 lexeme = L.lexeme sc
 
+symbol :: Text -> Parser Text
 symbol = L.symbol sc
 
 lineEnd :: Parser ()
---lineEnd = try ((newline *> pure ()) <|>  (eof >> pure ()))
 lineEnd = gaodbg "NL" $ recover $ choice [void (single '\n'), eof, void newline, void crlf, void eol]
-  where recover = withRecovery $ \e -> do
-                    registerParseError e
-                    some (anySingleBut '\n')
-                    void newline
-
+  where
+    recover = withRecovery $ \e -> do
+      registerParseError e
+      void $ some (anySingleBut '\n')
+      void newline
 
 signedNum :: Parser Float
---signedNum = lexeme $ toRealFloat <$> L.signed sc L.scientific
 signedNum = lexeme $ toRealFloat <$> L.signed sc cscientific
 
 gaoparser :: Parser GAOModel
 gaoparser = do
   GAOModel <$> deck
 
+deck :: Parser (Deck Expr)
 deck = do
   zip [1 ..] <$> (card `sepEndBy` lineEnd)
 
+card :: Parser (Card Expr)
 card = do
   choice
     [ cmCard,
       ceCard,
       symCard,
       gsymCard,
+      bndCard,
       gwCard,
-      geCard,
       frCard,
+      enCard,
+      geCard,
       ldCard,
       gnCard,
       exCard,
       ekCard,
       rpCard,
-      enCard,
       symCard,
       gsymCard,
       gaoCard,
@@ -82,56 +84,70 @@ cmCard,
   ceCard,
   symCard,
   gsymCard,
+  bndCard,
   gwCard,
   geCard,
   frCard,
   ldCard,
   gnCard,
   exCard,
+  enCard,
   ekCard,
   rpCard,
   gaoCard,
   otherCard ::
     Parser (Card Expr)
 cmCard = gaodbg "CM" $ do
-  lexeme $ symbol "CM"
+  void $ lexeme $ symbol "CM"
   t <- some printChar
   return $ Card $ CM $ T.pack t
 ceCard = gaodbg "CE" $ do
-  lexeme $ symbol "CE"
+  void $ lexeme $ symbol "CE"
   t <- some printChar
   return $ Card $ CE $ T.pack t
 symCard = gaodbg "symCard" $ do
-  lexeme $ symbol "SYM"
+  void $ lexeme $ symbol "SYM"
   i <- identifier
-  lexeme $ string ":="
+  void $ lexeme $ string ":="
   e <- lexeme $ gaodbg "Expr" expr
   return $ Card $ SYM i e
 gsymCard = gaodbg "gsymCard" $ do
-  lexeme $ symbol "GSYM"
+  void $ lexeme $ symbol "GSYM"
   i <- gaodbg "Ident" identifier
-  lexeme $ string ":="
+  void $ lexeme $ string ":="
   r <- lexeme $ gaodbg "Rnge" range
   return $ Card $ GSYM i r
 gaoCard = otherCard
 gwCard = do
-  lexeme $ symbol "GW"
+  void $ lexeme $ symbol "GW"
   ct <- lexeme L.decimal
-  sc <- lexeme L.decimal
+  segc <- lexeme L.decimal
   [spx, spy, spz] <- count 3 $ lexeme expr
   [epx, epy, epz] <- count 3 $ lexeme expr
   r <- lexeme expr
-  return $ Card $ GW ct sc (Point3 spx spy spz) (Point3 epx epy epz) (Radius r)
+  return $ Card $ GW ct segc (Point3 spx spy spz) (Point3 epx epy epz) (Radius r)
+bndCard = gaodbg "bndCard" $ do
+  void $ lexeme $ symbol "BND"
+  i <- T.pack <$> lexeme (some alphaNumChar)
+  w <- do
+    l <- toRealFloat <$> lexeme cscientific
+    h <- toRealFloat <$> lexeme cscientific
+    pure (l, h)
+  Card . BND . Band i w <$> L.decimal
 geCard = otherCard
-frCard = otherCard
+frCard = gaodbg "FR" $ do
+  void $ lexeme $ symbol "FR"
+  t <- some printChar
+  return $ Card $ FR $ T.pack t
 ldCard = otherCard
 gnCard = otherCard
 exCard = otherCard
 ekCard = otherCard
 rpCard = otherCard
-
-enCard = otherCard
-
+enCard = gaodbg "enCard" $ do
+  void $ lexeme $ symbol "EN"
+  void $ lexeme $ many printChar
+  return $ Card EN
 otherCard = gaodbg "Other" $ do
   ct <- lexeme $ some alphaNumChar
   rest <- lexeme $ many printChar
@@ -139,7 +155,7 @@ otherCard = gaodbg "Other" $ do
 
 identifier :: Parser Text
 identifier = do
-  t <- lexeme ((:) <$> letterChar <*> some (choice [alphaNumChar,single '_']))
+  t <- lexeme ((:) <$> letterChar <*> many (choice [alphaNumChar, single '_']))
   return $ T.pack t
 
 range :: Parser Range
@@ -147,13 +163,17 @@ range = do
   lexeme $
     between "[" "]" $ do
       l <- signedNum
-      lexeme $ string "..."
+      void $ lexeme $ string "..."
       u <- signedNum
       return (l, u)
 
+literal :: Parser Expr
 literal = Lit <$> lexeme signedNum
 
+variable :: Parser Expr
 variable = Var <$> identifier
+
+parens :: Parser a -> Parser a
 parens p = between (symbol "(") (symbol ")") $ lexeme p
 
 term :: Parser Expr
@@ -164,15 +184,17 @@ term = do
       literal
     ]
 
+expr :: Parser Expr
 expr = makeExprParser term optable
 
+optable :: [[Operator Parser Expr]]
 optable =
   [ [ prefix "SIN" (UnOp Sin),
       prefix "COS" (UnOp Cos),
       prefix "SQRT" (UnOp Sqrt),
       prefix "-" (UnOp Negate)
     ],
-    [ binary "^" (BiOp Exp)],
+    [binary "^" (BiOp Exp)],
     [ binary "*" (BiOp Mult),
       binary "/" (BiOp Div)
     ],
diff --git a/app/Genetics.hs b/app/Genetics.hs
index b836291..8c9e601 100644
--- a/app/Genetics.hs
+++ b/app/Genetics.hs
@@ -6,9 +6,14 @@ import Control.Monad.State
 import Data.List
 import qualified Data.Map as M
 import Data.Maybe
+import Display
 import Genotype
 import System.Random
 import Types
+import Utils
+
+survivorRatio :: Double
+survivorRatio = 0.3
 
 selectSurvivors :: GAO ()
 selectSurvivors = do
@@ -16,52 +21,78 @@ selectSurvivors = do
   let selector = if selectDistinct $ opts s then nub else id
       g = generation s
       gc = fromIntegral $ length g
-      sc = floor $ gc * 0.6
+      sc = floor $ gc * survivorRatio
       g' =
         take sc $
           sortBy
-            ( \i i' ->
-                let getScore x =
-                      if fromJust (score x) < 0
-                        then 100000
-                        else fromJust $ score x
-                 in compare (getScore i) (getScore i')
+            ( flip
+                ( \i i' ->
+                    let getScore ind =
+                          if score (optfun s) (fromJust (phenotype ind)) < 0
+                            then 0
+                            else score (optfun s) $ fromJust (phenotype ind)
+                     in compare (getScore i) (getScore i')
+                )
             )
-            $ selector g
-      dupmap = zipWith (\a b -> floor $ fromIntegral (length g' - a) / 5.0) [1 ..] g'
-      g'' = concat $ zipWith replicate dupmap g'
-  modify (\s -> s {generation = g''})
-  liftIO $ mapM (print . score) g''
-  pure ()
+            (selector g)
+  --    dupmap = zipWith (\a _ -> floor $ fromIntegral (length g' - a) / (3.0 :: Double)) [1 ..] g'
+  --    g'' = concat $ zipWith replicate dupmap g'
+  -- printGenerationSummary g''
+  -- modify (\u -> u {generation = g''})
+  printGenerationSummary g'
+  modify (\u -> u {generation = g'})
 
 applyGenOperations :: GAO ()
 applyGenOperations = do
   g <- gets generation
-  g' <- mapM genetics g
-  modify (\s -> s {generation = g <> g'})
+  g' <- concat <$> mapM genetics (nub g)
+  modify (\s -> s {generation = nub (g <> g')})
 
 -- TODO clamp,crossover
-genetics :: Individual -> GAO Individual
+genetics :: Individual -> GAO [Individual]
 genetics i = do
-  g <- gets generation
-  mutate i
-  where
-    crossover = id
-    mutate i = do
-      let Genotype gtm = genotype i
-      gtm' <-
-        liftIO $
-          mapM
-            ( \v -> do
-                mf <- randomRIO @Gene (-0.1, 0.1)
-                af <- randomRIO @Gene (0, 1.0)
-                pure $
-                  if af > 0.3
-                    then abs $ v + mf
-                    else v
-            )
-            gtm
-      return $ i {genotype = Genotype gtm', phenotype = Nothing, env = M.empty, score = Nothing}
+  r <- liftIO $ randomRIO (0.0, 1.0) :: GAO Float
+  if r < 0.8
+    then mutate i
+    else crossover i
+
+crossover :: Individual -> GAO [Individual]
+crossover ind = do
+  s <- get
+  mateidx <- liftIO $ randomRIO (1, length (generation s) - 1)
+  rcidx <- liftIO $ randomRIO (1, length (prototype s) - 1)
+  let mate = generation s !! mateidx
+      (ab, ba) = recombine rcidx ind mate
+  pure $ nub [ab, ba]
+
+recombine :: Int -> Individual -> Individual -> (Individual, Individual)
+recombine i ind mate =
+  let a = M.assocs $ let (Genotype g) = genotype ind in g
+      b = M.assocs $ let (Genotype g) = genotype mate in g
+      ha = take i a
+      hb = take i b
+      ta = drop i a
+      tb = drop i b
+      ab = M.fromList $ ha ++ tb
+      ba = M.fromList $ hb ++ ta
+   in (ind {genotype = Genotype ab}, mate {genotype = Genotype ba})
+
+mutate :: Individual -> GAO [Individual]
+mutate ind = do
+  let Genotype gtm = genotype ind
+  gtm' <-
+    liftIO $
+      mapM
+        ( \v -> do
+            mf <- randomRIO @Gene (-0.1, 0.1)
+            af <- randomRIO @Gene (0, 1.0)
+            pure $
+              if af > 0.3
+                then abs $ v + mf
+                else v
+        )
+        gtm
+  return $ pure $ ind {genotype = Genotype gtm', phenotype = Nothing, environment = M.empty}
 
 genNextGen :: GAO ()
 genNextGen = do
@@ -73,4 +104,4 @@ genNextGen = do
     if delta > 0
       then (generation s <>) <$> generateNindividuals delta
       else pure $ take ps $ generation s
-  modify (\s -> s {generation = gts})
+  modify (\u -> u {generation = gts})
diff --git a/app/Genotype.hs b/app/Genotype.hs
index 889ee75..cba8abe 100644
--- a/app/Genotype.hs
+++ b/app/Genotype.hs
@@ -2,22 +2,22 @@ module Genotype where
 
 import Control.Monad.State
 import qualified Data.Map as M
-import Data.Maybe (fromJust)
 import System.Random
 import Types
 
+genInitGenotypes :: GAO ()
 genInitGenotypes = do
   proto <- extractPrototypeFromModel
-  modify (\s -> s {prototype = proto})
+  bs <- extractBandsFromModel
+  modify (\s -> s {prototype = proto, bands = bs})
   s <- get
   let ps = popSize $ opts s
   gts <- generateNindividuals ps
-  modify (\s -> s {generation = gts})
+  modify (\e -> e {generation = gts})
 
 generateNindividuals :: Int -> GAO [Individual]
 generateNindividuals n = do
-  s <- get
-  fmap (\gt -> Individual gt Nothing M.empty Nothing) <$> replicateM n proto2geno
+  fmap (\gt -> Individual gt Nothing M.empty) <$> replicateM n proto2geno
 
 extractPrototypeFromModel :: GAO (Genotype Range)
 extractPrototypeFromModel = do
@@ -26,14 +26,20 @@ extractPrototypeFromModel = do
   pure $
     Genotype $
       M.fromList $
-        (\(i, Card (GSYM s r)) -> (s, r))
+        (\(_, Card (GSYM n r)) -> (n, r))
           <$> filter
-            ( \(i, Card ct) -> case ct of
-                GSYM s r -> True
+            ( \(_, Card ct) -> case ct of
+                GSYM _ _ -> True
                 _ -> False
             )
             cs
 
+extractBandsFromModel :: GAO [Band]
+extractBandsFromModel = do
+  s <- get
+  let (GAOModel cs) = gaomodel s
+  pure $ [b | c <- filter (\(Card ct) -> case ct of BND _ -> True; _ -> False) $ map snd cs, let (Card (BND b)) = c]
+
 proto2geno :: GAO (Genotype Gene)
 proto2geno = do
   s <- get
diff --git a/app/Main.hs b/app/Main.hs
index 409f81f..a46482d 100644
--- a/app/Main.hs
+++ b/app/Main.hs
@@ -2,26 +2,18 @@
 
 module Main where
 
-import Control.Concurrent
 import Control.Monad
 import Control.Monad.IO.Class
 import Control.Monad.State
 import Data.Char
-import Data.List
-import qualified Data.Map as M
-import Data.Maybe
-import qualified Data.Text.IO as T
-import System.INotify
-import System.IO
-import Text.Pretty.Simple
-
-import Types
-import Utils
-import Options
 import GAOParser
+import Genetics
 import Genotype
+import Options
+import Output
 import Phenotype
-import Genetics
+import System.IO
+import Types
 
 main :: IO ()
 main = do
@@ -34,6 +26,7 @@ optimizer = do
   genInitGenotypes
   goGAO
   outputResult
+  tidyUp
 
 goGAO :: GAO ()
 goGAO = do
@@ -44,6 +37,7 @@ goGAO = do
   genNextGen
   rerun
 
+rerun :: GAO ()
 rerun = do
   s <- get
   if done s
@@ -53,7 +47,7 @@ rerun = do
 askProceed :: GAO ()
 askProceed = do
   s <- get
-  a <- liftIO $ do 
+  a <- liftIO $ do
     hSetBuffering stdin NoBuffering
     putStrLn $ "Are you satisfied and want to [Q]uit or [P]roceed for another " ++ show (initGenCount $ opts s) ++ " generations ?"
     a <- liftIO getChar
@@ -62,23 +56,5 @@ askProceed = do
     return a
   case toLower a of
     'q' -> pure ()
-    'p' -> modify (\s -> s {genCount = genCount s + initGenCount (opts s)}) >> goGAO
+    'p' -> modify (\u -> u {genCount = genCount u + initGenCount (opts u)}) >> goGAO
     _ -> askProceed
-
-outputResult = do
-  s <- get
-  let survivors = nub $ filter (isJust . score) $ generation s
-  mapM_ toFile $ zip [1 ..] survivors
-  liftIO $ killThread $ fromJust $ xnec2c s
-  pPrint s
-
-toFile :: (Int, Individual) -> GAO ()
-toFile (n, i) = do
-  s <- get
-  let fn =
-        gaoFile (opts s) ++ "_" ++ show n
-          ++ "_[AVSVR:"
-          ++ show (fromJust $ score i)
-          ++ "].nec"
-      p = let Phenotype t = fromJust $ phenotype i in t
-  liftIO $ T.writeFile fn p
diff --git a/app/OptFunc.hs b/app/OptFunc.hs
new file mode 100644
index 0000000..2f77678
--- /dev/null
+++ b/app/OptFunc.hs
@@ -0,0 +1,26 @@
+module OptFunc where
+
+import Types
+
+optFunc :: OptimizingMode -> DirectiveMode -> (Fitness -> Float)
+optFunc VSWR dm = \(Fitness s _ f) -> swr s + fbr dm f
+optFunc GAIN dm = \(Fitness _ g f) -> gain g + fbr dm f
+optFunc VSWRGAIN dm = \(Fitness s g f) -> swr s + gain g + fbr dm f
+
+swrMaxThrs :: Float
+swrMaxThrs = 1.5
+
+swr :: Float -> Float
+swr s
+  | s < swrMaxThrs = (swrMaxThrs - s) * 100 + (100 / s)
+  | otherwise = 100 / s
+
+gain :: Float -> Float
+gain g = 30 * g
+
+fbr :: DirectiveMode -> Float -> Float
+fbr SYMMETRICAL 0 = 1000
+fbr SYMMETRICAL f
+  | f < 0.3 = 100 / 0.3 + (f * f)
+  | otherwise = 10 / f
+fbr DIRECTIVE f = f
diff --git a/app/Options.hs b/app/Options.hs
index 956ffdc..8e66973 100644
--- a/app/Options.hs
+++ b/app/Options.hs
@@ -1,3 +1,6 @@
+{-# LANGUAGE TemplateHaskell #-}
+{-# LANGUAGE LambdaCase #-}
+
 module Options
   ( GAOOpts (..),
     parseOptions,
@@ -6,10 +9,15 @@ where
 
 import Control.Monad.IO.Class
 import Control.Monad.State
-import Data.Semigroup ((<>))
+import Data.Version (showVersion)
+import Development.GitRev (gitHash)
+import OptFunc
 import Options.Applicative
+import Paths_xnec2c_gao (version)
 import Types
+import Utils
 
+gaoopts :: Parser GAOOpts
 gaoopts =
   GAOOpts
     <$> strOption
@@ -18,48 +26,85 @@ gaoopts =
           <> metavar "FILENAME"
           <> help "GAO Model to optimize"
       )
-    <*> ( length
-            <$> many
-              ( flag'
-                  ()
-                  ( short 'v'
-                      <> help "How verbose to optimize (can be specified multiple times)"
-                      <> showDefault
-                  )
-              )
+      <*> ( length
+              <$> many
+                ( flag'
+                    ()
+                    ( short 'v'
+                        <> help "How verbose to optimize (can be specified multiple times)"
+                        <> showDefault
+                    )
+                )
+          )
+      <*> switch
+        ( short 'd'
+            <> long "select-distinct"
+            <> help "Select only distinct individuals as survivors"
         )
-    <*> switch
-      ( short 'd'
-          <> long "select-distinct"
-          <> help "Select only distinct individuals as survivors"
-      )
-    <*> option
-      auto
-      ( long "population-size"
-          <> short 's'
-          <> help "How many individuals are in one generation"
-          <> showDefault
-          <> value 20
-          <> metavar "INT"
-      )
-    <*> option
-      auto
-      ( long "generation-count"
-          <> short 'c'
-          <> help "How many generations to run the optimizer"
-          <> showDefault
-          <> value 10
-          <> metavar "INT"
-      )
+      <*> option
+        auto
+        ( long "population-size"
+            <> short 's'
+            <> help "How many individuals are in one generation"
+            <> showDefault
+            <> value 20
+            <> metavar "INT"
+        )
+      <*> option
+        auto
+        ( long "generation-count"
+            <> short 'c'
+            <> help "How many generations to run the optimizer"
+            <> showDefault
+            <> value 10
+            <> metavar "INT"
+        )
+      <*> option
+        oMode
+        ( long "optimization-mode"
+            <> short 'o'
+            <> help "What should we optimize for: vswr, gain, vswr+gain"
+            <> showDefaultWith omodeShow
+            <> value VSWRGAIN
+            <> metavar "omode"
+        )
+      <*> option
+        dMode
+        ( long "directional-mode"
+            <> short 'y'
+            <> help "Are we optimizing a symmetrical or directive antenna: symmetrical, directive"
+            <> showDefaultWith dmodeShow
+            <> value SYMMETRICAL
+            <> metavar "dmode"
+        )
+
+oMode :: ReadM OptimizingMode
+oMode =
+  str >>= \case
+    "vswr" -> return VSWR
+    "gain" -> return GAIN
+    "vswr+gain" -> return VSWRGAIN
+    _ -> readerError "Accepted optimiziation modes are 'vswr', 'gain', and 'vswr+gain'."
+
+dMode :: ReadM DirectiveMode
+dMode =
+  str >>= \case
+    "symmetrical" -> return SYMMETRICAL
+    "directive" -> return DIRECTIVE
+    _ -> readerError "Accepted directional modes are 'symmetrical' or 'directive'."
 
 parseOptions :: GAO ()
 parseOptions = do
-  gopts <- liftIO $ execParser opts
-  modify (\s -> s {opts = gopts, genCount = initGenCount gopts})
+  gopts <- liftIO $ execParser lopts
+  modify (\s -> s {opts = gopts, genCount = initGenCount gopts, optfun = OF (optFunc (omode gopts) (dmode gopts))})
   where
-    opts =
+    versionOption =
+      infoOption
+        (concat [showVersion version, " ", $(gitHash)])
+        (long "version" <> help "Show version")
+    lopts =
       info
-        (gaoopts <**> helper)
+        (helper <*> versionOption <*> gaoopts)
         ( fullDesc
             <> progDesc "Run an optimizer for GAOModel FILENAME"
             <> header "xnec2c-gao - a genetic algorithm optimizer for your antenna model"
diff --git a/app/Output.hs b/app/Output.hs
new file mode 100644
index 0000000..d1906b0
--- /dev/null
+++ b/app/Output.hs
@@ -0,0 +1,83 @@
+{-# LANGUAGE OverloadedStrings #-}
+{-# LANGUAGE TypeApplications #-}
+
+module Output where
+
+import Control.Concurrent
+import Control.Monad.Loops
+import Control.Monad.State
+import qualified Data.ByteString.Char8 as B
+import Data.IORef
+import Data.List (nub)
+import Data.Maybe
+import Data.Text (Text)
+import qualified Data.Text as T
+import qualified Data.Text.IO as T
+import System.Directory (removeFile)
+import System.FilePath.Posix.ByteString
+import System.INotify
+import System.Posix.Directory.ByteString
+import System.Process (spawnCommand, system, terminateProcess, waitForProcess)
+import Types
+import Utils
+
+startXnec2c :: RawFilePath -> GAO ()
+startXnec2c necfile = do
+  s <- get
+  unless (isJust $ xnec2c s) $ do
+    let cmd = "xnec2c  --optimize -j2  -i " ++ B.unpack necfile ++ " > /dev/null 2>&1"
+    xnec <- liftIO $ spawnCommand cmd
+    modify (\u -> u {xnec2c = Just (XN xnec)})
+
+runWithXnec :: RawFilePath -> Text -> IO Fitness
+runWithXnec necfile bpheno = do
+  csvReady <- liftIO $ newIORef False
+  T.writeFile (B.unpack necfile) bpheno
+  let touchcmd = "touch " ++ B.unpack necfile ++ ".csv"
+  void $ system touchcmd
+  withINotify $ \notify -> do
+    d <- addWatch notify [Modify, Create] (necfile <.> "csv") (\_ -> writeIORef csvReady True)
+    untilM_ (threadDelay 300000) (readIORef csvReady)
+    threadDelay 300000
+    removeWatch d
+  csvData <- drop 1 . T.lines <$> T.readFile (decodeFilePath $ necfile <.> "csv")
+  let linecount = fromIntegral $ Prelude.length csvData
+  let bvswr = (1 / linecount) * sum ((\l -> read @Float $ T.unpack $ T.splitOn "," l !! 5) <$> csvData)
+  let bgain = (1 / linecount) * sum ((\l -> read @Float $ T.unpack $ T.splitOn "," l !! 10) <$> csvData)
+  let bfbr = (1 / linecount) * sum ((\l -> read @Float $ T.unpack $ T.splitOn "," l !! 16) <$> csvData)
+
+  return $ Fitness bvswr bgain bfbr
+
+tidyUp :: GAO ()
+tidyUp = do
+  s <- get
+  cwd <- liftIO getWorkingDirectory
+  let basefile = gaoFile $ opts s
+  let necfile = cwd </> takeBaseName (B.pack basefile) <.> ".run.nec"
+  let csvfile = necfile <.> "csv"
+  liftIO $ removeFile $ decodeFilePath necfile
+  liftIO $ removeFile $ decodeFilePath csvfile
+  liftIO $ do
+    when (isJust $ xnec2c s) $
+      do
+        let (XN p) = fromJust $ xnec2c s
+        terminateProcess p
+        void $ waitForProcess p
+
+outputResult :: GAO ()
+outputResult = do
+  s <- get
+  let survivors = nub $ filter (\i -> isJust (phenotype i) && (hasFitness . fromJust . phenotype) i) $ generation s
+  mapM_ toFile $ zip [1 ..] survivors
+
+toFile :: (Int, Individual) -> GAO ()
+toFile (n, i) = do
+  s <- get
+  let (Fitness vswr _ _) = getFitness $ fromJust $ phenotype i
+      fn =
+        gaoFile (opts s) ++ "_" ++ show n
+          ++ "_[AVSVR:"
+          ++ show vswr
+          ++ "].nec"
+      p = let t = fromJust $ phenotype i in outputPhenotype t
+  liftIO $ T.writeFile fn p
diff --git a/app/Phenotype.hs b/app/Phenotype.hs
index b807a9a..62f657e 100644
--- a/app/Phenotype.hs
+++ b/app/Phenotype.hs
@@ -1,16 +1,11 @@
 {-# LANGUAGE OverloadedStrings #-}
-{-# LANGUAGE TypeApplications #-}
 
 module Phenotype where
 
---import Utils
-
-import Control.Concurrent
 import Control.Monad.ListM
-import Control.Monad.Loops
 import Control.Monad.State
+import Data.Angle
 import qualified Data.ByteString.Char8 as B
-import Data.IORef
 import Data.List (intersperse)
 import Data.Map (Map)
 import qualified Data.Map as M
@@ -18,14 +13,14 @@ import Data.Maybe
 import Data.Text (Text)
 import qualified Data.Text as T
 import qualified Data.Text.IO as T
-import Data.Text.Lazy (toStrict)
+import Data.These
+import Display
+import Output
 import System.FilePath.Posix.ByteString
-import System.INotify
 import System.Posix.Directory.ByteString
-import System.Process
 import Text.Builder
-import Text.Pretty.Simple
 import Types
+import Utils
 
 type Env = Map Symbol Float
 
@@ -35,15 +30,13 @@ genPhenotypes = do
   pts <- mapM geno2pheno $ generation s
   put $ s {generation = pts}
 
---printState
-
 geno2pheno :: Individual -> GAO Individual
-geno2pheno (Individual g _ env Nothing) = do
+geno2pheno (Individual g Nothing env) = do
   s <- get
   let (GAOModel cards) = gaomodel s
   (env'', cards') <- mapAccumM (repack evalCard) env' cards :: GAO (Env, [(Int, Card Float)])
-  let d@(Phenotype x) = renderDeck cards'
-  return $ Individual g (Just d) env'' Nothing
+  d <- renderDeck (bands s) cards'
+  return $ Individual g (Just d) env''
   where
     env' = let Genotype e = g in env `M.union` e
     repack :: (Env -> Card Expr -> GAO (Env, Card Float)) -> Env -> (Int, Card Expr) -> GAO (Env, (Int, Card Float))
@@ -51,7 +44,7 @@ geno2pheno (Individual g _ env Nothing) = do
       e' <- fst <$> f e ce
       cf <- snd <$> f e ce
       pure (e', (i, cf))
-geno2pheno i@(Individual g _ env _) = pure i
+geno2pheno i = pure i
 
 evalCard :: Env -> Card Expr -> GAO (Env, Card Float)
 evalCard env (Card (SYM s e)) = do
@@ -70,43 +63,48 @@ evalCard env (Card (CM t)) = pure (env, Card (CM t))
 evalCard env (Card (CE t)) = pure (env, Card (CE t))
 evalCard env (Card (EX t)) = pure (env, Card (EX t))
 evalCard env (Card (EK t)) = pure (env, Card (EK t))
-evalCard env (Card FR) = pure (env, Card FR)
+evalCard env (Card (FR t)) = pure (env, Card (FR t))
 evalCard env (Card LD) = pure (env, Card LD)
 evalCard env (Card GN) = pure (env, Card GN)
 evalCard env (Card RP) = pure (env, Card RP)
+evalCard env (Card (GE g)) = pure (env, Card (GE g))
 evalCard env (Card EN) = pure (env, Card EN)
 evalCard env (Card GAOP) = pure (env, Card GAOP)
 evalCard env (Card (GSYM a b)) = pure (env, Card (GSYM a b))
+evalCard env (Card (BND b)) = pure (env, Card (BND b))
 evalCard env (Card (Other a b)) = pure (env, Card (Other a b))
 
 extend :: Env -> Symbol -> Float -> Env
 extend e s v = M.insert s v e
 
 eval :: Env -> Expr -> Float
-eval env (Lit x) = x
-eval env (Var x) = case M.lookup x env of
-  Nothing -> error $ "Undefined Variable" ++ show x
-  Just v -> v
+eval _ (Lit l) = l
+eval env (Var v) = case M.lookup v env of
+  Nothing -> error $ "Undefined Variable" ++ show v
+  Just var -> var
 eval env (BiOp Add e1 e2) = eval env e1 + eval env e2
 eval env (BiOp Sub e1 e2) = eval env e1 - eval env e2
 eval env (BiOp Mult e1 e2) = eval env e1 * eval env e2
 eval env (BiOp Div e1 e2) = eval env e1 / eval env e2
 eval env (BiOp Exp e1 e2) = eval env e1 ** eval env e2
 eval env (UnOp Negate e) = negate $ eval env e
-eval env (UnOp Sin e) = sin $ eval env e
-eval env (UnOp Cos e) = cos $ eval env e
+eval env (UnOp Sin e) = sine $ Degrees $ eval env e
+eval env (UnOp Cos e) = cosine $ Degrees $ eval env e
 eval env (UnOp Sqrt e) = sqrt $ eval env e
 
-renderDeck cs = Phenotype (T.concat (renderCard . snd <$> cs))
-
-tab = char '\t'
-
-nl = char '\n'
-
-renderCard :: Card Float -> Text
-renderCard (Card (CM t)) = run $ padr (string "CM") <> tab <> text t <> nl
-renderCard (Card (CE t)) = run $ padr (string "CE") <> tab <> text t <> nl
-renderCard (Card (GW ct sc (Point3 spx spy spz) (Point3 epx epy epz) (Radius r))) =
+renderDeck :: [Band] -> Deck Float -> GAO Phenotype
+renderDeck bs cs = do
+  let p = PhenotypeData (T.concat (renderCard Nothing . snd <$> cs)) None
+      ppb = M.fromList $ zip bs $ map (\b -> PhenotypeData (T.concat (renderCard (Just b) . snd <$> cs)) None) bs
+  pure $
+    if M.null ppb
+      then Phenotype $ This p
+      else Phenotype $ These p ppb
+
+renderCard :: Maybe Band -> Card Float -> Text
+renderCard _ (Card (CM t)) = run $ padr (string "CM") <> tab <> text t <> nl
+renderCard _ (Card (CE t)) = run $ padr (string "CE") <> tab <> text t <> nl
+renderCard _ (Card (GW ct sc (Point3 spx spy spz) (Point3 epx epy epz) (Radius r))) =
   run $
     padr (string "GW") <> tab <> padl (decimal ct) <> tab <> padl (decimal sc) <> tab
       <> toText spx
@@ -123,77 +121,107 @@ renderCard (Card (GW ct sc (Point3 spx spy spz) (Point3 epx epy epz) (Radius r))
       <> tab
       <> toText r
       <> nl
-renderCard (Card (Other t1 t2)) = run $ padr (text t1) <> tab <> text (T.concat $ intersperse "\t" $ T.words t2) <> nl
-renderCard (Card (SYM _ _)) = ""
-renderCard (Card (GSYM _ _)) = ""
-
-padr = padFromRight 3 ' '
-
-padl = padFromLeft 11 ' '
-
-toText :: Show a => a -> Builder
-toText a = padl $ text $ T.replace "." "," $ T.pack $ show a
+renderCard _ (Card (SYM _ _)) = ""
+renderCard _ (Card (GSYM _ _)) = ""
+renderCard _ (Card (BND _)) = ""
+renderCard Nothing (Card (FR t)) = run $ padr (string "FR") <> tab <> text (T.concat $ intersperse "\t" $ T.words t) <> nl
+renderCard (Just _) (Card (FR _)) = ""
+renderCard Nothing (Card EN) = run $ padr (string "EN") <> nl
+renderCard (Just b) (Card EN) =
+  let low = fst $ width b
+      high = snd $ width b
+      stps = steps b
+      delta = (high - low) / fromIntegral stps
+   in run $
+        padr (string "FR") <> tab
+          <> padl (decimal (0 :: Int))
+          <> tab
+          <> padl (decimal stps)
+          <> tab
+          <> decimal (0 :: Int)
+          <> tab
+          <> decimal (0 :: Int)
+          <> tab
+          <> toText low
+          <> tab
+          <> toText delta
+          <> tab
+          <> toText high
+          <> nl
+          <> padr (string "EN")
+          <> nl
+renderCard _ (Card (Other t1 t2)) = run $ padr (text t1) <> tab <> text (T.concat $ intersperse "\t" $ T.words t2) <> nl
+renderCard _ _ = error "catchall"
 
 evalPhenotypes :: GAO ()
 evalPhenotypes = do
   s <- get
+
   printGeneration
   let ptc = Prelude.length $ generation s
   g' <-
     mapM
-      ( \(i, p) -> do
+      ( \(idx, i) -> do
           liftIO $
-            putStrLn $
-              "\nRunning Phenotype number "
-                ++ show i
-                ++ " of "
-                ++ show ptc
-          p' <- runPhenotype p
-          liftIO $ putStrLn $ "AVG VSWR: " ++ show (fromJust $ score p')
-          return p'
+            T.putStrLn $
+              run $
+                nl <> tab <> string "Phenotype "
+                  <> decimal (idx :: Int)
+                  <> string " of "
+                  <> decimal ptc
+                  <> nl
+          i' <- runPhenotype i
+          liftIO $ printGenotype i'
+          liftIO $ T.putStr "\n"
+          liftIO $ case getPhenotype $ fromJust $ phenotype i' of
+            This (PhenotypeData _ f) -> T.putStrLn $ renderFitness "" f
+            That bpm -> mapM_ (\(Band bi _ _, PhenotypeData _ f) -> T.putStrLn $ renderFitness bi f) $ M.assocs bpm
+            These _ bpm -> mapM_ (\(Band bi _ _, PhenotypeData _ f) -> T.putStrLn $ renderFitness bi f) $ M.assocs bpm
+
+          liftIO $ T.putStr "\n"
+          liftIO $ T.putStrLn $ renderScore (optfun s) $ fromJust $ phenotype i'
+          liftIO $ T.putStr "\n"
+          modes <- renderOptModes
+          liftIO $ T.putStrLn modes
+          liftIO $ T.putStrLn "\n\n\n"
+
+          return i'
       )
       $ zip [1 ..] $ generation s
-  modify (\s -> s {generation = g', genNum = genNum s + 1, done = genNum s + 1 > genCount s})
+  modify (\u -> u {generation = g', genNum = genNum u + 1, done = genNum u + 1 > genCount u})
 
 runPhenotype :: Individual -> GAO Individual
-runPhenotype i@(Individual g _ env (Just _)) = do
-  printGenotype i
-  pure i
-runPhenotype i = do
-  s <- get
-  csvReady <- liftIO $ newIORef False
-  cwd <- liftIO getWorkingDirectory
-  let necfile = cwd </> takeBaseName (B.pack $ gaoFile $ opts s) <.> "nec"
-
-  liftIO $ T.writeFile (B.unpack necfile) $ let (Phenotype t) = fromJust $ phenotype i in t
-  printGenotype i
-  unless (isJust $ xnec2c s) $ do
-    let cmd = "sleep 2 && xnec2c  --optimize -j2  -i " ++ B.unpack necfile ++ " > /dev/null 2>&1"
-    --liftIO $ putStrLn cmd
-    tid <- liftIO $ forkIO $ do
-      void $ system cmd
-    modify (\s -> s {xnec2c = Just tid})
-  liftIO $ do
-    let touchcmd = "touch " ++ B.unpack necfile ++ ".csv"
-    void $ system touchcmd
-    withINotify $ \notify -> do
-      d <- addWatch notify [Modify, Create] (necfile <.> "csv") (\_ -> writeIORef csvReady True)
-      untilM_ (threadDelay 300000) (readIORef csvReady)
-      threadDelay 300000
-      removeWatch d
-  csvData <- liftIO $ drop 1 . T.lines <$> T.readFile (decodeFilePath $ necfile <.> "csv")
-  let linecount = fromIntegral $ Prelude.length csvData
-  let avgvswr = (1 / linecount) * sum ((\l -> read @Float $ T.unpack $ T.splitOn "," l !! 5) <$> csvData)
-
-  return $ i {score = Just avgvswr}
-
-printGenotype :: Individual -> GAO ()
-printGenotype i = liftIO $ do
-  putStr "<|"
-  mapM_ (\(k, v) -> putStr $ T.unpack k ++ ": " ++ show v ++ "|") $ M.assocs $ let Genotype g = genotype i in g
-  putStr ">\n"
-
-printGeneration :: GAO ()
-printGeneration = do
-  s <- get
-  liftIO $ putStrLn $ "Generation " ++ show (genNum s) ++ " of " ++ show (genCount s)
+runPhenotype i =
+  do
+    s <- get
+    cwd <- liftIO getWorkingDirectory
+    let necfile = cwd </> takeBaseName (B.pack $ gaoFile $ opts s) <.> "run.nec"
+    startXnec2c necfile
+
+    let pt@(Phenotype ps) = fromJust $ phenotype i
+    if not $ hasFitness pt
+      then case ps of
+        This p -> liftIO $ do
+          f <- runWithXnec necfile $ data_ p
+          pure $ i {phenotype = Just $ Phenotype $ This $ PhenotypeData (data_ p) f}
+        That bdm -> do
+          brm <-
+            liftIO $
+              mapM
+                ( \(b, PhenotypeData d _) -> do
+                    f' <- runWithXnec necfile d
+                    pure (b, PhenotypeData d f')
+                )
+                $ M.assocs bdm
+          pure $ i {phenotype = Just $ Phenotype $ That $ M.fromList brm}
+        These p bdm -> do
+          brm <-
+            liftIO $
+              mapM
+                ( \(b, PhenotypeData d _) -> do
+                    f' <- runWithXnec necfile d
+                    pure (b, PhenotypeData d f')
+                )
+                $ M.assocs bdm
+          pure $ i {phenotype = Just $ Phenotype $ These p $ M.fromList brm}
+      else pure i
diff --git a/app/Types.hs b/app/Types.hs
index 3dc1402..579efec 100644
--- a/app/Types.hs
+++ b/app/Types.hs
@@ -9,16 +9,14 @@
 module Types where
 
 import Control.Applicative
-import Control.Concurrent
 import Control.Monad.IO.Class
 import Control.Monad.Reader
 import Control.Monad.State
 import Data.Map (Map)
 import qualified Data.Map as M
 import Data.Text (Text)
-import qualified Data.Text as T
-import qualified Data.Text.IO as T
-import Text.Pretty.Simple
+import Data.These
+import System.Process
 
 newtype Genotype a = Genotype {getGenotype :: Map Symbol a}
   deriving (Eq, Show, Functor, Foldable, Traversable)
@@ -26,7 +24,23 @@ newtype Genotype a = Genotype {getGenotype :: Map Symbol a}
 emptyGenotype :: Genotype a
 emptyGenotype = Genotype M.empty
 
-newtype Phenotype = Phenotype {getPhenotype :: Text}
+data Fitness = None | Fitness Float Float Float
+  deriving (Eq, Show)
+
+instance Monoid Fitness where
+  mempty = None
+
+instance Semigroup Fitness where
+  (<>) None None = None
+  (<>) None a = a
+  (<>) a None = a
+  (<>) (Fitness a b c) (Fitness a' b' c') = Fitness (a + a') (b + b') (c + c')
+
+data PhenotypeData = PhenotypeData {data_ :: Text, fitness :: Fitness}
+  deriving (Eq, Show)
+
+newtype Phenotype = Phenotype
+  {getPhenotype :: These PhenotypeData (Map Band PhenotypeData)}
   deriving (Eq, Show)
 
 type Gene = Float
@@ -36,21 +50,39 @@ type Range = (Float, Float)
 data Individual = Individual
   { genotype :: Genotype Gene,
     phenotype :: Maybe Phenotype,
-    env :: Map Symbol Float,
-    score :: Maybe Float
+    environment :: Map Symbol Float
   }
   deriving (Eq, Show)
 
+data OptimizingMode = VSWR | GAIN | VSWRGAIN
+  deriving (Eq, Ord, Enum, Show, Read)
+
+data DirectiveMode = SYMMETRICAL | DIRECTIVE
+  deriving (Eq, Ord, Enum, Show, Read)
+
 data GAOOpts = GAOOpts
   { gaoFile :: String,
     verbosity :: Int,
     selectDistinct :: Bool,
     popSize :: Int,
-    initGenCount :: Int
+    initGenCount :: Int,
+    omode :: OptimizingMode,
+    dmode :: DirectiveMode
   }
   deriving (Show)
 
-defaultGAOOpts = GAOOpts "" 0 False 20 10
+defaultGAOOpts :: GAOOpts
+defaultGAOOpts = GAOOpts "" 0 False 20 10 VSWRGAIN SYMMETRICAL
+
+newtype OptFun = OF {runOf :: Fitness -> Float}
+
+instance Show OptFun where
+  show = const "OptFun"
+
+newtype Xnec2cProc = XN ProcessHandle
+
+instance Show Xnec2cProc where
+  show = const "Xnec2c Process"
 
 data GAOEnv = GAOEnv
   { done :: Bool,
@@ -60,11 +92,14 @@ data GAOEnv = GAOEnv
     genCount :: Int,
     opts :: GAOOpts,
     gaomodel :: GAOModel,
-    xnec2c :: Maybe ThreadId
+    bands :: [Band],
+    xnec2c :: Maybe Xnec2cProc,
+    optfun :: OptFun
   }
   deriving (Show)
 
-defaultGAOEnv = GAOEnv False emptyGenotype [] 1 0 defaultGAOOpts (GAOModel []) Nothing
+defaultGAOEnv :: GAOEnv
+defaultGAOEnv = GAOEnv False emptyGenotype [] 1 0 defaultGAOOpts (GAOModel []) [] Nothing (OF (const 0))
 
 newtype GAO a = GAO {runGAO :: GAOEnv -> IO (a, GAOEnv)}
   deriving
@@ -97,7 +132,7 @@ data CardType a
   | CE Text
   | GW CardTag SegmentCount (Point3 a) (Point3 a) (Radius a)
   | GE GroundType
-  | FR
+  | FR Text
   | LD
   | GN
   | EX ExType
@@ -107,9 +142,17 @@ data CardType a
   | SYM Symbol a
   | GSYM Symbol Range
   | GAOP
+  | BND Band
   | Other Text Text
   deriving (Show)
 
+data Band = Band
+  { ident :: Text,
+    width :: Range,
+    steps :: Int
+  }
+  deriving (Eq, Ord, Show)
+
 type Symbol = Text
 
 type CardTag = Int
diff --git a/app/Utils.hs b/app/Utils.hs
index 7e5248c..00466e4 100644
--- a/app/Utils.hs
+++ b/app/Utils.hs
@@ -1,4 +1,61 @@
 module Utils where
 
+import Control.Monad.State
+import qualified Data.Map as M
+import Data.Text (Text)
+import Data.These
+import Text.Pretty.Simple
 import Types
 
+hasFitness :: Phenotype -> Bool
+hasFitness (Phenotype (This (PhenotypeData _ (Fitness {})))) = True
+hasFitness (Phenotype (That bpm)) =
+  any
+    ( \(PhenotypeData _ f) -> case f of
+        Fitness {} -> True
+        _ -> False
+    )
+    bpm
+hasFitness (Phenotype (These _ bpm)) =
+  any
+    ( \(PhenotypeData _ f) -> case f of
+        Fitness {} -> True
+        _ -> False
+    )
+    bpm
+hasFitness _ = False
+
+getFitness :: Phenotype -> Fitness
+getFitness (Phenotype (This (PhenotypeData _ f@(Fitness {})))) = f
+getFitness (Phenotype (This (PhenotypeData _ None))) = None
+getFitness (Phenotype (That bmd)) =
+  let (Fitness s g fb) = foldr (<>) None $ fmap fitness bmd
+      l = fromIntegral $ M.size bmd
+   in Fitness (s / l) (g / l) (fb / l)
+getFitness (Phenotype (These _ bmd)) = getFitness (Phenotype (That bmd))
+
+outputPhenotype :: Phenotype -> Text
+outputPhenotype (Phenotype (This (PhenotypeData d _))) = d
+outputPhenotype (Phenotype (These (PhenotypeData d _) _)) = d
+outputPhenotype (Phenotype (That _)) = error "We won't output BNDs"
+
+score :: OptFun -> Phenotype -> Float
+score (OF h) (Phenotype (This (PhenotypeData _ f@(Fitness {})))) = h f
+score _ (Phenotype (This (PhenotypeData _ None))) = 0
+score (OF h) p = case getFitness p of
+  f@(Fitness {}) -> h f
+  None -> 0
+
+omodeShow :: OptimizingMode -> String
+omodeShow VSWR = "vswr"
+omodeShow GAIN = "gain"
+omodeShow VSWRGAIN = "vswr+gain"
+
+dmodeShow :: DirectiveMode -> String
+dmodeShow SYMMETRICAL = "symetrical"
+dmodeShow DIRECTIVE = "directive"
+
+ppState :: GAO ()
+ppState = do
+  s <- get
+  pPrint s
diff --git a/examples/2m_extended_yagi/2m_extended_yagi.gao b/examples/2m_extended_yagi/2m_extended_yagi.gao
new file mode 100644
index 0000000..fcefe2b
--- /dev/null
+++ b/examples/2m_extended_yagi/2m_extended_yagi.gao
@@ -0,0 +1,23 @@
+CM --- NEC2 Input File created or edited by xnec2c 3.5 ---
+CM 2m Extended (Double Zeppelin) Yagi
+CE --- End Comments ---
+GSYM boom_delta_radi := [-0,3...0,3]
+GSYM boom_delta_dire := [-0,3...0,3]
+SYM boom_pos_refl := 0,0
+SYM boom_pos_radi := boom_pos_refl + 0,4511 + boom_delta_radi
+SYM boom_pos_dire := boom_pos_radi + 0,8605 + boom_delta_dire
+GSYM delta_refl := [-0,3...0,3]
+GSYM delta_radi := [-0,3...0,3]
+GSYM delta_dire := [-0,3...0,3]
+SYM half_refl_len := 1,5137 + delta_refl
+SYM half_radi_len := 1,2656 + delta_radi
+SYM half_dire_len := 0,4386 + delta_dire
+SYM radius := 0,0075
+GW     1    73  half_refl_len  boom_pos_refl  0,00000E+00 (-half_refl_len)  boom_pos_refl  0,00000E+00  radius
+GW     2    61  half_radi_len  boom_pos_radi  0,00000E+00 (-half_radi_len)  boom_pos_radi  0,00000E+00  radius
+GW     3    21  half_dire_len  boom_pos_dire  0,00000E+00 (-half_dire_len)  boom_pos_dire  0,00000E+00  radius
+GE     0     0   0,00000E+00   0,00000E+00    0,00000E+00  0,00000E+00    0,00000E+00    0,00000E+00  0,00000E+00
+EX     0     2    31      0  1,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+RP     0    73    73      0  0,00000E+00  0,00000E+00  2,50000E+00  5,00000E+00  0,00000E+00  0,00000E+00
+FR     0    30     0     0  1,43800E+02  0,08200E+00  1,46200E+02  0,00000E+00  0,00000E+00  0,00000E+00
+EN     0     0     0      0  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
diff --git a/examples/2m_extended_yagi/2m_extended_yagi.nec b/examples/2m_extended_yagi/2m_extended_yagi.nec
new file mode 100644
index 0000000..9802d58
--- /dev/null
+++ b/examples/2m_extended_yagi/2m_extended_yagi.nec
@@ -0,0 +1,11 @@
+CM --- NEC2 Input File created or edited by xnec2c 3.5 ---
+CM 2m Extended (Double Zeppelin) Yagi
+CE --- End Comments ---
+GW     1    73  1.51370E+00  0.00000E+00  0.00000E+00 -1.51370E+00  0.00000E+00  0.00000E+00  7.50000E-03
+GW     2    61   1.26560E+00  4.51100E-01  0.00000E+00 -1.26560E+00  4.51100E-01  0.00000E+00  7.50000E-03
+GW     3    21   4.38600E-01  8.60500E-01  0.00000E+00 -4.38600E-01  8.60500E-01  0.00000E+00  7.50000E-03
+GE     0     0   0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
+EX     0     2    31      0  1.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
+RP     0    73    73      0  0.00000E+00  0.00000E+00  2.50000E+00  5.00000E+00  0.00000E+00  0.00000E+00
+FR     0    51     0      0  1.40000E+02  2.00000E-01  1.50000E+02  0.00000E+00  0.00000E+00  0.00000E+00
+EN     0     0     0      0  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
diff --git a/examples/awx/awx.gao b/examples/awx/awx.gao
new file mode 100644
index 0000000..1d2ee20
--- /dev/null
+++ b/examples/awx/awx.gao
@@ -0,0 +1,42 @@
+CM --- NEC2 Input File created or edited by xnec2c 4.4.12 ---
+CE --- End Comments ---
+GSYM gbase := [0,5...4,0]
+SYM afreq := 145800
+SYM bfreq := 436800
+SYM c := 300000
+SYM vk := 0,96
+SYM alambda := vk * c / afreq
+SYM blambda := vk * c / bfreq
+SYM alamq := scale * (alambda / 4)
+SYM blamq := scale * (blambda / 4)
+SYM base := gbase * alambda 
+SYM feeddist := 0,02
+SYM r := 0,0025
+GSYM scale := [0,7...1,5]
+BND 2m    144,000 146,000 20
+BND 70cm  430,000 440,000 20
+SYM tbase := base + (feeddist / 2)
+SYM bbase := base - (feeddist / 2)
+SYM dpart := alamq - blamq
+SYM vpart := blamq
+GSYM angle := [5...85]
+SYM dpartx := COS(angle) * dpart
+SYM dpartz := SIN(angle) * dpart
+GW     1     3   0            0               tbase               0            0            bbase        r
+GW     2     10  (-dpartx)    0               (tbase+dpartz+vpart)      (-dpartx)    0            (tbase+dpartz)     r
+GW     3     10  (-dpartx)    0               (tbase+dpartz)            0            0            (tbase)      r
+GW     4     10  (-dpartx)    0               (bbase-(dpartz+vpart))   (-dpartx)    0            (bbase-dpartz)    r
+GW     5     10  (-dpartx)    0               (bbase-dpartz)           0            0            (bbase)      r
+GW     6     10  (dpartx)    0               (tbase+dpartz+vpart)      (dpartx)     0            (tbase+dpartz)     r
+GW     7     10  (dpartx)    0               (tbase+dpartz)            0            0            (tbase)      r
+GW     8     10  (dpartx)    0               (bbase-(dpartz+vpart))   (dpartx)     0            (bbase-dpartz)    r
+GW     9     10  (dpartx)    0               (bbase-dpartz)           0            0            (bbase)      r
+GE    -1     0   0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+FR     0    50     0     0  1,43800E+02  4,89796E-02  1,46200E+02  0,00000E+00  0,00000E+00  0,00000E+00
+FR     0    50     0     0  4,30000E+02  2,04082E-01  4,40000E+02  0,00000E+00  0,00000E+00  0,00000E+00
+LD     5     0     0      0  2,48000E+07  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+GN     2     0     0      0  1,20000E+01  1,00000E-02  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+EX     0     1     2      0  1,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+EK     -1     0     0      0  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+RP     0    18    72   1000  0,00000E+00  0,00000E+00  5,00000E+00  5,00000E+00  1,00000E+01  0,00000E+00
+EN     0     0     0      0  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
diff --git a/examples/dipole/dipole.nec.csv b/examples/dipole/dipole.nec.csv
deleted file mode 100644
index b759e0b..0000000
--- a/examples/dipole/dipole.nec.csv
+++ /dev/null
@@ -1,31 +0,0 @@
-mhz,zreal,zimag,zmag,zphase,vswr,s11,s11_real,s11_imag,s11_ang,gain_max,gain_net,gain_max_theta,gain_max_phi,gain_viewer,gain_viewer_net,fb_ratio
-143.80000000000001,197.02391927059051,242.2166745945384,312.22963058681773,50.874399575279021,10.05026700519671,-1.7342275951692165,-1.7342275951692165,2.1686750756653379,14.30549248744938,5.5921222087738043,0.76704268184034063,80,0,-10.976111670040467,-15.801191196973932,0.15391721920644311
-143.88200000000001,197.57483013501039,242.80126669392587,313.0307790155037,50.863694490759293,10.072896791848951,-1.7303056269948336,-1.7303056269948345,2.1627572000659603,14.266455692182767,5.5825170727699467,0.74943578810817169,80,0,-10.980102170303041,-15.813183454964816,0.15443254741880552
-143.964,198.12646026624597,243.38534400159958,313.83199316259856,50.852882889880661,10.095497420277754,-1.7264064352981303,-1.726406435298131,2.1568661320262583,14.227595730849389,5.5727541764954891,0.73169577557136822,80,0,-10.984063886164531,-15.825122287088652,0.15494900432778683
-144.04599999999999,198.67880585401389,243.96891307261569,314.63327611905652,50.841966884590583,10.11806931784853,-1.722529758272622,-1.7225297582726209,2.1510017048339125,14.188911503744709,5.5628332999352681,0.71382216895621564,80,0,-10.987996712665101,-15.837007843644153,0.15546658778338721
-144.12799999999999,199.23186316531829,244.55198056677662,315.43463110358886,50.830948567047059,10.140612915612436,-1.7186753368969463,-1.7186753368969452,2.1451637523522034,14.150401914961082,5.5527542155850416,0.69581448597092876,80,0,-10.991900546404509,-15.84884027601862,0.15598529565018315
-144.20999999999998,199.78562854596254,245.13455324824599,316.23606146314569,50.819830009640839,10.163128648247541,-1.7148429149035076,-1.714842914903508,2.1393521090182515,14.112065872376151,5.5425166884223849,0.67767223729034143,80,0,-10.995775285541873,-15.860619736673916,0.15650512580810927
-144.29199999999997,200.34009842196599,245.71663798492003,317.03757067315422,50.808613265003068,10.185616953989996,-1.7110322387493233,-1.7110322387493233,2.1335666098424317,14.073902287648989,5.5321204758767788,0.65939492654317444,80,0,-10.999620829795589,-15.872346379129194,0.15702607615321845
-144.37399999999997,200.89526930116079,246.29824174840326,317.83916233834498,50.797300366056923,10.208078274594046,-1.7072430575824649,-1.7072430575824658,2.1278070904049979,14.035910076197837,5.5215653277994559,0.64098205028884259,80,0,-11.003437080442922,-15.884020357953535,0.15754814459841732
-144.45599999999996,201.45113777467523,246.87937161334841,318.64084019301885,50.785893326012364,10.23051305526268,-1.703475123214004,-1.7034751232140031,2.1220733868548254,13.998088157191786,5.5108509864329234,0.62243309800411861,80,0,-11.007223940319573,-15.895641828748378,0.15807132907433019
-144.53799999999995,202.00770051827857,247.4600347567816,319.44260810122069,50.774394138378874,10.252921744580219,-1.6997281900900698,-1.6997281900900698,2.1163653359092027,13.960435453549433,5.4999771863797413,0.60374755206827668,80,0,-11.010981313819187,-15.907210948130651,0.15859562752999207
-144.61999999999995,202.56495429406601,248.04023845801632,320.24447005758998,50.762804776997655,10.275304794470953,-1.6960020152600233,-1.6960020152600221,2.11068277484959,13.922950891910885,5.4889436545713322,0.58492488773902362,80,0,-11.014709106892674,-15.918727873724983,0.15912103793369212
-144.70199999999994,203.12289595189446,248.61999009813667,321.04643018772492,50.751127196059656,10.297662660140327,-1.6922963583480568,-1.6922963583480564,2.1050255415202095,13.885633402628457,5.4777501102359958,0.56596457313371484,80,0,-11.018407227047554,-15.930192764149835,0.15964755827376642
-144.78399999999993,203.68152243080866,249.19929715916689,321.8484927483035,50.739363330089432,10.319995800001577,-1.6886109815276602,-1.6886109815276586,2.0993934743271834,13.84848191976101,5.4663962648665052,0.5468660692150511,80,0,-11.022075583347132,-15.941605778998586,0.16017518655937715
-144.86599999999993,204.24083076050394,249.7781672236776,322.65066212757142,50.72751509397353,10.342304675624531,-1.6849456494928383,-1.6849456494928381,2.0937864122364602,13.811495381060235,5.4548818221875699,0.52762882976886605,80,0,-11.025714086409639,-15.952967078828344,0.16070392082127685
-144.94799999999992,204.80081806304028,250.35660797465135,323.45294284619484,50.715584382983344,10.364589751695204,-1.6813001294280503,-1.6813001294280503,2.0882041947691556,13.774672727939949,5.4432064781225176,0.50825230137731214,80,0,-11.029322648407284,-15.964276825152488,0.16123375911273868
-145.02999999999992,205.36148155398834,250.93462719422314,324.25533955688815,50.703573072749698,10.386851495927232,-1.6776741909864923,-1.677674190986493,2.0826466620034987,13.738012905488922,5.4366352030649887,0.49400120571367623,75,0,-11.032901183065182,-15.975535180416493,0.2431241000392827
-145.11199999999991,205.92281854428072,251.51223276390249,325.05785704564079,50.691483019306979,10.409090379039341,-1.6740676062580033,-1.6740676062580042,2.0771136545682669,13.701514862427571,5.4306938713855928,0.48040116904920094,75,0,-11.036449605660081,-15.986742307996472,0.24392692131135713
-145.1939999999999,206.48482644135953,252.0894326632791,325.86050023132952,50.679316059063886,10.431306874666801,-1.6704801497481832,-1.6704801497481825,2.0716050136445938,13.665177551119889,5.4246077903590812,0.46667725120317893,75,0,-11.039967833019412,-15.997898372175314,0.24473141402873644
-145.2759999999999,207.04750275083791,252.66623496983297,326.66327416651706,50.667074008827555,10.453501459322554,-1.666911598349833,-1.6669115983498326,2.0661205809615453,13.62899992754426,5.4183768419570981,0.45282908734140559,75,0,-11.043455783519789,-16.009003538135481,0.24553757546014926
-145.35799999999989,207.61084507805376,253.24264785839264,327.46618403791484,50.654758665803101,10.475674612342965,-1.6633617313173554,-1.6633617313173545,2.0606601987931676,13.592980951273983,5.4120009016050119,0.43885630674393195,75,0,-11.046913377085648,-16.020057971946727,0.24634540291394108
-145.43999999999988,208.17485112926943,253.81867960006309,328.26923516621684,50.642371807582457,10.497826815812401,-1.6598303302449327,-1.6598303302449331,2.0552237099591761,13.557119585481827,5.4054798381661451,0.42475853280746101,75,0,-11.050340535187759,-16.031061840546442,0.24715489373913252
-145.52199999999988,208.73951871341288,254.39433856205392,329.0724330069782,50.629915192159849,10.519958554526839,-1.6563171790388596,-1.6563171790388593,2.0498109578197163,13.521414796905459,5.3988135139253508,0.41053538303354387,75,0,-11.053737180841779,-16.042015311733586,0.24796604532685285
-145.60399999999987,209.30484574348262,254.96963320689363,329.87578315081106,50.617390557924729,10.542070315932055,-1.6528220638943707,-1.6528220638943703,2.0444217862737988,13.485865555837131,5.3920017845726163,0.39618646902558741,75,0,-11.057103238606439,-16.052918554153468,0.24877885511157247
-145.68599999999986,209.87083023788188,255.54457209151815,330.67929132344369,50.604799623650543,10.56416259005705,-1.6493447732735764,-1.6493447732735769,2.0390560397584165,13.450470836117853,5.3850444991857858,0.38171139648641894,75,0,-11.060438634582102,-16.063771737281471,0.24959332057226433
-145.76799999999986,210.43747032204442,256.11916386686318,331.4829633863593,50.592144088499964,10.586235869470071,-1.6458850978802215,-1.6458850978802215,2.0337135632444454,13.41522961511035,5.3779415002134385,0.36710976520830396,75,0,-11.063743296408759,-16.074575031413893,0.25040943923380887
-145.84999999999985,211.00476422981328,256.69341727704125,332.28680533696235,50.579425632017575,10.608290649217663,-1.6424428306375107,-1.6424428306375094,2.028394202235142,13.38014087368915,5.3706926234571926,0.35238116906809314,75,0,-11.067017153264178,-16.085328607653278,0.25122720866818737
-145.93199999999985,211.57271030492413,257.26734115860648,333.09082330888396,50.566645914118403,10.630327426767799,-1.6390177666656607,-1.639017766665662,2.023097802763592,13.345203596223763,5.3632976980533522,0.33752519602028741,75,0,-11.070260135862021,-16.096032637895085,0.25204662649581816
-146.01399999999984,212.14130700250587,257.8409444400184,333.89502357245772,50.553806575095976,10.652346701963342,-1.6356097032582193,-1.6356097032582197,2.017824211389657,13.310416770558533,5.3557565464547521,0.32254142808672359,75,0,-11.073472176449897,-16.106687294817924,0.25286769038683243
-146.09599999999983,212.7105528903995,258.41423614052042,334.69941253463247,50.540909235590739,10.674348976949659,-1.6322184398627595,-1.6322184398627586,2.0125732751992169,13.275779388007633,5.3480689844115714,0.30742944135395156,75,0,-11.076653208807194,-16.117292751864813,0.25369039806234456
-146.17799999999983,213.28044665072852,258.9872253696563,335.5039967395432,50.527955496594792,10.696334756131627,-1.6288437780573704,-1.6288437780573706,2.0073448418003594,13.241290443329953,5.3402348209525519,0.29218880596102625,75,0,-11.079803168242766,-16.127849183234289,0.25451474729578538
diff --git a/examples/jpole/jpole.nec.csv b/examples/jpole/jpole.nec.csv
deleted file mode 100644
index 9171a07..0000000
--- a/examples/jpole/jpole.nec.csv
+++ /dev/null
@@ -1,41 +0,0 @@
-mhz,zreal,zimag,zmag,zphase,vswr,s11,s11_real,s11_imag,s11_ang,gain_max,gain_net,gain_max_theta,gain_max_phi,gain_viewer,gain_viewer_net,fb_ratio
-143,6.1996645926405778,33.49929591412026,34.068147407850574,79.514986774307715,11.723854502937256,-1.4853556260085448,-1.4853556260085441,16.947435238129792,111.79240735554676,0.22984567985004306,-5.1512240925814927,55,0,-5.6944890716253385,-11.075558844056875,0.40127598869306086
-143.102564,6.1060350199676536,33.614760618948338,34.164832725103039,79.704623865794019,11.928003708761418,-1.4598126858230711,-1.4598126858230722,16.922424594684962,111.62742663717127,0.24107861667491962,-5.2032731756509625,55,0,-5.7200696296982212,-11.164421422024104,0.40031095663775745
-143.205128,6.0157390346033601,33.729450024801544,34.261712086643236,79.887472838008307,12.131748141122399,-1.4351834578959741,-1.4351834578959739,16.897525163153112,111.4631795193261,0.25906176950231358,-5.247555353491844,50,0,-5.7460723670248139,-11.252689490018971,0.41471303998699144
-143.307692,5.9281814547015204,33.843373806480194,34.358656637956855,80.064570558076326,12.335955824013835,-1.411320067904086,-1.4113200679040852,16.872753505657222,111.29977525273704,0.27926200912115795,-5.2888929869400405,50,0,-5.7721732701068849,-11.340328266168083,0.41280464782052667
-143.410256,5.8432650970633606,33.956524870031551,34.455613885746253,80.23612373559277,12.5405635693009,-1.3881944567740401,-1.3881944567740399,16.848116109678962,111.13725663156036,0.29921934708346509,-5.329745007937726,50,0,-5.7983717444898542,-11.427336099511045,0.41092861179745677
-143.51282,5.7608956424608415,34.068896869857539,34.55253467594401,80.402333047201722,12.745511031454331,-1.365779472212636,-1.3657794722126371,16.823618943849073,110.97566302739367,0.31893486752974509,-5.3701105013414931,50,0,-5.824666500909875,-11.513711869781114,0.40908434416704931
-143.61538400000001,5.6809817805857712,34.180485552241379,34.649374403852214,80.56339322801638,12.95074073683902,-1.344048860045558,-1.3440488600455593,16.799267157008703,110.81502840415762,0.33841016690002779,-5.4099891829744537,50,0,-5.8510556185944651,-11.599454968468947,0.40727131076752232
-143.71794800000001,5.6034383547061317,34.291289142431488,34.746093772463567,80.719488001649154,13.156190808534237,-1.3229779890122297,-1.3229779890122282,16.77506493131358,110.65538034910956,0.35764533089826389,-5.4493812756218647,50,0,-5.8775386716715294,-11.684565278191659,0.40548898598340355
-143.82051200000001,5.528180856005152,34.401306207499822,34.842655644466419,80.870798471587051,13.361806518805787,-1.3025425942373372,-1.3025425942373379,16.751016147823595,110.49674446335213,0.37664241386164665,-5.4882873177433318,50,0,-5.9041134072591772,-11.769043138864156,0.4037368336960015
-143.92307600000001,5.4551302072775716,34.510536948907223,34.939026976723937,81.01749567512519,13.567529431796249,-1.2827198645291624,-1.2827198645291624,16.727123979313717,110.33914167584275,0.39540210537915621,-5.5267083417146283,50,0,-5.9307788649748643,-11.852889312068648,0.40201432745343796
-144.02564000000001,5.3842110151966134,34.618983444144874,35.035178648925729,81.159743508462171,13.773300882512661,-1.2634880136422244,-1.2634880136422226,16.703390897677611,110.18258829223126,0.41392480102944473,-5.564645797403621,50,0,-5.9575342992782785,-11.936104897711344,0.40032097902116565
-144.12820400000001,5.3153486685662257,34.726647773448512,35.131083061152346,81.297702953407466,13.97906860320019,-1.2448256368151258,-1.2448256368151267,16.679819198873872,110.02709945762572,0.43221262176465325,-5.6021012991854722,50,0,-5.9843774223995956,-12.018691343349721,0.39865627506534329
-144.23076800000001,5.2484738047081274,34.833533964470014,35.226716647064819,81.431525371346055,14.184775718888869,-1.2267127027715794,-1.2267127027715812,16.656410454598365,109.87268553959277,0.45026602030452312,-5.6390769668903102,50,0,-6.0113074532653643,-12.100650440460198,0.39701974225281128
-144.33333200000001,5.1835195001049623,34.939646282976426,35.322057652230242,81.561356657708941,14.390366918921574,-1.2091299260926345,-1.209129926092634,16.633165994442834,109.71935530867852,0.46808588707112325,-5.6755750022762248,50,0,-6.0383233157584817,-12.181984205105831,0.39541087727623581
-144.43589600000001,5.1204262107337781,35.04499513342531,35.417092603450286,81.687330756828516,14.595775873039651,-1.1920597741297485,-1.1920597741297487,16.610085433661471,109.56710622693782,0.48566920419431531,-5.7115987298900839,50,0,-6.0654283773462883,-12.262696311430688,0.39382947673338797
-144.53846000000001,5.0591267994734954,35.149580288855297,35.511797454590152,81.809589738677374,14.800962195146731,-1.1754833556330722,-1.1754833556330728,16.587170927433551,109.41595251083929,0.50302156919243535,-5.7471497040775432,50,0,-6.0926160929705979,-12.342787366240577,0.39227478860968473
-144.64102400000002,4.999564960712962,35.253412569488106,35.606161652036164,81.92825942207223,15.005865810263103,-1.1593842094163291,-1.1593842094163282,16.564421716978945,109.26588915515013,0.52014118139906373,-5.7822317515603379,50,0,-6.1198885807767711,-12.422261513736172,0.39074657531636547
-144.74358800000002,4.9416700867785446,35.356487444888984,35.700158370617267,82.04348345742207,15.210473142370104,-1.1437430567254037,-1.1437430567254045,16.541840820297971,109.1169360678692,0.53703964913051161,-5.8168470214148869,50,0,-6.1472325169985407,-12.501119187543939,0.38924380282628979
-144.84615200000002,4.8853982264027929,35.458820962682019,35.793785772031399,82.155366260699637,15.414704693157576,-1.1285464873895439,-1.1285464873895446,16.519425819109721,108.96907729699902,0.55370980408619697,-5.8509996712970205,50,0,-6.1746557952877454,-12.579365270670962,0.38776643888687651
-144.94871600000002,4.8306955350659484,35.560420831297883,35.887033160342334,82.264024753953294,15.618514880065341,-1.1137791229261567,-1.113779122926156,16.49717637483862,108.82231060915275,0.57015378465783129,-5.884693151021672,50,0,-6.2021561054669379,-12.657003041146442,0.3863140520812135
-145.05128000000002,4.7775121137723602,35.661295161798854,35.97989153143989,82.369568357572618,15.821853204799964,-1.0994266218232458,-1.0994266218232469,16.475091902865262,108.67663214793504,0.58637216665317393,-5.9179311880317673,50,0,-6.2297326123115493,-12.73403596699649,0.38488621932495354
-145.15384400000002,4.7257999142752274,35.761452753793407,36.072353511956521,82.472102432022979,16.024670395513571,-1.0854751825661539,-1.0854751825661539,16.453171559676402,108.53203634917789,0.60236568693363113,-5.9507175961026206,50,0,-6.2573843773198732,-12.810467660356124,0.38348251872423572
-145.25640800000002,4.6755113793051608,35.860902558765012,36.164412603929449,82.571730345194126,16.226922716747779,-1.0719112377024091,-1.0719112377024114,16.431414400471414,108.3885169805792,0.61813629097777556,-5.9830562259040896,50,0,-6.2851092116608882,-12.886301728542755,0.38210255405316579
-145.35897200000002,4.6266022791191119,35.95965367547057,36.256063508176588,82.668549051447044,16.428562216845403,-1.0587220988670745,-1.0587220988670747,16.409819326965124,108.24606679738827,0.63368466702683757,-6.014951073803692,50,0,-6.3129061883377258,-12.961541929168256,0.38074592339145696
-145.46153600000002,4.579030136535609,36.05771547398394,36.347301993326575,82.762651656870332,16.629541586364869,-1.045895585916234,-1.0458955859162327,16.38838509749592,108.10467760908311,0.64901143205667378,-6.0464062727477552,50,0,-6.3407744930285688,-13.036192197832998,0.37941222404799196
-145.56410000000002,4.5327524573471125,36.155097356457283,36.438124398691144,82.854130180500093,16.829820353488426,-1.0334196232135242,-1.0334196232135247,16.367110421460708,107.96434090230314,0.6641187041792227,-6.0774259286374308,50,0,-6.3687117842035548,-13.110256417020208,0.37810108149856403
-145.66666400000003,4.4877278964631246,36.251808020327076,36.528527022251303,82.943073616151736,17.029360655896248,-1.0212825169841591,-1.0212825169841591,16.345994106446209,107.82504881138912,0.67900899964030637,-6.1080141815234086,50,0,-6.3967155998368881,-13.183738781000603,0.37681198542108263
-145.76922800000003,4.4439212224657627,36.347858369033112,36.61850957981185,83.029560824941171,17.228108932860742,-1.009474020266576,-1.0094740202665751,16.325034304265504,107.68678914492318,0.69368082862201519,-6.138175472826001,50,0,-6.4247868981623117,-13.256643199610327,0.37554470457125372
-145.87179200000003,4.401295673092827,36.44325877154472,36.708071503305057,83.113671964761423,17.426023080073808,-0.99798364200682887,-0.99798364200682954,16.304229290266228,107.54955052637327,0.70813530958927018,-6.1679141444046692,50,0,-6.4529241892163549,-13.328973643210293,0.3742989667966235
-145.97435600000003,4.3598163946792479,36.538018515218532,36.797211796734899,83.195483071618682,17.623059380528677,-0.98680143557653777,-0.98680143557653699,16.283577521159206,107.41332277555522,0.72237306215682873,-6.1972345396078516,50,0,-6.4811267995561872,-13.400734401320868,0.37307429634023459
-146.07692000000003,4.319446944628849,36.632147178205628,36.885929957006397,83.275071771174751,17.819186876212342,-0.97591717727958138,-0.9759171772795816,16.263077392081701,107.27809530611977,0.7363974011861153,-6.2261409581644696,50,0,-6.5093912607598883,-13.471929620110473,0.3718703522966646
-146.17948400000003,4.279902852460185,36.725650815859616,36.974193652802583,83.352895899978719,18.015420847458525,-0.96526493036824701,-0.96526493036824745,16.242732651916235,107.14389284727199,0.75045685316936062,-6.2546179471971248,50,0,-6.537469869665526,-13.542544670032012,0.37068905882145931
-146.28204800000003,4.2416578345605815,36.818547810108342,37.062071232335668,83.428246809174979,18.209627127975004,-0.95494948385635858,-0.95494948385635658,16.222530189977647,107.01062891542097,0.76405869843938501,-6.2827095581505539,50,0,-6.5658527686588242,-13.612621025248764,0.36952557646582035
-146.38461200000003,4.2044286633331609,36.910843542020636,37.14953016349714,83.501579501660999,18.402806018749139,-0.94490524948372157,-0.94490524948372068,16.202473897390643,106.87832917806108,0.77744989431446054,-6.3104005453597818,50,0,-6.5942941702157594,-13.682144609890003,0.36838177554677198
-146.48717600000003,4.1681866765447584,37.002548024700531,37.236572620084068,83.572956354688742,18.594918011389154,-0.93512406395350278,-0.93512406395350189,16.182561851558653,106.74698095293711,0.79063111921140472,-6.337695302461734,50,0,-6.6227930334782865,-13.751119455151425,0.36725735564220396
-146.58974000000003,4.1329037817160961,37.093671193234549,37.323201045204065,83.642438113251174,18.785926512834646,-0.92559796285008378,-0.92559796285008367,16.162792124063085,106.61657151938928,0.80360335460839205,-6.3645983754048716,50,0,-6.6513481115130526,-13.819549841526317,0.36615200099754885
-146.69230400000004,4.0985532545743837,37.18422286350534,37.409418182369741,83.710082694044402,18.975794151571826,-0.91631935716881197,-0.9163193571688113,16.143162777169092,106.487088095304,0.81636739816648185,-6.3911141636569262,50,0,-6.6799582503628754,-13.887439812186283,0.365065414787002
-146.79486800000004,4.0651080878163137,37.274212390282827,37.495226537274625,83.775947676637699,19.164490068769052,-0.90728066721762801,-0.9072806672176279,16.123671970712017,106.35851854218159,0.82892519241636231,-6.4172470903740688,50,0,-6.7086215900598392,-13.954793872850271,0.36399723903440062
-146.89743200000004,4.0325437234583124,37.363650031086529,37.580630017698155,83.840086400040406,19.351977686373768,-0.89847490604835112,-0.89847490604834934,16.104317611164056,106.23084905023254,0.84127731604031197,-6.4430015810155616,50,0,-6.7373371282297096,-14.021616025285583,0.36294719674340548
-146.99999600000004,4.0008355631852908,37.452544891500565,37.665631603007974,83.90255144163352,19.538225252122402,-0.88989515648475415,-0.88989515648475459,16.085097855223111,106.10406746026132,0.85342512786348457,-6.4683820900743854,50,0,-6.7661032319656993,-14.087910449903569,0.36191500121932624
diff --git a/examples/multibanddipole/multibanddipole.gao b/examples/multibanddipole/multibanddipole.gao
new file mode 100644
index 0000000..becd2bb
--- /dev/null
+++ b/examples/multibanddipole/multibanddipole.gao
@@ -0,0 +1,23 @@
+CM --- NEC2 Input File created or edited by xnec2c 4.4.12 ---
+CE --- End Comments ---
+BND 70cm  430,000 440,000 25
+BND 2m    144,000 146,000 25
+BND 6m    50,000 51,000 25
+GSYM gbase := [0,5...4,0]
+SYM freq := 50000
+SYM c := 300000
+SYM vk := 0,96
+SYM lambda := vk * c / freq
+SYM lamq := scale * (lambda / 4)
+SYM base := gbase * lambda 
+GSYM r := [0,002...0,01]
+GSYM scale := [0,7...1,5]
+GW     1    15   0,00000E+00-lamq  0,00000E+00  base       0,00000E+00+lamq  0,00000E+00  base    r
+GE    -1     0   0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+FR     0    60     0     0  4,90000E+01   6,86144      4,60000E+02  0,00000E+00  0,00000E+00  0,00000E+00
+LD     5     0     0      0  2,48000E+07  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+GN     2     0     0      0  1,20000E+01  1,00000E-02  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+EX     0     1     7      0  1,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+EK     -1     0     0      0  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
+RP     0    18    72   1000  0,00000E+00  0,00000E+00  5,00000E+00  5,00000E+00  1,00000E+01  0,00000E+00
+EN     0     0     0      0  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00  0,00000E+00
diff --git a/examples/serialpole/serialpole.nec.csv b/examples/serialpole/serialpole.nec.csv
deleted file mode 100644
index adc299c..0000000
--- a/examples/serialpole/serialpole.nec.csv
+++ /dev/null
@@ -1,31 +0,0 @@
-mhz,zreal,zimag,zmag,zphase,vswr,s11,s11_real,s11_imag,s11_ang,gain_max,gain_net,gain_max_theta,gain_max_phi,gain_viewer,gain_viewer_net,fb_ratio
-143.80000000000001,77.853780958491797,-31.479860374570897,83.97733514904661,-22.015745999666166,1.9378424337305558,-9.9179724086555598,-9.9179724086555581,-5.2551366147772018,-34.665090315639574,-25.698319813791926,-26.165105272928802,15,90,-999.99000000000001,-1000.4567854591369,2.2179531730725544e-14
-143.88200000000001,78.4324149067007,-31.087542140529639,84.368708562096955,-21.621438220672108,1.9360663121152317,-9.9291848851341964,-9.9291848851341964,-5.1463175466790574,-33.947273996064062,-25.725126251868321,-26.190641257363065,15,270,-999.99000000000001,-1000.4555150054947,8.6789471989795746e-15
-143.964,79.015746191051718,-30.697166745688822,84.76912287113349,-21.230864435383108,1.9347514062098283,-9.9375038444559003,-9.9375038444558985,-5.0374601112508905,-33.229204589450489,-25.751965449150216,-26.21654021378956,15,270,-999.99000000000001,-1000.4545747646393,1.3500584531745998e-14
-144.04599999999999,79.603804489374568,-30.308788438064408,85.178567408509039,-20.84410493743302,1.9338940430840379,-9.9429364022851097,-9.9429364022851097,-4.928684390167553,-32.511674205007267,-25.778836847643241,-26.242798684646214,15,270,-999.99000000000001,-1000.453961837003,4.5323390928004254e-14
-144.12799999999999,80.196619285218304,-29.922462569522821,85.597029802457186,-20.461237317779453,1.9334903394874148,-9.9454966893333285,-9.9454966893333268,-4.820109499586553,-31.795468582984515,-25.805739877670469,-26.269413146276971,15,270,-999.99000000000001,-1000.4536732686065,4.8216373327664331e-15
-144.20999999999998,80.794219828422456,-29.538244962681432,86.024495773929573,-20.082336044799536,1.9335361930653772,-9.9452058126018912,-9.9452058126018876,-4.7118530056205001,-31.081363238884563,-25.832673955719702,-26.296379999335386,15,270,-999.99000000000001,-1000.4537060436157,8.6789471989795746e-15
-144.29199999999997,81.39663512993306,-29.156192906913095,86.460949539665492,-19.707473064404311,1.9340273087798225,-9.9420915519215427,-9.9420915519215427,-4.6040304844614601,-30.370120561857313,-25.859638485885945,-26.323695587005659,15,270,-999.99000000000001,-1000.4540571011197,1.3500584531745998e-14
-144.37399999999997,82.003893941305932,-28.776364751898548,86.906373701077072,-19.336717527243604,1.9349591971360516,-9.9361881961983691,-9.9361881961983691,-4.496754995828331,-29.662486341337733,-25.886632859060668,-26.351356189928769,15,270,-999.99000000000001,-1000.4547233308681,1.3500584531745998e-14
-144.45599999999996,82.61602473564389,-28.398819538998449,87.360749163048752,-18.970135548429802,1.9363271753405507,-9.9275363246437269,-9.9275363246437252,-4.3901365823431346,-28.959186464721441,-25.913656452317362,-26.379358022591177,15,270,-999.99000000000001,-1000.4557015702738,1.3500584531745998e-14
-144.53799999999995,83.23305569726601,-28.023618030229716,87.824055521300096,-18.607790839164075,1.9381263962180979,-9.9161823792341934,-9.9161823792341934,-4.2842819227028279,-28.260924629540192,-25.940708630124256,-26.407697250633024,15,270,-999.99000000000001,-1000.4569886205088,1.7357894397959133e-14
-144.61999999999995,83.855014680648239,-27.650821239829604,88.296270602606668,-18.249743798903854,1.9403518303806448,-9.9021785033190817,-9.9021785033190799,-4.1792938042193404,-27.568378864114173,-25.967788741560348,-26.436369970050869,15,270,-999.99000000000001,-1000.4585812284905,1.3500584531745998e-14
-144.70199999999994,84.481929219280801,-27.280492049700488,88.777371053018584,-17.896052513698191,1.9429983076265158,-9.885581980735985,-9.8855819807359815,-4.0752709028232896,-26.882200076364072,-25.994896123104816,-26.465372233745747,15,270,-999.99000000000001,-1000.4604761106409,8.6789471989795746e-15
-144.78399999999993,85.113826478232454,-26.912693863145051,89.267331923487163,-17.546771938356681,1.9460605035900689,-9.8664549989216006,-9.8664549989215988,-3.9723073447167243,-26.203009162288289,-26.022030095738724,-26.494700032710664,15,270,-999.99000000000001,-1000.4626699369719,3.5680116262471476e-14
-144.86599999999993,85.750733256494712,-26.547491598889327,89.766127042552128,-17.201954517252933,1.9495329701740896,-9.8448641197457789,-9.8448641197457771,-3.8704925183283319,-25.531395765540385,-26.04918996688221,-26.524349312976018,15,270,-999.99000000000001,-1000.4651593460939,1.3500584531745998e-14
-144.94799999999992,86.392675937542904,-26.184950798109774,90.273728758419708,-16.861649660274711,1.9534101324000535,-9.8208799480990692,-9.8208799480990656,-3.7699107613187759,-24.867916212784781,-26.076375027834427,-26.554315963687596,15,270,-999.99000000000001,-1000.4679409358532,8.6789471989795746e-15
-145.02999999999992,87.039680500148421,-25.825138882235546,90.790108381115346,-16.525904524931928,1.9576863230642678,-9.7945765573426993,-9.7945765573426975,-3.6706412885687936,-24.213093038675925,-26.103584556497825,-26.58459583734578,15,270,-999.99000000000001,-1000.471011280848,4.4359063461450975e-14
-145.11199999999991,87.691772469369795,-25.46812359350271,91.315235739681754,-16.194763098149235,1.9623557696848082,-9.7660312062780328,-9.7660312062780292,-3.5727579022514417,-23.567413073372972,-26.130817814206239,-26.615184729788457,15,270,-999.99000000000001,-1000.4743669155822,1.3500584531745998e-14
-145.1939999999999,88.348976902165035,-25.113974596642979,91.84907968891207,-15.868267194745568,1.9674126339586444,-9.7353237323959565,-9.7353237323959547,-3.4763290441761385,-22.931327787823054,-26.158074047616623,-26.646078402875631,15,270,-999.99000000000001,-1000.478004355259,1.3500584531745998e-14
-145.2759999999999,89.011318362358125,-24.762762302103553,92.391607808477673,-15.546455777326404,1.9728510059343756,-9.7025362275330593,-9.7025362275330576,-3.3814176234876769,-22.305252157191582,-26.185352487186503,-26.677272571275672,15,270,-999.99000000000001,-1000.4819200840892,1.3500584531745998e-14
-145.35799999999989,89.678820897378785,-24.414558286115589,92.942786561681473,-15.229365246835815,1.9786649239616738,-9.6677525624016507,-9.6677525624016472,-3.2880810376794951,-21.689564799474574,-26.212652347449147,-26.708762909665754,15,270,-999.99000000000001,-1000.4861105622166,1.3500584531745998e-14
-145.43999999999988,90.35150801581635,-24.069435626122331,93.502581419403455,-14.91702967944839,1.9848483923047533,-9.6310579372642966,-9.6310579372642984,-3.1963712236947153,-21.08460831258262,-26.239972827243459,-26.740545058407886,15,270,-999.99000000000001,-1000.4905722311645,1.3500584531745998e-14
-145.52199999999988,91.029402640818731,-23.727468313904868,94.070956717425048,-14.609480520882903,1.9913953844803609,-9.5925385310316269,-9.5925385310316251,-3.106334654806624,-20.490689253763193,-26.267313107935543,-26.772614617163484,15,270,-999.99000000000001,-1000.4953015092279,1.8322221864512419e-14
-145.60399999999987,91.712527113738147,-23.388731668973698,94.647875827571582,-14.306746871553674,1.9982998618352981,-9.5522810747805522,-9.5522810747805522,-3.0180124631862211,-19.908078948107743,-26.294672354971386,-26.804967153118142,15,270,-999.99000000000001,-1000.5002947981468,1.7357894397959133e-14
-145.68599999999986,92.400903140587204,-23.053302116088748,95.233301106555373,-14.008855405065503,2.0055557814140186,-9.5103725065232716,-9.5103725065232716,-2.9314405333560312,-19.337014105012322,-26.322049716254416,-26.837598198648188,15,270,-999.99000000000001,-1000.5055484823938,1.3500584531745998e-14
-145.76799999999986,93.094551783632085,-22.721257659989899,95.827194060177959,-13.715830690511321,2.0131571135265878,-9.4668995887746696,-9.4668995887746679,-2.8466496860061961,-18.777698030705803,-26.349444323200654,-26.870503259065668,15,270,-999.99000000000001,-1000.511058935865,3.9537426128684593e-14
-145.84999999999985,93.793493412423246,-22.392677132590357,96.429515168700902,-13.427694805000469,2.021097841617399,-9.4219486484784696,-9.4219486484784696,-2.7636657659370369,-18.230301208355986,-26.376855288927697,-26.903677805093118,15,270,-999.99000000000001,-1000.5168225161655,8.6789471989795746e-15
-145.93199999999985,94.497747689159311,-22.06764100413185,97.040224122841224,-13.144467855187917,2.029371982828684,-9.3756052191542434,-9.3756052191542434,-2.6825099004352486,-17.694963002427599,-26.404281709419632,-26.937117283746456,15,270,-999.99000000000001,-1000.5228355743268,8.6789471989795746e-15
-146.01399999999984,95.207333529208029,-21.746230796504467,97.659279700378178,-12.866167689966458,2.0379735891646704,-9.3279538155080406,-9.3279538155080388,-2.6031986374411713,-17.171792570091714,-26.431722662336206,-26.970817113253354,15,90,-999.99000000000001,-1000.5290944509171,4.8216373327664331e-15
-146.09599999999983,95.922269069620867,-21.428529249280391,98.286639830915007,-12.592810051931929,2.0468967580509516,-9.2790776758499938,-9.279077675849992,-2.5257441700391254,-16.660870342058349,-26.459177206834973,-27.004772686258683,15,90,-999.99000000000001,-1000.5355954794237,8.6789471989795746e-15
-146.17799999999983,96.642571652200687,-21.114620614042174,98.922261696878309,-12.324408797778322,2.056135643652345,-9.2290585229677564,-9.2290585229677546,-2.4501545818995139,-16.162249641615777,-26.486644384400115,-27.038979375074614,15,270,-999.99000000000001,-1000.5423349906745,1.7357894397959133e-14
diff --git a/examples/simpleyagi/simple_yagi.nec.csv b/examples/simpleyagi/simple_yagi.nec.csv
deleted file mode 100644
index 880c471..0000000
--- a/examples/simpleyagi/simple_yagi.nec.csv
+++ /dev/null
@@ -1,31 +0,0 @@
-mhz,zreal,zimag,zmag,zphase,vswr,s11,s11_real,s11_imag,s11_ang,gain_max,gain_net,gain_max_theta,gain_max_phi,gain_viewer,gain_viewer_net,fb_ratio
-143.80000000000001,54.867747769279759,222.2854601065317,228.95697307582881,76.134577588678425,19.969442333041602,-0.87064634527070595,-0.87064634527070461,3.6386496045919063,24.002062442194433,7.6380313135185594,0.23055297854777468,15,85,-999.99000000000001,-1007.3974783349707,0.29661057228187371
-143.88200000000001,55.240744638014753,223.08734690030573,229.82494254362319,76.092205666937829,19.978461125715828,-0.87025265500854543,-0.87025265500854476,3.6246487314598546,23.909706797194524,7.6361101633855819,0.22685785849344153,15,85,-999.99000000000001,-1007.3992523048921,0.28994135868472981
-143.964,55.614833487577066,223.88916800509492,230.69323625512851,76.049885469567656,19.987561649051148,-0.86985575805235604,-0.86985575805235582,3.6107424054735819,23.817974824914803,7.6342381480375989,0.22319652716347438,15,85,-999.99000000000001,-1007.4010416208741,0.28335242323494442
-144.04599999999999,55.990011105615878,224.69092593843277,231.56185683889578,76.007617491662742,19.99674392532151,-0.86945566269952246,-0.86945566269952301,3.5969296585445001,23.726860139435974,7.6324147647407807,0.21956849660251887,15,85,-999.99000000000001,-1007.4028462681383,0.27684299974076343
-144.12799999999999,56.366274267186803,225.49262329615365,232.43080698509993,75.965402229904413,20.006007990203614,-0.86905237672941882,-0.86905237672941849,3.5832095349927884,23.636356436692182,7.6306395148374548,0.21597328045178532,15,85,-999.99000000000001,-1007.4046662343857,0.27041233088032196
-144.20999999999998,56.743619735754947,226.2942627534735,233.30008944673395,75.923240182404356,20.015353892420148,-0.86864590742119008,-0.86864590742118963,3.569581091338097,23.546457493090408,7.6289119036773743,0.21241039395542316,15,85,-999.99000000000001,-1007.4065015097219,0.26405966806924586
-144.29199999999997,57.122044264162518,227.09584706616252,234.16970704088695,75.881131848566838,20.024781693423442,-0.86823626156974687,-0.86823626156974709,3.5560433960933766,23.457157164150484,7.627231440548969,0.20887935395881829,15,85,-999.99000000000001,-1007.4083520865902,0.25778427132580639
-144.37399999999997,57.501544595665401,227.89737907145181,235.03966264978408,75.839077728922803,20.034291467029025,-0.86782344550391233,-0.86782344550391011,3.5425955295659071,23.368449383192605,7.625597638609956,0.20537967891607511,15,85,-999.99000000000001,-1007.4102179596939,0.25158540913658267
-144.45599999999996,57.882117464971124,228.69886168924461,235.90995922212556,75.797078324986316,20.043883299111982,-0.86740746510180489,-0.867407465101804,3.5292365836581112,23.280328160023391,7.6240100148245267,0.2019108888917005,15,85,-999.99000000000001,-1007.4120991259329,0.24546235833302255
-144.53799999999995,58.26375959926947,229.50029792300896,236.7805997741186,75.755134139098004,20.053557287268386,-0.86698832580771212,-0.86698832580771223,3.5159656616765744,23.192787579676121,7.6224680898953379,0.19847250556387053,15,85,-999.99000000000001,-1007.4139955843315,0.23941440395822297
-144.61999999999995,58.646467719398238,230.30169086103331,237.65158739089745,75.713245674262623,20.063313540506762,-0.86656603264762611,-0.86656603264762566,3.5027818781390625,23.105821801137729,7.620971388200994,0.19506405222725398,15,85,-999.99000000000001,-1007.4159073359738,0.23344083914362823
-144.70199999999994,59.030238540855152,231.10304367735299,238.5229252275906,75.671413434007761,20.073152178948078,-0.86614059024438061,-0.86614059024438061,3.4896843585906909,23.019425056136196,7.6195194377354243,0.19168505379600997,15,85,-999.99000000000001,-1007.4178343839394,0.22754096498868343
-144.78399999999993,59.41506877499824,231.90435963255186,239.39461651031695,75.629637922195016,20.083073333462355,-0.86571200283549243,-0.86571200283549354,3.4766722394239751,22.933591647953484,7.618111770043356,0.18833503681526054,15,85,-999.99000000000001,-1007.4197767332281,0.2217140904337665
-144.86599999999993,59.80095513013206,232.70564207503179,240.26666453761095,75.587919642890341,20.09307714542658,-0.86528027428572185,-0.86528027428572118,3.4637446676970312,22.848315950226375,7.6167479201635242,0.18501352945641969,15,85,-999.99000000000001,-1007.4217343907071,0.2159595321462276
-144.94799999999992,60.187894312783733,233.50689444199404,241.1390726816195,75.546259100177352,20.103163766393493,-0.86484540810345112,-0.8648454081034509,3.450900800957784,22.76359240578681,7.6154274265681794,0.18172006152626619,15,85,-999.99000000000001,-1007.4237073650419,0.21027661439980289
-145.02999999999992,60.575883028742979,234.30812026009761,242.01184438893134,75.5046567980107,20.113333357793401,-0.86440740745559641,-0.86440740745559586,3.438139807077186,22.679415525561783,7.6141498311076399,0.17845416447451523,15,85,-999.99000000000001,-1007.4256956666331,0.20466466896338564
-145.11199999999991,60.964917984394468,235.10932314676188,242.88498318210711,75.463113240044208,20.123586090662162,-0.86396627518135471,-0.86396627518135283,3.4254608640755055,22.59577988742743,7.612914678951471,0.17521537139280241,15,85,-999.99000000000001,-1007.4276993075587,0.19912303498264977
-145.1939999999999,61.354995887924858,235.91050681077925,243.75849266051429,75.421628929451842,20.133922145312198,-0.86352201380836169,-0.86352201380836158,3.4128631599614478,22.512680135147818,7.6117215185345746,0.1720032170297161,15,85,-999.99000000000001,-1007.4297183015049,0.1936510588712366
-145.2759999999999,61.746113450543504,236.71167505350905,244.63237650173261,75.380204368784021,20.144341711101518,-0.86307462556459436,-0.86307462556459447,3.4003458925669565,22.430110977285221,7.6105699015032586,0.16881723778649604,15,85,-999.99000000000001,-1007.4317526637168,0.18824809420132665
-145.35799999999989,62.138267387824861,237.51283176971131,245.50663846264524,75.338840059779315,20.154844986124026,-0.86262411239361825,-0.86262411239362002,3.3879082693888343,22.34806718615539,7.609459382661802,0.16565697172859295,15,85,-999.99000000000001,-1007.4338024109333,0.18291350159487327
-145.43999999999988,62.531454420886682,238.31398094820136,246.38128238032027,75.297536503207354,20.165432176932057,-0.86217047596840435,-0.86217047596840357,3.3755495074365784,22.266543596823809,7.6083895199199612,0.16252195859180052,15,85,-999.99000000000001,-1007.4358675613282,0.17764664861651469
-145.52199999999988,62.925671277789263,239.11512667285444,247.25631217329038,75.25629419868639,20.176103498261028,-0.86171371770501204,-0.86171371770501104,3.363268833077325,22.185535106082888,7.6073598742465727,0.1594117397939705,15,85,-999.99000000000001,-1007.4379481344527,0.17244690967801674
-145.60399999999987,63.320914694882596,239.9162731236261,248.1317318428398,75.215113644518425,20.18685917278334,-0.86125383877497641,-0.86125383877497641,3.3510654818838876,22.105036671449533,7.6063700096116555,0.15632585842946067,15,85,-999.99000000000001,-1007.4400441511822,0.16731366592024777
-145.68599999999986,63.717181418126714,240.717424577272,249.00754547399785,75.173995337516644,20.197699430835439,-0.86079084011880136,-0.86079084011880258,3.3389386984894145,22.025043310206453,7.60541949294395,0.15326385928414066,15,85,-999.99000000000001,-1007.4421556336598,0.16224630512420454
-145.76799999999986,64.114468204419154,241.51858540787296,249.88375723635079,75.132939772824173,20.20862451012496,-0.86032472246025549,-0.86032472246025482,3.32688773644696,21.945550098475795,7.6045078940820048,0.15022528884779618,15,85,-999.99000000000001,-1007.4442826052342,0.15724422161036083
-145.84999999999985,64.512771823007967,242.31976008805182,250.76037138554904,75.091947443758812,20.219634655530097,-0.85985548631668252,-0.85985548631668163,3.3149118580829122,21.866552170252298,7.6036347857275608,0.14720969531077976,15,85,-999.99000000000001,-1007.4464250904168,0.15230681614092564
-145.93199999999985,64.912089057022783,243.12095318943801,251.63739226412233,75.051018841591727,20.230730118763798,-0.85938313201503291,-0.85938313201503513,3.3030103343612187,21.788044716507688,7.6027997433947894,0.14421662858268025,15,85,-999.99000000000001,-1007.4485831148121,0.14743349581552195
-146.01399999999984,65.312416704644619,243.92216938360099,252.51482430266006,75.010154455437885,20.241911158219292,-0.85890765970011551,-0.85890765970011662,3.29118244474691,21.710022984290397,7.6020023453756709,0.1412456402918032,15,85,-999.99000000000001,-1007.4507567050839,0.14262367399541928
-146.09599999999983,65.71375158077943,244.72341344274383,253.39267202089434,74.96935477202544,20.253178038652003,-0.85842906934977636,-0.85842906934977747,3.279427477071704,21.632482275839049,7.6012421726870087,0.13829628379809833,15,85,-999.99000000000001,-1007.452945888889,0.13787677019470851
-146.17799999999983,66.116090518412548,245.52469024024379,254.2709400285581,74.928620275530818,20.264531030936737,-0.85794736078681288,-0.85794736078681388,3.2677447274058031,21.555417947736775,7.6005188090297491,0.13536811420351658,15,85,-999.99000000000001,-1007.4551506948262,0.13319220999374659
diff --git a/xnec2c-gao.cabal b/xnec2c-gao.cabal
index 9b579f7..dcc7cfb 100644
--- a/xnec2c-gao.cabal
+++ b/xnec2c-gao.cabal
@@ -1,6 +1,6 @@
 cabal-version:      2.4
 name:               xnec2c-gao
-version:            0.1.2.0
+version:            0.2.0
 
 -- A short (one-line) description of the package.
 synopsis: g(enetic)a(lgorithm)o(ptimizer) - An external genetic optimizer for xnec2c
@@ -45,6 +45,10 @@ executable xnec2c-gao
                  , Genotype
                  , Phenotype
                  , Genetics
+                 , OptFunc
+                 , Display
+                 , Output
+                 , Paths_xnec2c_gao
                  
 
     -- LANGUAGE extensions used by modules in this package.
@@ -67,9 +71,13 @@ executable xnec2c-gao
         process                     >= 1.6.13 && < 1.7,
         pretty-simple               >= 4.1.2 && < 4.2,
         random                      >= 1.2.1 && < 1.3,
+        these,
+        directory,
+        gitrev,
+        AC-Angle,
         text-builder                >= 0.6.7 && < 0.7,
 
     hs-source-dirs:   app
     default-language: Haskell2010
 
-    GHC-Options: -threaded 
+    GHC-Options: -threaded -Wincomplete-patterns -Wall