reading netcdf files (Julia)

read_net.jl

using NCDataFrame
using NCDatasets
using PlotlyJS
using Pkg
using Statistics
using NaNStatistics
import NaNStatistics ; na=NaNStatistics


ds = NCDataset("ECMWF_ERA-40_subset.nc","r")

@show keys(ds)


temp = collect(ds["p2t"])

n_temp = coalesce.(temp, NaN)

#nanmean(x) = mean(filter(!isnan,x),2)
A = [1 2; 3 4]
#ff = nanmean(collect(temp))
dd = na.nanmean(n_temp,dims=(2,3))[:]
times = ds["time"][:]

p = scatter(y=dd,x=times)

layout = Layout(

    title="Temperature Time Series",

    xaxis_title="Date",

    yaxis_title="Temperature(K)"

    )

plot(p,layout)