Skip to content

Add aerosol DA capabilities to the C96 atm-snow DA test#2323

Closed
CoryMartin-NOAA wants to merge 10 commits into
NOAA-EMC:developfrom
CoryMartin-NOAA:feature/gocart_inst_aod
Closed

Add aerosol DA capabilities to the C96 atm-snow DA test#2323
CoryMartin-NOAA wants to merge 10 commits into
NOAA-EMC:developfrom
CoryMartin-NOAA:feature/gocart_inst_aod

Conversation

@CoryMartin-NOAA
Copy link
Copy Markdown
Contributor

Description

This PR adds in 3DVar JEDI-based aerosol DA to the C96 atm-snow DA CI test.
Also included is a fix to copy the gocart inst_aod files that are hard coded to every 6 hours of the forecast. Thanks to @aerorahul for the help with that.

Things to note:
The B is 'identity', a future PR will likely add a realistic background error to the test.

This requires 1) new obs added to the dump, and 2) modified ICs for this test period (the gfs_ctrl and gfs_data.tile* files, the sfc* files are untouched). Both can be found on Hera at: /scratch1/NCEPDEV/da/Cory.R.Martin/testdata_aero or on Orion at: /work2/noaa/da/cmartin/work/feb2024/aero/testdata_aero

Resolves #2072
Partially resolves #1821

-->

Type of change

  • New feature (adds functionality)

Change characteristics

  • Is this a breaking change (a change in existing functionality)? NO
  • Does this change require a documentation update? NO

How has this been tested?

Example:

  • Cycled test on Orion
    -->

Checklist

  • Any dependent changes have been merged and published
  • My code follows the style guidelines of this project
  • I have performed a self-review of my own code
  • I have commented my code, particularly in hard-to-understand areas
  • My changes generate no new warnings
  • New and existing tests pass with my changes
  • I have made corresponding changes to the documentation if necessary

Comment thread ush/forecast_predet.sh

if [[ ! -d "${COM_CHEM_HISTORY}" ]]; then mkdir -p "${COM_CHEM_HISTORY}"; fi

GOCART_OUTPUT_FH=$(seq -s ' ' "${FHMIN}" "6" "${FHMAX}")

Check warning

Code scanning / shellcheck

GOCART_OUTPUT_FH appears unused. Verify use (or export if used externally).

GOCART_OUTPUT_FH appears unused. Verify use (or export if used externally).
Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

@aerorahul I thought I tested this but it seems it doesn't work...

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

oh, it's the ensemble fcst, somehow it's treated differently here?

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

I think once #2318 is addressed the line here as-is will be fine, but as it stands now, if you run in ATMA mode with an ensemble forecast FHMIN is 3 instead of 0, and then there isn't a gocart file at hour 3.

@CoryMartin-NOAA
Copy link
Copy Markdown
Contributor Author

This is in draft, because the ICs and obs need to be moved to the appropriate locations first

@KateFriedman-NOAA
Copy link
Copy Markdown
Contributor

KateFriedman-NOAA commented Feb 16, 2024

@CoryMartin-NOAA I have copied the provided obs into the GDA on Hera:

