From 565c8682e3c1aead0dfc7263e89eee4110336138 Mon Sep 17 00:00:00 2001 From: Rolf Reichle Date: Fri, 22 Dec 2023 13:56:29 -0500 Subject: [PATCH] fix bad no-data check (get_model_and_obs_clim_stats*.m) --- .../inputs/obs_scaling_params/get_model_and_obs_clim_stats.m | 5 ++--- .../get_model_and_obs_clim_stats_latlon_grid.m | 5 +++-- 2 files changed, 5 insertions(+), 5 deletions(-) diff --git a/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats.m b/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats.m index 9d0a7ebb..8d8b242e 100644 --- a/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats.m +++ b/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats.m @@ -419,10 +419,9 @@ ] = ... read_ObsFcstAna( fname ); - % remove tiles when there is no obs_fcst (obs_fcst == 0 in innov output when - % missing) + % remove tiles where obs_fcst is no-data (note: read_ObsFcstAna() returns NaN) - idx = find(obs_fcst == 0); + idx = isnan(obs_fcst); obs_assim( idx) = []; obs_species(idx) = []; diff --git a/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats_latlon_grid.m b/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats_latlon_grid.m index 0ff0e534..1d3a0c74 100644 --- a/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats_latlon_grid.m +++ b/src/Applications/LDAS_App/util/inputs/obs_scaling_params/get_model_and_obs_clim_stats_latlon_grid.m @@ -157,8 +157,9 @@ fclose(ifp); [date_time, obs_assim, obs_species, obs_tilenum, obs_lon, obs_lat, obs_obs, obs_obsvar, obs_fcst, obs_fcstvar, obs_ana, obs_anavar] = read_ObsFcstAna(fname); - % remove tiles when there is no obs_fcst (obs_fcst == 0 in innov output when missing) - idx = find(obs_fcst == 0); + % remove tiles where obs_fcst is no-data (note: read_ObsFcstAna() returns NaN) + + idx = isnan(obs_fcst); obs_assim( idx) = []; obs_species(idx) = [];