#!/usr/bin/env python
# coding: utf-8

# In[1]:


import numpy as np
from fmskill import ModelResult, PointObservation, Connector
import fmskill.metrics as mtr 
get_ipython().run_line_magic('matplotlib', 'inline')


# In[2]:


fn = '../tests/testdata/TS.dfs0'
mr = ModelResult(fn, 'OresundHD', item='Klagshamn: Surface elevation')
print(mr)
print(mr.dfs)


# In[3]:


fn = '../tests/testdata/smhi_2095_klagshamn.dfs0'
o1 = PointObservation(fn, x=366844.15, y=6154291.6, item=0, name='Klagshamn') 
print(o1)
print(o1.n_points, o1.start_time)


# In[4]:


con = Connector(o1, mr, validate=False)
con


# In[5]:


cc = con.extract()
cc


# In[6]:


cc["Klagshamn"].plot_timeseries();


# In[7]:


cmp = cc["Klagshamn"]
cmp.skill()


# In[8]:


mtr.rmse(cmp.obs, cmp.mod, unbiased=False)


# In[ ]: