search

openwfm / WRF-SFIRE

A coupled weather-fire forecasting model built on top of Weather Research and Forecasting (WRF). This is the original https://github.com/openwfm/wrf-fire transitioned to a fork of WRF and selected as ifire=1. Graphic log at https://repo.or.cz/git-browser/by-commit.html?r=WRF-SFIRE.git
https://wiki.openwfm.org
Other
42 stars 12 forks source link

Carry over WRF-Fire-Chem #16

Closed janmandel closed 4 years ago

janmandel commented 4 years ago

Check if WRF-SFIRE-Chem works after the transition to WRF 4 and fix if needed. There will be differences in chem since we took it in and changes in chem we made long time ago...

janmandel commented 4 years ago

test/em_fire/chem no longer works because of unknown flag mozcart_lbc_set_onoff in namelist.input namelist chem . Turns out I added the mozcart_lbc_set_onoff flag and some code in 92b5d3442 on April 26, 2014. Files changed: Registry/registry.chem | 1 + chem/chem_driver.F | 14 +++++++++++--- phys/module_fr_sfire_atm.F | 78 ++++++++++++++++++++++++++++++++++++++++++++++++------------------------------ test/em_fire/chem/namelist.input | 2 ++

Commit 92b5d3442 also adds some computation and prints of the sum of emissions over all MPI processes. In wrf-fire-filtered/master commit 92b5d3442 is commit bb8a1fb7d7a22e61b792577a2491ea2f2395accb but directory chem is not included in filtered/master. Files Registry/registry.chem and chem/chem_driver.F did not change after bb8a1fb7d7a22e61b792577a2491ea2f2395accb. In master and filtered/devel, 92b5d3442 is filtered as c703224aaad8efe9b832533f1458c27914ae0984 but contains only the two files outside of chem.

Found email:

Date: Sat, 26 Apr 2014 22:11:07 -0600 Subject: Re: Missing fluxes in WRF-SFIRE-CHEM - solved From: Jan Mandel jan.mandel@gmail.com To: Jan Mandel jan.mandel@gmail.com, ADAM KRZYSZTOF KOCHANSKI adam.kochanski@utah.edu, Jonathan Beezley jon.beezley@gmail.com Cc: Martin Vejmelka vejmelkam@gmail.com

Hi Adam,

the values of h2 and ch4 in the lowest layer get nuked in every step of chem by subroutine mozcart_lbc_set. I have added new switch mozcart_lbc_set_onoff in chem section in namelist input to bypass all calls to it. Please set mozcart_lbc_set_onoff=0, default is 1. I am not sure what this may break, only Georg Grell would know. We should ask. It does seem to take care of the problem. You might want to run with emiss_inpt_opt=3, the default is 1 which calls emissions_driver and I am not sure you want that, we have our own. BTW even before the code would bomb out on janus on segfault. No idea why. Your case is running on janus in /lustre/janus_scratch/jmandel/wrf-fire-chem/wrfv2_fire/test/em_chefire The updated code is in master, commit b41f70c81e19205ce1960f330f202027b851a463. I hope that helps.

Thanks,

Jan

janmandel commented 4 years ago

git cherry-pick -x bb8a1fb7d7a22e61b792577a2491ea2f2395accb resolved conflict in chem/chem_driver.F now at fa8c72ad7368d2c0b3ed92b08e3377722df4e111 branch develop-chem now ideal.exe wants wrf_chem_input_d01 did we have that before? stop in: main/ideal_em.F: CALL med_read_wrf_chem_input ( grid , config_flags) share/mediation_integrate.F:SUBROUTINE med_read_wrf_chem_input ( grid , config_flags ) reading namelists in: frame/module_configure.f90 SUBROUTINE initial_config

janmandel commented 4 years ago

