diff --git a/heave_correction.md b/heave_correction.md
index 16e9691..7870620 100644
--- a/heave_correction.md
+++ b/heave_correction.md
@@ -183,26 +183,17 @@ def heave_correction(moments, date):
 
 ### Quality control
 
-![uncorrected MDV](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_MDV_uncorrected.png)   
-Uncorrected Doppler Velocity  
+| Uncorrected Doppler Velocity                                 | Corrected Doppler Velocity                                   |
+| ------------------------------------------------------------ | ------------------------------------------------------------ |
+| ![uncorrected MDV](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_MDV_uncorrected.png) | ![corrected MDV](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_MDV_corrected.png) |
+| ![uncorrected MDV zoom](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_MDV_uncorrected.png) | ![corrected MDV zoom](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_MDV_corrected.png) |
 
-![uncorrected MDV zoom](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_MDV_uncorrected.png)  
-Zoom in of uncorrected Doppler Velocity
+  
 
-![heave correction](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_heave_correction.png)  
-Heave correction as calculated by function  
-
-![corrected MDV](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_MDV_corrected.png)  
-Corrected Doppler Velocity 
-
-![corrected MDV zoom](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_MDV_corrected.png)  
-Zoom in of corrected Doppler Velocity
-
-![difference between corrected and uncorrected](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_MDV_corrected-measured.png)    
-Difference between corrected and uncorrected Doppler Velocity
-
-![difference between corrected and uncorrected zoom](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_MDV_corrected-measured.png)      
-Zoom in on difference between corrected and uncorrected Doppler Velocity  
+| heave correction                                             | corrected - uncorrected                                      |
+| ------------------------------------------------------------ | ------------------------------------------------------------ |
+| ![heave correction](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_heave_correction.png) | ![difference between corrected and uncorrected](quality_control\20200205_0000_20200205_2359_3km_cloudnet_input_MDV_corrected-measured.png) |
+| ![20200205_0900_20200205_1100_3km_cloudnet_input_heave_correction](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_heave_correction.png) | ![difference between corrected and uncorrected zoom](quality_control\20200205_0900_20200205_1100_3km_cloudnet_input_MDV_corrected-measured.png) |
 
 ![heave components unnormalized](quality_control\tmp.png)  
 Heave components unnormalized, chirp 1,2 and 3 (from top to bottom)
diff --git a/heave_correction.py b/heave_correction.py
index eecc983..a8399d5 100644
--- a/heave_correction.py
+++ b/heave_correction.py
@@ -48,7 +48,7 @@
 ########################################################################################################################
 plot_path = "/projekt1/remsens/work/jroettenbacher/plots/heave_correction"
 plot_range = [0, 3000]
-# mdv['var_lims'] = [-7, 7]
+mdv['var_lims'] = [-7, 7]
 # uncorrected MDV
 fig, _ = pyLARDA.Transformations.plot_timeheight(mdv, rg_converter=False, title=True, range_interval=plot_range)
 name = f'{plot_path}/{begin_dt:%Y%m%d_%H%M}_{end_dt:%Y%m%d_%H%M}_{plot_range[1] / 1000:.0f}km_cloudnet_input'
@@ -59,7 +59,7 @@
 
 # uncorrected MDV zoom
 begin_dt = dt.datetime(2020, 2, 5, 9, 0, 0)
-end_dt = dt.datetime(2020, 2, 5, 11, 0, 0)
+end_dt = dt.datetime(2020, 2, 5, 10, 0, 0)
 fig, _ = pyLARDA.Transformations.plot_timeheight(mdv, rg_converter=False, title=True, time_interval=[begin_dt, end_dt],
                                                  range_interval=plot_range)
 name = f'{plot_path}/{begin_dt:%Y%m%d_%H%M}_{end_dt:%Y%m%d_%H%M}_{plot_range[1] / 1000:.0f}km_cloudnet_input'
@@ -69,7 +69,7 @@
 print(f'figure saved :: {fig_name}')
 
 Vel_cor = moments['Vel_cor']
-# Vel_cor['var_lims'] = [-7, 7]
+Vel_cor['var_lims'] = [-7, 7]
 # corrected MDV
 begin_dt = dt.datetime(2020, 2, 5, 0, 0, 5)
 end_dt = dt.datetime(2020, 2, 5, 23, 59, 55)
@@ -82,7 +82,7 @@
 
 # corrected MDV zoom
 begin_dt = dt.datetime(2020, 2, 5, 9, 0, 0)
-end_dt = dt.datetime(2020, 2, 5, 11, 0, 0)
+end_dt = dt.datetime(2020, 2, 5, 10, 0, 0)
 fig, _ = pyLARDA.Transformations.plot_timeheight(Vel_cor, rg_converter=False, title=True, range_interval=plot_range,
                                                  time_interval=[begin_dt, end_dt])
 name = f'{plot_path}/{begin_dt:%Y%m%d_%H%M}_{end_dt:%Y%m%d_%H%M}_{plot_range[1] / 1000:.0f}km_cloudnet_input'
@@ -92,6 +92,7 @@
 print(f'figure saved :: {fig_name}')
 
 cor_meas = moments['Vel_cor-Vel']
+cor_meas['var_lims'] = [-1, 1]
 # corrected MDV - measured MDV
 begin_dt = dt.datetime(2020, 2, 5, 0, 0, 5)
 end_dt = dt.datetime(2020, 2, 5, 23, 59, 55)
@@ -114,7 +115,8 @@
 print(f'figure saved :: {fig_name}')
 
 meas_cor = moments['Vel-Vel_cor']
-# corrected MDV - measured MDV
+meas_cor['var_lims'] = [-1, 1]
+# measured MDV - corrected MDV
 begin_dt = dt.datetime(2020, 2, 5, 0, 0, 5)
 end_dt = dt.datetime(2020, 2, 5, 23, 59, 55)
 fig, _ = pyLARDA.Transformations.plot_timeheight(meas_cor, rg_converter=False, title=True, range_interval=plot_range)
@@ -124,7 +126,7 @@
 plt.close()
 print(f'figure saved :: {fig_name}')
 
-# corrected MDV - measured MDV zoom
+# measured MDV - corrected MDV zoom
 begin_dt = dt.datetime(2020, 2, 5, 9, 0, 0)
 end_dt = dt.datetime(2020, 2, 5, 11, 0, 0)
 fig, _ = pyLARDA.Transformations.plot_timeheight(meas_cor, rg_converter=False, title=True, range_interval=plot_range,