diff --git a/FIGSEXAMPLES/p-value_map_39_202001-202012.png b/FIGSEXAMPLES/p-value_map_39_202001-202012.png
new file mode 100644
index 0000000000000000000000000000000000000000..412ae0410f699703cbb2d49202dae65e79c88d07
Binary files /dev/null and b/FIGSEXAMPLES/p-value_map_39_202001-202012.png differ
diff --git a/FIGSEXAMPLES/rank-hist-sumstats_39_202001-202012.png b/FIGSEXAMPLES/rank-hist-sumstats_39_202001-202012.png
new file mode 100644
index 0000000000000000000000000000000000000000..f9e85dd44739b883e70f3ba4a047da05cef8428a
Binary files /dev/null and b/FIGSEXAMPLES/rank-hist-sumstats_39_202001-202012.png differ
diff --git a/FIGSEXAMPLES/standard_plot_39_126736_202001-202012_python.png b/FIGSEXAMPLES/standard_plot_39_126736_202001-202012_python.png
new file mode 100644
index 0000000000000000000000000000000000000000..6d207aecf7415ce563abf9d648add11f172ed7d8
Binary files /dev/null and b/FIGSEXAMPLES/standard_plot_39_126736_202001-202012_python.png differ
diff --git a/README.md b/README.md
index 2d717e556197c7400df6202cc6b38dffcbfa74ff..c416d6aa869447a00039460375c5a20ff3673c97 100644
--- a/README.md
+++ b/README.md
@@ -14,10 +14,10 @@ The OBS application is currently being developed as a set of scripts (written in
 3. `synop_obs.sh` - bash-script to run reading and extraction (using `sql` with `odb_api` software) of selected meteorological variables at hourly time-slices for observations from ground-based synoptical stations at fixed at the surface geographical (latitude, longitude) locations over selected geographical domain (or over the globe) and saving (temporary) to dat-file. Note, that list of available  meteorological variables in open multi-year data ODB includes: 91 - total amount of clouds; 108 - sea level pressure; 80 - 1-hour precipitation; 58 - relative humidity; 999 - 10-minute precipitation intensity; 71 - snow depth; 39 - 2m temperature; 40 - dew point temperature; 62 - visibility; 111 - wind direction; 261 - maximum wind gust in last 10 minutes; 221 -surface wind speed.
 4. `synop_mod.sh` - bash-script to run reading, extraction and interpolation (using `cdo` operators) of modelled data for selected meteorological variables at the same hourly time-slices to geographical (latitude, longitude) locations of synoptical stations, and saving (temporary) to dat-file, and adding (using `import` with `odb_api`) such data to ODB. Interpolation is preceeded by constructing unstructured grid based on locations of synop stations and calculating weights for such grid. 
 5. `graph_mod_obs.py` - python-script to run an internal self-control in calculating differences between observed and modelled data for locations of synoptical stations 
-6. `synop_stats.sh` - bash-script to run (addtional `bash`, `fortran`, `python` scripts) calculating/producing
-- (i) rank histograms for all synop stations for 00, 06, 12 and 18 UTCs (with `produce_rank_histograms_all_stations.sh`, using `rank_histograms_one_station.f95` and `sql/import` with `oadb_api`), 
-- (ii) standard plots for each synop station (with `produce_standard_plots_all_stations.sh` using `plot_quantiles_rankhist.py` and  `sql` with `oadb_api`; see an example of such plot), and 
-- (iii) summary rank histograms for all stations (`summary_rank_histograms_all_stations.sh` using `rank_histogram_summary_statistics.f95`, `plot_p_values_map.py`, `plot_rank_hist_sum_all_stations.py` and `sql` with `oadb_api`; see examples of such plots). 
+6. `synop_stats.sh` - bash-script to run (addtional `bash`, `fortran`, `python` scripts) calculating/producing:
+- (i) rank histograms for all synop stations for 00, 06, 12 and 18 UTCs (with `produce_rank_histograms_all_stations.sh`, using `rank_histograms_one_station.f95` and `sql/import` with `odb_api`), 
+- (ii) standard plots for each synop station (with `produce_standard_plots_all_stations.sh` using `plot_quantiles_rankhist.py` and  `sql` with `oadb_api`; see an example of such plot ![Semantic description of image](FIGSEXAMPLES/standard_plot_39_126736_202001-202012_python.png "Image Title 1"), and 
+- (iii) summary rank histograms for all stations (`summary_rank_histograms_all_stations.sh` using `rank_histogram_summary_statistics.f95`, `plot_p_values_map.py`, `plot_rank_hist_sum_all_stations.py` and `sql` with `odb_api`; see examples of such plots ![Semantic description of image](FIGSEXAMPLES/p-value_map_39_202001-202012.png "Image Title 2") and ![Semantic description of image](FIGSEXAMPLES/rank-hist-sumstats_39_202001-202012.png "Image Title 3")). 
 - Note, needed input includes: list of synop stations with geographical coordinates, modelled data at stations coordinates, pre-computed quantiles as a function of time of year and rank histogram bootstrap mean square deviation (MSD) values. 
 
 # Disclaimers