chem_in_opt=0 in namelist.input 104f7f256a6171ec1849996f6d9d1b8001575b05 to prevent reading wrf_chem_input_d01 but now ideal.exe says error in namelist.fire but nothing wrong with namelist.fire cp ../hill/namelist.fire . b0e37768fe8ddb67d7240f2991946e6a888ba5ef to verify, still the same error but wrf-fire/master a14489984701a0ceb6c70e27148482a445572385 is fine

janmandel commented 4 years ago

test/em_fire/chem/namelist.input had ifire=2 so ideal.exe trying to read the wrong namelist fire fixed in a5df4fbab504cbf4a0aa38caec794fcc7d716d14 ideal.exe calls read_namelist_fire from phys/module_fr_fire_phys.F but with ifire=1, but for ifire=2 it should use phys/module_fr_sfire_phys.F

I changed dyn_em/module_initialize_fire.F and phys/module_fr_fire_phys.F in wrf-model#792 53d63b053fcb90586466596389e7cfa0ebd1e75b

Diagnosis: incompatibility when ifire=2 and ifire=2 share dyn_em/module_initialize_fire.F, which calls routine from ifire=1. I kind of expected that but not so soom.

Possible solutions:

  1. put separate code for ifire=1 and ifire=2 in dyn_em/module_initialize_fire.F and USE both phys/module_fr_fire_phys.F phys/module_fr_sfire_phys.F with different local names by =>
  2. divorce em_sfire from em_fire to generate a different ideal.exe #17

Bypassed already in edb91d99 - the guilty code is called if ifire=2 only

janmandel commented 4 years ago

To compile:

export WRF_CHEM=1
./configure chem
./compile em_fire

On cheyenne: qcmd -l walltime=4:00:00 -- ./compile -j 32 em_fire >& compile.log & https://www.openwfm.org/wiki/Coupling_with_WRF-Chem https://ruc.noaa.gov/wrf/wrf-chem/Users_guide.pdf https://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.2/

janmandel commented 4 years ago

Chem files commits: commit 92b5d3442bfd02ec03fd4bbc675b85ef4ec4c893 Author: jmandel jmandel@colibri.ucdenver.pvt Date: Sat Apr 26 19:02:53 2014 -0600

adding mozcart_lbc_set_onoff config flag
more chem diagnostics for Layer1, fixing test/em_fire/chem/namelist.input

:100644 100644 817b97aa1 d6c5cd026 M wrfv2_fire/Registry/registry.chem :100755 100755 3570946dd b336c8f6f M wrfv2_fire/chem/chem_driver.F :100644 100644 e74ab6902 7bc64101e M wrfv2_fire/phys/module_fr_sfire_atm.F :100644 100644 b9ce32a9d 03c773128 M wrfv2_fire/test/em_fire/chem/namelist.input

commit d505f86491fe868450c94bb4da6304e16ceedb2a Author: Jan Mandel Jan.Mandel@gmail.com Date: Wed Dec 26 13:04:40 2012 -0700

added sfc_full_init =.true. to initialize usgs, no need for make_usgs.*

:100644 000000 f47dfe2c3 000000000 D wrfv2_fire/test/em_fire/chem/make_usgs.m :100755 000000 b33485757 000000000 D wrfv2_fire/test/em_fire/chem/make_usgs.py :100644 100644 f498b1eaf b9ce32a9d M wrfv2_fire/test/em_fire/chem/namelist.input

commit 578534e46264b800f98d1d0f707f0099a20dcbe8 Author: Jan Mandel jmandel@frontend.cluster Date: Tue Dec 18 20:36:48 2012 -0700

in module_aerosols_sorgam.F:5156 also treat the case aa<0

:100644 100644 c5d49ca9a 644c3a0fd M wrfv2_fire/chem/module_aerosols_sorgam.F :100644 100644 41ef474b6 aa7e60d29 M wrfv2_fire/phys/module_fr_sfire_atm.F

commit c69fd04e63063288bba367fb3a77c7c7ad0b2daa Author: Jan Mandel Jan.Mandel@gmail.com Date: Tue Dec 18 20:09:41 2012 -0700

