Merge branch 'jinboxie/atm/e3smv3_oro_merge' into next (PR #6667)

Implementation of new orographic drag parameterization scheme to alleviate wind bias in E3SM

Add a new suite of orographic drag parameterization schemes into EAM.
It includes 4 components all combined in one module (i.e. subroutine gwdo2d).
The schemes include orographic gravity wave drag (oGWD, Xie et al.,2020),
flow-blocking drag (FBD, Xie et al.,2020), small-scale GWD (sGWD, Tsiringakis et al.,2017),
and turbulent-scale orographic form drag (TOFD, Beljaars et al., 2004).

The scheme is distributed by the interface named gwdo_gsd. The oGWD, FBD, and
sGWD are implemented in gw_drag.F90, while TOFD is implemented in clubb_intr.F90
to join the vertical diffusion process. Another change is in the startup_initialconds.F90
and gw_drag.F90 where the new topographic parameters (OA, OC, OL) are input into the model.

Although it is possible to turn each new scheme on/off individually, the only validated
config is to have them all on (and the original scheme off).

Technical guide for this scheme is available at
https://docs.e3sm.org/E3SM/EAM/tech-guide/orodrag/

[BFB]

Author: wlin7

components/eam/src/physics/crm/dummy_modules/clubb_intr.F90

@@ -0,0 +1,31 @@
+module clubb_intr
+!-------------------------------------------------------------------------------
+! Dummy module to override src/physics/cam/clubb_intr.F90
+!-------------------------------------------------------------------------------
+use shr_kind_mod,  only: r8=>shr_kind_r8
+public :: clubb_implements_cnst
+public :: clubb_init_cnst
+public :: clubb_readnl
+contains
+!===============================================================================
+function clubb_implements_cnst(name)
+  ! Return true if specified constituent is implemented
+  character(len=*), intent(in) :: name ! constituent name
+  logical :: clubb_implements_cnst     ! return value
+  clubb_implements_cnst = .false.
+end function clubb_implements_cnst
+!===============================================================================
+subroutine clubb_init_cnst(name, q, gcid)
+  ! Initialize the state if clubb_do_adv
+  character(len=*), intent(in)  :: name     ! constituent name
+  real(r8),         intent(out) :: q(:,:)   ! mass mixing ratio (gcol, plev)
+  integer,          intent(in)  :: gcid(:)  ! global column id
+  return
+end subroutine clubb_init_cnst
+!===============================================================================
+subroutine clubb_readnl(nlfile)
+  character(len=*), intent(in) :: nlfile  ! filepath for file containing namelist input
+  return
+end subroutine clubb_readnl
+!===============================================================================
+end module clubb_intr

components/eam/src/physics/crm/physpkg.F90

@@ -638,7 +638,7 @@ subroutine phys_init( phys_state, phys_tend, pbuf2d, cam_out )
   if (co2_transport()) call co2_init()
   call co2_diags_init(phys_state)
 
-  call gw_init()
+  call gw_init(pbuf2d)
 
   call rayleigh_friction_init()