Skip to content
GitLab
Projects Groups Topics Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • S s2dverification
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributor statistics
    • Graph
    • Compare revisions
  • Issues 65
    • Issues 65
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 7
    • Merge requests 7
  • Deployments
    • Deployments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Earth SciencesEarth Sciences
  • s2dverification
  • Issues
  • #257
Closed
Open
Issue created Feb 01, 2021 by Carlos Delgado Torres@cdelgadoMaintainer

Load(): Error in seq.default()... wrong sign in 'by' argument

Hi @nperez and @aho,

I'm trying to use Load() and I'm getting the following error:

Error in seq.default(day, min(days_in_month, (length(leadtimes) - jleadtime) *  : 
 wrong sign in 'by' argument

This error happens when trying to load data, for example, for February (leadtimemin = 32, leadtimemax = 59):

data_Feb <- Load(var = 'psl', 
                 exp = list(list(path = '/esarchive/exp/ecearth/a1st/cmorfiles/CMIP/EC-Earth-Consortium/EC-Earth3/historical/r7i1p1f1/day/psl/gr/v20200120/psl_*$YEAR$$MONTH$*.nc',nc_var_name = 'psl')),
                 obs = list(list(path = '/esarchive/recon/jma/jra55/daily_mean/psl_f6h/psl_*$YEAR$$MONTH$*.nc', nc_var_name = 'psl')),
                 sdates = paste0(1960:2010,'0101'), 
                 leadtimemin = 32, leadtimemax = 59,
                 latmin = 27, latmax = 81, 
                 lonmin = 274.5, lonmax = 45,
                 storefreq = 'daily', output = 'lonlat')

The same code but for January (leadtimemin = 1, leadtimemax = 31) works:

data_Jan <- Load(var = 'psl', 
                 exp = list(list(path = '/esarchive/exp/ecearth/a1st/cmorfiles/CMIP/EC-Earth-Consortium/EC-Earth3/historical/r7i1p1f1/day/psl/gr/v20200120/psl_*$YEAR$$MONTH$*.nc',nc_var_name = 'psl')),
                 obs = list(list(path = '/esarchive/recon/jma/jra55/daily_mean/psl_f6h/psl_*$YEAR$$MONTH$*.nc', nc_var_name = 'psl')),
                 sdates = paste0(1960:2010,'0101'), 
                 leadtimemin = 1, leadtimemax = 31,
                 latmin = 27, latmax = 81, 
                 lonmin = 274.5, lonmax = 45,
                 storefreq = 'daily', output = 'lonlat')

Also, it works for both January and February with an old version of the package (the one in develop-Regimes). Please tell me if you need more information.

Thank you very much,

Carlos

Assignee
Assign to
Time tracking