def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
# for CMIP5-ALL_med_9P_DA_fix_best.nc
ds['weights_ic_all'] = [
0.0799212297865281, 0.0991501104631161, 0.00355090416936244,
0.0175800526997022, 0.0175800526997022, 0.0175800526997022,
0.0175800526997022, 0.0175800526997022, 0.0175800526997022,
0.137847841006054, 0.0558673816564371, 0.0558673816564371,
0.0558673816564371, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.00140336196773973, 0.00140336196773973, 0.00140336196773973,
0.0767836545231081, 0.0848404634182867, 0.0922159550282929,
0.0252370776505070, 0.0252370776505070, 0.0252370776505070,
0.0252370776505070, 0.0252370776505070, 0.00744998626725123,
0.00744998626725123, 0.00744998626725123, 0.00744998626725123,
0.00744998626725123, 0.00744998626725123, 0.00744998626725123,
0.00744998626725123, 0.00744998626725123, 0.00744998626725123,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.000926992696255780,
0.000926992696255780, 0.000926992696255780, 0.00820181588058064,
0.00820181588058064, 0.00820181588058064, 0.00820181588058064,
0.00820181588058064, 0.0236051118821537, 0.0236051118821537,
0.0236051118821537, 0.0236051118821537, 0.0236051118821537,
0.0736937402123117, 0.00862980763323842, 0.00862980763323842,
0.00862980763323842, 0.0864302907878062, 0.0168519269380017,
0.00381643513449763, 0.0361237757994686, 0.0190798185657880,
0.0372471601611859, 0.0382458592020622, 0.0190798185657880,
0.0382458592020622, 0.0769799195984690, 0.0320453923875567,
0.0688547758192422, 0.0417020799052823, 0.0320453923875567,
0.0417020799052823, 0.111611631630745, 0.117738785961798,
0.0352548041590223, 0.0352548041590223, 0.0352548041590223,
0.0352548041590223, 0.00449714417985167, 0.00449714417985167,
0.00449714417985167, 0.00449714417985167, 0.113935386563146,
0.00455420633885610, 0.0379465364707411, 0.0470947210136297,
0.0251055960430957, 0.0251055960430957, 0.0251055960430957,
0.0405837017395739, 0.0405837017395739, 0.0405837017395739,
0.135838042630498, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.000836738614324269,
0.000836738614324269, 0.000836738614324269, 0.0362092735672884,
0.0414815108681919, 0.163988115194310, 0.132484175147473,
0.0630194025670830, 0.0996666304662157, 0.0342111705909201
]
ds['weights_ic_all'] = ds['weights_ic_all'] / ds['weights_ic_all'].sum()
ds['weights_mc_all'] = [
0.0447678350624111, 0.0447678350624111, 0.00355090416936244,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.00132831697012572, 0.00132831697012572, 0.00132831697012572,
0.0282044525521875, 0.0282044525521875, 0.0282044525521875,
0.0252370776505070, 0.0252370776505070, 0.0252370776505070,
0.0252370776505070, 0.0252370776505070, 0.00744998626725123,
0.00744998626725123, 0.00744998626725123, 0.00744998626725123,
0.00744998626725123, 0.00744998626725123, 0.00744998626725123,
0.00744998626725123, 0.00744998626725123, 0.00744998626725123,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.000833893268645939, 0.000833893268645939,
0.000833893268645939, 0.0236051118821537, 0.0236051118821537,
0.0236051118821537, 0.0236051118821537, 0.0236051118821537,
0.0736937402123117, 0.00862980763323842, 0.00862980763323842,
0.00862980763323842, 0.0118998503178117, 0.0118998503178117,
0.0118998503178117, 0.00516318216403558, 0.00516318216403558,
0.00516318216403558, 0.00516318216403558, 0.00516318216403558,
0.00516318216403558, 0.00516318216403558, 0.00516318216403558,
0.00516318216403558, 0.00516318216403558, 0.00516318216403558,
0.00516318216403558, 0.0220396467815806, 0.0220396467815806,
0.0220396467815806, 0.0220396467815806, 0.0220396467815806,
0.0220396467815806, 0.00529010832721190, 0.00529010832721190,
0.00529010832721190, 0.00529010832721190, 0.00529010832721190,
0.00529010832721190, 0.0124396648748893, 0.0124396648748893,
0.0124396648748893, 0.0124396648748893, 0.0124396648748893,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.000819940597404711,
0.000819940597404711, 0.000819940597404711, 0.0194226961088701,
0.0194226961088701, 0.0741180725854457, 0.0741180725854457,
0.0406715082583247, 0.0406715082583247, 0.0342111705909201
]
ds['weights_mc_all'] = ds['weights_mc_all'] / ds['weights_mc_all'].sum()
ds['weights_q_norm_all'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none_all'] = np.ones(288) / 288
ds['weights_ic'] = [
0.0814281063120076, 0.101360732375565, 0.00387158068646408,
0.0182109998786213, 0.0182109998786213, 0.0182109998786213,
0.0182109998786213, 0.0182109998786213, 0.0182109998786213,
0.140863919084632, 0.0574894076825984, 0.0574894076825984,
0.0574894076825984, 0.0782840853808204, 0.0873137993435310,
0.0953285417013194, 0.0257098534242143, 0.0257098534242143,
0.0257098534242143, 0.0257098534242143, 0.0257098534242143,
0.00765358393092044, 0.00765358393092044, 0.00765358393092044,
0.00765358393092044, 0.00765358393092044, 0.00765358393092044,
0.00765358393092044, 0.00765358393092044, 0.00765358393092044,
0.00765358393092044, 0.00833272338258814, 0.00833272338258814,
0.00833272338258814, 0.00833272338258814, 0.00833272338258814,
0.0228193533713553, 0.0228193533713553, 0.0228193533713553,
0.0228193533713553, 0.0228193533713553, 0.0228193533713553,
0.00880436359436853, 0.00880436359436853, 0.00880436359436853,
0.0897732071535642, 0.0182577607622864, 0.00443527060764555,
0.0378926428892257, 0.0198821618837288, 0.0386367304696207,
0.0395013119901435, 0.0198821618837288, 0.0395013119901435,
0.0788228878603702, 0.0332085558215958, 0.0715696432500609,
0.0434600839900184, 0.0332085558215958, 0.0434600839900184,
0.113046636510380, 0.119124118976484, 0.0358403976669652,
0.0358403976669652, 0.0358403976669652, 0.0358403976669652,
0.00455147377668773, 0.00455147377668773, 0.00455147377668773,
0.00455147377668773, 0.114734727530186, 0.00467592166118176,
0.0390015092086924, 0.0484132590434984, 0.0255324254078565,
0.0255324254078565, 0.0255324254078565, 0.0417708387761604,
0.0417708387761604, 0.0417708387761604, 0.137843157908220,
0.0368302891575653, 0.0423791381811358, 0.166587246898677,
0.136179964116823, 0.0658653676434850, 0.104283047403064,
0.0349836708545518
]
ds['weights_ic'] = ds['weights_ic'] / ds['weights_ic'].sum()
ds['weights_mc'] = [
0.0456972096718932, 0.0456972096718932, 0.00387158068646408,
0.0131386458385839, 0.0131386458385839, 0.0131386458385839,
0.0131386458385839, 0.0131386458385839, 0.0131386458385839,
0.0131386458385839, 0.0131386458385839, 0.0131386458385839,
0.0131386458385839, 0.0289918251584079, 0.0289918251584079,
0.0289918251584079, 0.0257098534242143, 0.0257098534242143,
0.0257098534242143, 0.0257098534242143, 0.0257098534242143,
0.00765358393092044, 0.00765358393092044, 0.00765358393092044,
0.00765358393092044, 0.00765358393092044, 0.00765358393092044,
0.00765358393092044, 0.00765358393092044, 0.00765358393092044,
0.00765358393092044, 0.00833272338258814, 0.00833272338258814,
0.00833272338258814, 0.00833272338258814, 0.00833272338258814,
0.0228193533713553, 0.0228193533713553, 0.0228193533713553,
0.0228193533713553, 0.0228193533713553, 0.0228193533713553,
0.00880436359436853, 0.00880436359436853, 0.00880436359436853,
0.0124962487248329, 0.0124962487248329, 0.0124962487248329,
0.00535507193896683, 0.00535507193896683, 0.00535507193896683,
0.00535507193896683, 0.00535507193896683, 0.00535507193896683,
0.00535507193896683, 0.00535507193896683, 0.00535507193896683,
0.00535507193896683, 0.00535507193896683, 0.00535507193896683,
0.0223782532821752, 0.0223782532821752, 0.0223782532821752,
0.0223782532821752, 0.0223782532821752, 0.0223782532821752,
0.00533983971162143, 0.00533983971162143, 0.00533983971162143,
0.00533983971162143, 0.00533983971162143, 0.00533983971162143,
0.0126882638769160, 0.0126882638769160, 0.0126882638769160,
0.0126882638769160, 0.0126882638769160, 0.0321112941808540,
0.0321112941808540, 0.0321112941808540, 0.0321112941808540,
0.0198023568346753, 0.0198023568346753, 0.0756918027538749,
0.0756918027538749, 0.0425371037616373, 0.0425371037616373,
0.0349836708545518
]
ds['weights_mc'] = ds['weights_mc'] / ds['weights_mc'].sum()
ds['weights_q_norm'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none'] = np.ones(88) / 88
return ds
def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
# for CMIP5-ALL_neu_9P_DA_fix.nc
ds['weights_mc_all'] = [
0.0369529903296525, 0.0369529903296525, 0.00789557571531341,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.000223622172762458, 0.000223622172762458, 0.000223622172762458,
0.00757568838398542, 0.00757568838398542, 0.00757568838398542,
0.0207288754670184, 0.0207288754670184, 0.0207288754670184,
0.0207288754670184, 0.0207288754670184, 3.14802067638986e-05,
3.14802067638986e-05, 3.14802067638986e-05, 3.14802067638986e-05,
3.14802067638986e-05, 3.14802067638986e-05, 3.14802067638986e-05,
3.14802067638986e-05, 3.14802067638986e-05, 3.14802067638986e-05,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.00100158618894435, 0.00100158618894435,
0.00100158618894435, 0.0370404668427630, 0.0370404668427630,
0.0370404668427630, 0.0370404668427630, 0.0370404668427630,
0.00186400702563852, 0.00527487256310576, 0.00527487256310576,
0.00527487256310576, 0.0220506854381394, 0.0220506854381394,
0.0220506854381394, 0.00129499009590567, 0.00129499009590567,
0.00129499009590567, 0.00129499009590567, 0.00129499009590567,
0.00129499009590567, 0.00129499009590567, 0.00129499009590567,
0.00129499009590567, 0.00129499009590567, 0.00129499009590567,
0.00129499009590567, 0.0115142107716664, 0.0115142107716664,
0.0115142107716664, 0.0115142107716664, 0.0115142107716664,
0.0115142107716664, 0.000152682252506608, 0.000152682252506608,
0.000152682252506608, 0.000152682252506608, 0.000152682252506608,
0.000152682252506608, 0.00944801694165041, 0.00944801694165041,
0.00944801694165041, 0.00944801694165041, 0.00944801694165041,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00135847908868839,
0.00135847908868839, 0.00135847908868839, 0.00400840036469699,
0.00400840036469699, 0.0155524889846518, 0.0155524889846518,
0.00335113762097806, 0.00335113762097806, 0.0174320531748485
]
ds['weights_mc_all'] = ds['weights_mc_all'] / ds['weights_mc_all'].sum()
ds['weights_ic_all'] = [
0.0782604641068231, 0.0695514972117871, 0.00789557571531341,
0.00201692317667272, 0.00201692317667272, 0.00201692317667272,
0.00201692317667272, 0.00201692317667272, 0.00201692317667272,
0.0101788053554796, 0.0152051205905696, 0.0152051205905696,
0.0152051205905696, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.000234162278765610, 0.000234162278765610, 0.000234162278765610,
0.00111052286441345, 0.0255645854618377, 0.0415060871296176,
0.0207288754670184, 0.0207288754670184, 0.0207288754670184,
0.0207288754670184, 0.0207288754670184, 3.14802067638986e-05,
3.14802067638986e-05, 3.14802067638986e-05, 3.14802067638986e-05,
3.14802067638986e-05, 3.14802067638986e-05, 3.14802067638986e-05,
3.14802067638986e-05, 3.14802067638986e-05, 3.14802067638986e-05,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00112260979563779,
0.00112260979563779, 0.00112260979563779, 0.00893094929848742,
0.00893094929848742, 0.00893094929848742, 0.00893094929848742,
0.00893094929848742, 0.0370404668427630, 0.0370404668427630,
0.0370404668427630, 0.0370404668427630, 0.0370404668427630,
0.00186400702563852, 0.00527487256310576, 0.00527487256310576,
0.00527487256310576, 0.144972913747555, 0.00780658495350885,
0.0456766702421899, 0.0113679426595628, 0.0146711522831160,
0.0176025662035585, 0.0100004382099991, 0.0146711522831160,
0.0100004382099991, 0.0141452724257955, 0.00497077094070862,
0.00490831069191025, 0.00497125902357383, 0.00497077094070862,
0.00497125902357383, 0.0526374597401102, 0.0130786797378713,
0.0217997155188756, 0.0217997155188756, 0.0217997155188756,
0.0217997155188756, 0.000130243631551600, 0.000130243631551600,
0.000130243631551600, 0.000130243631551600, 0.00341203435551315,
6.28629899129800e-07, 0.0227676785264735, 0.0297211268972951,
0.0204124020130546, 0.0204124020130546, 0.0204124020130546,
0.0554388846168218, 0.0554388846168218, 0.0554388846168218,
0.122616198420997, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00140717436632812,
0.00140717436632812, 0.00140717436632812, 0.00501284018497503,
0.0110207612738129, 0.0404851833511034, 0.0217247725875036,
0.00384832785928823, 0.00955622262462400, 0.0174320531748485
]
ds['weights_ic_all'] = ds['weights_ic_all'] / ds['weights_ic_all'].sum()
ds['weights_q_norm_all'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none_all'] = np.ones(288) / 288
# for CMIP5_neu_9P_DA_fix.nc
ds['weights_ic'] = [
0.0659090254416251, 0.0634127600564207, 0.00697737986343313,
0.00176778863539344, 0.00176778863539344, 0.00176778863539344,
0.00176778863539344, 0.00176778863539344, 0.00176778863539344,
0.00772898142491627, 0.0129845675654988, 0.0129845675654988,
0.0129845675654988, 0.000630392182257593, 0.0202620985586898,
0.0376089111046764, 0.0189469854054094, 0.0189469854054094,
0.0189469854054094, 0.0189469854054094, 0.0189469854054094,
2.28316788224882e-05, 2.28316788224882e-05, 2.28316788224882e-05,
2.28316788224882e-05, 2.28316788224882e-05, 2.28316788224882e-05,
2.28316788224882e-05, 2.28316788224882e-05, 2.28316788224882e-05,
2.28316788224882e-05, 0.00767861786838925, 0.00767861786838925,
0.00767861786838925, 0.00767861786838925, 0.00767861786838925,
0.0287720089061949, 0.0287720089061949, 0.0287720089061949,
0.0287720089061949, 0.0287720089061949, 0.0287720089061949,
0.00492280254659331, 0.00492280254659331, 0.00492280254659331,
0.130476413455570, 0.00534689955942724, 0.0434719757044674,
0.00934484197869562, 0.0136279210138440, 0.0153936563993022,
0.00930002838536717, 0.0136279210138440, 0.00930002838536717,
0.0123354884543809, 0.00389155184104562, 0.00377624535941361,
0.00414469645241206, 0.00389155184104562, 0.00414469645241206,
0.0474148394981459, 0.00973578470940712, 0.0198334791230916,
0.0198334791230916, 0.0198334791230916, 0.0198334791230916,
0.000103765234539357, 0.000103765234539357, 0.000103765234539357,
0.000103765234539357, 0.00302127295485722, 3.48734608226817e-07,
0.0220917117263817, 0.0289694064978044, 0.0188376185201033,
0.0188376185201033, 0.0188376185201033, 0.0507066467167848,
0.0507066467167848, 0.0507066467167848, 0.113252355767191,
0.00428106937538889, 0.0102313598944472, 0.0380835384079005,
0.0180762746181618, 0.00300924413926564, 0.00735635201056383,
0.0157031030646123
]
ds['weights_ic'] = ds['weights_ic'] / ds['weights_ic'].sum()
ds['weights_mc'] = [
0.0323304463745115, 0.0323304463745115, 0.00697737986343313,
0.00188230480388570, 0.00188230480388570, 0.00188230480388570,
0.00188230480388570, 0.00188230480388570, 0.00188230480388570,
0.00188230480388570, 0.00188230480388570, 0.00188230480388570,
0.00188230480388570, 0.00650015576062486, 0.00650015576062486,
0.00650015576062486, 0.0189469854054094, 0.0189469854054094,
0.0189469854054094, 0.0189469854054094, 0.0189469854054094,
2.28316788224882e-05, 2.28316788224882e-05, 2.28316788224882e-05,
2.28316788224882e-05, 2.28316788224882e-05, 2.28316788224882e-05,
2.28316788224882e-05, 2.28316788224882e-05, 2.28316788224882e-05,
2.28316788224882e-05, 0.00767861786838925, 0.00767861786838925,
0.00767861786838925, 0.00767861786838925, 0.00767861786838925,
0.0287720089061949, 0.0287720089061949, 0.0287720089061949,
0.0287720089061949, 0.0287720089061949, 0.0287720089061949,
0.00492280254659331, 0.00492280254659331, 0.00492280254659331,
0.0199216987466071, 0.0199216987466071, 0.0199216987466071,
0.00114379877057269, 0.00114379877057269, 0.00114379877057269,
0.00114379877057269, 0.00114379877057269, 0.00114379877057269,
0.00114379877057269, 0.00114379877057269, 0.00114379877057269,
0.00114379877057269, 0.00114379877057269, 0.00114379877057269,
0.0104023969493616, 0.0104023969493616, 0.0104023969493616,
0.0104023969493616, 0.0104023969493616, 0.0104023969493616,
0.000130051817835977, 0.000130051817835977, 0.000130051817835977,
0.000130051817835977, 0.000130051817835977, 0.000130051817835977,
0.00882398739620463, 0.00882398739620463, 0.00882398739620463,
0.00882398739620463, 0.00882398739620463, 0.0356007610136409,
0.0356007610136409, 0.0356007610136409, 0.0356007610136409,
0.00362810731745902, 0.00362810731745902, 0.0140399532565156,
0.0140399532565156, 0.00259139903745737, 0.00259139903745737,
0.0157031030646123
]
ds['weights_mc'] = ds['weights_mc'] / ds['weights_mc'].sum()
ds['weights_q_norm'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none'] = np.ones(88) / 88
return ds
def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
# for CMIP5-ALL_neu_9P_DA_fix_best.nc
ds['weights_mc_all'] = [
0.0679316202498926, 0.0679316202498926, 0.0217759358403415,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.000547510478316325, 0.000547510478316325, 0.000547510478316325,
0.0179448964438854, 0.0179448964438854, 0.0179448964438854,
0.0448157293910969, 0.0448157293910969, 0.0448157293910969,
0.0448157293910969, 0.0448157293910969, 0.000111141946048448,
0.000111141946048448, 0.000111141946048448, 0.000111141946048448,
0.000111141946048448, 0.000111141946048448, 0.000111141946048448,
0.000111141946048448, 0.000111141946048448, 0.000111141946048448,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.00175468998781131, 0.00175468998781131,
0.00175468998781131, 0.0416763531348045, 0.0416763531348045,
0.0416763531348045, 0.0416763531348045, 0.0416763531348045,
0.00514634292965141, 0.00997724508583466, 0.00997724508583466,
0.00997724508583466, 0.0342477908930083, 0.0342477908930083,
0.0342477908930083, 0.00346010607564650, 0.00346010607564650,
0.00346010607564650, 0.00346010607564650, 0.00346010607564650,
0.00346010607564650, 0.00346010607564650, 0.00346010607564650,
0.00346010607564650, 0.00346010607564650, 0.00346010607564650,
0.00346010607564650, 0.0239444169793284, 0.0239444169793284,
0.0239444169793284, 0.0239444169793284, 0.0239444169793284,
0.0239444169793284, 0.000530077196174071, 0.000530077196174071,
0.000530077196174071, 0.000530077196174071, 0.000530077196174071,
0.000530077196174071, 0.0173161039061022, 0.0173161039061022,
0.0173161039061022, 0.0173161039061022, 0.0173161039061022,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.00216707468632810,
0.00216707468632810, 0.00216707468632810, 0.0121434844643639,
0.0121434844643639, 0.0259996087793537, 0.0259996087793537,
0.0103515135399188, 0.0103515135399188, 0.0403290698971819
]
ds['weights_mc_all'] = ds['weights_mc_all'] / ds['weights_mc_all'].sum()
ds['weights_ic_all'] = [
0.170035690706966, 0.101690790292604, 0.0217759358403415,
0.00511133912710308, 0.00511133912710308, 0.00511133912710308,
0.00511133912710308, 0.00511133912710308, 0.00511133912710308,
0.0215375990716382, 0.0314063504684956, 0.0314063504684956,
0.0314063504684956, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.000593133904029985, 0.000593133904029985, 0.000593133904029985,
0.00344822058993842, 0.0645264080279494, 0.0935294393770811,
0.0448157293910969, 0.0448157293910969, 0.0448157293910969,
0.0448157293910969, 0.0448157293910969, 0.000111141946048448,
0.000111141946048448, 0.000111141946048448, 0.000111141946048448,
0.000111141946048448, 0.000111141946048448, 0.000111141946048448,
0.000111141946048448, 0.000111141946048448, 0.000111141946048448,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.00196080115042499,
0.00196080115042499, 0.00196080115042499, 0.0162373734826694,
0.0162373734826694, 0.0162373734826694, 0.0162373734826694,
0.0162373734826694, 0.0416763531348045, 0.0416763531348045,
0.0416763531348045, 0.0416763531348045, 0.0416763531348045,
0.00514634292965141, 0.00997724508583466, 0.00997724508583466,
0.00997724508583466, 0.188024263836110, 0.0223931979047178,
0.0978126562962474, 0.0328875020233990, 0.0368479903205261,
0.0478192925518620, 0.0260613543818759, 0.0368479903205261,
0.0260613543818759, 0.0379884041624291, 0.0143260627413111,
0.0149785563204718, 0.0139099725150202, 0.0143260627413111,
0.0139099725150202, 0.114792401935603, 0.0344509934976026,
0.0445472259889136, 0.0445472259889136, 0.0445472259889136,
0.0445472259889136, 0.000460478719524570, 0.000460478719524570,
0.000460478719524570, 0.000460478719524570, 0.0117123843113200,
2.73523855340493e-06, 0.0308679260616577, 0.0315811562188857,
0.0411615017080012, 0.0411615017080012, 0.0411615017080012,
0.0979945760658846, 0.0979945760658846, 0.0979945760658846,
0.230880111985416, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.00223262485107463,
0.00223262485107463, 0.00223262485107463, 0.0161064643939164,
0.0324674734635393, 0.0683569849667027, 0.0356414501507119,
0.0126647547696419, 0.0287412993900331, 0.0403290698971819
]
ds['weights_ic_all'] = ds['weights_ic_all'] / ds['weights_ic_all'].sum()
ds['weights_q_norm_all'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none_all'] = np.ones(288) / 288
# for CMIP5_neu_9P_DA_fix.nc
ds['weights_ic'] = [
0.147662908933192, 0.0928492811927518, 0.0195113343253419,
0.00455272244048226, 0.00455272244048226, 0.00455272244048226,
0.00455272244048226, 0.00455272244048226, 0.00455272244048226,
0.0168251587855511, 0.0274636695579662, 0.0274636695579662,
0.0274636695579662, 0.00205689897624903, 0.0527473806482939,
0.0852710545719979, 0.0415629168130336, 0.0415629168130336,
0.0415629168130336, 0.0415629168130336, 0.0415629168130336,
8.23198683529589e-05, 8.23198683529589e-05, 8.23198683529589e-05,
8.23198683529589e-05, 8.23198683529589e-05, 8.23198683529589e-05,
8.23198683529589e-05, 8.23198683529589e-05, 8.23198683529589e-05,
8.23198683529589e-05, 0.0142040710077886, 0.0142040710077886,
0.0142040710077886, 0.0142040710077886, 0.0142040710077886,
0.0325206196276115, 0.0325206196276115, 0.0325206196276115,
0.0325206196276115, 0.0325206196276115, 0.0325206196276115,
0.00933229305832333, 0.00933229305832333, 0.00933229305832333,
0.170960772905746, 0.0160334920388853, 0.0931546244882425,
0.0276970305249978, 0.0346009510959196, 0.0426135224765772,
0.0243500738898684, 0.0346009510959196, 0.0243500738898684,
0.0335297237504207, 0.0115629720798375, 0.0118243767337137,
0.0118098583213309, 0.0115629720798375, 0.0118098583213309,
0.104870579819529, 0.0264929760235278, 0.0409185449797546,
0.0409185449797546, 0.0409185449797546, 0.0409185449797546,
0.000373695855254702, 0.000373695855254702, 0.000373695855254702,
0.000373695855254702, 0.0105199157010856, 1.55800727166352e-06,
0.0301196287219131, 0.0308271880937737, 0.0381342234080748,
0.0381342234080748, 0.0381342234080748, 0.0910160512202272,
0.0910160512202272, 0.0910160512202272, 0.215686884225117,
0.0139126169941522, 0.0301059279560559, 0.0648262624725157,
0.0301867739551750, 0.0102306570072921, 0.0229493728221231,
0.0364483156107965
]
ds['weights_ic'] = ds['weights_ic'] / ds['weights_ic'].sum()
ds['weights_mc'] = [
0.0601280475314860, 0.0601280475314860, 0.0195113343253419,
0.00427896192664608, 0.00427896192664608, 0.00427896192664608,
0.00427896192664608, 0.00427896192664608, 0.00427896192664608,
0.00427896192664608, 0.00427896192664608, 0.00427896192664608,
0.00427896192664608, 0.0155639260218379, 0.0155639260218379,
0.0155639260218379, 0.0415629168130336, 0.0415629168130336,
0.0415629168130336, 0.0415629168130336, 0.0415629168130336,
8.23198683529589e-05, 8.23198683529589e-05, 8.23198683529589e-05,
8.23198683529589e-05, 8.23198683529589e-05, 8.23198683529589e-05,
8.23198683529589e-05, 8.23198683529589e-05, 8.23198683529589e-05,
8.23198683529589e-05, 0.0142040710077886, 0.0142040710077886,
0.0142040710077886, 0.0142040710077886, 0.0142040710077886,
0.0325206196276115, 0.0325206196276115, 0.0325206196276115,
0.0325206196276115, 0.0325206196276115, 0.0325206196276115,
0.00933229305832333, 0.00933229305832333, 0.00933229305832333,
0.0311276543814304, 0.0311276543814304, 0.0311276543814304,
0.00309000052106621, 0.00309000052106621, 0.00309000052106621,
0.00309000052106621, 0.00309000052106621, 0.00309000052106621,
0.00309000052106621, 0.00309000052106621, 0.00309000052106621,
0.00309000052106621, 0.00309000052106621, 0.00309000052106621,
0.0218350076533092, 0.0218350076533092, 0.0218350076533092,
0.0218350076533092, 0.0218350076533092, 0.0218350076533092,
0.000458350205345348, 0.000458350205345348, 0.000458350205345348,
0.000458350205345348, 0.000458350205345348, 0.000458350205345348,
0.0161661930995344, 0.0161661930995344, 0.0161661930995344,
0.0161661930995344, 0.0161661930995344, 0.0646769590754476,
0.0646769590754476, 0.0646769590754476, 0.0646769590754476,
0.0110046362375520, 0.0110046362375520, 0.0237532591069227,
0.0237532591069227, 0.00829500745735379, 0.00829500745735379,
0.0364483156107965
]
ds['weights_mc'] = ds['weights_mc'] / ds['weights_mc'].sum()
ds['weights_q_norm'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none'] = np.ones(88) / 88
return ds
def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
# for CMIP5-ALL_med_9P_v0.nc
ds['weights_ic_all'] = [
0.0306498936112816, 0.0391448485764449, 0.000412880542569758,
0.00107991659481872, 0.00107991659481872, 0.00107991659481872,
0.00107991659481872, 0.00107991659481872, 0.00107991659481872,
0.00979905036521699, 0.00637405217366092, 0.00637405217366092,
0.00637405217366092, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.000133208267797779, 0.000133208267797779, 0.000133208267797779,
0.00665379138213003, 0.0281684808990121, 0.0614801788762348,
0.0139972063232535, 0.0139972063232535, 0.0139972063232535,
0.0139972063232535, 0.0139972063232535, 0.000710470681476458,
0.000710470681476458, 0.000710470681476458, 0.000710470681476458,
0.000710470681476458, 0.000710470681476458, 0.000710470681476458,
0.000710470681476458, 0.000710470681476458, 0.000710470681476458,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.000147220282127516,
0.000147220282127516, 0.000147220282127516, 0.00125879282299671,
0.00125879282299671, 0.00125879282299671, 0.00125879282299671,
0.00125879282299671, 0.0170988328647253, 0.0170988328647253,
0.0170988328647253, 0.0170988328647253, 0.0170988328647253,
0.0289160324694173, 0.00117020740708971, 0.00117020740708971,
0.00117020740708971, 0.0194950592035346, 0.00199152893081816,
0.000136576952445064, 0.0154925503807138, 0.0100360852066653,
0.0219662147817345, 0.0178402980037795, 0.0100360852066653,
0.0178402980037795, 0.0264391154516917, 0.00931062472382337,
0.0225549594485357, 0.00894476758728697, 0.00931062472382337,
0.00894476758728697, 0.0306318927945253, 0.0491376700487571,
0.00881383469632304, 0.00881383469632304, 0.00881383469632304,
0.00881383469632304, 0.00444698964579276, 0.00444698964579276,
0.00444698964579276, 0.00444698964579276, 0.0398942009239250,
0.00816167221145302, 0.00469211374490190, 0.00552942006092372,
0.00352896874131218, 0.00352896874131218, 0.00352896874131218,
0.0107981283433813, 0.0107981283433813, 0.0107981283433813,
0.0360523995255677, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.000145514276845206,
0.000145514276845206, 0.000145514276845206, 0.0191301244960145,
0.0136540976648741, 0.0193373418176255, 0.00841189054061916,
0.00271082723842561, 0.0132330275876382, 0.0212700194457166
]
ds['weights_ic_all'] = ds['weights_ic_all'] / ds['weights_ic_all'].sum()
ds['weights_mc_all'] = [
0.0390667137981572, 0.0390667137981572, 0.000924418400802522,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.000273066666975135, 0.000273066666975135, 0.000273066666975135,
0.0239573620813699, 0.0239573620813699, 0.0239573620813699,
0.0313390284863299, 0.0313390284863299, 0.0313390284863299,
0.0313390284863299, 0.0313390284863299, 0.00159070748914398,
0.00159070748914398, 0.00159070748914398, 0.00159070748914398,
0.00159070748914398, 0.00159070748914398, 0.00159070748914398,
0.00159070748914398, 0.00159070748914398, 0.00159070748914398,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.000295704465229955, 0.000295704465229955,
0.000295704465229955, 0.0382834115505173, 0.0382834115505173,
0.0382834115505173, 0.0382834115505173, 0.0382834115505173,
0.0647415165814365, 0.00262003448536538, 0.00262003448536538,
0.00262003448536538, 0.00537923997887741, 0.00537923997887741,
0.00537923997887741, 0.00421326541590348, 0.00421326541590348,
0.00421326541590348, 0.00421326541590348, 0.00421326541590348,
0.00421326541590348, 0.00421326541590348, 0.00421326541590348,
0.00421326541590348, 0.00421326541590348, 0.00421326541590348,
0.00421326541590348, 0.0137316545499640, 0.0137316545499640,
0.0137316545499640, 0.0137316545499640, 0.0137316545499640,
0.0137316545499640, 0.00741388764351201, 0.00741388764351201,
0.00741388764351201, 0.00741388764351201, 0.00741388764351201,
0.00741388764351201, 0.00375984457450244, 0.00375984457450244,
0.00375984457450244, 0.00375984457450244, 0.00375984457450244,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.000328799719067376,
0.000328799719067376, 0.000328799719067376, 0.0183505488251501,
0.0183505488251501, 0.0155322777143053, 0.0155322777143053,
0.00892436871038025, 0.00892436871038025, 0.0476224847958918
]
ds['weights_mc_all'] = ds['weights_mc_all'] / ds['weights_mc_all'].sum()
ds['weights_q_norm_all'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none_all'] = np.ones(288) / 288
ds['weights_ic'] = [
0.0140339967740686, 0.0178654487311369, 0.000210302156713192,
0.000515809393998982, 0.000515809393998982, 0.000515809393998982,
0.000515809393998982, 0.000515809393998982, 0.000515809393998982,
0.00458976973766726, 0.00298228443567221, 0.00298228443567221,
0.00298228443567221, 0.00310444283674798, 0.0131891308489770,
0.0284861601330155, 0.00635297383168565, 0.00635297383168565,
0.00635297383168565, 0.00635297383168565, 0.00635297383168565,
0.000329624231026383, 0.000329624231026383, 0.000329624231026383,
0.000329624231026383, 0.000329624231026383, 0.000329624231026383,
0.000329624231026383, 0.000329624231026383, 0.000329624231026383,
0.000329624231026383, 0.000579360270382479, 0.000579360270382479,
0.000579360270382479, 0.000579360270382479, 0.000579360270382479,
0.00698767794336135, 0.00698767794336135, 0.00698767794336135,
0.00698767794336135, 0.00698767794336135, 0.00698767794336135,
0.000552791371608343, 0.000552791371608343, 0.000552791371608343,
0.00926287211927062, 0.00100137048089641, 7.65217240752083e-05,
0.00729264684946363, 0.00467439831770549, 0.0101799902432345,
0.00825362211846436, 0.00467439831770549, 0.00825362211846436,
0.0121236966293752, 0.00434459425161063, 0.0105525268947770,
0.00422843102259940, 0.00434459425161063, 0.00422843102259940,
0.0140000340672880, 0.0221904157052906, 0.00403707358817109,
0.00403707358817109, 0.00403707358817109, 0.00403707358817109,
0.00202973315111198, 0.00202973315111198, 0.00202973315111198,
0.00202973315111198, 0.0180901613072435, 0.00378795537800128,
0.00220108418241656, 0.00256204931982498, 0.00162613209617560,
0.00162613209617560, 0.00162613209617560, 0.00504960668064616,
0.00504960668064616, 0.00504960668064616, 0.0165779411890502,
0.00880913338401020, 0.00637829160577958, 0.00888171831567428,
0.00398178829276437, 0.00132879905976752, 0.00638567088192441,
0.00975534657615189
]
ds['weights_ic'] = ds['weights_ic'] / ds['weights_ic'].sum()
ds['weights_mc'] = [
0.00797486137630137, 0.00797486137630137, 0.000210302156713192,
0.000499994678426805, 0.000499994678426805, 0.000499994678426805,
0.000499994678426805, 0.000499994678426805, 0.000499994678426805,
0.000499994678426805, 0.000499994678426805, 0.000499994678426805,
0.000499994678426805, 0.00497552597986006, 0.00497552597986006,
0.00497552597986006, 0.00635297383168565, 0.00635297383168565,
0.00635297383168565, 0.00635297383168565, 0.00635297383168565,
0.000329624231026383, 0.000329624231026383, 0.000329624231026383,
0.000329624231026383, 0.000329624231026383, 0.000329624231026383,
0.000329624231026383, 0.000329624231026383, 0.000329624231026383,
0.000329624231026383, 0.000579360270382479, 0.000579360270382479,
0.000579360270382479, 0.000579360270382479, 0.000579360270382479,
0.00698767794336135, 0.00698767794336135, 0.00698767794336135,
0.00698767794336135, 0.00698767794336135, 0.00698767794336135,
0.000552791371608343, 0.000552791371608343, 0.000552791371608343,
0.00114897381380469, 0.00114897381380469, 0.00114897381380469,
0.000876062801794235, 0.000876062801794235, 0.000876062801794235,
0.000876062801794235, 0.000876062801794235, 0.000876062801794235,
0.000876062801794235, 0.000876062801794235, 0.000876062801794235,
0.000876062801794235, 0.000876062801794235, 0.000876062801794235,
0.00279954519953656, 0.00279954519953656, 0.00279954519953656,
0.00279954519953656, 0.00279954519953656, 0.00279954519953656,
0.00150982908619546, 0.00150982908619546, 0.00150982908619546,
0.00150982908619546, 0.00150982908619546, 0.00150982908619546,
0.000775932894712877, 0.000775932894712877, 0.000775932894712877,
0.000775932894712877, 0.000775932894712877, 0.00387652508217910,
0.00387652508217910, 0.00387652508217910, 0.00387652508217910,
0.00379685624744744, 0.00379685624744744, 0.00321587665210966,
0.00321587665210966, 0.00192861748542298, 0.00192861748542298,
0.00975534657615189
]
ds['weights_mc'] = ds['weights_mc'] / ds['weights_mc'].sum()
ds['weights_q_norm'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none'] = np.ones(88) / 88
return ds
def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
# for CMIP5-ALL_med_9P_DA_fix.nc
ds['weights_ic_all'] = [
0.0173645245053380, 0.0203569712595225, 0.000386623999800020,
0.00175385476816864, 0.00175385476816864, 0.00175385476816864,
0.00175385476816864, 0.00175385476816864, 0.00175385476816864,
0.0149877194350907, 0.00534133867487110, 0.00534133867487110,
0.00534133867487110, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.000150857296767254, 0.000150857296767254, 0.000150857296767254,
0.0111500739284606, 0.0173243017086583, 0.0254478697381020,
0.00660444157632810, 0.00660444157632810, 0.00660444157632810,
0.00660444157632810, 0.00660444157632810, 0.000789256962921410,
0.000789256962921410, 0.000789256962921410, 0.000789256962921410,
0.000789256962921410, 0.000789256962921410, 0.000789256962921410,
0.000789256962921410, 0.000789256962921410, 0.000789256962921410,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 7.81218954017782e-05,
7.81218954017782e-05, 7.81218954017782e-05, 0.000657141865269707,
0.000657141865269707, 0.000657141865269707, 0.000657141865269707,
0.000657141865269707, 0.00975293795439704, 0.00975293795439704,
0.00975293795439704, 0.00975293795439704, 0.00975293795439704,
0.0205828413526940, 0.000777714869769811, 0.000777714869769811,
0.000777714869769811, 0.00936654632862810, 0.00150766657362369,
0.000257303183878334, 0.0101293911099127, 0.00634237391463581,
0.0136328786842070, 0.0108653021490644, 0.00634237391463581,
0.0108653021490644, 0.0149985193452363, 0.00590432468557565,
0.0149412285347167, 0.00592697102067078, 0.00590432468557565,
0.00592697102067078, 0.0163786490319080, 0.0266036508275658,
0.00591207366188384, 0.00591207366188384, 0.00591207366188384,
0.00591207366188384, 0.00113746952131555, 0.00113746952131555,
0.00113746952131555, 0.00113746952131555, 0.0130571399126652,
0.00237580290248032, 0.00438123681715656, 0.00592322682981926,
0.00239185126410343, 0.00239185126410343, 0.00239185126410343,
0.00646243993086172, 0.00646243993086172, 0.00646243993086172,
0.0189238866516331, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 7.69126865451268e-05,
7.69126865451268e-05, 7.69126865451268e-05, 0.00807208094414755,
0.00748992176859595, 0.0293103380988083, 0.0200810011517527,
0.00433793743220949, 0.0119406692125648, 0.0108898285395486
]
ds['weights_ic_all'] = ds['weights_ic_all'] / ds['weights_ic_all'].sum()
ds['weights_mc_all'] = [
0.00943037394121511, 0.00943037394121511, 0.000386623999800020,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.000139817161411427, 0.000139817161411427, 0.000139817161411427,
0.00599136059724678, 0.00599136059724678, 0.00599136059724678,
0.00660444157632810, 0.00660444157632810, 0.00660444157632810,
0.00660444157632810, 0.00660444157632810, 0.000789256962921410,
0.000789256962921410, 0.000789256962921410, 0.000789256962921410,
0.000789256962921410, 0.000789256962921410, 0.000789256962921410,
0.000789256962921410, 0.000789256962921410, 0.000789256962921410,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 6.99944744251862e-05, 6.99944744251862e-05,
6.99944744251862e-05, 0.00975293795439704, 0.00975293795439704,
0.00975293795439704, 0.00975293795439704, 0.00975293795439704,
0.0205828413526940, 0.000777714869769811, 0.000777714869769811,
0.000777714869769811, 0.00123683512068112, 0.00123683512068112,
0.00123683512068112, 0.00117956878301291, 0.00117956878301291,
0.00117956878301291, 0.00117956878301291, 0.00117956878301291,
0.00117956878301291, 0.00117956878301291, 0.00117956878301291,
0.00117956878301291, 0.00117956878301291, 0.00117956878301291,
0.00117956878301291, 0.00382154106804487, 0.00382154106804487,
0.00382154106804487, 0.00382154106804487, 0.00382154106804487,
0.00382154106804487, 0.000934234865449844, 0.000934234865449844,
0.000934234865449844, 0.000934234865449844, 0.000934234865449844,
0.000934234865449844, 0.00127324500095627, 0.00127324500095627,
0.00127324500095627, 0.00127324500095627, 0.00127324500095627,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 7.82371219934039e-05,
7.82371219934039e-05, 7.82371219934039e-05, 0.00389050067818587,
0.00389050067818587, 0.0123478348126402, 0.0123478348126402,
0.00406965166119357, 0.00406965166119357, 0.0108898285395486
]
ds['weights_mc_all'] = ds['weights_mc_all'] / ds['weights_mc_all'].sum()
ds['weights_q_norm_all'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none_all'] = np.ones(288) / 288
ds['weights_ic'] = [
0.0174718770203454, 0.0205613863535534, 0.000412558961838729,
0.00181139252421693, 0.00181139252421693, 0.00181139252421693,
0.00181139252421693, 0.00181139252421693, 0.00181139252421693,
0.0152283273055325, 0.00550018051180113, 0.00550018051180113,
0.00550018051180113, 0.0112088503665843, 0.0176843558125814,
0.0260773435951545, 0.00665451624081320, 0.00665451624081320,
0.00665451624081320, 0.00665451624081320, 0.00665451624081320,
0.000801662917011347, 0.000801662917011347, 0.000801662917011347,
0.000801662917011347, 0.000801662917011347, 0.000801662917011347,
0.000801662917011347, 0.000801662917011347, 0.000801662917011347,
0.000801662917011347, 0.000656719634757235, 0.000656719634757235,
0.000656719634757235, 0.000656719634757235, 0.000656719634757235,
0.00899614071852682, 0.00899614071852682, 0.00899614071852682,
0.00899614071852682, 0.00899614071852682, 0.00899614071852682,
0.000781379408867015, 0.000781379408867015, 0.000781379408867015,
0.00972743527533972, 0.00163266888094045, 0.000301337471280230,
0.0105003069513024, 0.00651606493379239, 0.0139223845374914,
0.0111001101725287, 0.00651606493379239, 0.0111001101725287,
0.0152014391759463, 0.00607054179337712, 0.0154118121987557,
0.00616597283229141, 0.00607054179337712, 0.00616597283229141,
0.0164066070898815, 0.0265744029546302, 0.00595595344081403,
0.00595595344081403, 0.00595595344081403, 0.00595595344081403,
0.00112909631444207, 0.00112909631444207, 0.00112909631444207,
0.00112909631444207, 0.0130100732010882, 0.00237656896992004,
0.00442777010500610, 0.00598900416553578, 0.00239570964391535,
0.00239570964391535, 0.00239570964391535, 0.00662762048171131,
0.00662762048171131, 0.00662762048171131, 0.0190513881039158,
0.00808316367371425, 0.00755266722505562, 0.0294790423288641,
0.0205224435384989, 0.00453997994364382, 0.0125080241634372,
0.0109314202361373
]
ds['weights_ic'] = ds['weights_ic'] / ds['weights_ic'].sum()
ds['weights_mc'] = [
0.00950831584347471, 0.00950831584347471, 0.000412558961838729,
0.00129940082783552, 0.00129940082783552, 0.00129940082783552,
0.00129940082783552, 0.00129940082783552, 0.00129940082783552,
0.00129940082783552, 0.00129940082783552, 0.00129940082783552,
0.00129940082783552, 0.00610783886381335, 0.00610783886381335,
0.00610783886381335, 0.00665451624081320, 0.00665451624081320,
0.00665451624081320, 0.00665451624081320, 0.00665451624081320,
0.000801662917011347, 0.000801662917011347, 0.000801662917011347,
0.000801662917011347, 0.000801662917011347, 0.000801662917011347,
0.000801662917011347, 0.000801662917011347, 0.000801662917011347,
0.000801662917011347, 0.000656719634757235, 0.000656719634757235,
0.000656719634757235, 0.000656719634757235, 0.000656719634757235,
0.00899614071852682, 0.00899614071852682, 0.00899614071852682,
0.00899614071852682, 0.00899614071852682, 0.00899614071852682,
0.000781379408867015, 0.000781379408867015, 0.000781379408867015,
0.00129571573639560, 0.00129571573639560, 0.00129571573639560,
0.00121143542910732, 0.00121143542910732, 0.00121143542910732,
0.00121143542910732, 0.00121143542910732, 0.00121143542910732,
0.00121143542910732, 0.00121143542910732, 0.00121143542910732,
0.00121143542910732, 0.00121143542910732, 0.00121143542910732,
0.00384100736382045, 0.00384100736382045, 0.00384100736382045,
0.00384100736382045, 0.00384100736382045, 0.00384100736382045,
0.000929227311168926, 0.000929227311168926, 0.000929227311168926,
0.000929227311168926, 0.000929227311168926, 0.000929227311168926,
0.00127912644263120, 0.00127912644263120, 0.00127912644263120,
0.00127912644263120, 0.00127912644263120, 0.00491874827745735,
0.00491874827745735, 0.00491874827745735, 0.00491874827745735,
0.00390895772469247, 0.00390895772469247, 0.0125003714668407,
0.0125003714668407, 0.00426200102677025, 0.00426200102677025,
0.0109314202361373
]
ds['weights_mc'] = ds['weights_mc'] / ds['weights_mc'].sum()
ds['weights_q_norm'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none'] = np.ones(88) / 88
return ds
def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
return ds
def read_data(filename, path):
ds = xr.open_dataset(os.path.join(path, filename))
ds = area_weighted_mean(ds)
# for CMIP5-ALL_neu_9P_v2.nc
ds['weights_mc_all'] = [
0.0519684495104962, 0.0519684495104962, 0.00528928315039839,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.000196826232546333, 0.000196826232546333, 0.000196826232546333,
0.00843590211664639, 0.00843590211664639, 0.00843590211664639,
0.0211853533713132, 0.0211853533713132, 0.0211853533713132,
0.0211853533713132, 0.0211853533713132, 2.43842579015591e-05,
2.43842579015591e-05, 2.43842579015591e-05, 2.43842579015591e-05,
2.43842579015591e-05, 2.43842579015591e-05, 2.43842579015591e-05,
2.43842579015591e-05, 2.43842579015591e-05, 2.43842579015591e-05,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.00121271885331108, 0.00121271885331108,
0.00121271885331108, 0.0461537219125901, 0.0461537219125901,
0.0461537219125901, 0.0461537219125901, 0.0461537219125901,
0.00214294419723561, 0.00324148233535707, 0.00324148233535707,
0.00324148233535707, 0.0223349242532901, 0.0223349242532901,
0.0223349242532901, 0.00156503148476069, 0.00156503148476069,
0.00156503148476069, 0.00156503148476069, 0.00156503148476069,
0.00156503148476069, 0.00156503148476069, 0.00156503148476069,
0.00156503148476069, 0.00156503148476069, 0.00156503148476069,
0.00156503148476069, 0.0149216035052486, 0.0149216035052486,
0.0149216035052486, 0.0149216035052486, 0.0149216035052486,
0.0149216035052486, 0.000345992901257310, 0.000345992901257310,
0.000345992901257310, 0.000345992901257310, 0.000345992901257310,
0.000345992901257310, 0.00883044830394363, 0.00883044830394363,
0.00883044830394363, 0.00883044830394363, 0.00883044830394363,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00157976836762813,
0.00157976836762813, 0.00157976836762813, 0.00331303452460497,
0.00331303452460497, 0.0122832682662121, 0.0122832682662121,
0.00266064244524378, 0.00266064244524378, 0.0157976191507417
]
ds['weights_mc_all'] = ds['weights_mc_all'] / ds['weights_mc_all'].sum()
ds['weights_ic_all'] = [
0.0506873152345471, 0.0500870001010942, 0.00256417063198165,
0.000863803896267789, 0.000863803896267789, 0.000863803896267789,
0.000863803896267789, 0.000863803896267789, 0.000863803896267789,
0.00374726110951423, 0.00849275823070520, 0.00849275823070520,
0.00849275823070520, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
9.28911879448266e-05, 9.28911879448266e-05, 9.28911879448266e-05,
0.000414496911640593, 0.0140805986517069, 0.0223113726517877,
0.0102703635631197, 0.0102703635631197, 0.0102703635631197,
0.0102703635631197, 0.0102703635631197, 1.18211478221080e-05,
1.18211478221080e-05, 1.18211478221080e-05, 1.18211478221080e-05,
1.18211478221080e-05, 1.18211478221080e-05, 1.18211478221080e-05,
1.18211478221080e-05, 1.18211478221080e-05, 1.18211478221080e-05,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.000654579491439552,
0.000654579491439552, 0.000654579491439552, 0.00567906046413463,
0.00567906046413463, 0.00567906046413463, 0.00567906046413463,
0.00567906046413463, 0.0223746800690747, 0.0223746800690747,
0.0223746800690747, 0.0223746800690747, 0.0223746800690747,
0.00103886943093851, 0.00157142538451243, 0.00157142538451243,
0.00157142538451243, 0.0812793305187681, 0.00226701554333645,
0.0139025892442928, 0.00753333263494764, 0.00816363349905153,
0.0117666962195444, 0.00570982496298650, 0.00816363349905153,
0.00570982496298650, 0.00935050783219849, 0.00291145854301123,
0.00326982068608221, 0.00254838861676683, 0.00291145854301123,
0.00254838861676683, 0.0341249993798806, 0.00830690565725499,
0.0136240223451796, 0.0136240223451796, 0.0136240223451796,
0.0136240223451796, 0.000128454343999261, 0.000128454343999261,
0.000128454343999261, 0.000128454343999261, 0.00398206036946370,
1.04114920868242e-06, 0.00515866393160315, 0.00683931954433429,
0.0105582186780247, 0.0105582186780247, 0.0105582186780247,
0.0297544144936968, 0.0297544144936968, 0.0297544144936968,
0.0656230909358085, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.000795001700145025,
0.000795001700145025, 0.000795001700145025, 0.00204225983692364,
0.00438219171099055, 0.0147223657308445, 0.00909666218457424,
0.00107845081009328, 0.00408091844756054, 0.00765846522671256
]
ds['weights_ic_all'] = ds['weights_ic_all'] / ds['weights_ic_all'].sum()
ds['weights_q_norm_all'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none_all'] = np.ones(288) / 288
# for CMIP5_neu_9P_v2.nc
ds['weights_ic'] = [
0.0569924821317768, 0.0615528850278657, 0.00287684839617446,
0.000989291623655752, 0.000989291623655752, 0.000989291623655752,
0.000989291623655752, 0.000989291623655752, 0.000989291623655752,
0.00362645200930820, 0.00957453823734612, 0.00957453823734612,
0.00957453823734612, 0.000301128471870280, 0.0149254076527316,
0.0268280356898401, 0.0125544850039963, 0.0125544850039963,
0.0125544850039963, 0.0125544850039963, 0.0125544850039963,
1.09775298295311e-05, 1.09775298295311e-05, 1.09775298295311e-05,
1.09775298295311e-05, 1.09775298295311e-05, 1.09775298295311e-05,
1.09775298295311e-05, 1.09775298295311e-05, 1.09775298295311e-05,
1.09775298295311e-05, 0.00649887889305261, 0.00649887889305261,
0.00649887889305261, 0.00649887889305261, 0.00649887889305261,
0.0233506495019221, 0.0233506495019221, 0.0233506495019221,
0.0233506495019221, 0.0233506495019221, 0.0233506495019221,
0.00191938091945592, 0.00191938091945592, 0.00191938091945592,
0.0982006580148391, 0.00198923648617189, 0.0172702473315458,
0.00810434428732132, 0.0100558218242553, 0.0135510743387476,
0.00705077718919593, 0.0100558218242553, 0.00705077718919593,
0.0108517344791539, 0.00301334534638672, 0.00333058850947048,
0.00279290654674263, 0.00301334534638672, 0.00279290654674263,
0.0412561815761165, 0.00815363175620175, 0.0166144320893954,
0.0166144320893954, 0.0166144320893954, 0.0166144320893954,
0.000132585689409592, 0.000132585689409592, 0.000132585689409592,
0.000132585689409592, 0.00460635814363639, 7.40443205286654e-07,
0.00650119361831786, 0.00854252486032772, 0.0128658720505988,
0.0128658720505988, 0.0128658720505988, 0.0364649310179668,
0.0364649310179668, 0.0364649310179668, 0.0808748419014684,
0.00228307904917499, 0.00529688231491314, 0.0180994187970468,
0.00983288695810432, 0.00106673597388643, 0.00410206611882681,
0.00910730781858947
]
ds['weights_ic'] = ds['weights_ic'] / ds['weights_ic'].sum()
ds['weights_mc'] = [
0.0148181708949553, 0.0148181708949553, 0.00287684839617446,
0.000125411794597030, 0.000125411794597030, 0.000125411794597030,
0.000125411794597030, 0.000125411794597030, 0.000125411794597030,
0.000125411794597030, 0.000125411794597030, 0.000125411794597030,
0.000125411794597030, 0.00155757673386822, 0.00155757673386822,
0.00155757673386822, 0.00251089700079925, 0.00251089700079925,
0.00251089700079925, 0.00251089700079925, 0.00251089700079925,
1.09775298295311e-06, 1.09775298295311e-06, 1.09775298295311e-06,
1.09775298295311e-06, 1.09775298295311e-06, 1.09775298295311e-06,
1.09775298295311e-06, 1.09775298295311e-06, 1.09775298295311e-06,
1.09775298295311e-06, 0.00129977577861052, 0.00129977577861052,
0.00129977577861052, 0.00129977577861052, 0.00129977577861052,
0.00467012990038442, 0.00467012990038442, 0.00467012990038442,
0.00467012990038442, 0.00467012990038442, 0.00467012990038442,
0.000639793639818639, 0.000639793639818639, 0.000639793639818639,
0.00435037562342803, 0.00435037562342803, 0.00435037562342803,
7.37784405329951e-05, 7.37784405329951e-05, 7.37784405329951e-05,
7.37784405329951e-05, 7.37784405329951e-05, 7.37784405329951e-05,
7.37784405329951e-05, 7.37784405329951e-05, 7.37784405329951e-05,
7.37784405329951e-05, 7.37784405329951e-05, 7.37784405329951e-05,
0.00145944780908632, 0.00145944780908632, 0.00145944780908632,
0.00145944780908632, 0.00145944780908632, 0.00145944780908632,
3.11503223027553e-05, 3.11503223027553e-05, 3.11503223027553e-05,
3.11503223027553e-05, 3.11503223027553e-05, 3.11503223027553e-05,
0.00104669253547228, 0.00104669253547228, 0.00104669253547228,
0.00104669253547228, 0.00104669253547228, 0.00639155032911203,
0.00639155032911203, 0.00639155032911203, 0.00639155032911203,
0.000947495170511015, 0.000947495170511015, 0.00349153821939390,
0.00349153821939390, 0.000646100261589155, 0.000646100261589155,
0.00910730781858947
]
ds['weights_mc'] = ds['weights_mc'] / ds['weights_mc'].sum()
ds['weights_q_norm'] = ds['weights_q'] / ds['weights_q'].sum()
ds['weights_none'] = np.ones(88) / 88
return ds