fix in module_aerosols_sorgam.F:5156 when no positive root exists treat like complex roots
from Jon:
Specifically, the fatal error at line 1205 in module_aerosols_sorgam.F is caught.
The message here is "Ionic strength below zero...negative concentrations'".  As far
as I can tell none of the quantities coming in indicate negative concentrations.
I traced the problem back to module_aerosols_sorgam.F:5156 where some quadratic polynomial
is being solved for the nitrate concentration.   There seems to be an implicit
assumption that at least one root of the polynomial is positive, but that is not the
case here.  I'll let you know if I find out more.

:100644 100644 72996f25c c5d49ca9a M wrfv2_fire/chem/module_aerosols_sorgam.F

commit b9f8db3105365f0b03cd88da21c9314441a4a683 Author: Jonathan Beezley jon.beezley@gmail.com Date: Sat Apr 21 17:46:39 2012 -0600

adding python script to get chem example working

:000000 100755 000000000 b33485757 A wrfv2_fire/test/em_fire/chem/make_usgs.py

commit e7b52ab1baef2565ae1743a25bc1ed7f7f2fc109 Author: Jonathan Beezley jon.beezley@gmail.com Date: Sat Apr 21 16:11:54 2012 -0600

putting fire emissions into tracer array

:100644 100644 f1bf27eb1 f3ee782ac M wrfv2_fire/phys/module_fr_sfire_atm.F :100644 100644 c3bbfd5c7 1636e39f4 M wrfv2_fire/phys/module_fr_sfire_driver_wrf.F :000000 100644 000000000 0f91581b8 A wrfv2_fire/test/em_fire/chem/input_sounding :100644 100644 d38054f13 5e5daa3d0 M wrfv2_fire/test/em_fire/chem/namelist.input

commit 2fd240869e798d48dad969051f390cdf56c90373 Author: Jonathan Beezley jon.beezley@gmail.com Date: Sat Apr 21 11:51:51 2012 -0600

tracers showing up in atmosphere

:100644 100644 94b1d7e3b f1bf27eb1 M wrfv2_fire/phys/module_fr_sfire_atm.F :100644 100644 81c58ea51 d38054f13 M wrfv2_fire/test/em_fire/chem/namelist.input

commit 299981cb47d67eb2ef66c4ec7a285ad98a4a7dd5 Author: Jonathan Beezley jon.beezley@gmail.com Date: Sat Apr 21 10:39:54 2012 -0600

work in progress for chem/fire integration

:100644 000000 fdfd7ff59 000000000 D wrfv2_fire/Registry/.registry.chem.swp :100644 100644 5648dfe43 94b1d7e3b M wrfv2_fire/phys/module_fr_sfire_atm.F :100644 100644 d2d45ed91 c3bbfd5c7 M wrfv2_fire/phys/module_fr_sfire_driver_wrf.F :000000 100644 000000000 837a890b4 A wrfv2_fire/test/em_fire/chem/LANDUSE.TBL :000000 100644 000000000 4f670771e A wrfv2_fire/test/em_fire/chem/namelist.fire :000000 100644 000000000 81c58ea51 A wrfv2_fire/test/em_fire/chem/namelist.input

commit 0969a14d659d202c74025327da4f608da8a9b9de (tag: wrf-3.4) Author: Jonathan Beezley jon.beezley@gmail.com Date: Fri Apr 20 17:31:34 2012 -0600

r5744 | weiwang | 2012-04-06 18:08:27 -0600 (Fri, 06 Apr 2012)
janmandel commented 4 years ago

