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

Change Rayleigh damping settings to improve model instability and acoustic wave time step for C768 to improve computational efficiency #1612

Closed
XiaqiongZhou-NOAA opened this issue May 16, 2023 · 19 comments
Labels
feature New feature or request

Comments

@XiaqiongZhou-NOAA
Copy link

XiaqiongZhou-NOAA commented May 16, 2023

Description
This request is to add "fast_tau_w_sec" in dynamics input namelist (&fv_core_nml). This parameter controls implicit Rayleigh damping of vertical velocity introduced to UFS with PR#1720.
The suggested default value:
fast_tau_w_sec=0.2
tau=0
rf_cutoff=10
n_split=4

n_split is changed from 6 to 4 to save computational resources
and meantime the model remains stable.

The results from the corresponding sensitivity tests are here: https://docs.google.com/presentation/d/1n4DUr_4KqbJ5t04n2kgZDpLCJV_KpvXS/edit?usp=sharing&ouid=117079466243404562293&rtpof=true&sd=true

The impact of the change of Rayleigh damping on C384L127:
https://docs.google.com/presentation/d/1KMAuJ4H5uczhE5dZPrTiMhapktVMtGwKIa2ahGJYLO4/edit?usp=sharing

Requirements
HR2 can run with this new option.

Acceptance Criteria (Definition of Done)
The modification of the parameters is determined based on the verification of sensitivity experiments.

**Related to
ufs-community/ufs-weather-model#1824

@XiaqiongZhou-NOAA XiaqiongZhou-NOAA added the feature New feature or request label May 16, 2023
@aerorahul aerorahul added this to the HR2 milestone May 16, 2023
@junwang-noaa
Copy link
Contributor

junwang-noaa commented May 16, 2023

@XiaqiongZhou-NOAA Is the fast_tau_w_sec=0.2 used by all the resolutions (C48-C768)? Would you please create an issue in ufs weather model repo so that the UFS RT will use the same value in the namelist? I saw the default value for fast_tau_w_sec is 0.
@SadeghTabas-NOAA FYI

@XiaqiongZhou-NOAA
Copy link
Author

I would use this for HR2 at first and need more experiments to confirm.

@JessicaMeixner-NOAA
Copy link
Contributor

@XiaqiongZhou-NOAA the settings used in HR2 should match whatever settings we're using in the ufs-weather-model cpld_control_gfsv17 test so if it's not already there, so when you're ready to confirm what settings we should use, they should be updated in both that test in ufs-weather-model and here in the g-w.

@XiaqiongZhou-NOAA XiaqiongZhou-NOAA changed the title Change Rayleigh damping settings Change Rayleigh damping settings to improve model instability and acoustic wave time step for C768 to improve computational efficiency Jul 6, 2023
@yangfanglin
Copy link
Contributor

yangfanglin commented Jul 6, 2023

@XiaqiongZhou-NOAA Kate, please add a link to the Google Slides showing the test and evaluation of implicit Rayleigh damping you have done.

This PR contains two parts. The credit of the work goes to Kate (Xiaqiong) Zhou, Wei Li, Kevin Viner and our GFDL collaborators.

@JessicaMeixner-NOAA
Copy link
Contributor

@XiaqiongZhou-NOAA I removed the HR2 milestone from this issue as we don't have a corresponding update to ufs-weather-model yet.

@yangfanglin
Copy link
Contributor

yangfanglin commented Jul 19, 2023 via email

@JessicaMeixner-NOAA
Copy link
Contributor

By model changes, maybe I should be more specific by saying "the coefficient updates need to be updated in the GFSv17 ufs-weather-model regression test". Does that make more sense? If it's C768 setting only, then there's not a corresponding regtest. Otherwise, we try to keep things in sync.

@yangfanglin
Copy link
Contributor

yangfanglin commented Jul 19, 2023 via email

@XiaqiongZhou-NOAA
Copy link
Author

I only created an issue in the UFS-weather model not an PR and assume that the code manager of the ufs-weather-model will take care that since it is only a namelist change:
ufs-community/ufs-weather-model#1824

@JessicaMeixner-NOAA
Copy link
Contributor

The namelists for the regression tests are version controlled on github and therefore require a PR for any changes. Could you please submit one to the ufs-weather-model and then after that's merged we can also update the global-workflow?

@SadeghTabas-NOAA
Copy link

Dear @JessicaMeixner-NOAA and @XiaqiongZhou-NOAA, I'm taking care of the ufswm namelists to be matched with G-W side. Currently I'm updating the ufswm C$resolution tests to GFSv17. I'll take care of this in that PR or another one. Thanks

@yangfanglin
Copy link
Contributor

@SadeghTabas-NOAA, When can PR#1833 " ufswm C$resolution tests to GFSv17" be completed ? We'd like to include the new Rayleigh damping and new time-stepping for acoustic waves in the upcoming HR experiments sooner than later. Thanks.

@aerorahul
Copy link
Contributor

@yangfanglin
If this is a priority, then perhaps it is best if @XiaqiongZhou-NOAA could open a PR directly to global-workflow and if necessary to ufs-weather-model with a RT.
I would like to see how this fits with the configuration of other resolutions in the global-workflow.

@SadeghTabas-NOAA
Copy link

@SadeghTabas-NOAA, When can PR#1833 " ufswm C$resolution tests to GFSv17" be completed ? We'd like to include the new Rayleigh damping and new time-stepping for acoustic waves in the upcoming HR experiments sooner than later. Thanks.

@yangfanglin Sorry I missed this massage, I've updated the tests C48 and C96 however for the tests C192 and C384 we need inputs and ICs of the common time that are not currently available, so I don't know how long it may take to merge PR#1833. As @aerorahul mentioned, if this is a priority please open a new PR. Thanks and sorry for inconvenience.

@yangfanglin
Copy link
Contributor

Since HR2 has started to run, let's try to get this update into HR3. After the HR2 tags for workflow and model are created, @XiaqiongZhou-NOAA please create a PR to update model regression tests and another PR to update workflow.

@WalterKolczynski-NOAA WalterKolczynski-NOAA removed this from the HR2 milestone Sep 12, 2023
@WalterKolczynski-NOAA
Copy link
Contributor

@XiaqiongZhou-NOAA @yangfanglin

Is this something that you're still looking to change in workflow? I haven't seen a PR for it.

@JessicaMeixner-NOAA
Copy link
Contributor

I think the issue just wasn't linked. Here's a recent input.nml: /scratch1/NCEPDEV/climate/Jessica.Meixner/forlisa.input.nml (from a branch but no input.nml parts were updated). I see fast_tau_w_sec=0.2 updated. But it'd be good to confirm everything else is as expected.

@XiaqiongZhou-NOAA
Copy link
Author

XiaqiongZhou-NOAA commented Jan 16, 2024 via email

@WalterKolczynski-NOAA
Copy link
Contributor

Thanks @XiaqiongZhou-NOAA @JessicaMeixner-NOAA

Completed in #1987

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
feature New feature or request
Projects
None yet
Development

No branches or pull requests

7 participants