Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

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

  • namelistinput.txt
    6.4 KB · Views: 65
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
 
Top