updates related to MYNN-EDMF (#135)

* adding new variable (Sm3D) stability function for momentum * update submodule pointer for ccpp-physics * Point to gsl/develop for ccpp/physics Co-authored-by: samuel.trahan <[email protected]>
NOAA-GSL · Mar 24, 2022 · 036f0b3 · 036f0b3
1 parent ae7feac
commit 036f0b3
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 1 deletion.
diff --git a/ccpp/data/GFS_typedefs.F90 b/ccpp/data/GFS_typedefs.F90
@@ -1465,6 +1465,7 @@ module GFS_typedefs
     real (kind=kind_phys), pointer :: QI_BL      (:,:)   => null()  !
     real (kind=kind_phys), pointer :: el_pbl     (:,:)   => null()  !
     real (kind=kind_phys), pointer :: Sh3D       (:,:)   => null()  !
+    real (kind=kind_phys), pointer :: Sm3D       (:,:)   => null()  !
     real (kind=kind_phys), pointer :: qke        (:,:)   => null()  !
     real (kind=kind_phys), pointer :: tsq        (:,:)   => null()  !
     real (kind=kind_phys), pointer :: qsq        (:,:)   => null()  !
@@ -6241,6 +6242,7 @@ subroutine tbd_create (Tbd, IM, Model)
        allocate (Tbd%qi_bl     (IM,Model%levs))
        allocate (Tbd%el_pbl    (IM,Model%levs))
        allocate (Tbd%sh3d      (IM,Model%levs))
+       allocate (Tbd%sm3d      (IM,Model%levs))
        allocate (Tbd%qke       (IM,Model%levs))
        allocate (Tbd%tsq       (IM,Model%levs))
        allocate (Tbd%qsq       (IM,Model%levs))
@@ -6251,6 +6253,7 @@ subroutine tbd_create (Tbd, IM, Model)
        Tbd%qi_bl         = clear_val
        Tbd%el_pbl        = clear_val
        Tbd%sh3d          = clear_val
+       Tbd%sm3d          = clear_val
        Tbd%qke           = zero
        Tbd%tsq           = clear_val
        Tbd%qsq           = clear_val

diff --git a/ccpp/data/GFS_typedefs.meta b/ccpp/data/GFS_typedefs.meta
@@ -6268,6 +6268,14 @@
   type = real
   kind = kind_phys
   active = (flag_for_mellor_yamada_nakanishi_niino_pbl_scheme)
+[Sm3D]
+  standard_name = stability_function_for_momentum
+  long_name = stability function for momentum
+  units = none
+  dimensions = (horizontal_loop_extent,vertical_layer_dimension)
+  type = real
+  kind = kind_phys
+  active = (flag_for_mellor_yamada_nakanishi_niino_pbl_scheme)
 [qke]
   standard_name = nonadvected_turbulent_kinetic_energy_multiplied_by_2
   long_name = 2 x tke at mass points

diff --git a/ccpp/physics b/ccpp/physics
+1 −16		physics/GFS_rrtmg_pre.F90
+2 −9		physics/GFS_rrtmgp_gfdlmp_pre.F90
+37 −37		physics/module_MYNNPBL_wrapper.F90
+8 −0		physics/module_MYNNPBL_wrapper.meta
+199 −302		physics/module_bl_mynn.F90