Skip to content

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')
Edited by aho