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.

Sign up for GitHub

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

GSI changes to enable the use of historical or reprocessed SATWND and other observations for Reanalysis #837

Merged
merged 19 commits into from
Feb 13, 2025
+337 −541
Merged

GSI changes to enable the use of historical or reprocessed SATWND and other observations for Reanalysis #837

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
19 commits
Select commit Hold shift + click to select a range
5b0cfde
add index for ompsnmnc in TypeManager and fine tune do_qc flag in rea…
jack-woollen Dec 12, 2023
66a8cd4
add satwnd types 005041-42-43 and 005001-002-003 to the read_obs satw…
jack-woollen Dec 14, 2023
7b87b4a
look for CMCM in older satwnd data
jack-woollen Dec 15, 2023
c3c1249
ra changes for bc spinup and lza estimation where need be
jack-woollen Jan 9, 2024
a591ca1
add ompsnpnc ozone type
jack-woollen Feb 16, 2024
d4d79f9
add index for ompsnmnc in TypeManager and fine tune do_qc flag in rea…
jack-woollen Dec 12, 2023
4e53bbf
add satwnd types 005041-42-43 and 005001-002-003 to the read_obs satw…
jack-woollen Dec 14, 2023
dc5dd02
look for CMCM in older satwnd data
jack-woollen Dec 15, 2023
4431f3e
ra changes for bc spinup and lza estimation where need be
jack-woollen Jan 9, 2024
1170fae
add ompsnpnc ozone type
jack-woollen Feb 16, 2024
b5d0c0d
merging GSI develop into fork
May 30, 2024
6b474b4
J.Derber fix for indxRF array bound in satellite processing
jack-woollen Jul 5, 2024
8c20b7c
refactored read_satwnd.f90 to include historical satellites for RA an…
Dec 25, 2024
023ff01
Merge branch 'NOAA-EMC:develop' into develop
jack-woollen Dec 25, 2024
87ec85a
check for CMCM for each time through the data
Dec 30, 2024
3ca5b9e
Merge branch 'NOAA-EMC:develop' into develop
jack-woollen Jan 30, 2025
4f98338
Merge branch 'NOAA-EMC:develop' into develop
jack-woollen Feb 10, 2025
98debf7
add docblock for sattabin
Feb 11, 2025
316f98d
some coding cleanups
Feb 13, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions src/gsi/gsi_obOperTypeManager.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -289,6 +289,8 @@ function dtype2index_(dtype) result(index_)
case("omieff" ); index_= iobOper_oz
case("tomseff" ); index_= iobOper_oz
case("ompsnmeff"); index_= iobOper_oz
case("ompsnmnc" ); index_= iobOper_oz
case("ompsnpnc" ); index_= iobOper_oz

case("o3l" ,"[o3loper]" ); index_= iobOper_o3l
case("o3lev" ); index_= iobOper_o3l
Expand Down
7 changes: 6 additions & 1 deletion src/gsi/radinfo.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1254,6 +1254,7 @@ subroutine radinfo_write(pe_out)
! !USES:

use mpimod, only: mype
use constants, only: r10000
implicit none

integer(i_kind),optional, intent(in) :: pe_out
Expand All @@ -1274,7 +1275,11 @@ subroutine radinfo_write(pe_out)
rewind lunout
do jch=1,jpch_rad
do i=1,npred
varx(i)=varA(i,jch)
if (inew_rad(jch) .or. abs(ostats(jch)) .le. tiny(ostats(1))) then
varx(i) = r10000
else
varx(i)=varA(i,jch)
endif
end do
write(lunout,'(I5,1x,A20,1x,I5,e15.7/2(4x,10e15.7/))') jch,nusis(jch),&
nuchan(jch),ostats(jch),(varx(ip),ip=1,npred)
Expand Down
2 changes: 1 addition & 1 deletion src/gsi/read_bufrtovs.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -701,7 +701,7 @@ subroutine read_bufrtovs(mype,val_tovs,ithin,isfcalc,&

panglr=(start+real(ifovmod-1,r_kind)*step)*deg2rad
lzaest = asin(rato*sin(panglr))
if( msu .or. hirs2 .or. ssu)then
if( msu .or. hirs2 .or. hirs3 .or. ssu .or. bfr2bhdr(1) > 1.e5)then
lza = lzaest
else
lza = bfr2bhdr(1)*deg2rad ! local zenith angle
Expand Down
2 changes: 2 additions & 0 deletions src/gsi/read_obs.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -459,6 +459,8 @@ subroutine read_obs_check (lexist,filename,jsatid,dtype,minuse,nread)
trim(subset) == 'NC005090' .or. trim(subset) == 'NC005091' .or.&
trim(subset) == 'NC005067' .or. trim(subset) == 'NC005068' .or. trim(subset) == 'NC005069' .or.&
trim(subset) == 'NC005047' .or. trim(subset) == 'NC005048' .or. trim(subset) == 'NC005049' .or.&
trim(subset) == 'NC005041' .or. trim(subset) == 'NC005042' .or. trim(subset) == 'NC005043' .or.&
trim(subset) == 'NC005001' .or. trim(subset) == 'NC005002' .or. trim(subset) == 'NC005003' .or.&
trim(subset) == 'NC005081' .or. &
trim(subset) == 'NC005072' ) then
lexist = .true.
Expand Down
Loading