CESM2 Noleap Calendar

[Z-Richard]

With the latest version I will encounter an indexing error with cftime.NoLeap using CESM2 output. The error can be traced back to L155 at read data.py, e.g.,

DATA['data'] = DS.sel({variable_names[2]:str(dt_to_use)},method='nearest')[variable_names[3]].values

which can be fixed by

DATA['data'] = DS.sel({variable_names[2]:dt_to_use},method='nearest')[variable_names[3]].values

[Z-Richard]

When running mask output on CESM2 data, I also need to change the code for function get_lpo_grid_points to:

def get_lpo_grid_points(lp_object_id, mask_times, lp_objects_dir, lp_objects_fn_format):

    nnnn = int(str(int(lp_object_id))[-4:])

    dt_this0 = lpt.helpers.get_objid_datetime(lp_object_id)
    dt_this = cftime.datetime(dt_this0.year,dt_this0.month,dt_this0.day,dt_this0.hour,calendar=mask_times[0].calendar) #,calendar=mask_times[0].calendar)

I find that I have to add the line calendar=mask_times[0].calendar back in.