koppen_io.fsx

#I @"C:/Users/username/.nuget/packages/fsharp.data/3.3.2/lib/net45"
#r @"FSharp.Data.dll"
#load @"koppen.fs"
open Koppen
open System.IO
open FSharp.Data


type KoppenFile = CsvProvider<"C:/Users/username/Documents/koppen_empirical/100_135.csv",
                              HasHeaders = true,
                              Schema = "float, string, float, float">
let file = KoppenFile.Load("C:/Users/username/Documents/koppen_empirical/100_135.csv")

let imax = 360
let jmax = 180

let readKoppen imax jmax (timer: System.Diagnostics.Stopwatch) name =

    let rec readKoppenCol i imax jmax (timer: System.Diagnostics.Stopwatch) =
        let rec readKoppenRow i j jmax (timer: System.Diagnostics.Stopwatch) =
            let fname = String.concat "" ["C:/Users/username/Documents/koppen_"; name; "/"; string i; "_"; string j; ".csv"]
            if not(File.Exists fname) then
                if j = jmax then [None] else None :: (readKoppenRow i (j+1) jmax timer)
            else
                let file = KoppenFile.Load(fname)

                let daysInFebruary = (file.Rows |> Seq.head).DaysInFebruary * 1.0<day>
                let hemisphere = (file.Rows |> Seq.head).Hemisphere |> function | "Northern" -> Koppen.Northern | "Southern" -> Koppen.Southern | _ -> failwith "fail"
                let temperature =   file.Rows |> Seq.map ((fun x -> x.Temperature)   >> float >> (*) 1.0<C>) |> Seq.toArray
                let precipitation = file.Rows |> Seq.map ((fun x -> x.Precipitation) >> float >> (*) 1.0<mm>) |> Seq.toArray

                let clim = Koppen.Climate (daysInFebruary, hemisphere, temperature, precipitation)
                timer.Start()
                let zone = Koppen.Zones clim
                timer.Stop()
                if j = jmax then [Some zone] else Some zone :: (readKoppenRow i (j+1) jmax timer)
        let row = readKoppenRow i 1 jmax timer
        if i = imax then [row] else row :: (readKoppenCol (i+1) imax jmax timer)

    readKoppenCol 1 imax jmax timer

let writeKoppen name x =
    x
    |> List.map (fun x -> x |> List.map (fun (y: Koppen.Zone Option) -> if y.IsNone then "NA\t" else y |> (Option.get >> string >> fun z -> String.concat "" [z; "\t"])))
    |> List.map (fun x -> x |> List.fold (fun acc elem -> String.concat "" [acc; elem]) "" )
    |> fun x -> File.WriteAllLines (String.concat "" ["C:/Users/username/Documents/koppen_"; name; "/koppen_"; name; ".out"], x)


let koppenTimer = System.Diagnostics.Stopwatch.StartNew()
let grid1 = readKoppen imax jmax koppenTimer "ESM4hist"
let grid2 = readKoppen imax jmax koppenTimer "ESM4hist2"
printfn "%i ms\n" koppenTimer.ElapsedMilliseconds
//writeKoppen grid1
//writeKoppen grid2

let costf = fun x ->
    let sum = List.sumBy float x
    x |> List.map (fun y -> 1000.0 / (float y) / sum)

let paths = Koppen.PathDictionary costf grid1

let asfdfa = Koppen.Distance costf grid1 Koppen.EF Koppen.Aw paths
let asdffa = Koppen.Distance costf grid1 Koppen.Dwc Koppen.Aw paths
let fdsfas = Koppen.Distance costf grid1 Koppen.Cfa Koppen.Aw paths
let fsdssf = Koppen.Distance costf grid1 Koppen.BSh Koppen.Aw paths
let asdfds = Koppen.Distance costf grid1 Koppen.Aw Koppen.Aw paths

let timer = System.Diagnostics.Stopwatch.StartNew()
timer.Start()
let dists = Koppen.DistanceDictionary costf grid1
timer.Stop()
printfn "%i ms\n" timer.ElapsedMilliseconds


let timer = System.Diagnostics.Stopwatch.StartNew()
timer.Start()
let diffs = Koppen.Difference costf grid1 grid2
timer.Stop()
printfn "%i ms\n" timer.ElapsedMilliseconds

let writeDiffs x =
    x
    |> List.map (fun x -> x |> List.map (fun (y: float Option) -> if y.IsNone then "NA," else y |> (Option.get >> string >> fun z -> String.concat "" [z; ","])))
    |> List.map (fun x -> x |> List.fold (fun acc elem -> String.concat "" [acc; elem]) "" )
    |> fun x -> File.WriteAllLines ("C:/Users/username/Documents/koppen_difference.out", x)

writeDiffs diffs