Skip to content

GitHub Issue NOAA-EMC/GSI#538 Add options to tune weights of mixed ensemble for FV3-LAM EnVar#541

Merged
ShunLiu-NOAA merged 1 commit into
NOAA-EMC:developfrom
shoyokota:feature/PR_NOAA-EMC_EnVar-MixEn
Mar 17, 2023
Merged

GitHub Issue NOAA-EMC/GSI#538 Add options to tune weights of mixed ensemble for FV3-LAM EnVar#541
ShunLiu-NOAA merged 1 commit into
NOAA-EMC:developfrom
shoyokota:feature/PR_NOAA-EMC_EnVar-MixEn

Conversation

@shoyokota
Copy link
Copy Markdown
Collaborator

@shoyokota shoyokota commented Feb 21, 2023

This PR adds options (weight_ens_gfs and weight_ens_fv3sar) to tune weights of mixed ensemble for FV3-LAM EnVar (#538). Regression tests for global 3dvar/4denvar/4dvar are not completed yet, but for other tests, issues are not found except for "failed the scalability test" and "exceeded maximum allowable hardware memory limit" on Orion.

Fixes #538

@shoyokota
Copy link
Copy Markdown
Collaborator Author

Is it possible to ask @TingLei-NOAA and @lbi2018 to review this PR? I think they are well-qualified because they are major contributers of #451.

Comment thread src/gsi/gsimod.F90
write(6,*)'n_ens_gfs and n_ens_fv3sar won"t be used if not regional_ensemble_option==5'
endif

weight_ens_gfs=one
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Here, the namelist parameters like l_both_fv3sar_gfs_ens had been readin as well as weight_ens_gfs and weight_ens_fv3sar, right? then those assignments of the later two parameters (weight_ens??) will override the readin values for them. So, they will always be one.
Do I misunderstand something here?

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

In the case of l_both_fv3sar_gfs_ens=.false., only either global or regional ensemble is used, so weight_ens_gfs and weight_ens_fv3sar should be one (only either is used). To force that, I added this part. (If ensemble B should be tuned in this case, please don't set weight_ens_gfs or weight_ens_fv3sar but set beta_e0.) In the case of l_both_fv3sar_gfs_ens=.true., weight_ens_gfs and weight_ens_fv3sar are not overwritten and the parameters in the namelist are used.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

@shoyokota Thanks for your explanation . That does make sense.
This question is resolved

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

@shoyokota Could you please run a regression test and document the test result here.

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

@ShunLiu-NOAA I tried regression tests on Orion, but global 3dvar/4denvar/4dvar tests were failed because accesses to some files belonging to group "rstprod" are not permitted. I requested access to the NOAA RDHPCS Program, but it has not been accepted yet after two weeks. Even after it is accepted, I can run tests only on Orion. Could I ask someone to run the regression tests on my behalf?

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

@TingLei-NOAA Could you please run a regression test on HERA or WCOSS2?

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

@ShunLiu-NOAA Sure, I will give an update here when I have finished the regression tests.

Copy link
Copy Markdown
Contributor

@lbi2018 lbi2018 left a comment

Choose a reason for hiding this comment

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

Approved.

@ShunLiu-NOAA ShunLiu-NOAA self-requested a review March 9, 2023 14:38
@TingLei-NOAA
Copy link
Copy Markdown
Contributor

@ShunLiu-NOAA On hera, two failed tests: 2:global_4dvar;7:rrfs_3denvar_glbens. Seems GSI/fix/anavinfo.rrfs is missing for 7 while GSI/fix/global_anavinfo_qlqi.l127.txt is missing for 2 in this PR. It seems I still use an old fix dir. But "git submodule update" gave the error with message
"error: Server does not allow request for unadvertised object 0be26971f834fe9b1d5b118e1e0ffed53facf671
Fetched in submodule path 'fix', but it did not contain 0be26971f834fe9b1d5b118e1e0ffed53facf671. Direct fetching of that commit failed"
Are anything that I missed for the current GSI downloading?

@TingLei-NOAA
Copy link
Copy Markdown
Contributor

@ShunLiu-NOAA , seems I should first run "git submodule sync" . After that, "fix " dir was updated . The regressiong tests of those two failed ones is ongoing.

@ShunLiu-NOAA
Copy link
Copy Markdown
Contributor

I don't have any regression test on Hera recently. @hu5970 Could you work with Ting on this regression test?

@TingLei-NOAA
Copy link
Copy Markdown
Contributor

@ShunLiu-NOAA An update on my regression tests. Thanks to Russ, the current GSI regression tests procedure works fine for me on hera. But a few re-run of the ctest always gave me a couple of failed tests (different for different runs) for their wall clock time exceeding the time limits (even I changed all of them to 30:00). I will try more required time while I am wondering what the culprit is .

@RussTreadon-NOAA
Copy link
Copy Markdown
Contributor

@TingLei-NOAA , the most important aspect of the ctests is reproducibility independent of task and thread count. A few of the checks in the current suite of ctests need to be revised. While the test indicates Failed it is not a fatal fail. For example, slow I/O on Orion can lead to an anomalously large wall time. This is a system fileset issue, not a gsi or enkf code issue.

@TingLei-NOAA
Copy link
Copy Markdown
Contributor

@RussTreadon-NOAA Thank you for your clarification/reminding. Yes. the job killed for exceeding time would cause the regression test become idle. Since, such tests with this issue were different for different runs of the whole ctest (they succeeded in different runs), I would conclude that the current PR passed the current regression tests . @ShunLiu-NOAA .

@ShunLiu-NOAA ShunLiu-NOAA merged commit 89ca542 into NOAA-EMC:develop Mar 17, 2023
TingLei-daprediction referenced this pull request in TingLei-daprediction/GSI-old May 23, 2023
…semble for FV3-LAM EnVar (#541)

This PR adds options (weight_ens_gfs and weight_ens_fv3sar) to tune
weights of mixed ensemble for FV3-LAM EnVar
(NOAA-EMC/GSI#538). Regression tests for
global 3dvar/4denvar/4dvar are not completed yet, but for other tests,
issues are not found except for "failed the scalability test" and
"exceeded maximum allowable hardware memory limit" on Orion.

Fixes #538

Co-authored-by: Sho Yokota <syokota@Orion-login-1.HPC.MsState.Edu>
xyzemc pushed a commit to xyzemc/GSI_develop-v16-tms that referenced this pull request Feb 27, 2026
…le for FV3-LAM EnVar (NOAA-EMC#541)

This PR adds options (weight_ens_gfs and weight_ens_fv3sar) to tune
weights of mixed ensemble for FV3-LAM EnVar
(NOAA-EMC#538). Regression tests for
global 3dvar/4denvar/4dvar are not completed yet, but for other tests,
issues are not found except for "failed the scalability test" and
"exceeded maximum allowable hardware memory limit" on Orion.

Fixes NOAA-EMC#538

Co-authored-by: Sho Yokota <syokota@Orion-login-1.HPC.MsState.Edu>
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.

Add options to tune weights of mixed ensemble for FV3-LAM EnVar

5 participants