Deng shallow convection option fails to work

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

sekluzia

Member
Dear Colleges,

I am trying to use the new shcu_physics=5 option of the WRF v4.1. The real.exe shows the following error:

Ntasks in X 1 , ntasks in Y 1
--- ERROR: Deng shallow convection can only work with MYJ or MYNN PBL
--- ERROR: Fix shcu_physics or bl_pbl_physics in namelist.input.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1935
NOTE: 1 namelist settings are wrong. Please check and reset these options
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1935
NOTE: 1 namelist settings are wrong. Please check and reset these options
-------------------------------------------

However, you can see in my namelist.input file (attached) that I set bl_pbl_physics=2 (MYJ PBL). Could you please help me to understand the problem?

I have two further questions regarding the use of Deng shallow convection option:
1) Would you recommend using this option for my nested domain which has convection-permitting resolution (3 km). Or it is better to turn off shallow convection for the nest, i.e. shcu_physics=0.
2) Is it recommended to use Deng shallow convection with the Grell-Freitas ensemble cumulus scheme (cu_physics = 3) for my parent domain?

Best regards,
Artur
 

Attachments

Artur,

This seems to me like a bug in module_check_a_mundo.F. Try this:

Go into share directory and in file module_check_a_mundo.F replace OR with AND at line 362:
Code:
(model_config_rec % bl_pbl_physics(i) /= myjpblscheme .OR. &
change into
Code:
(model_config_rec % bl_pbl_physics(i) /= myjpblscheme .AND. &

Recompile and try again, I believe it will work if my logic is correct.

Ivan
 
Ivan,
Good catch. Thanks! I'll put in a pull request to modify the code in our repository for the v4.1.1 bug-fix release.

Artur,
Hopefully Ivan's suggestion will correct the problem for you. To answer your other 2 questions:
1) Yes, a 3km domain would be a good resolution to use a shallow cumulus scheme.
2) Yes, it should be okay to use with the GF cumulus scheme.

Kelly
 
Back
Top