Skip to content

make usenewgfsberror configurable; add l_obsprvdiag #369

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
guoqing-noaa merged 2 commits into from
May 3, 2022
Merged

make usenewgfsberror configurable; add l_obsprvdiag #369

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
b1969df
make usenewgfsberror configurable; add l_obsprvdiag
guoqing-noaa May 3, 2022
200c980
add anavinfo.rtma
guoqing-noaa May 3, 2022
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
60 changes: 60 additions & 0 deletions fix/gsi/anavinfo.rtma
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,60 @@
met_guess::
!var level crtm_use desc orig_name
ps 1 -1 surface_pressure ps
z 1 -1 geopotential_height phis
u 65 2 zonal_wind u
v 65 2 meridional_wind v
tv 65 2 virtual_temperature tv
q 65 2 specific_humidity sphu
oz 65 2 ozone ozone
delp 65 -1 fv3 del pressure delp
t2m 1 -1 2-m_T t2m
q2m 1 -1 2-m_Q q2m
::

state_derivatives::
!var level src
ps 1 met_guess
u 65 met_guess
v 65 met_guess
tv 65 met_guess
q 65 met_guess
oz 65 met_guess
cw 65 met_guess
prse 66 met_guess
::

state_tendencies::
!var levels source
u 65 met_guess
v 65 met_guess
tv 65 met_guess
q 65 met_guess
oz 65 met_guess
prse 66 met_guess
::

state_vector::
!var level itracer source funcof
u 65 0 met_guess u
v 65 0 met_guess v
tv 65 0 met_guess tv
tsen 65 0 met_guess tv,q
q 65 1 met_guess q
oz 65 1 met_guess oz
prse 66 0 met_guess prse
ps 1 0 met_guess prse
sst 1 0 met_guess sst
::

control_vector::
!var level itracer as/tsfc_sdv an_amp0 source funcof
sf 65 0 1.00 -1.0 state u,v
vp 65 0 1.00 -1.0 state u,v
ps 1 0 0.50 -1.0 state prse
t 65 0 0.70 -1.0 state tv
q 65 1 0.70 -1.0 state q
sst 1 0 1.00 -1.0 state sst
stl 1 0 1.00 -1.0 motley sst
sti 1 0 1.00 -1.0 motley sst
::
3 changes: 2 additions & 1 deletion fix/gsi/gsiparm.anl.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ gsi_namelist="
l4densvar=.false.,nhr_obsbin=3,
use_gfs_nemsio=.false.,use_gfs_ncio=.true.,reset_bad_radbc=.true.,
netcdf_diag=.false.,
l_obsprvdiag=${l_obsprvdiag},
/
&GRIDOPTS
fv3_regional=.true.,grid_ratio_fv3_regional=${grid_ratio_fv3},nvege_type=20,
Expand All @@ -25,7 +26,7 @@ gsi_namelist="
vs=${bkgerr_vs},
hzscl=${bkgerr_hzscl},
bw=0.,fstat=.true.,
usenewgfsberror=.true.,
usenewgfsberror=${usenewgfsberror},
/
&ANBKGERR
anisotropic=.false.,
Expand Down
18 changes: 11 additions & 7 deletions ush/config.sh.RTMA_CONUS_3km
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@ VERBOSE="TRUE"
RUN_ENVIR="nco"

PREDEF_GRID_NAME=RRFS_CONUS_3km
ADDNL_OUTPUT_GRIDS=( )
ADDNL_OUTPUT_GRIDS=()

TILE_LABELS="CONUS REGIONS"
TILE_SETS="full NE,NC,NW,SE,SC,SW"
Expand All @@ -59,6 +59,7 @@ RADARREFL_TIMELEVEL=(0)
#FH_DFI_RADAR="0.0,0.25,0.5"
#DO_SOIL_ADJUST="TRUE"
#DO_RADDA="FALSE"
#DO_BUFRSND="FALSE"
#
#-------------------------------------------------------------------------------------
# GSI Namelist parameters configurable across differnt applications
Expand All @@ -71,11 +72,13 @@ RADARREFL_TIMELEVEL=(0)
#-------------------------------------------------------------------------------------
#
# &SETUP and &BKGERR
l_obsprvdiag=.true.
diag_radardbz=.false.
write_diag_2=.true.
bkgerr_vs=0.125 #1.0 is default setting
#bkgerr_hzscl=0.373,0.746,1.5 #default setting, no trailing ,
bkgerr_hzscl=0.046625,0.09325,0.1875 #hzscl/8
usenewgfsberror=.false.

