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.

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

create a 13km conus domain 3dVar run using static background error covariance from Bump #5

Closed
TingLei-NOAA opened this issue Nov 9, 2023 · 9 comments
Assignees

Comments

@TingLei-NOAA
Copy link
Contributor

TingLei-NOAA commented Nov 9, 2023

On top of existing 13km conus domain pure Ens3Dvar run, we want to add in the bump based static background covariance component.
Hence, we will first create a static 3Dvar.
This issue will document how it is done and provide a place for comments/discussions.

@HuiLiu-NOAA
Copy link

HuiLiu-NOAA commented Nov 13, 2023 via email

@TingLei-NOAA
Copy link
Contributor Author

TingLei-NOAA commented Nov 13, 2023 via email

@TingLei-NOAA
Copy link
Contributor Author

TingLei-NOAA commented Nov 13, 2023

Bump official documentation from JCSDA: https://jointcenterforsatellitedataassimilation-jedi-docs.readthedocs-hosted.com/en/latest/inside/jedi-components/saber/BUMP.html.
Bump doc from the Bump User Training at 2023 : https://drive.google.com/file/d/1AOeECej4-45wpTuAa4O5zXcQ07-_Wdcc/view?usp=sharing and https://drive.google.com/file/d/1A3KjQfQHu_hyQUcK4539WQu2Dawp5w_D/view?usp=sharing (thanks to Benjemin and Anna)
The current test fv3jedi_test_tier1_3dvar_lam_cmaq could be used as a good example .

@andytangborn
Copy link

I think what you are looking for is the yaml that computes stddev, right? I use the following:

/scratch1/NCEPDEV/da/Andrew.Tangborn/JEDI/gdas_app_july28_2023/work_aero_c96/staticb_var_aero.yaml for aerosols fv3-jedi (so it might be a bit different from what you need). This step is after the staticb_prep_YYYYMMDD.yaml step.

I have examples of staticb_prep in the same directory. It is a bit confusing, but it creates two kinds of files: var and m4 (probably moment files), that are needed for computing the stddev files. The hardest part is that you need to have 1. An ensemble of background files at each time and 2. Three different times. I'm not running an ensemble, so I got around this by using 3-4 different experiments that have different physics or emissions, etc and call that my ensemble.

Have a look at my staticb_prep yaml files in that directory, and also notice the Data directories below them. You'll need to construct something very similar.

@HuiLiu-NOAA
Copy link

HuiLiu-NOAA commented Nov 14, 2023 via email

@TingLei-daprediction
Copy link
Contributor

Another helpful link on the static B training with bump https://drive.google.com/file/d/1jiflsG-p5g-_g-AsBuvhIxdKQTkJPX-9/view?usp=sharing.

@TingLei-NOAA
Copy link
Contributor Author

TingLei-NOAA commented Nov 15, 2023

Following the regional aero 's procedure as @andytangborn listed/explained (Thanks again), I have finished a set of runs for our 13km conus test case . Those runs succeed to completion and without apparent problems. Further verification /validation , especially by comparison with GSI is to be done in the future.
In the following , scripts to run each step are listed, all information needed to repeat those runs could be found through tracing those scripts (and yaml pointed in them).
On hera in /scratch2/NCEPDEV/fv3-cam/Ting.Lei/dr-emc-regional_jedi/dr-rrfs_rundir/hyb-jedianl_tm05_00_gsidiag_obs_new_jedi,

  1. d-static-prep.sh ( to run the prep step to get raw estimates of stddev and other moments).
  2. d-static-sdtdev.sh ( to get the calibrated estimates of stddev ( only one cycle raw data (generated in the previous step ) is used).
  3. d-3dvar.sh ( 3dvar run)
  4. d-hyb3dvar.sh ( hybrid 3dvar run) .

@HuiLiu-NOAA
Copy link

HuiLiu-NOAA commented Nov 15, 2023 via email

@TingLei-NOAA
Copy link
Contributor Author

It is completed and, a rrfs ctest for the final hybrid analysis has been added to RDAS rrfs-test

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

No branches or pull requests

4 participants