diff --git a/fix b/fix index 3c983506e..a3142927a 160000 --- a/fix +++ b/fix @@ -1 +1 @@ -Subproject commit 3c983506eb4973cafb4a73e00e09d35b75971dba +Subproject commit a3142927a7504994bf8f90bf11b8781e046b13dc diff --git a/src/gsi/read_obs.F90 b/src/gsi/read_obs.F90 index 8312a7bdd..4821b721f 100644 --- a/src/gsi/read_obs.F90 +++ b/src/gsi/read_obs.F90 @@ -408,7 +408,8 @@ subroutine read_obs_check (lexist,filename,jsatid,dtype,minuse,nread) (said == 42) .or. (said == 43) .or. (said == 722) .or. & (said == 723).or. (said == 265).or. (said == 266) .or. & (said == 267).or. (said == 268).or. (said == 269) .or. & - (said == 803).or. (said == 804).or. (said == 66)) then + (said == 803).or. (said == 804).or. (said == 66) .or. & + (said == 768)) then lexist=.true. exit gpsloop end if diff --git a/src/gsi/setupbend.f90 b/src/gsi/setupbend.f90 index c4e382095..87eb2628f 100644 --- a/src/gsi/setupbend.f90 +++ b/src/gsi/setupbend.f90 @@ -111,6 +111,7 @@ subroutine setupbend(obsLL,odiagLL, & ! 2024-12-04 Li - remove the QC check for rejecting MetOp data <8 km ! 2024-12-04 Li - add GRACE-FO (803&804) data ! 2024-12-04 Li - add new obs error model by Chris Riedel +! 2025-12-08 Li - add PlanetiQ YAM-8 (768) data ! ! input argument list: ! lunin - unit from which to read observations @@ -653,7 +654,7 @@ subroutine setupbend(obsLL,odiagLL, & rdiagbuf(18,i) = trefges ! temperature at obs location (Kelvin) if monotone grid rdiagbuf(21,i) = qrefges ! specific humidity at obs location (kg/kg) if monotone grid commdat=.false. - if (data(isatid,i)>=265 .and. data(isatid,i)<=269) commdat=.true. + if ( (data(isatid,i)>=265 .and. data(isatid,i)<=269) .or. (data(isatid,i)==768) ) commdat=.true. if (.not. qcfail(i)) then ! not SR ratio_errors(i) = data(ier,i) @@ -855,7 +856,7 @@ subroutine setupbend(obsLL,odiagLL, & end do loopoverobs1 ! end of loop over observations - write(6,'("setupbend: Number of obs considered and accepted " 2I10)') nobs, count(mask=muse .neqv. .false.) +! write(6,'("setupbend: Number of obs considered and accepted " 2I10)') nobs, count(mask=muse .neqv. .false.) if (nobs_out>=1) then write(6,*)'WARNING GPSRO:',nobs_out,'obs outside integration grid. Increase nsig_ext to',&