Running test/em_fire/chem with git commit 1bf1e5f2ed6eb544a75658f217cc3673427c52c3

 *************************************
 Parent domain
 ids,ide,jds,jde            1          42           1          42
 ims,ime,jms,jme           -4          47          -4          47
 ips,ipe,jps,jpe            1          42           1          42
 *************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,              354031592  bytes allocated
  wrf main: calling open_r_dataset for wrfinput
  med_initialdata_input: calling input_input
   Input data is acceptable to use: wrfinput_d01
  mminlu = 'USGS'
   NOTE:  Ideal cases always use hypsometric_opt=1, regardless of namelist setting
 CURRENT DATE          = 0001-01-01_00:00:00
 SIMULATION START DATE = 0001-01-01_00:00:00
Timing for processing wrfinput file (stream 0) for domain        1:    0.76908 elapsed seconds
Max map factor in domain 1 =  0.00. Scale the dt in the model accordingly.
  WRF TILE   1 IS      1 IE     42 JS      1 JE     42
  set_tiles3: NUMBER OF TILES =   1
INPUT LandUse = "USGS"
 LANDUSE TYPE = "USGS" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =            0  SNOW CATEGORY =            0
*********************************************************************
*             PROGRAM:WRF-Chem V4.2 MODEL
*                                                                   *
*    PLEASE REPORT ANY BUGS TO WRF-Chem HELP at                     *
*                                                                   *
*              wrfchemhelp.gsd@noaa.gov                             *
*                                                                   *
*********************************************************************
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     502
 ERROR: CHEM_INIT: USGS land use map should have 24 or more catagories.
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     502
 ERROR: CHEM_INIT: USGS land use map should have 24 or more catagories.
-------------------------------------------
janmandel commented 4 years ago

With chem_opt=0 SFIRE:crash: add_fire_emissions: a chem pointer is out of bounds - pointers 1:2, chem array size 1 With chem_opt=300 runs no crash With chem_opt=2 TBD

janmandel commented 4 years ago

chem_opt=0 and tracer_opt=1 creates one tracer called smoke with value 2 chem_opt=0 and tracer_opt=0 does not call WRF-Chem at all, no banner

janmandel commented 4 years ago

Ready for testing on real with chem_opt=300

janmandel commented 4 years ago

code compiled with WRF-Chem: chem_opt=300 tracer_opt= 1 num_chem= 20 num_tracer= 2 max chem= 20 max tracer= 1 chem_opt= 0 tracer_opt= 1 num_chem= 1 num_tracer= 2 max chem= 2 max tracer= 1, smoke=2 chem_opt=300 tracer_opt= 0 num_chem= 20 num_tracer= 1 max chem= 20 max tracer= 1, smoke=1 chem_opt=0 and tracer_opt=0 we are not calling fire_add_emissions

janmandel commented 4 years ago

chem_opt=300 tracer_opt= 1 num_chem= 20 num_tracer= 2 max chem= 20 max tracer= 2 chem_opt=300 tracer_opt= 2 num_chem= 20 num_tracer= 9 max chem= 20 chem_opt=300 tracer_opt= 3 num_chem= 20 num_tracer= 10 max chem= 20 chem_opt=300 tracer_opt= 4 crash: read_emissions_table: increase max_tracer, have num_tracer=11 chem_opt=300 tracer_opt= 5 num_chem= 20 num_tracer= 1 max chem= 20 max tracer= 1 (= no tracers)

janmandel commented 4 years ago

From 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
janmandel commented 4 years ago

Setup in ~/p/em_rxcadre_ref Code and runs:

janmandel commented 4 years ago

