diff --git a/parm/metplus_config/stats/global_det/atmos_grid2grid/GridStat_fcstGLOBAL_DET_obsGHRSST_OSPO_DailyAvg.conf b/parm/metplus_config/stats/global_det/atmos_grid2grid/GridStat_fcstGLOBAL_DET_obsGHRSST_OSPO_DailyAvg.conf index 101881b5c9..7a19b29166 100644 --- a/parm/metplus_config/stats/global_det/atmos_grid2grid/GridStat_fcstGLOBAL_DET_obsGHRSST_OSPO_DailyAvg.conf +++ b/parm/metplus_config/stats/global_det/atmos_grid2grid/GridStat_fcstGLOBAL_DET_obsGHRSST_OSPO_DailyAvg.conf @@ -31,9 +31,10 @@ LEAD_SEQ = {ENV[fhr_list]} #### Variables FCST_VAR1_NAME = FCST_TMP_Z0_DAILYAVG FCST_VAR1_LEVELS = "(*,*)" -FCST_VAR1_OPTIONS = set_attr_name = "SST_DAILYAVG"; +FCST_VAR1_OPTIONS = set_attr_name = "SST_DAILYAVG"; cnt_thresh = [ >= 271.15 ]; cnt_logic = INTERSECTION; OBS_VAR1_NAME = analysed_sst OBS_VAR1_LEVELS = "(0,*,*)" +OBS_VAR1_OPTIONS = cnt_thresh = [ >= 271.15 ]; cnt_logic = INTERSECTION; #### GridStat GRID_STAT_CONFIG_FILE = {PARM_BASE}/met_config/GridStatConfig_wrapped GRID_STAT_ONCE_PER_FIELD = False diff --git a/ush/global_det/global_det_atmos_plots_grid2grid_create_job_scripts.py b/ush/global_det/global_det_atmos_plots_grid2grid_create_job_scripts.py index be7654f308..ee7c7d62df 100644 --- a/ush/global_det/global_det_atmos_plots_grid2grid_create_job_scripts.py +++ b/ush/global_det/global_det_atmos_plots_grid2grid_create_job_scripts.py @@ -464,8 +464,8 @@ #### sst for sst_job in list(filter_stats_jobs_dict['sst'].keys()): filter_stats_jobs_dict['sst'][sst_job]['grid'] = 'G004' - filter_stats_jobs_dict['sst'][sst_job]['fcst_var_dict']['threshs'] = ['NA'] - filter_stats_jobs_dict['sst'][sst_job]['obs_var_dict']['threshs'] = ['NA'] + filter_stats_jobs_dict['sst'][sst_job]['fcst_var_dict']['threshs'] = ['ge271.15&&'] + filter_stats_jobs_dict['sst'][sst_job]['obs_var_dict']['threshs'] = ['ge271.15'] filter_stats_jobs_dict['sst'][sst_job]['interps'] = ['NEAREST/1'] if JOB_GROUP == 'filter_stats': JOB_GROUP_dict = filter_stats_jobs_dict diff --git a/ush/global_det/global_det_atmos_plots_specs.py b/ush/global_det/global_det_atmos_plots_specs.py index 7c9872d8cb..153d45ca9a 100644 --- a/ush/global_det/global_det_atmos_plots_specs.py +++ b/ush/global_det/global_det_atmos_plots_specs.py @@ -683,29 +683,13 @@ def get_plot_title(self, plot_info_dict, date_info_dict, units): and units == '10^6_km^2': units = 'x'+units.replace('_', ' ') plot_title = plot_title+' '+'('+units+')' - if var_thresh_for_title != 'NA': + if var_thresh_for_title != 'NA'\ + and plot_info_dict['fcst_var_name'] not in ['SST_DAILYAVG']: plot_title = plot_title+', '+var_thresh_for_title+' '+units - thresh_value = float(plot_info_dict['fcst_var_thresh'][2:]) if plot_info_dict['fcst_var_name'] == 'APCP': + thresh_value = float(plot_info_dict['fcst_var_thresh'][2:]) thresh_in = round(thresh_value*0.0393701, 3) plot_title = plot_title+' ('+str(thresh_in)+' in)' - elif plot_info_dict['fcst_var_name'] in ['SNOD_A24', 'WEASD_A24']: - thresh_in = round(thresh_value*39.3701,3) - plot_title = plot_title+' ('+str(thresh_in)+' in)' - elif plot_info_dict['fcst_var_name'] == 'DPT': - thresh_F = round((((thresh_value-273.15)*9)/5)+32) - plot_title = plot_title+' ('+str(thresh_F)+' F)' - elif plot_info_dict['fcst_var_name'] == 'HGT' \ - and plot_info_dict['fcst_var_level'] == 'CEILING': - thresh_kft = round(thresh_value/304.8,1) - if int(thresh_kft) == thresh_kft: - thresh_kft = int(thresh_kft) - plot_title = plot_title+' ('+str(thresh_kft)+' kft)' - elif plot_info_dict['fcst_var_name'] == 'VIS': - thresh_mile = round(thresh_value * 0.000621371,1) - if int(thresh_mile) == thresh_mile: - thresh_mile = int(thresh_mile) - plot_title = plot_title+' ('+str(thresh_mile)+' mile)' if plot_info_dict['interp_method'] == 'NBRHD_SQUARE': plot_title = (plot_title+' ' +'Neighborhood Points: ' @@ -747,9 +731,10 @@ def get_savefig_name(self, image_dir, plot_info_dict, date_info_dict): thresh_symbol, thresh_letter = gda_util.format_thresh( plot_info_dict['fcst_var_thresh'] ) - metric_savefig_name = ( - metric_savefig_name+'_' - +thresh_letter.replace('.','p') + if plot_info_dict['fcst_var_name'] not in ['SST_DAILYAVG']: + metric_savefig_name = ( + metric_savefig_name+'_' + +thresh_letter.replace('.','p') ) parameter_savefig_name = plot_info_dict['fcst_var_name'] if plot_info_dict['fcst_var_name'] == 'HGT_DECOMP':