Skip to content

Commit

Permalink
Update the QC for the enhanced high-resolution GOES-R mesoscale float…
Browse files Browse the repository at this point in the history
…er AMVs (#724)

**Description**
This is an update of the QC process for the enhanced high-resolution
GOES-R mesoscale floater AMVs to resolve the issue #713.
As mentioned in the issue, enhanced AMV data are derived from IR band,
which turns out to be not reliable in the mid-layer of atmosphere.
Adding additional QC process to remove data in those layers are
necessary and proved to enhance the Hurricane intensity forecast.

**Type of change**

Please delete options that are not relevant.

- [ ] Bug fix (non-breaking change which fixes an issue)
- [x] New feature (non-breaking change which adds functionality)
- [ ] Breaking change (fix or feature that would cause existing
functionality to not work as expected)
- [ ] This change requires a documentation update

**How Has This Been Tested?**
This change has been tested with GSI regression test on Hera.
  
**Checklist**

- [x] My code follows the style guidelines of this project
- [x] I have performed a self-review of my own code
- [x] I have commented my code, particularly in hard-to-understand areas
- [x] New and existing tests pass with my changes
- [x] Any dependent changes have been merged and published
  • Loading branch information
JingCheng-NOAA authored Apr 1, 2024
1 parent b2fc4fd commit db477e3
Showing 1 changed file with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion src/gsi/setupw.f90
Original file line number Diff line number Diff line change
Expand Up @@ -945,7 +945,11 @@ subroutine setupw(obsLL,odiagLL,lunin,mype,bwork,awork,nele,nobs,is,conv_diagsav
else
if( presw >950.0_r_kind) error =zero ! screen data beloww 950mb
endif
if(itype ==242 .or. itype ==243 ) then ! visible winds from JMA and EUMETSAT
if(itype == 241 ) then
if( presw >399.0_r_kind .and. presw <601.0_r_kind) then !CIMISS(enhanced AMV) winds
error=zero ! no data between400-600mb
endif
else if(itype ==242 .or. itype ==243 ) then ! visible winds from JMA and EUMETSAT
if(presw <700.0_r_kind) error=zero ! no visible winds above 700mb
else if(itype ==245 ) then
if( presw >399.0_r_kind .and. presw <801.0_r_kind) then !GOES IR winds
Expand Down

0 comments on commit db477e3

Please sign in to comment.