Skip to content

Commit

Permalink
Bugfix in EnKF assim of conventional q obs (NOAA-EMC#545)
Browse files Browse the repository at this point in the history
**Description**

Added code to divide hx_modens by qsat, for consistency with treatment
of q as q/qsat.

Fixes issue NOAA-EMC#544
  • Loading branch information
ClaraDraper-NOAA authored and TingLei-daprediction committed May 23, 2023
1 parent 87ff05e commit 54c2f52
Showing 1 changed file with 3 additions and 0 deletions.
3 changes: 3 additions & 0 deletions src/enkf/readconvobs.f90
Original file line number Diff line number Diff line change
Expand Up @@ -789,6 +789,9 @@ subroutine get_convobs_data_nc(obspath, datestring, nobs_max, nobs_maxdiag, &
x_obs(nob) = x_obs(nob) /Forecast_Saturation_Spec_Hum(i)
hx_mean(nob) = hx_mean(nob) /Forecast_Saturation_Spec_Hum(i)
hx_mean_nobc(nob) = hx_mean_nobc(nob) /Forecast_Saturation_Spec_Hum(i)
if (neigv>0) then
hx_modens(:,nob) = hx_modens(:,nob)/ Forecast_Saturation_Spec_Hum(i)
endif
endif

! for wind, also read v-component
Expand Down

0 comments on commit 54c2f52

Please sign in to comment.