[role.glopara@hfe02 dump]$ pwd
/scratch1/NCEPDEV/global/glopara/dump
[role.glopara@hfe02 dump]$ ll g*2021122*/*/*viirs*                                                                                                               
-rw-r--r-- 1 role.glopara global 967425 Feb 16 16:13 gdas.20211220/18/gdas.t18z.viirs_n20.2021122018.nc4
-rw-r--r-- 1 role.glopara global 967425 Feb 16 16:13 gdas.20211220/18/gdas.t18z.viirs_npp.2021122018.nc4
-rw-r--r-- 1 role.glopara global 917601 Feb 16 16:13 gdas.20211221/00/gdas.t00z.viirs_n20.2021122100.nc4
-rw-r--r-- 1 role.glopara global 985169 Feb 16 16:13 gdas.20211221/00/gdas.t00z.viirs_npp.2021122100.nc4
-rw-r--r-- 1 role.glopara global 967425 Feb 16 16:13 gfs.20211220/18/gfs.t18z.viirs_n20.2021122018.nc4
-rw-r--r-- 1 role.glopara global 967425 Feb 16 16:13 gfs.20211220/18/gfs.t18z.viirs_npp.2021122018.nc4
-rw-r--r-- 1 role.glopara global 917601 Feb 16 16:13 gfs.20211221/00/gfs.t00z.viirs_n20.2021122100.nc4
-rw-r--r-- 1 role.glopara global 985169 Feb 16 16:13 gfs.20211221/00/gfs.t00z.viirs_npp.2021122100.nc4

Let me know if you need these synced to Orion or WCOSS2 GDAs.

@CoryMartin-NOAA
Copy link
Copy Markdown
Contributor Author

@KateFriedman-NOAA orion would be great, WCOSS is not necessary at this point (can't even build there!)

@KateFriedman-NOAA
Copy link
Copy Markdown
Contributor

@KateFriedman-NOAA orion would be great, WCOSS is not necessary at this point (can't even build there!)

@WalterKolczynski-NOAA Please sync the 20211220 and 20211221 GDA from Hera to Orion. Thanks!

@KateFriedman-NOAA
Copy link
Copy Markdown
Contributor

  1. modified ICs for this test period (the gfs_ctrl and gfs_data.tile* files, the sfc* files are untouched).

@CoryMartin-NOAA Are these a new set of ICs or they should overwrite existing ones?

@CoryMartin-NOAA
Copy link
Copy Markdown
Contributor Author

I believe they should overwrite the C96 gdas ICs but we should double check with @aerorahul if they should be overwritten or supplemented (for now)

@CoryMartin-NOAA CoryMartin-NOAA marked this pull request as ready for review February 29, 2024 19:18
@WalterKolczynski-NOAA
Copy link
Copy Markdown
Contributor

Dump sync'd to Orion

@KateFriedman-NOAA
Copy link
Copy Markdown
Contributor

@CoryMartin-NOAA Getting back to this, apologies for not completing the ICs part. Question, I don't see an enkf folder in your ICs but our current C96C48 set has enkf ICs too:

[role.glopara@hfe07 ICSDIR]$ ll /scratch1/NCEPDEV/da/Cory.R.Martin/testdata_aero/ICs/
total 4
drwxr-sr-x 4 Cory.R.Martin da 4096 Feb 16 16:12 gdas.20211220
[role.glopara@hfe07 ICSDIR]$ ll C96C48 
total 8
drwxr-sr-x 4 role.glopara global 4096 Jan  5 16:33 enkfgdas.20211220
drwxr-sr-x 4 role.glopara global 4096 Jan  5 16:40 gdas.20211220

Am I overwriting just the gdas.PDY folder with your ICs?

Also, once I pull your ICs in are you testing or shall we test with CI?

@CoryMartin-NOAA
Copy link
Copy Markdown
Contributor Author

@KateFriedman-NOAA we won't have aerosols in the EnKF, so it will just be the GDAS directory. Once pulling those in, please test with CI, it should work, but I can debug if things fail. Let me merge in develop here first since it's been a month.

@KateFriedman-NOAA
Copy link
Copy Markdown
Contributor

KateFriedman-NOAA commented Mar 11, 2024

@CoryMartin-NOAA Sounds good. I've pulled your gdas.PDY ICs into the C96C48 set on Hera:

[role.glopara@hfe07 C96C48]$ pwd
/scratch1/NCEPDEV/global/glopara/data/ICSDIR/C96C48
[role.glopara@hfe07 C96C48]$ ll
total 12
drwxr-sr-x 4 role.glopara global 4096 Jan  5 16:33 enkfgdas.20211220
drwxr-sr-x 4 role.glopara da     4096 Feb 16 16:12 gdas.20211220
drwxr-sr-x 4 role.glopara global 4096 Jan  5 16:40 gdas.20211220_SAVE_ORIG

I saved the prior ones just in case.

Confirm this meets your request and we'll sync these ICs to the other platforms and can fire off CI tests.

@CoryMartin-NOAA
Copy link
Copy Markdown
Contributor Author

@KateFriedman-NOAA yes, please do, thanks!

@KateFriedman-NOAA
Copy link
Copy Markdown
Contributor

I have synced the ICs to both WCOSS2s. Jet is down for maintenance. @WalterKolczynski-NOAA please sync the updated C96C48 ICs to Orion, thanks!

@aerorahul
Copy link
Copy Markdown
Contributor

@CoryMartin-NOAA
I am merging this PR into #2387 and will address the issues from jinja templates in the gdas app.

@aerorahul aerorahul closed this Mar 12, 2024
@aerorahul aerorahul mentioned this pull request Mar 12, 2024
7 tasks
@CoryMartin-NOAA CoryMartin-NOAA deleted the feature/gocart_inst_aod branch June 14, 2024 18:11
bbakernoaa pushed a commit to bbakernoaa/global-workflow that referenced this pull request Mar 19, 2026
…s_cubed_sphere to fix a typo when printing out the nest grid corner locations NOAA-EMC#2323+Bump urllib3 from 2.2.1 to 2.2.2 in /doc/UsersGuide NOAA-EMC#2329 + Orion changes for switch to new OS - Rocky-9 NOAA-EMC#2336 (NOAA-EMC#2324)

  * CICE - Adds restart write capability at the end of a run for CICE
  * FV3 - 
    * atmos_cubed_sphere - Update submodule FV3/atmos_cubed_sphere to fix a typo when printing out the nest grid corner locations
    * doc/UsersGuide/requirements.txt: Bumps [urllib3](https://github.com/urllib3/urllib3) from 2.2.1 to 2.2.2.
    * Orion Spack stack modulefile update: Rocky-9
    * update detect_machine.sh: hercules
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.

GOCART AOD files not copied properly Test cases for JEDI DA configurations

5 participants