Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

add WRFDA satwnd.bufr ingest #1617

Merged
merged 1 commit into from
Dec 24, 2021

Conversation

jamiebresch
Copy link
Contributor

@jamiebresch jamiebresch commented Dec 23, 2021

TYPE: new feature

KEYWORDS: WRFDA, satwnd.bufr

SOURCE: Jamie Bresch (NCAR)

DESCRIPTION OF CHANGES:

When ob_format=1 (bufr), WRFDA reads NCEP's prepbufr file and processes satellite winds (AMV) contained in prepbufr
as ob type geoamv. This PR adds the capability to read a separate satwnd bufr file and decode additional AMVs (as
ob type polaramv) that are not included in prepbufr.

The use_satwnd_bufr default is true, when ob_format=1 and satwnd.bufr does not exist in the working directory.
There is no impact except for a warning message of "satwnd.bufr does not exist".

Data sources:
https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gdas.YYYYMMDD/HH/atmos/gdas.tHHz.satwnd.tm00.bufr_d
cheyenne:/gpfs/fs1/collections/rda/data/ds351.0/bufr/YYYY/gdas.satwnd.tHHz.YYYYMMDD.bufr

Thinning for satwnd.bufr (done in da_read_obs_bufr_satwnd.inc) is implemented differently than all other ob types.
Superobbing is available when thin_conv_opt (polaramv) is set to 3 (2-D thinning boxes) or 4 (3-D thinning boxes).

The reading of an external obs error table is moved one level up, so the table values can be used by both
da_read_obs_bufr.inc and da_read_obs_bufr_satwnd.inc

Piggybacking on top of this PR are some other README.namelist update to fix a few typos pointed out by Yi-Chuan Lo
(Central Weather Bureau, Taiwan) and adding a descrption for recently added gpsref options.

LIST OF MODIFIED FILES:
M Registry/registry.var
M var/README.namelist
M var/build/depend.txt
M var/da/da_control/da_control.f90
M var/da/da_obs_io/da_obs_io.f90
M var/da/da_obs_io/da_read_obs_bufr.inc
A var/da/da_obs_io/da_read_obs_bufr_satwnd.inc
M var/da/da_radiance/gsi_thinning.f90
M var/da/da_setup_structures/da_setup_obs_structures.inc
M var/da/da_setup_structures/da_setup_obs_structures_bufr.inc
M var/da/da_setup_structures/da_setup_structures.f90

TESTS CONDUCTED:

  1. A few tests were conducted to check that this version of the code runs.
  2. Jenkins tests are all passing.

RELEASE NOTE: WRFDA can ingest NCEP's satwnd.bufr (gdas.satwnd.tHHz.YYYYMMDD.bufr) to assimilate more AMVs that are not included in NCEP's prepbufr files. More flexible thinning/superobbing options (thin_conv_opt) and obs error specification (uv_error_opt) are implemented for using satwnd.bufr. See WRFDA/var/README.namelist for details.

Copy link
Contributor

@liujake liujake left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Great to have the capability for ingesting satwnd.bufr.

@liujake liujake requested a review from davegill December 23, 2021 18:23
@davegill davegill merged commit 1436973 into wrf-model:develop Dec 24, 2021
@jamiebresch jamiebresch deleted the satwnd_bufr_ingest branch January 13, 2022 14:29
vlakshmanan-scala pushed a commit to scala-computing/WRF that referenced this pull request Apr 4, 2024
TYPE: new feature

KEYWORDS: WRFDA, satwnd.bufr

SOURCE: Jamie Bresch (NCAR)

DESCRIPTION OF CHANGES:

When ob_format=1 (bufr), WRFDA reads NCEP's prepbufr file and processes satellite winds (AMV) contained in prepbufr 
as ob type geoamv. This PR adds the capability to read a separate satwnd bufr file and decode additional AMVs (as 
ob type polaramv) that are not included in prepbufr.

The `use_satwnd_bufr` default is true, when `ob_format=1` and `satwnd.bufr` does not exist in the working directory.
There is no impact except for a warning message of "satwnd.bufr does not exist".

Data sources:
https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gdas.YYYYMMDD/HH/atmos/gdas.tHHz.satwnd.tm00.bufr_d
cheyenne:/gpfs/fs1/collections/rda/data/ds351.0/bufr/YYYY/gdas.satwnd.tHHz.YYYYMMDD.bufr

Thinning for satwnd.bufr (done in da_read_obs_bufr_satwnd.inc) is implemented differently than all other ob types.
Superobbing is available when `thin_conv_opt` (polaramv) is set to 3 (2-D thinning boxes) or 4 (3-D thinning boxes).

The reading of an external obs error table is moved one level up, so the table values can be used by both 
da_read_obs_bufr.inc and da_read_obs_bufr_satwnd.inc

Piggybacking on top of this PR are some other README.namelist update to fix a few typos pointed out by Yi-Chuan Lo 
(Central Weather Bureau, Taiwan) and adding a descrption for recently added gpsref options.

LIST OF MODIFIED FILES:
M       Registry/registry.var
M       var/README.namelist
M       var/build/depend.txt
M       var/da/da_control/da_control.f90
M       var/da/da_obs_io/da_obs_io.f90
M       var/da/da_obs_io/da_read_obs_bufr.inc
A       var/da/da_obs_io/da_read_obs_bufr_satwnd.inc
M       var/da/da_radiance/gsi_thinning.f90
M       var/da/da_setup_structures/da_setup_obs_structures.inc
M       var/da/da_setup_structures/da_setup_obs_structures_bufr.inc
M       var/da/da_setup_structures/da_setup_structures.f90

TESTS CONDUCTED:
1. A few tests were conducted to check that this version of the code runs.
2. Jenkins tests are all passing.

RELEASE NOTE: WRFDA can ingest NCEP's satwnd.bufr (gdas.satwnd.tHHz.YYYYMMDD.bufr) to assimilate more AMVs that are not included in NCEP's prepbufr files. More flexible thinning/superobbing options (thin_conv_opt) and obs error specification (uv_error_opt) are implemented for using satwnd.bufr. See WRFDA/var/README.namelist for details.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants