diff --git a/Registry/Registry.EM_COMMON b/Registry/Registry.EM_COMMON index f2811ea6bf..f7fe4632d0 100644 --- a/Registry/Registry.EM_COMMON +++ b/Registry/Registry.EM_COMMON @@ -3086,6 +3086,7 @@ package drip sf_surf_irr_alloc==2 - state:irrigat package sprinkler sf_surf_irr_alloc==3 - state:irrigation,irr_rand_field +package lesscheme bl_pbl_physics==0 - - package ysuscheme bl_pbl_physics==1 - - package myjpblscheme bl_pbl_physics==2 - state:tke_pbl,el_pbl package gfsscheme bl_pbl_physics==3 - - diff --git a/share/module_check_a_mundo.F b/share/module_check_a_mundo.F index 8acdd95e8f..44b227b5b3 100644 --- a/share/module_check_a_mundo.F +++ b/share/module_check_a_mundo.F @@ -406,6 +406,35 @@ END FUNCTION bep_bem_ngr_u END IF ENDDO + +!----------------------------------------------------------------------- +! Check that LES PBL is only paired with acceptable other PBL options. +! Currently, problems occur with any CG PBL option that has a packaged +! scalar component: MYNN2, MYNN3, EEPS. This test is also if a user +! chooses to not run a PBL scheme on a finer domain, wbut use a PBL +! parameterized scheme on a coarser domain (obviously, just for testing +! purposes). +!----------------------------------------------------------------------- + DO i = 1, model_config_rec % max_dom + IF ( .NOT. model_config_rec % grid_allowed(i) ) CYCLE + IF ( model_config_rec % bl_pbl_physics(i) .EQ. LESscheme ) THEN + exists = .TRUE. + END IF + END DO + IF ( ( exists ) .AND. & + ( ( model_config_rec % bl_pbl_physics(1) .EQ. MYNNPBLSCHEME2 ) .OR. & + ( model_config_rec % bl_pbl_physics(1) .EQ. MYNNPBLSCHEME3 ) .OR. & + ( model_config_rec % bl_pbl_physics(1) .EQ. EEPSSCHEME ) ) ) THEN + WRITE(wrf_err_message,fmt='(a,i2)') '--- ERROR: LES PBL on fine grid does not work with CG PBL option ',model_config_rec % bl_pbl_physics(1) + CALL wrf_message ( TRIM( wrf_err_message ) ) + wrf_err_message = ' Fix bl_pbl_physics in namelist.input: choose a CG PBL option without any scalar components' + CALL wrf_message ( TRIM( wrf_err_message ) ) + wrf_err_message = ' Alternatively, remove all of the packaged variables from the CG PBL selection' + CALL wrf_message ( TRIM( wrf_err_message ) ) + count_fatal_error = count_fatal_error + 1 + END IF + + !----------------------------------------------------------------------- ! Urban physics set up. If the run-time option for use_wudapt_lcz = 0, ! then the number of urban classes is 3. Else, if the use_wudapt_lcz = 1,