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

Some questions about LCZ;WARNING, FRC_URB2D = 0 BUT IVGTYP IS URBAN WARNING, THE URBAN FRACTION WILL BE READ FROM URBPARM.TBL

likeyou511

New member
I use the 'cglc_modis_lcz_global'.
I want to use the 'URBPARM_LCZ.TBL' , but why does this warning appear?
 

Attachments

  • namelist.input
    3.8 KB · Views: 4
  • namelist.wps
    881 bytes · Views: 4
physics should be the same for all the domain, try to turn on same scheme for all your domains. In your namelist.input file you have used different schemes for domain(parent and child) , make it the same.

it would look like this ,

use_wudapt_lcz = 1,
mp_physics = 2, 2, 2,
cu_physics = 3, 0, 0,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
bl_pbl_physics = 11, 11, 11,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
radt = 12.5, 12.5, 12.5,
bldt = 0, 0, 0,
cudt = 5, 5, 5,
icloud = 1,
num_land_cat = 61,
num_soil_layers = 4,
sf_urban_physics = 1, 1, 1,
pxlsm_smois_init = 0, 0, 0.
 
Last edited:
physics should be the same for all the domain, try to turn on same scheme for all your domains. In your namelist.input file you have used different schemes for domain(parent and child) , make it the same.

it would look like this ,

use_wudapt_lcz = 1,
mp_physics = 2, 2, 2,
cu_physics = 3, 0, 0,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
bl_pbl_physics = 11, 11, 11,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
radt = 12.5, 12.5, 12.5,
bldt = 0, 0, 0,
cudt = 5, 5, 5,
icloud = 1,
num_land_cat = 61,
num_soil_layers = 4,
sf_urban_physics = 1, 1, 1,
pxlsm_smois_init = 0, 0, 0.
But about the cu_physics,shouldn't it be used when the distance is less than 10 kilometers.
And Is the physics plan related to my warning?
 
It depends on what you intend to do? what is your purpose of simulation? For dx ≥ 10 km .need cumulus scheme. About IVGTYP , if you see the module , module_sf_urban.F, it acts as conditional pointer for your land category.(for example :
IF( IVGTYP(I,J) == 31) THEN
UTYPE_URB2D(I,J) = 3 ! low-‐intensity residential,,,,,
....
and your urban physics should be the same for all domain. that is what works for me. check it maybe it works . (or maybe check without LCZ option if it is giving any warnings)

 
Last edited:
It depends on what you intend to do? what is your purpose of simulation? For dx ≥ 10 km .need cumulus scheme. About IVGTYP , if you see the module , module_sf_urban.F, it acts as conditional pointer for your land category.(for example :
IF( IVGTYP(I,J) == 31) THEN
UTYPE_URB2D(I,J) = 3 ! low-‐intensity residential,,,,,
....
and your urban physics should be the same for all domain. that is what works for me. check it maybe it works . (or maybe check without LCZ option if it is giving any warnings)

Thank you. I have set Urban physics to the same settings and I found that the warning does not appear.
But I don't really understand IVGTYP. Is there any document that describes this?
 
Glad it worked. Your are welcome. WUDAPT uses LCZ classification system that defines land cover types ( low-intenstity, high residetial,commercial, industrial .etc). IVGTYP(I,J) is used to identify this land cover type.

" " .
you can find some basic concept here but this document is a bit old. Starting from WPS4.5 , LCZ has been injected to WPS that lets you use the "CGLC-MODIS-LCZ" as geo data resolution. If you want to learn how it changed now then you can see the documentation here.


The classes used to be 31-41 before ,and one had to manually inject the LCZ by defining your classes (It has a separate procedure, you can find it in WUDAPT official website). Now, the classes are changed from 51-61 statrting from WPS-4.5 and that is why we use num of land cat = 61. IVGTYP is simply pointing to this classification , so that it can identify what is what and deploy the parameters accordingly.
if you want to learn more about how IVGTYP is working in the program, you can simply go to WRF/physics directory and study the code written for module_sf_urban.F.
 
Glad it worked. Your are welcome. WUDAPT uses LCZ classification system that defines land cover types ( low-intenstity, high residetial,commercial, industrial .etc). IVGTYP(I,J) is used to identify this land cover type.

" " .
you can find some basic concept here but this document is a bit old. Starting from WPS4.5 , LCZ has been injected to WPS that lets you use the "CGLC-MODIS-LCZ" as geo data resolution. If you want to learn how it changed now then you can see the documentation here.


The classes used to be 31-41 before ,and one had to manually inject the LCZ by defining your classes (It has a separate procedure, you can find it in WUDAPT official website). Now, the classes are changed from 51-61 statrting from WPS-4.5 and that is why we use num of land cat = 61. IVGTYP is simply pointing to this classification , so that it can identify what is what and deploy the parameters accordingly.
if you want to learn more about how IVGTYP is working in the program, you can simply go to WRF/physics directory and study the code written for module_sf_urban.F.
Thank you very much.Now I know.
 
It depends on what you intend to do? what is your purpose of simulation? For dx ≥ 10 km .need cumulus scheme. About IVGTYP , if you see the module , module_sf_urban.F, it acts as conditional pointer for your land category.(for example :
IF( IVGTYP(I,J) == 31) THEN
UTYPE_URB2D(I,J) = 3 ! low-‐intensity residential,,,,,
....
and your urban physics should be the same for all domain. that is what works for me. check it maybe it works . (or maybe check without LCZ option if it is giving any warnings)

I now encounter such a warning again.I want to know what does this FRC_URB2D mean?
How can I make sure that FRC_URB2D is not equal to 0.
I am using 'CGLC-MODIS-LCZ' data.
 

Attachments

  • namelist.wps
    905 bytes · Views: 1
Last edited:

FRC_URB2D is urban fraction at each grid cell. It is derived from static input data and used later when urban physics is activated. IVGTYP is land use type. In some cases FRC_URB2D = 0 at a grid cell of urban landuse type. In this case WRF will use values of FRC_URB2D from URBPARM.TBL to make it consistent with IVGTYP.

 
Top