Test of parallel consistency on test/em_fire/chem: small differences in chem fields only. Starting as rounding error, then grow. Example after 51 frames in wrfout: ncdiff('chem-4x4/wrfout_d01_0001-01-01_00:00:00','chem-2x3/wrfout_d01_0001-01-01_00:00:00') **Finished comparing files chem-4x4/wrfout_d01_0001-01-01_00:00:00 chem-2x3/wrfout_d01_0001-01-01_00:00:00 ACGRDFLX relative diff 0 p-value 0 ACHFX relative diff 0 p-value 0 ACLHF relative diff 0 p-value 0 ACSNOM relative diff 0 p-value 0 ALBBCK relative diff 0 p-value 0 ALBEDO relative diff 0 p-value 0 AOD2D_OUT relative diff 0 p-value 0 AOD_OUT relative diff 0 p-value 0 AREA2D relative diff 0 p-value 0 ATOP2D_OUT relative diff 0 p-value 0 AVG_FUEL_FRAC relative diff 8.58903e-05 p-value 0.975608 BBB relative diff 0 p-value 0 BC1 relative diff 0.00362439 p-value 1 BC2 relative diff 0.00362377 p-value 1 C1F relative diff 0 p-value 0 C1H relative diff 0 p-value 0 C2F relative diff 0 p-value 0 C2H relative diff 0 p-value 0 C3F relative diff 0 p-value 0 C3H relative diff 0 p-value 0 C4F relative diff 0 p-value 0 C4H relative diff 0 p-value 0 CANHFX relative diff 0 p-value 0 CANQFX relative diff 0 p-value 0 CANWAT relative diff 0 p-value 0 CAN_TOP relative diff 0 p-value 0 CF1 relative diff 0 p-value 0 CF2 relative diff 0 p-value 0 CF3 relative diff 0 p-value 0 CFN relative diff 0 p-value 0 CFN1 relative diff 0 p-value 0 CLAT relative diff 0 p-value 0 CLDFRA relative diff 0 p-value 0 CLDFRA2 relative diff 0 p-value 0 CN2O5 relative diff 0 p-value 0 COSALPHA relative diff 0 p-value 0 COSZEN relative diff 0 p-value 0 CUF relative diff 0 p-value 0 CVF relative diff 0 p-value 0 DMS_0 relative diff 0 p-value 0 DN relative diff 0 p-value 0 DNW relative diff 0 p-value 0 DRYDEPVEL relative diff 0 p-value 0 DRY_DEP_LEN relative diff 0 p-value 0 DUST_1 relative diff 0.00362412 p-value 1 DUST_2 relative diff 0.00362412 p-value 1 DUST_3 relative diff 0.00362412 p-value 1 DUST_4 relative diff 0.00362412 p-value 1 DUST_5 relative diff 0.00362412 p-value 1 DX2D relative diff 0 p-value 0 DZDXF relative diff 0 p-value 0 DZDYF relative diff 0 p-value 0 DZS relative diff 0 p-value 0 E relative diff 0 p-value 0 EBIO_API relative diff 0 p-value 0 EBIO_ISO relative diff 0 p-value 0 EMISS relative diff 0 p-value 0 EVAPPROD relative diff 0 p-value 0 F relative diff 0 p-value 0 FCANHFX relative diff 0 p-value 0 FCANQFX relative diff 0 p-value 0 FGIP relative diff 0 p-value 0 FGRNHFX relative diff 0.00536993 p-value -1 FGRNQFX relative diff 0.00536992 p-value -1 FIRE_AREA relative diff 0.00548531 p-value -1 FIRE_HFX relative diff 0 p-value 0 FLINEINT relative diff 0.206871 p-value 1 FLINEINT2 relative diff 0.206871 p-value 1 FMC_G relative diff 0 p-value 0 FMC_GC relative diff 0 p-value 0 FNDWI relative diff 0 p-value 0 FNM relative diff 0 p-value 0 FNP relative diff 0 p-value 0 FUEL_FRAC relative diff 0.0010162 p-value 1 FUEL_FRAC_BURNT relative diff 0.00536994 p-value -1 FUEL_TIME relative diff 0 p-value 0 FWH relative diff 0 p-value 0 FXLAT relative diff 0 p-value 0 FXLONG relative diff 0 p-value 0 FZ0 relative diff 0 p-value 0 F_INT relative diff 0 p-value 0 F_LINEINT relative diff 7.5567e-06 p-value -1 F_LINEINT2 relative diff 7.58448e-06 p-value -1 F_ROS relative diff 7.54535e-06 p-value -1 F_ROS0 relative diff 0 p-value 0 F_ROS11 relative diff 7.82876e-06 p-value -1 F_ROS12 relative diff 8.48361e-06 p-value -1 F_ROS13 relative diff 1.12203e-05 p-value -1 F_ROS21 relative diff 7.9591e-06 p-value -1 F_ROS23 relative diff 0 p-value 0 F_ROS31 relative diff 1.87215e-05 p-value 1 F_ROS32 relative diff 0 p-value 0 F_ROS33 relative diff 0 p-value 0 F_ROSX relative diff 8.26802e-06 p-value 1 F_ROSY relative diff 8.15066e-06 p-value 1 GAMN2O5 relative diff 0 p-value 0 GLW relative diff 0 p-value 0 GRAUPELNC relative diff 0 p-value 0 GRDFLX relative diff 0 p-value 0 GRNHFX relative diff 0.000779961 p-value -0.984227 GRNQFX relative diff 0.000780289 p-value -0.984262 HAILNC relative diff 0 p-value 0 HFX relative diff 2.18506e-05 p-value 1 HFX_FORCE relative diff 0 p-value 0 HFX_FORCE_TEND relative diff 0 p-value 0 HGT relative diff 0 p-value 0 ICN_DIAG relative diff 0 p-value 0 ISCHAP relative diff 0 p-value 0 ISEEDARRAY_SPP_CONV relative diff 0 p-value 0 ISEEDARRAY_SPP_LSM relative diff 0 p-value 0 ISEEDARRAY_SPP_PBL relative diff 0 p-value 0 ISEEDARR_RAND_PERTURB relative diff 0 p-value 0 ISEEDARR_SKEBS relative diff 0 p-value 0 ISEEDARR_SPPT relative diff 0 p-value 0 ISLTYP relative diff 0 p-value 0 ITIMESTEP relative diff 0 p-value 0 IVGTYP relative diff 0 p-value 0 KN2O5 relative diff 0 p-value 0 LAI relative diff 0 p-value 0 LAI_VEGMASK relative diff 0 p-value 0 LAKEMASK relative diff 0 p-value 0 LANDMASK relative diff 0 p-value 0 LFN relative diff 0.000427716 p-value 1 LH relative diff 0 p-value 0 LH_FORCE relative diff 0 p-value 0 LH_FORCE_TEND relative diff 0 p-value 0 LU_INDEX relative diff 0 p-value 0 MAPFAC_M relative diff 0 p-value 0 MAPFAC_MX relative diff 0 p-value 0 MAPFAC_MY relative diff 0 p-value 0 MAPFAC_U relative diff 0 p-value 0 MAPFAC_UX relative diff 0 p-value 0 MAPFAC_UY relative diff 0 p-value 0 MAPFAC_V relative diff 0 p-value 0 MAPFAC_VX relative diff 0 p-value 0 MAPFAC_VY relative diff 0 p-value 0 MAX_MSTFX relative diff 0 p-value 0 MAX_MSTFY relative diff 0 p-value 0 MF_VX_INV relative diff 0 p-value 0 MU relative diff 2.75733e-05 p-value 1 MUB relative diff 0 p-value 0 NC_DIAG relative diff 0 p-value 0 NDVI relative diff 0 p-value 0 NDWI relative diff 0 p-value 0 NEST_POS relative diff 0 p-value 0 NFUEL_CAT relative diff 0 p-value 0 NOAHRES relative diff 0 p-value 0 OC1 relative diff 0.00362439 p-value 1 OC2 relative diff 0.00362081 p-value 1 OLR relative diff 0 p-value 0 P relative diff 0.000253781 p-value 0.999999 P00 relative diff 0 p-value 0 P10 relative diff 0.00362412 p-value 1 P25 relative diff 0.00362403 p-value 1 PB relative diff 0 p-value 0 PBLH relative diff 0 p-value 0 PC relative diff 0 p-value 0 PCB relative diff 0 p-value 0 PH relative diff 7.3812e-05 p-value 1 PHB relative diff 0 p-value 0 PHISC relative diff 0 p-value 0 PHIWC relative diff 0 p-value 0 PHOTR204 relative diff 0 p-value 0 PM10 relative diff 0.00368471 p-value 1 PM2_5_DRY relative diff 0.00368465 p-value 1 PSFC relative diff 1.56718e-07 p-value 0.99865 PV relative diff 0 p-value 0 P_HYD relative diff 1.57532e-07 p-value 1 P_STRAT relative diff 0 p-value 0 P_TOP relative diff 0 p-value 0 Q2 relative diff 3.13488e-06 p-value -1 QFX relative diff 0 p-value 0 QVAPOR relative diff 1.07416e-05 p-value -0.689399 RAINC relative diff 0 p-value 0 RAINNC relative diff 0 p-value 0 RAINPROD relative diff 0 p-value 0 RAINSH relative diff 0 p-value 0 RDN relative diff 0 p-value 0 RDNW relative diff 0 p-value 0 RDX relative diff 0 p-value 0 RDY relative diff 0 p-value 0 RESM relative diff 0 p-value 0 ROUGH_COR relative diff 0 p-value 0 RQVFRTEN relative diff 0.000781598 p-value -0.99537 RTHFRTEN relative diff 0.000781168 p-value -0.995368 R_0 relative diff 0 p-value 0 SAC relative diff 0 p-value 0 SAVE_TOPO_FROM_REAL relative diff 0 p-value 0 SEAICE relative diff 0 p-value 0 SEAS_1 relative diff 0.00362412 p-value 1 SEAS_2 relative diff 0.00362412 p-value 1 SEAS_3 relative diff 0.00362412 p-value 1 SEAS_4 relative diff 0.00362412 p-value 1 SFROFF relative diff 0 p-value 0 SH2O relative diff 0 p-value 0 SHDMAX relative diff 0 p-value 0 SHDMIN relative diff 0 p-value 0 SINALPHA relative diff 0 p-value 0 SMOIS relative diff 0 p-value 0 SMOIS_COR relative diff 0 p-value 0 SNOALB relative diff 0 p-value 0 SNOW relative diff 0 p-value 0 SNOWC relative diff 0 p-value 0 SNOWH relative diff 0 p-value 0 SNOWNC relative diff 0 p-value 0 SNU relative diff 0 p-value 0 SR relative diff 0 p-value 0 SST relative diff 0 p-value 0 SSTSK relative diff 0 p-value 0 SST_INPUT relative diff 0 p-value 0 SWDOWN relative diff 0 p-value 0 SWNORM relative diff 0 p-value 0 T relative diff 0.000198042 p-value -0.999647 T00 relative diff 0 p-value 0 T2 relative diff 2.40056e-06 p-value 0.782484 TH2 relative diff 2.39755e-06 p-value 0.503594 THIS_IS_AN_IDEAL_RUN relative diff 0 p-value 0 THM relative diff 0.000153512 p-value -0.995398 TIGN_G relative diff 2.05477e-05 p-value 1 TIGN_IN relative diff 0 p-value 0 TISO relative diff 0 p-value 0 TLP relative diff 0 p-value 0 TLP_STRAT relative diff 0 p-value 0 TMN relative diff 0 p-value 0 TSK relative diff 0 p-value 0 TSK_FORCE relative diff 0 p-value 0 TSK_FORCE_TEND relative diff 0 p-value 0 TSLB relative diff 0 p-value 0 Times relative diff 0 p-value 0 U relative diff 8.32584e-06 p-value 1 U10 relative diff 1.99208e-05 p-value 1 UAH relative diff 0 p-value 0 UDROFF relative diff 0 p-value 0 UF relative diff 6.76522e-06 p-value 1 UNIT_FXLAT relative diff 0 p-value 0 UNIT_FXLONG relative diff 0 p-value 0 UST relative diff 2.21928e-05 p-value -1 UST_T relative diff 0 p-value 0 V relative diff 1.17025e-05 p-value -1 V10 relative diff 2.53694e-05 p-value 0.961354 VAH relative diff 0 p-value 0 VAR relative diff 0 p-value 0 VAR_SSO relative diff 0 p-value 0 VEGFRA relative diff 0 p-value 0 VF relative diff 6.91935e-06 p-value 1 W relative diff 1.63983e-05 p-value 1 WZ0 relative diff 0 p-value 0 XICEM relative diff 0 p-value 0 XLAND relative diff 0 p-value 0 XLAT relative diff 0 p-value 0 XLAT_U relative diff 0 p-value 0 XLAT_V relative diff 0 p-value 0 XLONG relative diff 0 p-value 0 XLONG_U relative diff 0 p-value 0 XLONG_V relative diff 0 p-value 0 XTIME relative diff 0 p-value 0 YCLNO2 relative diff 0 p-value 0 ZETATOP relative diff 0 p-value 0 ZNU relative diff 0 p-value 0 ZNW relative diff 0 p-value 0 ZS relative diff 0 p-value 0 ZSF relative diff 0 p-value 0 dms relative diff 0.00362428 p-value 1 msa relative diff 0.00362428 p-value 1 so2 relative diff 0.00362404 p-value 1 sulf relative diff 0.0036242 p-value 1

