diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 07025cb2c42a05f9cb606df95f395151adfa7414..9e17a93708aa6d99c9028b88cb6a4388b7b3634f 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1232,7 +1232,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd opd_thin_threshold = 0.1 closeness_top = 100.0 # meters - max_cld_depth = 1500.0 + max_cld_depth = 2000.0 max_altitude = 5000.0 max_cld_altitude = 5000.0 @@ -1257,8 +1257,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd # if not (is_day(solzen[i,]) or is_night(solzen[i,])): # continue - if not (icing_alt[i] < max_altitude): - continue + # if not (icing_alt[i] < max_altitude): + # continue keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) @@ -1285,7 +1285,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd cld_tmp = cld_top_temp[i, ].flatten() med_cld_tmp = np.median(cld_tmp[keep]) - if med_cld_tmp > 274.5: + if med_cld_tmp > 275.0: continue if icing_alt[i] < med_cld_hgt - max_cld_depth: @@ -1333,8 +1333,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ # if not (is_day(solzen[i,]) or is_night(solzen[i,])): # continue - if not (icing_alt[i] < max_altitude): - continue + # if not (icing_alt[i] < max_altitude): + # continue keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) @@ -1361,7 +1361,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ cld_tmp = cld_top_temp[i, ].flatten() med_cld_tmp = np.median(cld_tmp[keep]) - if med_cld_tmp > 274.5: + if med_cld_tmp > 275.0: continue if icing_alt[i] < med_cld_hgt - max_cld_depth: