def read_primary_dataset(variq, dataset, lat_bounds=lat_bounds, lon_bounds=lon_bounds): data, lats, lons = df.readFiles(variq, dataset, monthlychoice) datar, lats, lons = df.getRegion(data, lats, lons, lat_bounds, lon_bounds) print('\nOur dataset: ', dataset, ' is shaped', data.shape) return datar, lats, lons
def read_obs_dataset(variq,dataset_obs,lat_bounds=lat_bounds,lon_bounds=lon_bounds): data_obs,lats_obs,lons_obs = df.readFiles(variq,dataset_obs,monthlychoice) data_obs,lats_obs,lons_obs = df.getRegion(data_obs,lats_obs,lons_obs, lat_bounds,lon_bounds) if dataset_obs == '20CRv3': if monthlychoice == 'DJF': year20cr = np.arange(1837,2015+1) else: year20cr = np.arange(1836,2015+1) year_obsall = np.arange(yearsall[sis].min(),yearsall[sis].max()+1,1) yearqq = np.where((year20cr >= year_obsall.min()) & (year20cr <= year_obsall.max()))[0] data_obs = data_obs[yearqq,:,:] print('our OBS dataset: ',dataset_obs,' is shaped',data_obs.shape) return data_obs,lats_obs,lons_obs