Error in wrf.exe when using LCZ maps via w2w

[Tarmactuna]

I am trying to use w2w to incorporate LCZ information into my simulations with WRF v4.6. I run the w2w package, and it works. I can generate the geo_em.dx.nc files correctly. The WPS and real.exe execution is also working fine. However, wrf.exe crashes almost immediately with the following error:

-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: LINE: 2419
Warning: too many input landuse types
-------------------------------------------

WPS, real.exe, and wrf.exe are running normally without LCZ.

How can I fix this? Thank you!

 

[Ming Chen]

What is w2w tool?
The global 100-m LCZ dataset has been introduced (WPSv4.5). Can you follow the normal WPS steps to process this data?

 

[Tarmactuna]

Hello, Mr. Chen.

Thank you so much for your answer. I incorporated the LCZs into my WRF simulation using the dataset you mentioned.

I am a WRF newbie and have spent an embarrassing amount of time trying to understand how to run the model. I am an undergraduate researcher and I want to measure outdoor thermal comfort/heat exposure at a local scale for my research project. I am still confused about what model to use and the differences between them (Noah LSM, SLUCM, MLUCM, MLUCM (BEP+BEM)). Could you provide me with some guidance on this? That is, which model could be appropriate for this, and where to find more resources about it (possibly a user's guide)?

Thank you so much for your time and help. It is greatly appreciated!

 

[Ming Chen]

What variables do you use to "measure outdoor thermal comfort/heat exposure "?

 

[Tarmactuna]

What variables do you use to "measure outdoor thermal comfort/heat exposure "?

2m air temperature. Ideally, this would take into account factors including but not limited to anthropogenic heat, urban geometry, materials, and wind flow. I want to get an accurate representation of "how hot" it feels outside.

 

[Ming Chen]

Both WRF and MPAS can well simulate surface conditions like 2-m temperature and humidity and 10-m winds, etc. While both models work fine, WRF has entered maintenance stage and will no longer be developed. In contrast, MPAS is the new generation of NWM model developed at NCAR-MMM lab. MPAS can run globally and regionally and has been experiencing robust development in recent years.
You can choose either WRF or MPAS for your work, and personally I would think MPAS could be a better option because it represents the cutting-edge modeling system.
if you run WRF, then Noah LSM with SLUCM, MLUCM, MLUCM (BEP+BEM) all can be a good options.
For WRF, please read the WRF User's Guide here
The MPAS User Guide can be found here.
Let me know if you have more questions. Thanks.

 

[Tarmactuna]

Dear Ming Chen,

Thank you for your detailed response and recommendations. Given my time constraints and growing familiarity with the WRF model, I have decided to continue using WRF for my research project. I have been simulating the 2m air temperature in New York City using the WRF model with the MLUCM incorporating the LCZs and the BEP model.

I believe I have successfully run the model and obtained results that seem accurate. However, I want to ensure that I am not missing any crucial parameters or steps. I have been following the instructions from the technical details provided in this ResearchGate document, which were designed for WRF V3.8. Since I am using WRF V4.6, I also referred to an example from the NCAR forum.

I have attached my current namelist.wps and namelist.input files to this email. Initially, I encountered an error while running the real executable with the following message:
module_io_quilt_old.F 2931 T
------ ERROR while reading namelist domains ------
Maybe here?: ! sfcp_to_sfcp = .true.,
Maybe here?: diff_6th_opt = 1, 1, 1
------ ERROR while reading namelist physics ------
Maybe here?: pxlsm_smois_init = 0,
Maybe here?: num_urban_layer = 360, 360, 360, ! Calculated value for all domains
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 11680
ERRORS while reading one or more namelists from namelist.input.
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0


I resolved this error by removing the lines:
diff_6th_opt = 1, 1, 1
num_urban_layer = 360, 360, 360, ! Calculated value for all domains

After doing this, the WRF executable ran successfully.

I would greatly appreciate your feedback on my namelist.wps and namelist.input files, particularly regarding any parameters I should add, delete, or modify to enhance the accuracy of my simulations. Please let me know if you need more information about my simulation setup to provide more specific advice.

Thank you very much for your time and assistance.

Best regards,

 

[Ming Chen]

Hi,
Your namelist.input looks fine to me except that:
(1) ra_sw_physics = 4, 4, 4,
(2) sf_urban_physics = 1, 1, 1 or sf_urban_physics = 2, 2, 2
(3) pxlsm_smois_init = 0, 0, 0
(4) delete num_urban_layer = 360, 360, 360,
(5) delete the following two options from &domains
diff_6th_opt = 1, 1, 1
diff_6th_factor = 0.12, 0.12, 0.12

 

[Tarmactuna]

Hello, Mr. Chen.

Thank you so much for your answer and I apologize for the late reply. Your help has been beneficial and I am immensely grateful for it.

I will apply the suggested modifications. I have a question on the second suggestion:

(2) sf_urban_physics = 1, 1, 1 or sf_urban_physics = 2, 2, 2

Would it not be computationally efficient to set it like:
sf_urban_physics = 1, 1, 2, ! Use SLUCM for the outer domains and BEP for the innermost domain

I am just curious as to whether this would work.

Thanks again for the help!

 

[Ming Chen]

It is recommended to run same physics (except cumulus scheme) over all domains. This is because when information is exchanged between parent and child domains, different physics over the domains may cause trouble.

 

[Tarmactuna]

Hi, Mr. Chen.

Apologies for the late reply. Thank you so much for your help!