janmandel commented 4 years ago

WRF-SFIRE based on WRF v4.2 is ready both with and without compilation with WRF-Chem. I fixed few bugs and improved WRF-Chem support, when compiled witjh WRF-Chem, you can also add tracer_opt=0,1,and 2. Tested on ideal examples including test/em_fire/chem, compiled with WRF_CHEM=1 and on real case em_rxcadre from Derek, without WRF_CHEM.

I would like to test on a real Chem case but I do not know how to turn em_rxcadre into a chem simulation. I tried to add chem section to the namelist from elsewhere but do not know what I am doing.

Ideal cases run fine but slightly different results with different number of processors in chem variables are creeping in, starting at at rounding error level. It would be great to test original WRF-Chem from wrf-model/WRF without any modification if it does that too.

The code is in https://github.com/openwfm/WRF-SFIRE master, at the moment at c178451c67bd8498fa4437f825ec01942f2bc597 Not closing https://github.com/openwfm/WRF-SFIRE/issues/16 yet until we can run some real chem case.

Any help with testing appreciated.

janmandel commented 4 years ago

Testing in progress if the results on real are similar as 3.4 (can't expect to be identical too much has changed since 3.4), but so far so good. Also #35 turned out to be false alarm.

adamk0 commented 3 years ago

Yes, it does. But you need to link the proper namelist.fire_emissions file compatible with RADM2 (namelist.fire_emissions.RADM2).

Adam Kochanski

Assistant Professor Department of Meteorology and Climate Science Wildfire Interdisciplinary Research Center San Jose State University

From: Zongrun @.> Reply-To: openwfm/WRF-SFIRE @.> Date: Thursday, September 16, 2021 at 7:29 AM To: openwfm/WRF-SFIRE @.> Cc: Subscribed @.> Subject: Re: [openwfm/WRF-SFIRE] Carry over WRF-Fire-Chem (#16)

Hi Dr. Mandel, does wrf-chem supports RAMD2? It seems like the module crashed: crash: read_emissions_table: increase max_chem How can I fix it? GOCART works well while RAMD2 cannot run. The version I used is: git clone git://github.com/openwfm/WRF-SFIRE.git

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/openwfm/WRF-SFIRE/issues/16#issuecomment-920954184, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AA6K3WPLNPOFY45ENVRM2HTUCH5MHANCNFSM4NNYXI4A. Triage notifications on the go with GitHub Mobile for iOShttps://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Androidhttps://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.