Skip to content

Feature/rrfs dev1: add retro and radar tten functions#108

Merged
hu5970 merged 8 commits into
NOAA-GSL:feature/RRFS_dev1from
hu5970:feature/RRFS_dev1
May 18, 2021
Merged

Feature/rrfs dev1: add retro and radar tten functions#108
hu5970 merged 8 commits into
NOAA-GSL:feature/RRFS_dev1from
hu5970:feature/RRFS_dev1

Conversation

@hu5970
Copy link
Copy Markdown

@hu5970 hu5970 commented May 17, 2021

Add functions to run retro with staged data on disk.
Add namelist option for radar tten.

</task>

{%- endif %}
{%- if machine in ["JET", "HERA"] %}
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

A general question, do retros need the archive task?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I am not sure. But I never archived my retros runs and it is OK so far. If we need to archive retro for paper, we should only archive grib2 files under com.

Comment thread ush/config.sh.RRFS_dev1
if [[ $DO_RETRO == "true" ]] ; then

if [[ $MACHINE == "jet" ]] ; then
EXTRN_MDL_SOURCE_BASEDIR_ICS=/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/hrrr/conus/wrfnat/grib2
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do we have data staging scripts available? Or most users will get retro data from /mnt/lfs4/BMC/wrfruc/Ruifang.Li/data?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, we only needs one copy of the data in the disk and all retros should use this one copy of the data to save the disk space.

Comment thread ush/config.sh.RRFS_dev1 Outdated
FIX_CRTM=/home/rtrr/FIX_EXEC_MODULE/crtm/CRTM_v2.3.0
OBSPATH_NSSLMOSIAC=/public/data/radar/nssl/mrms/conus
RADARREFL_TIMELEVEL=(0 15 30 45)
#FH_DFI_RADAR=(0.0,0.25,0.5)
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do you mean FH_DFI_RADAR=(0.0 0.25 0.5) instead of FH_DFI_RADAR=(0.0,0.25,0.5) (the latter is a one element array, i.e. "0.0,0.25,0.5" ) ?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It has to include comma as it will be used to fill in namelist variable.

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Okay, I see. In this case, to avoid confusion, would it be enough to just let FH_DFI_RADAR="0.0,0.25,0.5" and hence no need to deal with an array which contains only one item?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

That doesn't not work when python generates the namelist for FV3.

Comment thread ush/config_defaults.sh Outdated
LAYOUT_X=""
LAYOUT_Y=""
BLOCKSIZE=""
FH_DFI_RADAR=(-20000000000)
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is it possible to leave FH_DFI_RADAR as an empty array if expected to run without radar tten?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We cannot set empty value to namelist parameters. This is default value for this namelist variable.

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

suggest changing to FH_DFI_RADAR="-20000000000" so as to match another suggested change FH_DFI_RADAR="0.0,0.25,0.5"

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Does not work that way.

Comment thread ush/set_extrn_mdl_params.sh Outdated
#
# for retro, the external boundary could come from other location.
#
if [[ "${DO_RETRO}" = "true" && ! -z "${EXTRN_MDL_SOURCE_BASEDIR_ICS}" ]]; then
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

To make a consistent style, suggest to use double equal sign here: if [[ "${DO_RETRO}" == "true" && ! -z "${EXTRN_MDL_SOURCE_BASEDIR_ICS}" ]]; then

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done

Comment thread ush/set_extrn_mdl_params.sh Outdated
fi


if [[ "${DO_RETRO}" = "true" && ! -z "${EXTRN_MDL_SOURCE_BASEDIR_LBCS}" ]]; then
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Similarly, suggest to use double equal sign in this if statement

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done

Comment thread ush/config_defaults.sh
#
# FH_DFI_RADAR (hour):
# The forecast period of using radar tten.
#
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

When does the period start and end? Is this 0-4 if FH_DFI_RADAR=4? Should it be an array? If so, how many entries, and what do they correspond to?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Added more details.


FV3_HRRR:
fv_core_nml:
dz_min: 6.0
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is this an intended change for tten?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is bug fix from Tanya.

Comment thread ush/config.sh.RRFS_dev1
OBSPATH_NSSLMOSIAC=/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/reflectivity
LIGHTNING_ROOT=/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/lightning/vaisala/netcdf/
ENKF_FCST=/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/enkf/atm
fi
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Please remove user-specific paths.

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is the location we stage the retro data now. Will be used by all other users to run the same retro.

Comment thread ush/config.sh.RRFS_dev1 Outdated
FIX_CRTM=/home/rtrr/FIX_EXEC_MODULE/crtm/CRTM_v2.3.0
OBSPATH_NSSLMOSIAC=/public/data/radar/nssl/mrms/conus
RADARREFL_TIMELEVEL=(0 15 30 45)
#FH_DFI_RADAR=(0.0,0.25,0.5)
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should this setting be on? If not, please remove commented code.

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Some will be on and some will be off. It will be used in 3km CONUS if the parallel test in the next couple of week show positive impact.

Copy link
Copy Markdown
Collaborator

@christinaholtNOAA christinaholtNOAA left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good. Just a few comments/questions.