# &HYBRID_ENSEMBLE
readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_v" setting
Expand All @@ -95,7 +98,7 @@ i_use_2mT4B=0
#-------------------------------------------------------------------------------------
#
BERROR_FN="berror.rtma" #under $FIX_GSI
ANAVINFO_FN="anavinfo.fv3lam_hrrr"
ANAVINFO_FN="anavinfo.rtma"
CONVINFO_FN="convinfo.rtma"
OBERROR_FN="errtable.rtma"
HYBENSINFO_FN="hybens_info.rtma"
Expand Down Expand Up @@ -178,6 +181,7 @@ NCORES_PER_NODE=40 #kjet
FG_ROOTDIR="/lfs4/BMC/nrtrr/NCO_dirs/nwges/para/RRFS_CONUS" #only needed by RTMA
DT_ATMOS=1
#adjust to low priority queue and possible preemption, slurm socket error, etc
MAXTRIES_RUN_PREPSTART="3"
MAXTRIES_RUN_FCST="3"
MAXTRIES_ANAL_GSI="3"
MAXTRIES_RUN_POST="3"
Expand Down Expand Up @@ -213,7 +217,6 @@ MODEL="RTMA_CONUS" #used by NCL
RUN="RTMA_CONUS"
COMINgfs=""
OBSPATH="/lfs4/BMC/wrfruc/gge/rtma/wcoss_obs"
BERROR_FN="rrfs_glb_berror.l127y194.f77"

STMP="/lfs4/BMC/nrtrr/NCO_dirs/stmp" # Path to directory STMP that mostly contains input files.
PTMP="/lfs4/BMC/nrtrr/NCO_dirs/ptmp" # Path to directory STMP that mostly contains input files.
Expand All @@ -230,14 +233,15 @@ if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
PARTITION_ANALYSIS=""
PARTITION_GRAPHICS=""
PARTITION_WGRIB2=""
DO_BUFRSND="FALSE"

OBSPATH_NSSLMOSIAC="/gpfs/dell1/nco/ops/dcom/prod/ldmdata/obs/upperair/mrms/conus/MergedReflectivityQC"
AIRCRAFT_REJECT="/gpfs/dell2/emc/modeling/noscrub/Edward.Colon/FixData/obsuselist/amdar_reject_lists"
SFCOBS_USELIST="/gpfs/dell2/emc/modeling/noscrub/Edward.Colon/FixData/obsuselist/mesonet_uselists"

WTIME_RUN_ANAL="00:30:00"
NNODES_RUN_ANAL="70"
PPN_RUN_ANAL="7"
NNODES_RUN_ANAL="28"
PPN_RUN_ANAL="14"

LAYOUT_X=42
LAYOUT_Y=40
Expand All @@ -259,6 +263,6 @@ if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
NWGES="/gpfs/dell2/ptmp/${USER}/nwges"
fi
# Wait longer for rrfsdas
if [[ "$MACHINE" == "jet" ]] && [[ $regional_ensemble_option -eq 5 ]] ; then
#if [[ "$MACHINE" == "jet" ]] && [[ $regional_ensemble_option -eq 5 ]] ; then
# START_TIME_ANAL="5:00:00"
fi
#fi
4 changes: 3 additions & 1 deletion ush/config.sh.RTMA_NA_3km
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,11 +75,13 @@ DO_DACYCLE="TRUE"
#-------------------------------------------------------------------------------------
#
# &SETUP and &BKGERR
l_obsprvdiag=.true.
diag_radardbz=.false.
write_diag_2=.true.
bkgerr_vs=0.125 #1.0 is default setting
#bkgerr_hzscl=0.373,0.746,1.5 #default setting, no trailing ,
bkgerr_hzscl=0.046625,0.09325,0.1875 #hzscl/8
usenewgfsberror=.false.

# &HYBRID_ENSEMBLE
readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_v" setting
Expand All @@ -99,7 +101,7 @@ i_use_2mT4B=0
#-------------------------------------------------------------------------------------
#
BERROR_FN="berror.rtma" #under $FIX_GSI
ANAVINFO_FN="anavinfo.fv3lam_hrrr"
ANAVINFO_FN="anavinfo.rtma"
CONVINFO_FN="convinfo.rtma"
OBERROR_FN="errtable.rtma"
HYBENSINFO_FN="hybens_info.rtma"
Expand Down
2 changes: 2 additions & 0 deletions ush/config_defaults.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -581,10 +581,12 @@ ARCHIVE_CYCLEDEF="00 01 01 01 2100 *"
# (need to follow FORTRAN namelist convetion)
#-------------------------------------------------------------------------------------
# &SETUP and &BKGERR
l_obsprvdiag=.false.
diag_radardbz=.true.
write_diag_2=.false.
bkgerr_vs=1.0
bkgerr_hzscl=0.7,1.4,2.80 #no trailing ,
usenewgfsberror=.true.

# &HYBRID_ENSEMBLE
readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_v" setting
Expand Down