* In ./src/core_atmosphere/physics:

  -> in ./physics_mmm, added the first version of the CCPP-compliant parameterization
     of the YSU planetary boundary layer (PBL; bl_ysu.F). modified Makefile accordingly.

  -> in ./physics_wrf, modified module_bl_ysu.F to call bl_ysu.F. modified Makefile
     accordingly.

  -> in mpas_atmphys_driver_pbl.F, change the call to subroutine ysu to reflect
     changes made in module_bl_ysu.F, mainly adding an error flag and an error
     message for CCPP compliance. modifed Makefile to compile the physics_mmm
     sub-directory.

Author: ldfowler58

src/core_atmosphere/physics/mpas_atmphys_driver_pbl.F

@@ -60,7 +60,7 @@ module mpas_atmphys_driver_pbl
 ! * for the mynn parameterization, change the definition of dx_p to match that used in other physics
 !   parameterizations.
 !   Laura D. Fowler (laura@ucar.edu) / 2016-10-18.
-! * updated the call to subroutine ysu in comjunction with updating module_bl_ysu.F from WRF version 3.6.1 to
+! * updated the call to subroutine ysu in conjunction with updating module_bl_ysu.F from WRF version 3.6.1 to
 !   WRF version 3.8.1
 !   Laura D. Fowler (laura@ucar.edu) / 2016-10-27.
 ! * since we removed the local variable pbl_scheme from mpas_atmphys_vars.F, now defines pbl_scheme as a pointer
@@ -69,6 +69,10 @@ module mpas_atmphys_driver_pbl
 ! * after updating module_bl_ysu.F to WRF version 4.0.3, corrected call to subroutine ysu to output diagnostics of
 !   exchange coefficients exch_h and exch_m.
 !   Laura D. Fowler (laura@ucar.edu) / 2019-03-12.
+! * updated the call to subroutine ysu after updating the YSU PBL scheme to that in WRF 4.4.1. added the flags
+!   errmsg and errflg in the call to subroutine ysu for compliance with the CCPP framework. also removed local
+!   variable regime_p which is no longer needed in the call to subroutine ysu.
+!   Laura D. Fowler (laura@ucar.edu) / 2023-05-15.
 
 
  contains
@@ -124,7 +128,6 @@ subroutine allocate_pbl(configs)
        if(.not.allocated(ctopo2_p)) allocate(ctopo2_p(ims:ime,jms:jme)      )
        if(.not.allocated(psih_p)  ) allocate(psih_p(ims:ime,jms:jme)        )
        if(.not.allocated(psim_p)  ) allocate(psim_p(ims:ime,jms:jme)        )
-       if(.not.allocated(regime_p)) allocate(regime_p(ims:ime,jms:jme)      )
        if(.not.allocated(u10_p)   ) allocate(u10_p(ims:ime,jms:jme)         )
        if(.not.allocated(v10_p)   ) allocate(v10_p(ims:ime,jms:jme)         )
        if(.not.allocated(exch_p)  ) allocate(exch_p(ims:ime,kms:kme,jms:jme))
@@ -212,7 +215,6 @@ subroutine deallocate_pbl(configs)
        if(allocated(ctopo2_p)) deallocate(ctopo2_p)
        if(allocated(psih_p)  ) deallocate(psih_p  )
        if(allocated(psim_p)  ) deallocate(psim_p  )
-       if(allocated(regime_p)) deallocate(regime_p)
        if(allocated(u10_p)   ) deallocate(u10_p   )
        if(allocated(v10_p)   ) deallocate(v10_p   )
        if(allocated(exch_p)  ) deallocate(exch_p  )