@guoqing-noaa
Copy link
Copy Markdown
Collaborator

@hu5970 Thanks for the update. I suggest changing FH_DFI_RADAR=(0.0,0.25,0.5) to FH_DFI_RADAR="0.0,0.25,0.5" to avoid possible confusions. In the bash script, FH_DFI_RADAR is not an array of 3 or 5 items, but an array of one string "0.0,0.25,0.5".

@hu5970
Copy link
Copy Markdown
Author

hu5970 commented May 18, 2021

@guoqing-noaa string does not work. That is why I ended up with an array.

Copy link
Copy Markdown
Collaborator

@guoqing-noaa guoqing-noaa left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@hu5970 Sorry that I forgot to note that we need to make a minor update at ush/generate_FV3LAM_wflow.sh to allow to use FH_DFI_RADAR as a string. I just tested it and it worked. See my comment at ush/generate_FV3LAM_wflow.sh for the suggested change. I think a string is still preferred over an array for this situation. Thanks!

Comment thread ush/generate_FV3LAM_wflow.sh Outdated
}"
}
'gfs_physics_nml': {
'fh_dfi_radar': [${FH_DFI_RADAR[@]}],
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@hu5970 To accommodate changing FH_DFI_RADAR to a string, line 845 will need to be updated as:

-    'fh_dfi_radar': [${FH_DFI_RADAR[@]}],
+    'fh_dfi_radar': [${FH_DFI_RADAR}],

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

String works now and it is better solution. I have committed the changes. Thanks.

Copy link
Copy Markdown
Collaborator

@guoqing-noaa guoqing-noaa left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for making the updates for my picky comments. I will go ahead to approve the PR but hope the comments for FH_DFI_RADAR in ush/config_defaults.sh can be updated before the final merge. Thanks!

Comment thread ush/config_defaults.sh Outdated
# FH_DFI_RADAR:
# The forecast period of using radar tten. Array of 5:
# (0.0,0.25,0.5,0.75,1.0) to use radar tten in the first hour of the
# forecast with reading new observed tten every 15 minutes.
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggest to change the comments here to something as follows to adapt to the change of FH_DFI_RADAR to a string

# FH_DFI_RADAR:
# the forecast hour to use radar tten, this is used  to set the fh_dfi_radar parameter in input.nml*
# e.g. FH_DFI_RADAR="0.0,0.25,0.5,0.75,1.0" will set fh_dfi_radar = 0.0,0.25,0.5,0.75,1.0 in input.nml* and
# it tells the model to read and apply radar tten at the 0, 15, 30, 45, 60 minutes of forecasts respectively.
#

@hu5970 hu5970 merged commit 0ff99ef into NOAA-GSL:feature/RRFS_dev1 May 18, 2021
hu5970 added a commit that referenced this pull request May 29, 2021
* Feature/rrfs dev1: add retrospective function and changes for running on WCOSS (#104)

* Cold start RRFS_dev1 3km CONUS hourly cycles from 09/21Z HRRR analysis.
* Changes for running on RRFS cycles on WCOSS.

* Copy addtional grid from UPP run directory to comout space. (#105)

* Copy addtional grid from UPP run directory to comout space.

* Improvements for 3-km NA (#103)

* Extending walltime for "run_post" in 3-km NA

* Updated settings for layout_x, layout_y to improve run time.

* Increase RRFS_AK time step to match CONUS (#110)

* Feature/rrfs dev1: add retro and radar tten functions (#108)

* Add functions to run retro with staged data on disk.
* Add namelist option for radar tten.
* Bug fix for using GFS on JET for 3km NA.

* Feature/rrfs dev1: delete radar tten option (#112)

* delete tten option in namelist to match current USD model.

* Turn off cloud analysis for future debug. (#113)

* Change the "overwrite" option to "upgrade" (#111)

* Change $EXPDIR backup from inside to outside (#117)

* RRFS_dev1: Add python maps (#82)

* Updating job script to work with maps.

* Workflow for Python graphics.

* Use a metatask to manage map subsets.

Move skewt out of metatask.

* Sorting out plotting the wgrib2-produced domains.

* Updating configure scripts for appropriate domains.

* Adding domains for 3DRTMA

* Fixes for undefined variables and logic changes.

* RRFS_dev1: Turn on radar tten option and use "upgrade" to deploy. (#121)

* Add and turn on radar tten option
* Change PREEXISTING_DIR_METHOD from "rename" to "upgrade" for all RRFS_dev runs.

* Add configure options to decide how to clean the old cycles. (#123)

Change the reference date for clean from current run time to current cycle time

Co-authored-by: Jaymes Kenyon <jaymes.kenyon@noaa.gov>
Co-authored-by: Trevor Alcott <69324222+TrevorAlcott-NOAA@users.noreply.github.com>
Co-authored-by: Guoqing Ge <guoqing.ge@noaa.gov>
Co-authored-by: Christina Holt <56881914+christinaholtNOAA@users.noreply.github.com>
EdwardColon-NOAA pushed a commit to EdwardColon-NOAA/regional_workflow that referenced this pull request Sep 24, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants