Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Enable tracer_opt=3 #46

Open
janmandel opened this issue Feb 5, 2022 · 2 comments
Open

Enable tracer_opt=3 #46

janmandel opened this issue Feb 5, 2022 · 2 comments
Assignees
Labels
active On going work

Comments

@janmandel
Copy link
Member

See also #22 - reconcile fire registry entries

@janmandel
Copy link
Member Author

janmandel commented Feb 5, 2022

tracer_opt in the registry

As of 0a634ac
in Registry/registry.chem:

# tracer options. (smoke) only has one tracer for smoke (has to run with biomassburning).
#                 (test1) has lateral boundary tracer(1,2), stratospheric tracer(5,6),
#                     boundary layer tracer(7,8), surface tracer (3,4)
#                 (test2) has surface tracer replaced by smoke tracer
package   tracer_smoke  tracer_opt==1       -             tracer:smoke
package   tracer_test1  tracer_opt==2       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
package   tracer_test2  tracer_opt==3       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
package   tracer_test3  tracer_opt==4       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8,tr18_1,tr18_2

Anywhere in the registry:

 grep -n tracer_opt Registry/*
Registry/Registry.EM:37:package   tracer_test1  tracer_opt==2       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
Registry/Registry.EM_COMMON:2704:rconfig   integer tracer_opt              namelist,dynamics     max_domains    0       rh    "tracer_opt"          ""      ""
Registry/Registry.NMM:57:#package   tracer_option_1     ntracers==4                -             tracers:t1,t2,t3,t4
Registry/Registry.NMM:84:package tracer_option_2 ntracers==4 - szj:szj1,szj2,szj3,szj4;s1z:s1z1,s1z2,s1z3,s1z4;spz:spz1,spz2,spz3,spz4;tcs:tcs1,tcs2,tcs3,tcs4
Registry/Registry.tladj:38:package   tracer_test1  tracer_opt==2       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
Registry/Registry.wrfvar:33:package   tracer_test1  tracer_opt==2       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
Registry/registry.chem:4125:package   tracer_smoke  tracer_opt==1       -             tracer:smoke
Registry/registry.chem:4126:package   tracer_test1  tracer_opt==2       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
Registry/registry.chem:4127:package   tracer_test2  tracer_opt==3       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8
Registry/registry.chem:4128:package   tracer_test3  tracer_opt==4       -             tracer:tr17_1,tr17_2,tr17_3,tr17_4,tr17_5,tr17_6,tr17_7,tr17_8,tr18_1,tr18_2
Registry/registry.fire:347:rconfig  real     fire_tracer_smoke         namelist,fire         max_domains            0.02       -     "parts per unit of burned fuel becoming smoke (tracer_opt=3)" "g_smoke/kg_air"
Registry/registry.fire:355:package  tracer_fire  tracer_opt==3       -             tracer:fire_smoke

ifire=2 - the original WRF-Fire modified at NCAR

ifire=2 is using tracer_opt=3 only:

module_fr_fire_driver.F:346:     if(fire_ifun.eq.6 .AND. config_flags%tracer_opt.eq.3)then
module_fr_fire_driver.F:347:       call add_fire_tracer_emissions(config_flags%tracer_opt,grid%dt,grid%dx,grid%dy, &

in subroutine add_fire_tracer_emissions as

avgw = 1.0/(ir*jr)
...
emis=avgw*fire_tracer_smoke*burnt_area_dt(i_f,j_f)*fgip(i_f,j_f)*1000/(rho(i,kts,j)*dz8w(i,kts,j)) ! g_smoke/kg_air
tracer(i,kts,j,p_fire_smoke)=tracer(i,kts,j,p_fire_smoke)+emis

where in module_fr_fire_phys.F

real,pointer,dimension(:,:):: fgip                 ! init mass of surface fuel (kg/m^2)

and fire_tracer_smoke is namelist.input flag defined in registry.fire:

rconfig  real fire_tracer_smoke         namelist,fire         max_domains            0.02       -     "parts per unit of burned fuel becoming smoke (tracer_opt=3)" "g_smoke/kg_air"

ifire=1 - the original WRF-Fire continued as WRF-SFIRE

WRF-SFIRE is using tracer_opt=2 and sets all 8 tracers to the same value, but the postprocessing in wrfxpy uses only tr17_1. There is also tracer_opt=1 in the registry, which generates only one tracer smoke, but only if compiled with CHEM.

@liy25
Copy link

liy25 commented Apr 19, 2023

can i design the PM as a tracer?

@Fergui Fergui changed the title enable tracer_opt=3 Enable tracer_opt=3 Jun 25, 2023
@Fergui Fergui added the active On going work label Jun 25, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
active On going work
Projects
None yet
Development

No branches or pull requests

3 participants