diff --git a/video_prediction_tools/postprocess/get_metrics_joint_dom.ipynb b/video_prediction_tools/postprocess/get_metrics_joint_dom.ipynb
index 6c7274bd16d25684531558433a894cbc5c39ccf0..e54043e10d1caf51352f75c11ed4cbde656d5aee 100644
--- a/video_prediction_tools/postprocess/get_metrics_joint_dom.ipynb
+++ b/video_prediction_tools/postprocess/get_metrics_joint_dom.ipynb
@@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "code",
- "execution_count": 8,
+ "execution_count": null,
"id": "d636ee12-e299-485f-b84d-6f35c05fa766",
"metadata": {},
"outputs": [],
@@ -33,7 +33,7 @@
},
{
"cell_type": "code",
- "execution_count": 62,
+ "execution_count": null,
"id": "440b15fa-ecd4-4bb4-9100-ede5abb2b04f",
"metadata": {},
"outputs": [],
@@ -58,20 +58,10 @@
},
{
"cell_type": "code",
- "execution_count": 63,
+ "execution_count": null,
"id": "e6726da3-d774-4eda-89d6-e315a865bb99",
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "File '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_0_999.nc' does not exist. \n",
- " Start reading data with sample index between 0 and 999 from '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/'...\n",
- "Data loaded successfully. Save merged data to '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_0_999.nc'.\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"def get_fname_ndigits(indir, prefix, suffix, n, patt=\"[0-9]\"):\n",
" flist = []\n",
@@ -108,23 +98,7 @@
"execution_count": null,
"id": "54f4aa3e-3a39-496e-ae97-65f79d9cd598",
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "File '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_1000_1999.nc' does not exist. Start reading data with sample index between 1000 and 1999 from '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/'...\n",
- "Data loaded successfully. Save merged data to '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_1000_1999.nc'.\n",
- "File '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_2000_2999.nc' does not exist. Start reading data with sample index between 2000 and 2999 from '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/'...\n",
- "Data loaded successfully. Save merged data to '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_2000_2999.nc'.\n",
- "File '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_3000_3999.nc' does not exist. Start reading data with sample index between 3000 and 3999 from '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/'...\n",
- "Data loaded successfully. Save merged data to '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_3000_3999.nc'.\n",
- "File '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_4000_4999.nc' does not exist. Start reading data with sample index between 4000 and 4999 from '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/'...\n",
- "Data loaded successfully. Save merged data to '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_4000_4999.nc'.\n",
- "File '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/vfp_savp_forecasts_sample_ind_5000_5999.nc' does not exist. Start reading data with sample index between 5000 and 5999 from '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-92x56-3840N0000E-2t_tcc_t_850/savp/20210901T090059_gong1_savp_cv12/'...\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"for i in np.arange(1, 9):\n",
" outfile = os.path.join(indir, \"vfp_{0}_forecasts_sample_ind_{1:d}000_{1:d}999.nc\".format(model, i))\n",
@@ -150,34 +124,10 @@
},
{
"cell_type": "code",
- "execution_count": 37,
+ "execution_count": null,
"id": "92f15edf-c23f-4803-b3c5-618305194de5",
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_0_999.nc'\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_1000_1999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_2000_2999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_3000_3999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_4000_4999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_5000_5999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_6000_6999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_7000_7999.nc'\n",
- "Start merging\n",
- "Processing file '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/../vfp_savp_forecasts_sample_ind_8000_8999.nc'\n",
- "Start merging\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"outfile_all = os.path.join(indir, \"vfp_{0}_forecasts_all.nc\".format(model))\n",
"\n",
@@ -186,7 +136,7 @@
" print(\"netCDF-file with all forecasts '{0}' does not exist yet. Start merging and sorting all precursor files.\".format(outfile))\n",
" all_files = sorted(glob.glob(os.path.join(indir, \"vfp_{0}_forecasts_sample_ind_*.nc\".format(model))))\n",
" \n",
- " if len(all_files == 0):\n",
+ " if len(all_files) == 0:\n",
" raise FileNotFoundError(\"Could not find any precursor files.\")\n",
"\n",
" for i, f in enumerate(all_files):\n",
@@ -208,51 +158,22 @@
]
},
{
- "cell_type": "code",
- "execution_count": null,
- "id": "a75d8205-4841-440f-865f-cbf3e26da29b",
+ "cell_type": "markdown",
+ "id": "0fcf1cb1-ba0d-4262-8e23-12ba44b6e2d0",
"metadata": {},
- "outputs": [],
"source": [
"Now, we slice the dataset to the domain of interest (defined by `lonlatbox`)."
]
},
{
"cell_type": "code",
- "execution_count": 45,
+ "execution_count": null,
"id": "ede23e56-5be8-48be-b584-0eb8741acbf3",
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "<xarray.Dataset>\n",
- "Dimensions: (fcst_hour: 12, in_hour: 12, init_time: 8440, lat: 44, lon: 72)\n",
- "Coordinates:\n",
- " * init_time (init_time) datetime64[ns] 2019-01-01T11:00:00 .....\n",
- " * in_hour (in_hour) int64 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0\n",
- " * lat (lat) float64 53.1 52.8 52.5 52.2 ... 40.8 40.5 40.2\n",
- " * lon (lon) float64 3.0 3.3 3.6 3.9 ... 23.7 24.0 24.3\n",
- " * fcst_hour (fcst_hour) int64 1 2 3 4 5 6 7 8 9 10 11 12\n",
- "Data variables:\n",
- " 2t_in (init_time, in_hour, lat, lon) float32 281.9 ... ...\n",
- " tcc_in (init_time, in_hour, lat, lon) float32 0.7001 ......\n",
- " t_850_in (init_time, in_hour, lat, lon) float32 278.9 ... ...\n",
- " 2t_ref (init_time, fcst_hour, lat, lon) float32 281.8 .....\n",
- " tcc_ref (init_time, fcst_hour, lat, lon) float32 0.7272 ....\n",
- " t_850_ref (init_time, fcst_hour, lat, lon) float32 272.4 .....\n",
- " 2t_savp_fcst (init_time, fcst_hour, lat, lon) float32 282.0 .....\n",
- " tcc_savp_fcst (init_time, fcst_hour, lat, lon) float32 0.4467 ....\n",
- " t_850_savp_fcst (init_time, fcst_hour, lat, lon) float32 273.0 .....\n",
- " 2t_persistence_fcst (init_time, fcst_hour, lat, lon) float64 281.9 .....\n",
- " tcc_persistence_fcst (init_time, fcst_hour, lat, lon) float64 0.7888 ....\n",
- " t_850_persistence_fcst (init_time, fcst_hour, lat, lon) float64 277.5 .....\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
- "all_fcst_sl = all_fcst.sel({\"lon\": slice(lonlatbox[0], lonlatbox[1]), \"lat\": slice(lonlatbox[3], lonlatbox[2])}) "
+ "all_fcst_sl = all_fcst.sel({\"lon\": slice(lonlatbox[0], lonlatbox[1]), \"lat\": slice(lonlatbox[3], lonlatbox[2])}) \n",
+ "print(all_fcst_sl)"
]
},
{
@@ -273,7 +194,7 @@
},
{
"cell_type": "code",
- "execution_count": 46,
+ "execution_count": null,
"id": "c2b70b80-6b86-4674-b051-6a23aaa821ea",
"metadata": {},
"outputs": [],
@@ -295,7 +216,7 @@
},
{
"cell_type": "code",
- "execution_count": 48,
+ "execution_count": null,
"id": "b49db031-126c-44b1-b649-4f70587fac89",
"metadata": {},
"outputs": [],
@@ -324,377 +245,10 @@
},
{
"cell_type": "code",
- "execution_count": 53,
+ "execution_count": null,
"id": "5090e71c-f20f-43e6-94f6-71cbd0b6006d",
"metadata": {},
- "outputs": [
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- " 10%|█ | 100/1000 [00:00<00:00, 994.19it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1001.31it/s]\n",
- " 10%|█ | 102/1000 [00:00<00:00, 1013.59it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1007.59it/s]\n",
- " 21%|██ | 206/1000 [00:00<00:00, 1013.05it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1014.17it/s]\n",
- " 10%|█ | 104/1000 [00:00<00:00, 1030.59it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1023.93it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1013.22it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1011.14it/s]\n",
- " 10%|█ | 104/1000 [00:00<00:00, 1032.73it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1027.81it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1013.01it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1014.60it/s]\n",
- " 21%|██ | 208/1000 [00:00<00:00, 1022.36it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1024.22it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1009.07it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1008.34it/s]\n",
- " 21%|██ | 206/1000 [00:00<00:00, 1015.52it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1024.34it/s]\n",
- " 10%|█ | 101/1000 [00:00<00:00, 1006.28it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1009.03it/s]\n",
- " 21%|██ | 206/1000 [00:00<00:00, 1009.63it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1018.56it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1005.11it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1010.13it/s]\n",
- " 10%|█ | 103/1000 [00:00<00:00, 1021.97it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1022.39it/s]\n",
- " 20%|██ | 201/1000 [00:00<00:00, 1001.34it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1008.24it/s]\n",
- " 20%|██ | 205/1000 [00:00<00:00, 1014.66it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1020.29it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1007.14it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1001.87it/s]\n",
- " 21%|██ | 206/1000 [00:00<00:00, 1014.55it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1019.30it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1009.72it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:01<00:00, 980.91it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 978.29it/s] "
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:01<00:00, 957.25it/s]\n",
- " 10%|█ | 100/1000 [00:00<00:00, 992.80it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1001.39it/s]\n",
- " 10%|█ | 104/1000 [00:00<00:00, 1033.97it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1020.41it/s]\n",
- " 20%|██ | 204/1000 [00:00<00:00, 1013.23it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1006.80it/s]\n",
- " 21%|██ | 206/1000 [00:00<00:00, 1012.80it/s]"
- ]
- },
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "%perform_block_bootstrap_metric: Start block bootstrapping...\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "100%|██████████| 1000/1000 [00:00<00:00, 1018.64it/s]\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"for i, fh in enumerate(fcst_hours):\n",
" mse_model_curr = mse_model_all.sel(fcst_hour=fh)\n",
@@ -707,25 +261,20 @@
},
{
"cell_type": "code",
- "execution_count": 56,
- "id": "61457486-27b3-41eb-bd0c-3c5d3f6ba079",
+ "execution_count": null,
+ "id": "d526324d-5d19-4193-8208-e609d9c65205",
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "<xarray.DataArray (fcst_hour: 12)>\n",
- "array([0.18470440804958344, 0.43441563844680786, 0.6988919973373413,\n",
- " 0.9624305963516235, 1.2138677835464478, 1.4549846649169922,\n",
- " 1.6799184083938599, 1.8890923261642456, 2.0873427391052246,\n",
- " 2.275247573852539, 2.45733642578125, 2.6395061016082764],\n",
- " dtype=object)\n",
- "Coordinates:\n",
- " * fcst_hour (fcst_hour) int64 1 2 3 4 5 6 7 8 9 10 11 12\n"
- ]
- }
- ],
+ "outputs": [],
+ "source": [
+ "print(mse_model_fcst)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "b42cd738-b966-4b24-ad13-351d9b88f9e8",
+ "metadata": {},
+ "outputs": [],
"source": [
"print(mse_model_fcst)"
]
@@ -740,28 +289,10 @@
},
{
"cell_type": "code",
- "execution_count": 61,
+ "execution_count": null,
"id": "7cd455ee-4749-46dd-8095-9e43744a1563",
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Save evaluation metrics to '/p/project/deepacf/deeprain/video_prediction_shared_folder/results/era5-Y2007-2019M01to12-80x48-3960N0180E-2t_tcc_t_850/savp/20210713T143850_gong1_savp_t2opt_3vars/evaluation_metrics_72x44.nc'\n",
- "<xarray.Dataset>\n",
- "Dimensions: (fcst_hour: 12, iboot: 1000)\n",
- "Coordinates:\n",
- " * fcst_hour (fcst_hour) int64 1 2 3 4 5 ... 9 10 11 12\n",
- " * iboot (iboot) int64 0 1 2 3 4 ... 996 997 998 999\n",
- "Data variables:\n",
- " 2t_savp_mse_avg (fcst_hour) object 0.1847 0.4344 ... 2.64\n",
- " 2t_savp_mse_bootstrapped (fcst_hour, iboot) object 0.1863 ... 2.704\n",
- " 2t_persistence_mse_avg (fcst_hour) object 7.498 7.497 ... 7.48\n",
- " 2t_persistence_mse_bootstrapped (fcst_hour, iboot) object 7.768 ... 7.556\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"# create Dataset and save to netCDF-file\n",
"ds_mse = xr.Dataset({\"2t_{0}_mse_avg\".format(model): mse_model_fcst, \"2t_{0}_mse_bootstrapped\".format(model): mse_model_fcst_boot, \n",