import numpy as np
import xarray as xr
import sys
import matplotlib.pyplot as plt

var="Evap"


g="/bdd/MEDI/workspaces/polcher/LIAISEoutputs/GFDL/WFDE5_CRU_GPCC/TS_DA/GFDL_WFDE5_CRU_GPCC_19890101_20131231_1D_"+var+".nc"
c="/bdd/MEDI/workspaces/polcher/LIAISEoutputs/CNRM/WFDE5_CRU_GPCC/TS_DA/CNRM_WFDE5_CRU_GPCC_19890101_20131231_1D_"+var+".nc"

ModelVarNames={"ORCH":{"time_counter":"time"}, \
               "JULES":{"time":"time","longitude":"lon","latitude":"lat"}, \
               "SURFEX":{"time_counter":"time"}, \
               "ecLand":{"time":"time"}, \
               "LM4":{"time":"time"}}

gfdl=xr.open_dataset(g).rename(name_dict=ModelVarNames["LM4"]).mean(dim="time")
cnrm=xr.open_dataset(c).rename(name_dict=ModelVarNames["SURFEX"]).mean(dim="time")

(gfdl[var]*86400-cnrm[var]*86400).plot()
plt.show()