Bug when loading two datasets on same grid but with lons shifted.
Hi,
This script crashes, there's a bug.
#!/usr/bin/env Rscript
setwd("~/Desktop/Experiments/b02p/bias")
library(s2dverification)
# This script plots the bias of the 'b02p' sea ice concentration compared with the NSIDC data
#900 es 75 años multiplicados por 12 meses
# NOTE -> When the dataset is changed and corrected, 'ice' won't be loaded anymore, but 'sic' (and the path will change)
b02p <- list(path = '/home/Earth/rcruzgar/Desktop/Experiments/b02p/first_interpolation/ice_b02p_$START_DATE$_fc0_*_regular_drown.nc')
nsidc <- list(path = '/esnas/scratch/pabretonniere/data_cfunas_cleaning/ice_nsidc/ice_nsidc_$YEAR$$MONTH$.nc', nc_var_name = 'ice')
mask <- list(path = '/esnas/autosubmit/con_files/mask_hadisst.nc')
# March 1979 - March 2013
if (file.exists('sic_mar.sav')){
load('sic_mar.sav')
} else {
sic_mar <- Load('sic', NULL, list(nsidc), '19501101', output = 'lonlat', leadtimemin = 341, leadtimemax = 749,
sampleperiod = 12, latmin = 47, grid = 'r360x180')
save(sic_mar, file = 'sic_mar.sav')