WRF-Fire-Chem for real case is crashed without nesting about wrf-sfire

[janmandel]

@Porkking
Can you please provide the inputs and complete file rsl.error.0000?
Also, the files provided at the google drive link are for one test case with d01 only, while you mention several cases, with multiple domains.
You mention "the information is hidden when the fire model is in a daughter domain". I think you mean you do not see a print the total emissions. A low value of fire_print_msg in that domain would do that.

[zli867]

Hi Dr. Mandel,
For the one domain case:
rsl.error.0000.txt

For the nesting case, you may correct. I set fire_print_msg=1. I could increase the value and have a try.
I upload the nesting case namelist and rsl.error.0000:
rsl.error.0000_nesting.txt
namelist_nesting.txt
I can also provide all data to reproduce a nesting WRF-Fire-Chem case: https://drive.google.com/drive/folders/1fTeuU1ZfY2y0KzHiDXHLg9tWQCHY6xq9?usp=sharing

Thanks,
Zongrun

[janmandel]

@Porkking from the commit line in the rsl.error.0000.txt file, I see the source code was changed and I do not know what it is, I am sorry. Also, the printout rsl.error.0000.txt shows more species than my coupling code, it does not look like chem_opt=112 in namelist.fire_emisions. That's why the "increase max_chem" error message. You can try to do what it asks, or make sure you have chem_opt=112 in namelist.input and the code has no changes against github. It seems there are more species than when I wrote the coupling years ago. I'll have to review the chemical species and models supported soon anyway #51.

The case in the google drive link is not nested so I am assuming that's where the rsl.error.0000.txt is from.

It would be best to keep all files for each run together in a separate directory, otherwise I am not sure what is what.

In the namelist_nesting.txt, there is missing fire_print_msg for domain 2 so it is taken as 0 by default.

[janmandel]

@Porkking can I reproduce the case that gave you the crash: read_emissions_table: increase max_chem error please? There is not much I can do without that. Thanks, Jan

[zli867]

Hi Dr. Mandel,
I think this directory provide all inputs you may need for running the case:
https://drive.google.com/drive/folders/1fTeuU1ZfY2y0KzHiDXHLg9tWQCHY6xq9?usp=sharing.
wrfchemi_d*: anthropogenic emission
wrfchemi_gocart*: GOCART background
wrfbiochemi*: biogenic emission (MEGAN)
wrf_season* and exo_*: required for MOZCART
wrfinput: initial condition (I generate it from mother domain by ndown)
wrfbdy: boundary condition (I generate it from mother domain by ndown)
I also provide the all namelist I used. If there is anything I have left, please let me know.

Thanks,
Zongrun

[janmandel]

@Porkking could you please add tohe directory the file rsl.error.0000 from a run with exactly these files to the directory? I need to see what you see. Among other things, it tells me which version (git commit) you used. The file rsl.error.0000.txt you have provided earlier indicates a problem, it shows at the top
git commit b7ae876 1 file changed, 529 insertions(+)
and I do not know what exactly changed. If I see a commit hash with no change then I know what you have.
Thanks, Jan

[janmandel]

@Porkking I am sorry, it's not all files. After copying everything from ./run, it will still crash in WRF (not sfire) with
Possibly missing file for = auxinput5

[zli867]

Hi Dr. Mandel,
"git commit b7ae876 1 file changed, 529 insertions(+)"
Do you mean I have change the code for the model? If so, I may change the optical threshold for WRF (I cannot remember clearly whether I change it or not for this directory). I change it since optical value exceed the threshold when PM2.5 is really high and it kills the WRF job.

"it will still crash in WRF (not sfire) with Possibly missing file for = auxinput5"
I am not sure why it happens. auxinput5 is for anthropogenic emission. I use hourly emission and the name is wrfchemi_d*. I think I have uploaded all the emission files from 04-06_00:00 to 04-07_00:00 (total 25 files). By the way, it should be WRF-Chem with KPP option. I can also upload a environment variable setting when I install the model: https://drive.google.com/file/d/1LTAzKuWxEiwQXeNxXn3wP1U5J0sScurL/view?usp=sharing

[janmandel]

@Porkking I need to be able to run the code with the same files you have, including the same source code, so that I can help you. Can you rerun your case with WRF-SFIRE without any changes and provide all input files and rsl.error.0000 so that I see what you see, please?

[janmandel]

@Porkking Also, you can commit your changes to the source code and provide the modified git repository with your files. I'd like to see the same commit hash in git log and in rsl.error.0000 with no changes.

[zli867]

@janmandel After checking by the git diff, I think I have not change the optical threshold in the version I currently used. It seems the 1 file changed is the namelist.input. Anyway, I will build the WRF-Fire without any changes and rerun my case. I will update when it is done.
git diff.txt

[zli867]

@janmandel it seems you have merged current WRF-SFire to WRFv4.4 while I met problem on installing by using my previous setting. I have posted the issue in the WRF forum and what I can do is stick with WRFv4.2 and do the test. I will try to download all files I have provided before and try to run it with WRFv4.2 again.

[zli867]

@janmandel Dr. Mandel, you are correct. I have seen the same problem when I use the files provided in google drive to run WRF-Chem. I need to reproduce this work in WRFv4.2. I will share the data/results after reproducing.

[zli867]

@janmandel Hi Dr. Mandel, I have uploaded the input data for one domain case of WRF-Chem: https://drive.google.com/file/d/18gE-BmBcKC163cabJ-dPNlxDCrLgdaF1/view?usp=sharing. You can try the current data to set up the WRF-Chem simulation. I guess if you switch ifire from 0 to 1, you will find the error:
crash: read_emissions_table: increase max_chem.
Google drive changes the input files if I upload all files directly. So, I compressed the input data. At least, I can run the WRF-Chem model by downloading the .tar.gz file. Hope this time the WRF-Chem can work successfully. If there is any problems, please let me know. Thanks, Zongrun.

[janmandel]

@Porkking sorry about the delay, the semester started. I'll look into that asap.