202305 Release of GFDL_atmos_cubed_sphere (#273)

This is the 202305 public release.  This release is the work of the GFDL FV3 development team.
* Merge branch 'TC20_test_case' into 'main'

Introducing idealized TC test case with SHiELD physics

See merge request fv3team/atmos_cubed_sphere!91

* Merge branch 'user/lnz/shield2022_gfdlmp' into 'main'

remove the abandoned file fv_diagnostics.F90.epv

See merge request fv3team/atmos_cubed_sphere!97

* Merge branch 'domain_deg_fix' into 'main'

bug fix to prevent model crash when using domain_deg = 0. in debug mode

See merge request fv3team/atmos_cubed_sphere!99

* Merge branch 'user/lnz/shield2022_gfdlmp' into 'main'

Fix to avoid reproducibility issue.

See merge request fv3team/atmos_cubed_sphere!100

* Merge branch 'user/lnz/shield2022_gfdlmp' into 'main'

A number of updates

See merge request fv3team/atmos_cubed_sphere!101

* Merge branch 'main' into 'main'

New test cases

See merge request fv3team/atmos_cubed_sphere!102

* Merge branch 'tau_w_202210' into 'main'

Improved Rayleigh Damping on w

See merge request fv3team/atmos_cubed_sphere!103

* Merge branch 'user/lnz/shield2022_gfdlmp' into 'main'

Update namelist reading code to avoid model crash because of the absence of naemlist.

See merge request fv3team/atmos_cubed_sphere!104

* Merge branch 'smag2d_202211' into 'main'

Experimental 2D Smagorinsky damping and tau_w

See merge request fv3team/atmos_cubed_sphere!107

* Merge branch 'user/lnz/shield2022_gfdlmp' into 'main'

Add the option to disable intermediate physics.

See merge request fv3team/atmos_cubed_sphere!108

* Merge branch 'user/lnz/shield2022_gfdlmp' into 'main'

Add the options to sub-cycling condensation evaporation, control the time scale of evaporation, and delay condensation and evaporation.

See merge request fv3team/atmos_cubed_sphere!109

* Merge branch 'user/lnz/shield2023' into 'main'

Remove grid size in energy and mass calculation

See merge request fv3team/atmos_cubed_sphere!110

* Merge branch 'user/lnz/shield2023' into 'main'

2023/03 Jan-Huey Chen

See merge request fv3team/atmos_cubed_sphere!112

* Merge branch 'lmh_public_release_202205' into 'main'

FV3 Solver updates 202305

See merge request fv3team/atmos_cubed_sphere!114

* gnu updates

* Merge branch 'user/lnz/shield2023' into 'main'

Pass the namelist variables from the dycore to the physics during the initialization

See merge request fv3team/atmos_cubed_sphere!117

* Merge branch 'main_mayrelease_smag_rollback' into 'main'

Rolling back smag damping

See merge request fv3team/atmos_cubed_sphere!118

* Merge branch 'lmh_revised_mapz' into 'main'

Revised vertical remapping operators

See merge request fv3team/atmos_cubed_sphere!116

* Merge branch 'main_upstream_mayrelease_dudz' into 'main'

Removed dudz and dvdz arrays that are not currently used.

See merge request fv3team/atmos_cubed_sphere!119

* Merge branch 'dp+nest' into 'main'

Add nest to DP cartesian config

See merge request fv3team/atmos_cubed_sphere!121

* Merge branch 'solonest' into 'main'

fix nesting in solo mode and add a new idealized test case for multiple nests

See merge request fv3team/atmos_cubed_sphere!120

* Merge branch 'fixsquare' into 'main'

fix square_domain logic for one tile grids

See merge request fv3team/atmos_cubed_sphere!122

* updating release notes.

---------

Co-authored-by: Lucas Harris <[email protected]>
Co-authored-by: Linjiong Zhou <[email protected]>
Co-authored-by: Joseph Mouallem <[email protected]>

GFDL_tools/fv_ada_nudge.F90

@@ -41,7 +41,6 @@ module fv_ada_nudge_mod
  use fms_mod,           only: write_version_number, check_nml_error
  use mpp_mod,           only: mpp_error, FATAL, stdlog, get_unit, mpp_pe, input_nml_file
  use mpp_mod,           only: mpp_root_pe, stdout ! snz
- use mpp_mod,           only: mpp_clock_id, mpp_clock_begin, mpp_clock_end
  use mpp_mod,           only: CLOCK_COMPONENT, CLOCK_SUBCOMPONENT, CLOCK_MODULE, CLOCK_ROUTINE
  use mpp_domains_mod,   only: mpp_update_domains, domain2d
  use mpp_domains_mod,   only: mpp_get_data_domain ! snz
@@ -233,7 +232,6 @@ module fv_ada_nudge_mod
   integer :: id_u_adj, id_v_adj, id_t_adj, id_q_adj, id_ps_adj ! snz
   integer :: id_u_a, id_v_a, id_t_a, id_q_a, id_ps_a ! snz
   integer :: id_u_da, id_v_da, id_t_da, id_q_da, id_ps_da ! snz
-  integer :: id_ada
 
   type(FmsNetcdfDomainFile_t) :: ada_driver_restart ! snz
   character(len=*), parameter :: restart_file="INPUT/ada_driver.res.nc" ! snz
@@ -439,7 +437,6 @@ subroutine fv_ada_nudge ( Time, dt, npx, npy, npz, ps_dt, u_dt, v_dt, t_dt, q_dt
 ! *t_obs* is virtual temperature
 #ifdef ENABLE_ADA ! snz
 
- call mpp_clock_begin(id_ada)
 !  call get_time (time, seconds, days)
 
   if (mod(seconds, 21600) == 0) then
@@ -662,7 +659,6 @@ subroutine fv_ada_nudge ( Time, dt, npx, npy, npz, ps_dt, u_dt, v_dt, t_dt, q_dt
     call mpp_update_domains(pt(:,:,:), domain, complete=.true.)
   end if
 
-  call mpp_clock_end(id_ada)
 
 #endif ! snz
 
@@ -3801,9 +3797,9 @@ subroutine del2_scalar(qdt, cd, kmd, nmax, bd,  npx, npy, gridstruct, domain)
          enddo
       enddo
    enddo
-                     call timing_on('COMM_TOTAL')
+   call timing_on('COMM_TOTAL')
    call mpp_update_domains(q, domain, complete=.true.)
-                     call timing_off('COMM_TOTAL')
+   call timing_off('COMM_TOTAL')
 
    do n=1,ntimes
 

README.md

@@ -1,6 +1,6 @@
 # GFDL_atmos_cubed_sphere
 
-The source contained herein reflects the 202210 release of the Finite Volume Cubed-Sphere Dynamical Core (FV3) from GFDL
+The source contained herein reflects the 202305 release of the Finite Volume Cubed-Sphere Dynamical Core (FV3) from GFDL
 
 The GFDL Microphysics is also available within this repository.
 

RELEASE.md

@@ -1,3 +1,45 @@
+# RELEASE NOTES for FV3 202305: Summary
+FV3-202305-public --- May 2023
+Lucas Harris, GFDL [email protected]
+
+This version has been tested with SHiELD physics release 202305
+and with FMS release 2023.01 from https://github.com/NOAA-GFDL/FMS
+
+This release includes the following:
+- Revised Vertical Remapping Operators (Lucas)
+  - kord=10 reverted back to AM4 version.
+  - Post-AM4 version of kord=10 is now kord=12.
+  - do_am4_remap no longer does anything and is deprecated.
+  - New strictly-monotone operators kord=14, 15 for improving tracer correlations, and kord=13 without subgrid limiting.
+  - kord <= 7 now deprecated; may be removed in a future release.
+- New Test Cases: (Joseph, Kun, Lucas)
+  - Idealized TC test case with SHiELD physics
+  - Zurita-Gotor et al. 2022 Held-Suarez variant
+  - New Stable Boundary Layer (Beale at al.) doubly-periodic test case
+- New nesting updates: (Joseph)
+  - Enable nesting in solo core and add a new idealized test case (58)
+  - Enable adding multiple nests in doubly-periodic test cases using absolute coordinates
+- Additional idealized capability (Linjiong, Kun, Lucas)
+  - Added namelist variable is_ideal_case, which must be used for runs starting (or re-starting) from idealized states.
+  - Begin saving the initial wind fields (u0 and v0) to the restart files
+- GFDL MP and Integrated Physics (Linjiong):
+  - Added options to sub-cycling condensation evaporation (nconds), control timescale or evaporation (do_evap_timescale), and delay condensation and evaporation (delay_cond_evap)
+  - Removed unused 3d microphysics diagnostics to save time and memory
+  - Optimized the mpp domain updates for fast physics
+  - Update gfdl_mp_nml reading code to avoid model crash for absent gfdl_mp_nml
+  - Added an option (do_intermediate_phys) to disable intermediate phys
+  - Removed grid size in GFDL MP energy and mass calculation
+  - Updates to use dry_cp instead of moist_cp in a hydrostatic case
+- Added a function to use O3 data from IFS ICs (Jan-Huey)
+  - Namelist parameter: “use_gfsO3” with the default value = “false”
+  - This function only works when ecmwf_ic = T
+  - If the IFS IC does not include O3 data, or the run would like to use GFS O3 with other IFS ICs, set use_gfsO3 = T
+- Solver Updates (Lucas)
+  - Revised semi-implicit solver to partially linearize vertical sound wave propagation about the hydrostatic state. This removes a specific instability causing deep “columnar” modes in the vertical velocity field due to the equation for the pressure perturbation being updated partially forward-in-time. This removes the spurious modes, reduces vertical velocities, and makes the solver slightly more stable.
+  - MPI bug fix for tracer diffusion
+  - Fast Rayleigh Damping on w controlled by fast_tau_w_sec.
+
+
 # RELEASE NOTES for FV3 202210: Summary
 FV3-202210-public --- October 2022
 Lucas Harris, GFDL [email protected]