Use of sf_lake_physics

[GabrieleBentiv]

I am trying to run a WRF simulazione over the Alpine region, which works correctly unless I try to activate sf_lake_physics. In that case, it stops after a few seconds, with the last line in the output file being:
the total number of lake grid is : 0
(the result is the same with Noah and NoahMP)
What I did was using cglc_modis_lcz with 61 land uses and in the resulting wrfinput files I edited the LU_INDEX value of the water cells corresponding to lakes, by replacing their value of 17 with 21 (both domains have been edited).
I am linking the wrfinput files, the namelist, the ncl code I used and the rsl.out.0000 files (too big to be attached here): Proton Drive

The following is a screenshot of how the LU_INDEX field appears (as viewed with ncview):

 

[kwerner]

Hi,
Apologies for the long delay in response while our team tended to time-sensitive obligations. Thank you for your patience.

Since it's been a while, can you let me know if you're still experiencing this issue? If so, were you able to run wrf.exe to completion if you used cglc_modis_lcz as the geog data resolution and if :

1) you didn't turn on sf_lake_physics
2) you don't modify the wrfinput* files

 

[GabrieleBentiv]

Yes, I am curring running it correctly without the edited wrfinput files and the lake physics.

 

[kwerner]

What are the results of these two tests:

1) Using modified wrfinput* files, but turning sf_lake_physics OFF
2) Using UNMODIFIED wrfinput* files, with sf_lake_physics ON

And again, just to verify - when it's running correctly, you are still using cglc_modis_lcz for the geog data in wps?

 

[GabrieleBentiv]

The results in the two cases are:
1) The model runs with no issues, as if I didn't edit the file
2) Real is completed successfully but then, when running WRF, it stops immediately and the last lines in the rsl.out.0000 are the same as in the attachment provided in the initial question.

 

[kwerner]

Again, apologies for the delay. We've had back to back tutorials/workshops this summer and haven't had much opportunity to tend to forum inquiries.

I notice in the namelist.wps file, you have the following setting

geog_data_res = 'cglc_modis_lcz','cglc_modis_lcz',

What happens if you set that to

geog_data_res = 'cglc_modis_lcz+default','cglc_modis_lcz+default',

Does it make any difference when using lake physics?

 

[GabrieleBentiv]

In this case the operations fail at the real.exe step, with the attached rsl.out.0000 file

----- ERROR: sf_lake_physics = 1 AND flag_lake_depth = 0
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 347
Either modify the namelist settings, or rebuild the geogrid/metgrid data

 

[kwerner]

Thanks for trying that. Upon further research, it looks like you may need to use 'cglc_modis_lcz+modis_landuse_20class_15s' instead of '...+default." See Ming's post in this thread, which explains that. She also provides some references that may be useful.