correction bug in process_traj.py

meteoswiss-mdr · Feb 23, 2017 · bab61f5 · bab61f5
1 parent 0d88c35
commit bab61f5
Show file tree

Hide file tree

Showing 2 changed files with 50 additions and 13 deletions.
diff --git a/config/processing/mals_pay_PHD_ASR_traj_prod.txt b/config/processing/mals_pay_PHD_ASR_traj_prod.txt
@@ -4,10 +4,11 @@
 
 # List of datasets to generate.
 # The detailed specification of each dataset is given below.
-dataSetList STRARR 4   
+dataSetList STRARR 5
     reflectivity_traj_lb
     rainrate_traj_lb
     windshear_traj_lb    
+    temperature_traj
     temperature_traj_lb
 
 
@@ -124,9 +125,9 @@ reflectivity_traj_lb STRUCT 14
    products      STRUCT 1
       LOWBEAM_ANTENNA STRUCT 1
          type  STRING   PLOT_AND_WRITE
-#         ymin  FLOAT    0.0
-#         ymax  FLOAT    4.0
-#
+         ymin  FLOAT    0.0
+         ymax  FLOAT    4.0
+
 #rainrate_traj STRUCT 3
 #   type          STRING TRAJ_ATPLANE
 #   datatype      STRING PROC:RR
@@ -135,7 +136,7 @@ reflectivity_traj_lb STRUCT 14
 #         type  STRING   PLOT_AND_WRITE
 #         ymin  FLOAT    0.0
 #         ymax  FLOAT    4.0
-#
+
 rainrate_traj_lb STRUCT 14
    type          STRING TRAJ_ANTENNA_PATTERN
    datatype      STRING CFRADIAL:RR,RR_hydro,SAVEVOL_RR_hydro
@@ -162,7 +163,7 @@ rainrate_traj_lb STRUCT 14
 #   products STRUCT 1
 #      AT_PLANE STRUCT 1
 #         type  STRING PLOT_AND_WRITE
-#
+
 windshear_traj_lb STRUCT 14
    type     STRING TRAJ_ANTENNA_PATTERN
    datatype STRING CFRADIAL:windshear_v,windshear_v,SAVEVOL_windshear_v
@@ -183,12 +184,12 @@ windshear_traj_lb STRUCT 14
          ymin  FLOAT    -20.0
          ymax  FLOAT     20.0
 
-# temperature_traj STRUCT 3
-#    type     STRING TRAJ_ATPLANE
-#    datatype STRING RAD4ALPCOSMO:TEMP
-#    products STRUCT 1
-#       AT_PLANE STRUCT 1
-#          type  STRING PLOT_AND_WRITE
+temperature_traj STRUCT 3
+    type     STRING TRAJ_ATPLANE
+    datatype STRING RAD4ALPCOSMO:TEMP
+    products STRUCT 1
+       AT_PLANE STRUCT 1
+          type  STRING PLOT_AND_WRITE
 
 temperature_traj_lb STRUCT 14
    type     STRING TRAJ_ANTENNA_PATTERN
@@ -209,3 +210,39 @@ temperature_traj_lb STRUCT 14
          type  STRING   PLOT_AND_WRITE
 #         ymin  FLOAT    -20.0
 #         ymax  FLOAT     20.0
+
+#dBZ_traj STRUCT 3
+#    type     STRING TRAJ_ATPLANE
+#    datatype STRING RAD4ALP:dBZ
+#    products STRUCT 1
+#       AT_PLANE STRUCT 1
+#          type  STRING PLOT_AND_WRITE
+#          
+#dBZv_traj STRUCT 3
+#    type     STRING TRAJ_ATPLANE
+#    datatype STRING RAD4ALP:dBZv
+#    products STRUCT 1
+#       AT_PLANE STRUCT 1
+#          type  STRING PLOT_AND_WRITE
+#          
+#ZDR_traj STRUCT 3
+#    type     STRING TRAJ_ATPLANE
+#    datatype STRING RAD4ALP:ZDR
+#    products STRUCT 1
+#       AT_PLANE STRUCT 1
+#          type  STRING PLOT_AND_WRITE
+#          
+#uRhoHV_traj STRUCT 3
+#    type     STRING TRAJ_ATPLANE
+#    datatype STRING RAD4ALP:uRhoHV
+#    products STRUCT 1
+#       AT_PLANE STRUCT 1
+#          type  STRING PLOT_AND_WRITE
+#          
+#uPhiDP_traj STRUCT 3
+#    type     STRING TRAJ_ATPLANE
+#    datatype STRING RAD4ALP:uPhiDP
+#    products STRUCT 1
+#       AT_PLANE STRUCT 1
+#          type  STRING PLOT_AND_WRITE
+
diff --git a/src/pyrad_proc/pyrad/proc/process_traj.py b/src/pyrad_proc/pyrad/proc/process_traj.py
@@ -684,7 +684,7 @@ def process_traj_antenna_pattern(procstatus, dscfg, radar_list=None,
 
             # flag regions with colocated usable data in radar_sel
             gate_coloc_radar_sel = intersection(
-                radar_sel, r_radar, h_tol=alt_tol, latlon_tol=latlon_tol
+                radar_sel, r_radar, h_tol=alt_tol, latlon_tol=latlon_tol,
                 vol_d_tol=None, vismin=None, hmin=None, hmax=max_altitude,
                 rmin=None, rmax=None, elmin=None, elmax=None, azmin=None,
                 azmax=None, visib_field=None,