Skip to content

Support for high-resolution WACCM-X #1317

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

Open
wants to merge 7 commits into
base: cam_development
Choose a base branch
from

Conversation

fvitt
Copy link

@fvitt fvitt commented Jun 2, 2025

This PR sets default namelist settings for high-resolution (ne120pg3) WACCM-X.

Closes #1308

fvitt added 6 commits May 22, 2025 11:16
	modified:   bld/build-namelist
	modified:   bld/configure
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	modified:   cime_config/config_pes.xml
	modified:   bld/build-namelist
	modified:   bld/namelist_files/namelist_definition.xml
	modified:   src/dynamics/se/dycore/control_mod.F90
	modified:   src/dynamics/se/dycore/prim_advance_mod.F90
	modified:   src/dynamics/se/dyn_comp.F90
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	modified:   cime_config/testdefs/testlist_cam.xml
	modified:   cime_config/testdefs/testmods_dirs/cam/outfrq1h_wcmxhires/user_nl_clm

	new file:   cime_config/testdefs/testmods_dirs/cam/outfrq1h_wcmxhires/shell_commands
	new file:   cime_config/testdefs/testmods_dirs/cam/outfrq1h_wcmxhires/user_nl_cam
	new file:   cime_config/testdefs/testmods_dirs/cam/outfrq1h_wcmxhires/user_nl_clm
        modified:   bld/namelist_files/namelist_defaults_cam.xml
@fvitt fvitt added this to the CESM3 milestone Jun 2, 2025
@fvitt fvitt self-assigned this Jun 2, 2025
@fvitt fvitt added the BFB bit for bit tag label Jun 2, 2025
@cacraigucar cacraigucar self-requested a review June 3, 2025 17:25
Comment on lines 690 to +691
enddo
if (min_temperature>0._r8) then
Copy link
Collaborator

Choose a reason for hiding this comment

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

I was checking that being outside the loop over ic did not impact this new code. It isn't readily apparent to me what the loop over ic is actually doing, since the variable ic is not used in this subroutine. I assume that there must be something which is being advanced in biharmonic_wk_dp3d.

Copy link
Author

Choose a reason for hiding this comment

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

My change in within the substepping loop (within the ic loop). I assume the temperature floor should be applied each substep.

@PeterHjortLauritzen had some comments in the issue on where the temperature floor should be applied. I would like his input on this.

Copy link
Collaborator

Choose a reason for hiding this comment

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

I would follow DOE's implementation where the limiter is called in the last stage of the RK time-stepping:

https://github.com/E3SM-Project/E3SM/blob/28628343e7b5b25b9958e1ab601ef2b4e21a29ce/components/homme/src/theta-l/share/prim_advance_mod.F90#L326

That said, I don't know if it matters ... what we do know is that this seems to be an instability in this version of spetral-elements so putting it in the time-stepping of the inviscid equations of motion seems like a good place ...

More details from DOE here:

E3SM-Project/E3SM#5089

@@ -293,6 +297,8 @@ subroutine dyn_readnl(NLFileName)
call MPI_bcast(se_molecular_diff, 1, mpi_real8, masterprocid, mpicom, ierr)
call MPI_bcast(se_pgf_formulation, 1, mpi_integer, masterprocid, mpicom, ierr)
call MPI_bcast(se_dribble_in_rsplit_loop, 1, mpi_integer, masterprocid, mpicom, ierr)
call MPI_bcast(se_min_temperature, 1, mpi_real8, masterprocid, mpicom, ierr)
Copy link
Collaborator

Choose a reason for hiding this comment

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

@nusbaume - Probably next time we are in the SE dycore code, we should add checks for the MPI calls

group="dyn_se_inparm" valid_values="" >
Used by SE dycore to apply floor to temperature.....

NEEDS TO BE FILLED IN
Copy link
Collaborator

Choose a reason for hiding this comment

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

Just a reminder this does need to be filled in

Copy link
Author

Choose a reason for hiding this comment

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

done

@@ -212,6 +212,7 @@
<ncdata hgrid="ne120np4" nlev="30" aquaplanet="1" >atm/cam/inic/se/ape_cam5_ne120np4_L30_c170419.nc</ncdata>
<ncdata hgrid="ne120np4" nlev="32" >atm/cam/inic/se/F2000climo_ne120pg3_mt13_01-01-00000_c200420.nc</ncdata>
<ncdata hgrid="ne120np4" nlev="32" aquaplanet="1" >atm/cam/inic/se/ape_cam6_ne120np4_L32_c170908.nc</ncdata>
<ncdata hgrid="ne120np4" nlev="273" waccmx="1" >atm/waccm/ic/fx2000_ne120pg3L273.001.cam.i.0002-01-01-00000.noBRCL_c250522.nc</ncdata>
Copy link
Collaborator

Choose a reason for hiding this comment

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

Just a reminder to deal with this file after it is discussed in the CSEG meeting

	modified:   bld/namelist_files/namelist_definition.xml
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
BFB bit for bit tag
Projects
Status: No status
Development

Successfully merging this pull request may close these issues.

3 participants