Merge branch 'develop' into feature/msienkie/CrIS-SNPP_LWfail

GEOS-ESM · Sep 11, 2023 · e5dfe2c · e5dfe2c
2 parents 0b65faf + e88ffae
commit e5dfe2c
Show file tree

Hide file tree

Showing 6 changed files with 83 additions and 65 deletions.
diff --git a/GEOSaana_GridComp/GSI_GridComp/etc/gmao_global_tlapmean.txt b/GEOSaana_GridComp/GSI_GridComp/etc/gmao_global_tlapmean.txt
@@ -3039,6 +3039,9 @@
  avhrr3_n18               3   0.884430E-02
  avhrr3_n18               4   0.115395E-01
  avhrr3_n18               5   0.161449E-01
+ avhrr3_n19               3   0.884430E-02
+ avhrr3_n19               4   0.115395E-01
+ avhrr3_n19               5   0.161449E-01
  avhrr3_metop-b           3   0.870211E-02
  avhrr3_metop-b           4   0.101765E-01
  avhrr3_metop-b           5   0.136881E-01

diff --git a/GEOSaana_GridComp/GSI_GridComp/read_prepbufr.f90 b/GEOSaana_GridComp/GSI_GridComp/read_prepbufr.f90
@@ -143,6 +143,7 @@ subroutine read_prepbufr(nread,ndata,nodata,infile,obstype,lunout,twindin,sis,&
 !   2020-01-29  Sienkiewicz - allow obstypes marked as passive in convinfo to be thinned
 !   2020-10-27  sienkiewicz - update for BUFR drifting buoys, T29=564
 !   2022-09-21  Sienkiewicz - Add BUFR ship subtypes (524, 525) to ship definition
+!   2023-05-22  Sienkiewicz - fix for ship obs with zero obs height where POB different from PMO
 !
 
 !   input argument list:
@@ -2268,6 +2269,18 @@ subroutine read_prepbufr(nread,ndata,nodata,infile,obstype,lunout,twindin,sis,&
 
 !             Surface pressure
               else if(psob) then
+!
+! Substitute PMO for POB for ship obs if surface height is zero and POB and PMO disagree
+! (obs needs actual barometer height but that's not available)
+                 if ( kx == 180 .and. obsdat(4,k) == zero ) then
+                    it29=nint(hdr(8))
+                    if (it29 >= 522 .and. it29 <= 525 .and. .not. ibfms(obsdat(13,k))) then
+                       if (obsdat(13,k) - obsdat(1,k) > one_tenth) then
+                          plevs(k)=one_tenth*obsdat(13,k)
+                          dlnpob=log(plevs(k))  ! ln(pressure in cb)
+                       endif
+                    endif
+                 endif
 
                  poe=obserr(1,k)*one_tenth                  ! convert from mb to cb
                  if (inflate_error) poe=poe*r1_2