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

Does anthropogenic heat in SLUCM include industry or traffic?

htan2013

Member
Dear all,

I am wondering whether the AH in SLUCM includes industry and traffic. If not, what does the AH include in WRF then? Does anyone know where the code is?
Also, in the URBPARM.TBL. The AH is 20/50/90 for Low-dens Res, Hi-dens Res, and Commercial. Are those numbers based on some specific city as an example?

Thanks in advance,
Haochen
 
Also, I have a question about the 2-m temperature in UCM.
For the non-urban part (no buildings etc), the 2-m temperature is 2-m air temp above the ground.
But when green roofs are applied over buildings, does the 2-m temperature represent the 2-m air temp above the roof surface or it is still the 2-m temperature on the building's ground?

Thank you,
Haochen
 
Haochen,
I will try to answer your questions, but I am not an expert in urban module and hopefully you can find more information in the literature.
AH in SLUCM is dependent on urban land use category (i.e., Low-dens Res, Hi-dens Res, Commercial). It should include AH from domestic, industry and traffic. However, note that the AH values are highly case-dependent. You need to adjust the values based on your case.
The AH values are read into WRF by the code phys/module_sf_urban.F, and the value is later passed to urban module for calculation of other variables. For example, in phys/module_sf_urban.F, you can find the following piece of code that uses AH:

Code:
1453 !-------------------------------------------------------------------------------
1454 ! Total Fluxes from Urban Canopy
1455 !-------------------------------------------------------------------------------
1456 !===Yang, 2014/10/08, cal. ah. alh. green roof===
1457  if(groption==1) then
1458    if(ahoption==1) then
1459      FLXTH  = ((1.-FGR)*R*FLXTHR + FGR*R*FLXTHGR + W*FLXTHB + RW*FLXTHG)+ AH/RHOO/CPP
1460    else
1461      FLXTH  = ((1.-FGR)*R*FLXTHR + FGR*R*FLXTHGR + W*FLXTHB + RW*FLXTHG)
1462    endif
1463    if(alhoption==1) then
1464      FLXHUM  = ((1.-FGR)*R*FLXHUMR + FGR*R*FLXHUMGR + W*FLXHUMB + RW*FLXHUMG)+ ALH/RHOO/ELL
1465    else
1466      FLXHUM  = ((1.-FGR)*R*FLXHUMR + FGR*R*FLXHUMGR + W*FLXHUMB + RW*FLXHUMG)
1467    endif
 
Last edited:
For 2-m temperature, I believe it is T at 2-m above the ground level. This is a diagnostic variable.
 
Thank you so much, Ming!
I looked at several works of literature before. One of them states that the AH in SLUCM includes traffic and industry but others don't.
I know that the AH values in URBPARM.tbl are diurnal profiles identical for all urban grid points, but I didn't know whether that value includes traffic or industry.
 
My understanding is that AH contains domestic, industrial and traffic heating in urban area. The concept of AH in WRF is quite general and its application is simple. Specifically, the values of AH could vary greatly in different regions. Based on this. no matter whether or not AH in the original WRF code includes industrial and traffic heating, you can always change its value for your case and run sensitivity experiments, then identify what values are the best.
 
Thank you so much, Ming!
I looked at several works of literature before. One of them states that the AH in SLUCM includes traffic and industry but others don't.
I know that the AH values in URBPARM.tbl are diurnal profiles identical for all urban grid points, but I didn't know whether that value includes traffic or industry.
Hello, sorry to bother you, I remember that I seem to ask you a question before, I want to ask you, I am studying the content of urban meteorology, I found several variables related to air conditioning, such as SF_AC_URB3D, CM_AC_URB3D, LF_AC_URB3D, and a variable E_PV_URB3D of photovoltaic conversion electric energy are all 0. I think there is no problem, but I still can’t output the value normally. Do you have any good suggestions? Attach my namelist.input, wps and URBPARM.TBL
 

Attachments

  • namelist.input
    4.4 KB · Views: 5
  • namelist.wps
    800 bytes · Views: 2
  • URBPARM .TBL.txt
    11.8 KB · Views: 3
