scripts/CD_WP3/look_at_data.R

In NorskRegnesentral/SpatGEVBMA: Hierarchical spatial generalized extreme value (GEV) modeling with Bayesian Model Averaging (BMA)

rcp45=nc_open("/nr/project/stat/ClimDesign/WP3/Data/fromOskar/inference/rcp45/1991/MAP.nc")
tas1=ncvar_get(rcp45,"MAP")

image.plot(tas1)


rcp45=nc_open("/nr/project/stat/ClimDesign/WP3/Data/fromOskar/prediction/rcp45/2071/MAP.nc")
tas2=ncvar_get(rcp45,"MAP")

diff=tas2-tas1

diff=diff[is.na(diff)==FALSE]

hist(diff[1:10000])


senorge=nc_open("/nr/project/stat/ClimDesign/WP3/Data/fromOskar/inference/senorge/1991/MSP.nc")
var=ncvar_get(senorge,"MSP")
image.plot(var)

rcp=nc_open("/nr/project/stat/ClimDesign/WP3/Data/fromOskar/inference/rcp45/1991/MSP.nc")
var=ncvar_get(rcp,"MSP")
image.plot(var)

senorge=nc_open("/nr/project/stat/ClimDesign/WP3/Data/fromOskar/inference/senorge/1991/JJAtemp.nc")
var=ncvar_get(senorge,"tas")
image.plot(var)

rcp=nc_open("/nr/project/stat/ClimDesign/WP3/Data/fromOskar/inference/rcp45/1991/JJAtemp.nc")
var=ncvar_get(rcp,"tas")
image.plot(var)