Please help: Urban Canopy.

[Dipson11]

I have three questions.

1.My wrfout shows Lu_index from 1 to 56, I have use_wudapt_lcz =1 and num_land_cat=61 because now LCZ works for 51 to 61 according to the documentation of cglc modis 100m dataset In my geogrid I have LU_index: 1 TO 56 as well . Does it mean that my area doesnot of category from 57-61? .Am I doing anything wrong? I have attached my namelist.input files, and visuals.

2.I want to obtain variable tc_urb2d and qc_urb2d in wrfout , these variable are stated in urban canopy model documentation as urban state parameters but it is under noah mosaic scheme in registry.em_common, should I enable both sf_urban_physic and Noah Mosaic?. I cannot see urban state variables in my wrfout from sf_urban_physics=1.

3. When I run wrf.exe , I donot see how it is processing, it is generating the files as wrfout but no any progress movements in the screen(I have attached the piture named as Problem1.png). Also it doesnot go inside the third domain which is the innermost domain, it is simulating only for d01 and d02 (Attachment:wrfout.png). I have 160available processor but if I increase the number of processors it gives error.(Attachment: tail.error.png), I am using only 4 processor for this simulation..Could you please check whether my domain looks normal and let me know where I am lacking.

Thank you very much.

 

[Ming Chen]

Please see the answers below:

I have three questions.

1.My wrfout shows Lu_index from 1 to 56, I have use_wudapt_lcz =1 and num_land_cat=61 because now LCZ works for 51 to 61 according to the documentation of cglc modis 100m dataset In my geogrid I have LU_index: 1 TO 56 as well . Does it mean that my area doesnot of category from 57-61? .Am I doing anything wrong? I have attached my namelist.input files, and visuals.

This indicates that there is no landuse type of 57-61 in your domain. This is normal and you don't need to worry.

2.I want to obtain variable tc_urb2d and qc_urb2d in wrfout , these variable are stated in urban canopy model documentation as urban state parameters but it is under noah mosaic scheme in registry.em_common, should I enable both sf_urban_physic and Noah Mosaic?. I cannot see urban state variables in my wrfout from sf_urban_physics=1.

Please add 'h' to the 8th column of the variable definition, for example, change the line from:

state    real   TC_URB2D         ij     misc        1         -     rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm)     "TC_URB"              "URBAN CANOPY TEMPERATURE"           "K"

To:

state    real   TC_URB2D         ij     misc        1         -     rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm)     "TC_URB"              "URBAN CANOPY TEMPERATURE"           "K"

Then recompile WRF. With the above changes, you will have the variables included in your wrfout files.

3. When I run wrf.exe , I donot see how it is processing, it is generating the files as wrfout but no any progress movements in the screen(I have attached the piture named as Problem1.png). Also it doesnot go inside the third domain which is the innermost domain, it is simulating only for d01 and d02 (Attachment:wrfout.png). I have 160available processor but if I increase the number of processors it gives error.(Attachment: tail.error.png), I am using only 4 processor for this simulation..Could you please check whether my domain looks normal and let me know where I am lacking.

Thank you very much.

You should monitor the progress of the model run by looking at the RSL files in your working directory.
If the case works fine, D03 should be activated and involved in the run.

 

[Dipson11]

Please see the answers below:

This indicates that there is no landuse type of 57-61 in your domain. This is normal and you don't need to worry.

Please add 'h' to the 8th column of the variable definition, for example, change the line from:

state    real   TC_URB2D         ij     misc        1         -     rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm)     "TC_URB"              "URBAN CANOPY TEMPERATURE"           "K"

To:

state    real   TC_URB2D         ij     misc        1         -     rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm)     "TC_URB"              "URBAN CANOPY TEMPERATURE"           "K"

Then recompile WRF. With the above changes, you will have the variables included in your wrfout files.




You should monitor the progres sof the model run by looking at the RSL files in your working directory.
If the case works fine, D03 should be activated and involved in the run.

Hi chen,

I used the instruction in documentation to add the variables without compiling, I made a text file myoutfields.txt which has
+:h:24:TC_URB2D,QC_URB2D

In my namelist.input
iofields_filename = 'myoutfields.txt' , 'myoutfields.txt' , 'myoutfields.txt',
!io_form_auxhist24=2
auxhist24_interval = 20, 20, 20,

But, it gives me output,
Assume Noah LSM input
d02 2018-01-01_00:00:00 forcing artificial silty clay loam at 131 points, out of 35640
d02 2018-01-01_00:00:00 Timing for processing 1 s.
W A R N I N G : Unable to modify mask for tc_urb2d. Variable not found. File: myoutfields.txt at line 1
W A R N I N G : Unable to modify mask for qc_urb2d. Variable not found. File: myoutfields.txt at line 1
W A R N I N G : Unable to modify mask for tc_urb2d. Variable not found. File: myoutfields.txt at line 1
W A R N I N G : Unable to modify mask for qc_urb2d. Variable not found. File: myoutfields.txt at line 1
d02 2018-01-01_00:00:00 Timing for output 0 s.
d02 2018-01-01_00:00:00 Timing for loop # 1 = 2 s.
d01 2018-01-06_00:00:00 real_em: SUCCESS COMPLETE REAL_EM INIT

Is it not possible from auxhist method as mentioned in the documention of Registy?

Thank you for your support.

 

[Ming Chen]

The runtime I/O option only works for wrf.exe. It has no effects when you run real.exe.

 

[dipson]

The runtime I/O option only works for wrf.exe. It has no effects when you run real.exe.

Thanks ming, it works now!