The three variables (SF_AC_URB3D, CM_AC_URB3D, LF_AC_URB3D) are only available when you run with bep_bem scheme, i..e, sf_urban_physics= 3
You will also need to modify Registry.EM_COMMON. Below is an example. Change the line:
from
state real SF_AC_URB3D ij misc 1 - r "SF_AC_URB3D" "SENSIBLE HEAT FLUX FROM THE AIR COND." "W m{-2}"
to
state real SF_AC_URB3D ij misc 1 - rh "SF_AC_URB3D" "SENSIBLE HEAT FLUX FROM THE AIR COND." "W m{-2}"

Note the 8th column is changed from 'r' to 'rh'.


Please do the same for the other two variables.

Then you need to recompile WRF by
./clean -a
./configure
./compile em_real
 
Haochen,
I will try to answer your questions, but I am not an expert in urban module and hopefully you can find more information in the literature.
AH in SLUCM is dependent on urban land use category (i.e., Low-dens Res, Hi-dens Res, Commercial). It should include AH from domestic, industry and traffic. However, note that the AH values are highly case-dependent. You need to adjust the values based on your case.
The AH values are read into WRF by the code phys/module_sf_urban.F, and the value is later passed to urban module for calculation of other variables. For example, in phys/module_sf_urban.F, you can find the following piece of code that uses AH:


1453 !-------------------------------------------------------------------------------

1454 ! Total Fluxes from Urban Canopy


1455 !-------------------------------------------------------------------------------


1456 !===Yang, 2014/10/08, cal. ah. alh. green roof===


1457 if(groption==1) then


1458 if(ahoption==1) then


1459 FLXTH = ((1.-FGR)*R*FLXTHR + FGR*R*FLXTHGR + W*FLXTHB + RW*FLXTHG)+ AH/RHOO/CPP


1460 else


1461 FLXTH = ((1.-FGR)*R*FLXTHR + FGR*R*FLXTHGR + W*FLXTHB + RW*FLXTHG)


1462 endif


1463 if(alhoption==1) then


1464 FLXHUM = ((1.-FGR)*R*FLXHUMR + FGR*R*FLXHUMGR + W*FLXHUMB + RW*FLXHUMG)+ ALH/RHOO/ELL


1465 else


1466 FLXHUM = ((1.-FGR)*R*FLXHUMR + FGR*R*FLXHUMGR + W*FLXHUMB + RW*FLXHUMG)


1467 endif
Hello, if I want to output the heat emissions from air conditioning in the BEM while also outputting hourly traffic heat emissions at the base of the building, how should I modify the BEM-related code and set the relevant variables in the URBPARM.TBL file
 
@ Fenlin,
I don't have an immediate answer to your question. I will talk to our expert and hopefully he knows how to address your issue.
Note that it is hard for us to take care of individual requests, especially when the request involves great deals of coding and code changes. We appreciate users can dig into the codes, read relevant papers, and find answers to their specific questions.
 
sions from air conditioning in the BEM while also outputting hourly traffic heat emissions at the base of the building, how s
Hello, if I want to output the heat emissions from air conditioning in the BEM while also outputting hourly traffic heat emissions at the base of the building, how should I modify the BEM-related code and set the relevant variables in the URBPARM.TBL file
I believe there are outputs for sensible/latent heat from AC called SF_AC_URB3D, and LF_AC_URB3D.
 
Hello, I know that BEM can output the sensible and latent heat of air conditioning, but it cannot output other types of anthropogenic heat. The single-layer urban canopy model can input hourly anthropogenic heat values, which are set by ourselves. However, the multi-layer urban canopy model cannot input hourly anthropogenic heat. BEM can only input parameters for air conditioning, and then the model outputs the anthropogenic heat from air conditioning. Therefore, I want to set a fixed anthropogenic heat input at the bottom of BEM, but I need to modify the relevant code and I do not know how to modify it
 
@ Ming Chen
Hello, Professor Chen. I would like to ask how to replace the fixed artificial thermal parameters (20/50/90) of the UCM model with gridded artificial thermal data. Do I need to modify the source code or set the parameters in the namelist_input file? How should I specifically modify the source code or set the namelist_input parameters? I look forward to your reply. Thank you.
 
Top