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

Add soil increments to deterministic member #2773

Open
ClaraDraper-NOAA opened this issue Jul 17, 2024 · 3 comments · May be fixed by #3295
Open

Add soil increments to deterministic member #2773

ClaraDraper-NOAA opened this issue Jul 17, 2024 · 3 comments · May be fixed by #3295
Labels
feature New feature or request JEDI Feature development to support JEDI-based DA

Comments

@ClaraDraper-NOAA
Copy link
Contributor

What new functionality do you need?

The soil analysis (soil moisture, soil temperature) is done using the EnKF only, and we do not have a varitional soil analysis available for the deterministic member. Instead, we will apply the (upscaled) ensemble mean EnKF increment to the deterministic member.

This will require adding code for:

  1. calculation of the ensemble mean increment [add to enkfgdaseupd job]
  2. re-gridding the (Gaussian, ens res) mean increment to the deterministic native model resolution [in gdassfcanl?]

Initially, will use pre-existing sfc update executable (from UFS_UTILS) in gdassfcanl to add the deterministic increment (later, replace with IAU). This will require adding a dependency on enkfgdaseupd to gdassfcanl.

What are the requirements for the new functionality?

If GSI_SOILANAL option is selected, the ensemble mean soil increment is added to the deterministic member.

Acceptance Criteria

If GSI_SOILANAL option is selected, the ensemble mean soil increment is added to the deterministic member. If not, output should be zero diff.

Suggest a solution (optional)

No response

@ClaraDraper-NOAA ClaraDraper-NOAA added feature New feature or request triage Issues that are triage labels Jul 17, 2024
@ClaraDraper-NOAA
Copy link
Contributor Author

@aerorahul @WalterKolczynski-NOAA Do either of you know why we're calculating the ensemble mean increment using a purpose-designed executable, rather than using the option in the EnKF to output the ensemble mean increment? (asking, as I need to add calculation of the ensemble mean increment for the soil analysis, and am wondering if I can just use the EnKF output option).

@ClaraDraper-NOAA
Copy link
Contributor Author

@aerorahul @WalterKolczynski-NOAA Do either of you know why we're calculating the ensemble mean increment using a purpose-designed executable, rather than using the option in the EnKF to output the ensemble mean increment? (asking, as I need to add calculation of the ensemble mean increment for the soil analysis, and am wondering if I can just use the EnKF output option).

I figured it out. The code to write out increments from within the enkf requires the background (for conversion to Tv, and things like that). The background mean file is not available. Also it's non-linear, so using the mean background to write out the mean increment wouldn't give the same answer as calculating the mean of the ensemble output increments.

@WalterKolczynski-NOAA WalterKolczynski-NOAA removed the triage Issues that are triage label Jul 23, 2024
@ClaraDraper-NOAA
Copy link
Contributor Author

Update: I'll begin on this PR soon, and will likely merge in the capability to add the soil increments as IAU.

@WalterKolczynski-NOAA WalterKolczynski-NOAA added this to the GFS v17 milestone Jan 27, 2025
@WalterKolczynski-NOAA WalterKolczynski-NOAA added the JEDI Feature development to support JEDI-based DA label Jan 27, 2025
@RussTreadon-NOAA RussTreadon-NOAA removed this from the GFS v17 milestone Feb 8, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
feature New feature or request JEDI Feature development to support JEDI-based DA
Projects
None yet
Development

Successfully merging a pull request may close this issue.

3 participants