@@ -274,7 +276,7 @@ subroutine pbl_from_MPAS(configs,mesh,sfc_input,diag_physics,tend_physics,its,it
 
 !local pointers for YSU scheme:
  logical,pointer:: config_ysu_pblmix
- real(kind=RKIND),dimension(:),pointer:: br,fh,fm,regime,u10,v10
+ real(kind=RKIND),dimension(:),pointer:: br,fh,fm,u10,v10
  real(kind=RKIND),dimension(:,:),pointer:: rthratenlw,rthratensw
 
 !local pointers for MYNN scheme:
@@ -325,7 +327,6 @@ subroutine pbl_from_MPAS(configs,mesh,sfc_input,diag_physics,tend_physics,its,it
        call mpas_pool_get_array(diag_physics,'br'    ,br    )
        call mpas_pool_get_array(diag_physics,'fm'    ,fm    )
        call mpas_pool_get_array(diag_physics,'fh'    ,fh    )
-       call mpas_pool_get_array(diag_physics,'regime',regime)
        call mpas_pool_get_array(diag_physics,'u10'   ,u10   )
        call mpas_pool_get_array(diag_physics,'v10'   ,v10   )
 
@@ -341,7 +342,6 @@ subroutine pbl_from_MPAS(configs,mesh,sfc_input,diag_physics,tend_physics,its,it
           br_p(i,j)     = br(i)
           psim_p(i,j)   = fm(i)
           psih_p(i,j)   = fh(i)
-          regime_p(i,j) = regime(i)
           u10_p(i,j)    = u10(i)
           v10_p(i,j)    = v10(i)
           !initialization for YSU PBL scheme:
@@ -598,10 +598,18 @@ subroutine driver_pbl(itimestep,configs,mesh,sfc_input,diag_physics,tend_physics
  integer:: initflag
  integer:: i,k,j
 
+!CCPP-compliant flags:
+ character(len=StrKIND):: errmsg
+ integer:: errflg
+
 !-----------------------------------------------------------------------------------------------------------------
 !call mpas_log_write('')
 !call mpas_log_write('--- enter subroutine driver_pbl:')
 
+!initialization of CCPP-compliant flags:
+ errmsg = ' '
+ errflg = 0
+
  call mpas_pool_get_config(configs,'config_do_restart',config_do_restart)
  call mpas_pool_get_config(configs,'config_pbl_scheme',pbl_scheme       )
 
@@ -614,14 +622,14 @@ subroutine driver_pbl(itimestep,configs,mesh,sfc_input,diag_physics,tend_physics
  pbl_select: select case (trim(pbl_scheme))
 
     case("bl_ysu")
-       call mpas_timer_start('YSU')
+       call mpas_timer_start('bl_ysu')
        call ysu ( &
                  p3d      = pres_hyd_p , p3di     = pres2_hyd_p , psfc     = psfc_p     , &
-                 th3d     = th_p       , t3d      = t_p         , dz8w     = dz_p       , &
-                 pi3d     = pi_p       , u3d      = u_p         , v3d      = v_p        , &
-                 qv3d     = qv_p       , qc3d     = qc_p        , qi3d     = qi_p       , &
-                 rublten  = rublten_p  , rvblten  = rvblten_p   , rthblten = rthblten_p , &
-                 rqvblten = rqvblten_p , rqcblten = rqcblten_p  , rqiblten = rqiblten_p , & 
+                 t3d      = t_p        , dz8w     = dz_p        , pi3d     = pi_p       , &
+                 u3d      = u_p        , v3d      = v_p         , qv3d     = qv_p       , &
+                 qc3d     = qc_p       , qi3d     = qi_p        , rublten  = rublten_p  , &
+                 rvblten  = rvblten_p  , rthblten = rthblten_p  , rqvblten = rqvblten_p , &
+                 rqcblten = rqcblten_p , rqiblten = rqiblten_p  , flag_qc  = f_qc       , &
                  flag_qi  = f_qi       , cp       = cp          , g        = gravity    , &
                  rovcp    = rcp        , rd       = R_d         , rovg     = rdg        , & 
                  ep1      = ep_1       , ep2      = ep_2        , karman   = karman     , &
@@ -633,13 +641,14 @@ subroutine driver_pbl(itimestep,configs,mesh,sfc_input,diag_physics,tend_physics
                  exch_m   = kzm_p      , wstar    = wstar_p     , delta    = delta_p    , &
                  uoce     = uoce_p     , voce     = voce_p      , rthraten = rthraten_p , &
                  u10      = u10_p      , v10      = v10_p       , ctopo    = ctopo_p    , &
-                 ctopo2   = ctopo2_p   , regime   = regime_p    ,                         &
+                 ctopo2   = ctopo2_p   , flag_bep = flag_bep    , idiff    = idiff      , &
                  ysu_topdown_pblmix = ysu_pblmix ,                                        &
+                 errmsg   = errmsg     , errflg   = errflg      ,                         &
                  ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde  , &
                  ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme  , &
                  its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte    &
                 )
-       call mpas_timer_stop('YSU')
+       call mpas_timer_stop('bl_ysu')
 
     case("bl_mynn")
        call mpas_timer_start('MYNN_pbl')

src/core_atmosphere/physics/mpas_atmphys_vars.F

@@ -114,6 +114,9 @@ module mpas_atmphys_vars
 !   input of "observed" 2D seaice albedos. In conjunction with this update, we also change the initialization of
 !   albsi from albbck to seaice_albedo_default.
 !   Laura D. Fowler (laura@ucar.edu) / 2022-05-10.
+! * added the local parameters flag_bep and idiff in the call to subroutine ysu to update the YSU PBL scheme to
+!   that of WRF version 4.4.1.
+!   Laura D. Fowler (laura@ucar.edu) / 2023-05-15.
 
 
 !=================================================================================================================
@@ -348,6 +351,14 @@ module mpas_atmphys_vars
 !... variables and arrays related to parameterization of pbl:
 !=================================================================================================================
 
+ logical,parameter:: &
+    flag_bep = .false. !flag to use BEP/BEP+BEM for use in the YSU PBL scheme (with urban physics). since we do
+                       !not run urban physics, flag_bep is always set to false.
+
+ integer,parameter:: &
+    idiff = 0          !BEP/BEM+BEM diffusion flag for use in the YSU PBL scheme (with urban physics). since we
+                       !do not run urban physics, idiff is set to zero.
+
  integer:: ysu_pblmix
 
  integer,dimension(:,:),allocatable:: &

src/core_atmosphere/physics/physics_mmm/Makefile

@@ -6,7 +6,8 @@ dummy:
 	echo "****** compiling physics_mmm ******"
 
 OBJS = \
-	bl_gwdo.o
+	bl_gwdo.o \
+	bl_ysu.o
 
 physics_mmm: $(OBJS)
 	ar -ru ./../libphys.a $(OBJS)