Talk:Smoke and coupling with WRF-Chem - Revision history

Jmandel: Jmandel moved page Talk:Coupling with WRF-Chem to Talk:Smoke and coupling with WRF-Chem: Reflect better content of the page

2025-11-21T23:51:48Z

Jmandel moved page Talk:Coupling with WRF-Chem to Talk:Smoke and coupling with WRF-Chem: Reflect better content of the page

← Older revision	Revision as of 23:51, 21 November 2025
(No difference)

Jmandel: /* ch4 and h2 nuked */

2014-04-27T03:08:51Z

ch4 and h2 nuked

← Older revision		Revision as of 03:08, 27 April 2014
Line 36:		Line 36:

	==ch4 and h2 nuked==		==ch4 and h2 nuked==
	Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt>0]. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is emiss_inpt_opt=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. Not sure about call to [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_ltng_cpm.F#L211 adjust_ktop] which no longer gets called. [[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)		Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt>0]. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is emiss_inpt_opt=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. Not sure about call to [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_ltng_cpm.F#L211 adjust_ktop] which no longer gets called.
			[[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)

			Call mozcart_lbc_set is in few other places in [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F chem/chem_driver.F], added mozcart_lbc_set_onoff to the namelist.
			[[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 03:08, 27 April 2014 (UTC)

Jmandel: /* ch4 and h2 nuked */

2014-04-27T00:08:56Z

ch4 and h2 nuked

← Older revision		Revision as of 00:08, 27 April 2014
Line 36:		Line 36:

	==ch4 and h2 nuked==		==ch4 and h2 nuked==
	Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt>0]. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is ~~emiss_inp_opt~~=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. Not sure about call to [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_ltng_cpm.F#L211 adjust_ktop] which no longer gets called. [[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)		Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt>0]. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is emiss_inpt_opt=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. Not sure about call to [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_ltng_cpm.F#L211 adjust_ktop] which no longer gets called. [[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)

Jmandel: /* ch4 and h2 nuked */

2014-04-27T00:05:09Z

ch4 and h2 nuked

← Older revision		Revision as of 00:05, 27 April 2014
Line 36:		Line 36:

	==ch4 and h2 nuked==		==ch4 and h2 nuked==
	Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt]>0. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is emiss_inp_opt=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. [[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)		Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt>0]. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is emiss_inp_opt=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. Not sure about call to [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_ltng_cpm.F#L211 adjust_ktop] which no longer gets called. [[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)

Jmandel: /* Development notes */

2014-04-27T00:02:53Z

Development notes

← Older revision		Revision as of 00:02, 27 April 2014
Line 1:		Line 1:
	==~~Development~~ notes==		==Initial development notes==
	As of {{WRF-Fire-commit\|e7b52ab1baef2565ae1743a25bc1ed7f7f2fc109\|April 21, 2012}},		As of {{WRF-Fire-commit\|e7b52ab1baef2565ae1743a25bc1ed7f7f2fc109\|April 21, 2012}},
	the tracer is simply copied from the heat flux at the surface on the atmospheric mesh (<tt>GRNHFX</tt>).		the tracer is simply copied from the heat flux at the surface on the atmospheric mesh (<tt>GRNHFX</tt>).
Line 34:		Line 34:

	From Jon's text April 21 2012. [[User:Jmandel\|Jmandel]] 16:34, 23 December 2012 (UTC)		From Jon's text April 21 2012. [[User:Jmandel\|Jmandel]] 16:34, 23 December 2012 (UTC)

			==ch4 and h2 nuked==
			Even if we set ch4 and h2 in they get nuked in WRF-Chem in subroutine [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/module_input_chem_data.F#L3540 mozcart_lbc_set] called from [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L885 chem_driver] because [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/chem/chem_driver.F#L799 emiss_inpt_opt]>0. But the [https://github.com/jbeezley/wrf-fire/blob/e9051dc5ba0c59e1d7e4d59b47965d891c090142/wrfv2_fire/Registry/registry.chem#L2393 default] is emiss_inp_opt=1. This selection also calls emissions_driver, which seems unnecessary because we provide our own emissions. Specify emiss_inp_opt=0 in namelist.input. [[User:Jmandel\|Jmandel]] ([[User talk:Jmandel\|talk]]) 00:02, 27 April 2014 (UTC)

Jmandel: Created page with '==Development notes== As of {{WRF-Fire-commit|e7b52ab1baef2565ae1743a25bc1ed7f7f2fc109|April 21, 2012}}, the tracer is simply copied from the heat flux at the surface on the atmo…'

2012-12-23T16:34:35Z

Created page with '==Development notes== As of {{WRF-Fire-commit|e7b52ab1baef2565ae1743a25bc1ed7f7f2fc109|April 21, 2012}}, the tracer is simply copied from the heat flux at the surface on the atmo…'

New page

==Development notes==
As of {{WRF-Fire-commit|e7b52ab1baef2565ae1743a25bc1ed7f7f2fc109|April 21, 2012}},
the tracer is simply copied from the heat flux at the surface on the atmospheric mesh (<tt>GRNHFX</tt>).
A call to <tt>fire_emission</tt> has been added to <tt>module_fr_sfire_driver_wrf.F</tt> just after
calculation of heat flux tendencies. This subroutine has been added to the fire atmospheric module
and is preprocessed out of the source when WRF-Chem is not being compiled.

WRF-Chem contains a multitude of arrays representing concentrations of chemical species. The typical
method of injecting data into the simulation involves generating emission input data files from standard
sources. These sources are read into the code using auxiliary input streams and interpolated into
the species arrays in the <tt>emissions</tt> module. For the coupling with SFire, turning off file
input in the namelist seems to initialize the emission arrays to zero. This eliminates the need to generate
useless emission files.

WRF-Chem stores its state in several large arrays called <tt>*_emis</tt>, <tt>chem</tt>, and <tt>tracer</tt>. These arrays are four dimensional indexed like (x,z,y,s), where s is the chemical
species. The <tt>chem</tt>
and <tt>tracer</tt> arrays are located on the standard atmospheric grid, while the <tt>*_emis</tt> arrays
have an alternate vertical grid indexed to <tt>config%kemit</tt>.
An internal table keeps track of where each species is stored in these arrays. The domain structure
contains scalars such as <tt>p_smoke</tt> and <tt>p_e_so4</tt> that give the species index for
each type. The species quantities are specific to different chemistry options (<tt>chem_opt</tt>)
from the registry. We should eventually test that a compatible namelist option has been chosen to
avoid errors. The dynamics section of the namelist has an option for <tt>tracer_opt</tt>. This
controls what tracers are generated for Chem. For example, <tt>trace_opt = 1</tt> creates a single
tracer called <tt>smoke</tt>, see the registry for package definitions.

For the passive tracer scheme (<tt>chem_opt = 14</tt>), it appears to be sufficient just to add
the fire emission directly into the tracer array. If we want to fully utilize WRF-Chem in the future,
we may need to do something more complicated. The [http://ruc.noaa.gov/wrf/WG11/Users_guide.pdf User's guide] mentions that there is code to handle wildfire emissions data from MODIS satellites. This may be
helpful either directly or as a guide in how to inject the chemical species from the SFire model.

According to some comments in Registry.EM, basic tracer support exists in WRF without Chem. It is possible
we may be able to add smoke the WRF even without compile for Chem. I have not explored this yet, though.

From Jon's text April 21 2012. [[User:Jmandel|Jmandel]] 16:34, 23 December 2012 (UTC)