Skip to content

Fix for mistake in useflags #437

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

Merged
rtodling merged 2 commits into from
Jun 8, 2026
Merged

Fix for mistake in useflags #437

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
9bdc750
Fix for mistake in useflags
rtodling May 21, 2026
5b50a28
missing atms_n21
rtodling Jun 8, 2026
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
22 changes: 12 additions & 10 deletions src/Applications/GEOSdas_App/jedi/jedi_useflags.csh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -87,15 +87,17 @@ foreach fn ( $lst )
set this = `echo $fn | cut -d. -f1`
set use = (`grep $this satinfo | cut -c26-30`)
/bin/cp $inpobs/$this.yaml $outobs/$this.yaml
mksi_flags.py --input $outobs/$this.yaml \
--output $outobs/$this.yaml \
--keyword ">>>use_channels_${this}<<<" \
--values $use
set cld = (`grep $this satinfo | cut -c72-77`)
if ( ! $status ) then
mksi_flags.py --input $outobs/$this.yaml \
--output $outobs/$this.yaml \
--keyword ">>>clddet_channels_${this}<<<" \
--values $cld
if ( "$use" != "" ) then
mksi_flags.py --input $outobs/$this.yaml \
--output $outobs/$this.yaml \
--keyword ">>>use_channels_${this}<<<" \
--values $use
set cld = (`grep $this satinfo | cut -c72-77`)
if ( ! $status ) then
mksi_flags.py --input $outobs/$this.yaml \
--output $outobs/$this.yaml \
--keyword ">>>clddet_channels_${this}<<<" \
--values $cld
endif
endif
end
1 change: 1 addition & 0 deletions src/Applications/GEOSdas_App/jedi/setup_aanajedi.pl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -407,6 +407,7 @@ sub init {
amsua_n15.yaml
amsua_n19.yaml
atms_n20.yaml
atms_n21.yaml
atms_npp.yaml
avhrr3_metop-b.yaml
avhrr3_metop-c.yaml
Expand Down
41 changes: 30 additions & 11 deletions src/Applications/GEOSdas_App/jedi/util/obsstats.csh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,8 +3,11 @@
setenv FVROOT /home/dao_ops/GEOSadas-5_43_0/GEOSadas/install-SLES15
set path = ( . $FVROOT/bin $path )

setenv DRYRUN #echo
setenv DRYRUN # echo
setenv OUTFIGS $TMP/SwellExperiments/Figs
#setenv XTRA "--xGSI --common"
setenv XTRA
setenv TOTRAIN 0

set humKX = (`echorc.x -rc ostats.rc -ncol 2 specific_humidity`)
set humNM = (`echorc.x -rc ostats.rc -ncol 1 specific_humidity`)
Expand All @@ -24,25 +27,29 @@ set ozNM = (`echorc.x -rc ostats.rc -ncol 1 ozone`)
set stwKX = (`echorc.x -rc ostats.rc -ncol 2 satwind`)
set stwNM = (`echorc.x -rc ostats.rc -ncol 1 satwind`)

set humKX = ()
set acTKX = ()
set radKX = ()
set gpsKX = ()
#set humKX = ()
#set acTKX = ()
#set radKX = ()
#set gpsKX = ()
#set ozKX = ()
set stwKX = ()
#set stwKX = ()

set expidGSI = x0054
set expidGSI = null
set expidGSI = x0053RPY
set expidJEDI = j4drpy
set expidJEDI = null
set expidJEDI = j4drp1

#set nymd = 20260113
set nymd = 20260120
set nhms = 000000
foreach nymd ( 20260115 )
foreach nhms ( 000000 060000 120000 180000 )

set ODSarch = $DAD/archive/544
set ODSarch = $DAD/archive/543
set IODAarch = $DAD/archive/JEDI/543

if ( ! -d $OUTFIGS ) mkdir -p $OUTFIGS

# GEOS-GSI experiment output
if ( $expidGSI != "null" ) then

Expand Down Expand Up @@ -141,11 +148,17 @@ if ( $expidGSI != "null" ) then
@ ic++
end
wait
if ( $TOTRAIN ) then
scp $OUTFIGS/$expidGSI.*.${nymd}_${hh}z.png train:/san_agcm/geos5/$expidGSI/obs/images
endif

endif # GEOS-GSI

# GEOS-JEDI experiment output

set radKX = (`echorc.x -rc ostats.rc -ncol 2 jedi_radiance`)
set radNM = (`echorc.x -rc ostats.rc -ncol 1 jedi_radiance`)

if ( $expidJEDI != "null" ) then

set jedi_date = (`tick $nymd $nhms -10800`)
Expand Down Expand Up @@ -214,8 +227,8 @@ if ( $expidJEDI != "null" ) then
set jhh = `echo $jnhms | cut -c1-2`
set hh = `echo $nhms | cut -c1-2`
set instr = $radNM[$ic]
$DRYRUN ~/src/python/JEDI/OBS/ioda_prs.binned.py --obtype radiance --satid $kx \

$DRYRUN ~/src/python/JEDI/OBS/ioda_prs.binned.py --obtype radiance --satid $kx $XTRA \
--fig $OUTFIGS/$expidJEDI.${instr}.${nymd}_${hh}z.png \
--tarname $IODAarch/$expidJEDI/jedi/obs/Y$jyyyy/M$jmm/$expidJEDI.jedi_hofx.${jnymd}_${jhh}z.tar \
${instr}.${jnymd}T${jhh}0000Z.nc4 &
Expand Down Expand Up @@ -261,4 +274,10 @@ if ( $expidJEDI != "null" ) then
end
wait

if ( $TOTRAIN ) then
scp $OUTFIGS/$expidJEDI.*.${nymd}_${hh}z.png train:/san_agcm/geos5/$expidJEDI/obs/images
endif
endif # GEOS-JEDI

end # nhms
end # nymd
21 changes: 21 additions & 0 deletions src/Applications/GEOSdas_App/jedi/util/ostats.rc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,6 +9,7 @@ radiance::
iasi_metop-c 877
atms_npp 900
atms_n20 920
atms_n21 921
ssmis_f17 717
avhrr_metop-b 626
avhrr_metop-c 627
Expand All @@ -18,6 +19,26 @@ radiance::
amsr2_gcom-w1 550
::

jedi_radiance::
amsua_metop-b 326
amsua_metop-c 327
mhs_metop-b 826
mhs_metop-c 827
airs_aqua 49
iasi_metop-b 876
iasi_metop-c 877
atms_npp 900
atms_n20 920
atms_n21 921
ssmis_f17 717
avhrr3_metop-b 626
avhrr3_metop-c 627
gmi_gpm 706
cris-fsr_n20 960
#cris-fsr_n21 961
amsr2_gcom-w1 550
::

specific_humidity::
sondeQ 120
::
Expand Down