anchor.diff:      let obs2 = Table.MapToColumn ["hours since 01.01.2015  1:00:00"] "t" (fun t -> t+1.0) obs
   
       let predictorsT = Table.Filter ["t"] (fun t -> t > 40.0 ) ch_data
  -    let observationsT = Table.Filter ["t"] (fun t -> t > 40.0+168.0 ) obs2 //first 168 hours is burnin to enable first pulse reach the Earth
  +    let observationsT = Table.Filter ["t"] (fun t -> t > 600.0 ) obs2
   
       let time_start = 0.0
       let time_step = 0.2
  
positive.code:     spaceAxis.Append(spaceAxisData)
  
      ds.Commit()        
      
      0 // return an integer exit code
  
negative.code: 
  
      let wind = getWind predictorsT S_0 V_k P_w V_b
      
      let a = Table.Map ["t";"velocity"] (fun (t:float) (v:float) -> (t,v)) observationsT |> Array.ofSeq
  
      let current_lglk p =
          let t,v = p
          let avgVelocity = WindAvgAForTime wind t
          let predAtEarth = avgVelocity D
          let mu = V_b + predAtEarth //base level plus pulses        
          if System.Double.IsNaN mu then
              improbable
          else
              let distribution = Statistics.Normal(mu,Sigma)
              let lglk = log_pdf distribution (v * 60.0 * 60.0)
              lglk
              
      let res = Array.Parallel.map current_lglk a |> Seq.sum         
      printf "."
      res
          
  
  [<EntryPoint>]
  let main argv =         
  