Skip to content

Commit

Permalink
fix bad no-data check (get_model_and_obs_clim_stats*.m)
Browse files Browse the repository at this point in the history
  • Loading branch information
gmao-rreichle committed Dec 22, 2023
1 parent 4113033 commit 565c868
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 5 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -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) = [];
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -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) = [];
Expand Down

0 comments on commit 565c868

Please sign in to comment.