Getting very large rsl files
500 views
Skip to first unread message
elham baranizadeh
Aug 14, 2018, 3:53:58 PM8/14/18
to wrf-chem-bio_emiss
Hi all,
I'm running wrf-chem with MEGAN and NEI emissions. However, I'm getting very large rsl files (>100GB after a few days are simulated). I only post this in this group beacuse the rsl tail is as follows (so something to do with MEGAN or radiation scheme):
calculate MEGAN emissions at ktau, gmtp, tmidh = 3032 12.0000000 12.1260004
Timing for main: time 2016-05-26_12:07:40 on domain 2: 0.51121 elapsed seconds
calculate MEGAN emissions at ktau, gmtp, tmidh = 3033 12.0000000 12.1300001
Timing for main: time 2016-05-26_12:07:55 on domain 2: 0.49534 elapsed seconds
calculate MEGAN emissions at ktau, gmtp, tmidh = 3034 12.0000000 12.1339998
Timing for main: time 2016-05-26_12:08:09 on domain 2: 0.51136 elapsed seconds
calculate MEGAN emissions at ktau, gmtp, tmidh = 3035 12.0000000 12.1379995
Timing for main: time 2016-05-26_12:08:24 on domain 2: 0.51153 elapsed seconds
Timing for main: time 2016-05-26_12:08:24 on domain 1: 3.37041 elapsed seconds
calculate MEGAN emissions at ktau, gmtp, tmidh = 608 12.0000000 12.1499996
Please know that the model still keeps running but with these large rsl files I run out of space. Any idea why this is?
Best,
Elham
Stacy Walters
Aug 14, 2018, 5:18:30 PM8/14/18
to elham baranizadeh, wrf-chem-bio_emiss
Elham,
What is the value of the namelist.input variable "debug_level"?
Stacy
elham baranizadeh
Aug 14, 2018, 5:28:54 PM8/14/18
to Stacy Walters, wrf-chem-bio_emiss
Stacy,
That was my first guess too- but it's set to 0. Then I looked at the rsl files and I think it’s all the “Diagnostics” output that’s filling up your rsl files:
WARNING: Large total sw optical depth of ******** at point i,j,nb= 50 7 7
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1****************************************
2****************************************
3****************************************
4****************************************
5****************************************
6****************************************
7****************************************
8****************************************
9****************************************
10****************************************
11****************************************
12****************************************
13****************************************
14****************************************
15****************************************
16****************************************
17****************************************
18****************************************
19****************************************
20****************************************
21****************************************
22****************************************
23****************************************
24 150.25 108.64 59.84 21.60 16.86
25 132.07 95.60 52.71 19.04 14.87
26 124.29 90.00 49.64 17.93 14.01
27 49.37 35.75 19.72 7.12 5.57
28 0.04 0.03 0.01 0.01 0.00
29 0.00 0.00 0.00 0.00 0.00
30 0.00 0.00 0.00 0.00 0.00
Then, I just realized that my wrfbiochemi files have been created from slightly a different wrfinput_d0<n> files, and tried to create them again but I'm getting this error:
Reading megan2 bio emiss file TAS.nc
No such file or directory
any idea what's the error for?
-Elham
--
You received this message because you are subscribed to the Google Groups "wrf-chem-bio_emiss" group.
To unsubscribe from this group and stop receiving emails from it, send an email to wrf-chem-bio_em...@ucar.edu.
Stacy Walters
Aug 14, 2018, 5:42:07 PM8/14/18
to elham baranizadeh, wrf-chem-bio_emiss
Elham,
When running the megain_bio_emiss utility you have to download
basic input data files first. I'm not sure but believe they should be
on the ACOM wrfchem utilities website.
Stacy
To unsubscribe from this group and stop receiving emails from it, send an email to wrf-chem-bio_emiss+unsub...@ucar.edu.
Ankan Sarkar
Feb 6, 2023, 2:23:06 AM2/6/23
to wrf-chem-bio_emiss, st...@ucar.edu, wrf-chem-bio_emiss, elham.ba...@gmail.com, Gabriele Pfister
Dear all,
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1****************************************
2****************************************
3****************************************
4****************************************
5****************************************
6****************************************
7****************************************
8****************************************
9****************************************
10****************************************
11****************************************
12****************************************
13****************************************
14****************************************
15****************************************
16****************************************
17****************************************
18****************************************
19****************************************
20****************************************
21****************************************
22****************************************
23****************************************
I am also having the same issue regarding the large rsl files after running WRF-Chem with the MOZART-MOSAIC option (chem_opt=202). Similar warnings are also printed in my rsl.error files, like Elham shared in the previous thread. The warning messages are as follows:
-------------------------
-------------------------
WARNING: Large total sw optical depth of ******** at point i,j,nb= 27 42 5
-------------------------
WARNING: Large total sw optical depth of ******** at point i,j,nb= 27 42 5
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1****************************************
2****************************************
3****************************************
4****************************************
5****************************************
6****************************************
7****************************************
8****************************************
9****************************************
10****************************************
11****************************************
12****************************************
13****************************************
14****************************************
15****************************************
16****************************************
17****************************************
18****************************************
19****************************************
20****************************************
21****************************************
22****************************************
23****************************************
24****************************************
25****************************************
26****************************************
27****************************************
28****************************************
29****************************************
30****************************************
31****************************************
32****************************************
33****************************************
34****************************************
Diagnostics 2: k, gaer300, gaer400, gaer600, gaer999
1 0.82 0.82 0.81 0.80
2 0.82 0.82 0.81 0.80
3 0.82 0.82 0.81 0.80
4 0.82 0.82 0.81 0.80
5 0.82 0.82 0.81 0.80
6 0.82 0.82 0.81 0.80
7 0.82 0.82 0.81 0.80
8 0.82 0.82 0.81 0.80
9 0.82 0.82 0.81 0.80
10 0.82 0.82 0.81 0.80
11 0.82 0.82 0.81 0.80
12 0.82 0.82 0.81 0.80
13 0.82 0.82 0.81 0.80
14 0.82 0.82 0.81 0.80
15 0.82 0.82 0.81 0.80
16 0.81 0.82 0.80 0.77
17 0.81 0.82 0.80 0.77
18 0.81 0.82 0.80 0.77
19 0.81 0.81 0.79 0.75
20 0.83 0.81 0.81 0.78
21 0.82 0.81 0.79 0.76
22 0.82 0.81 0.81 0.77
23 0.82 0.82 0.81 0.80
24 0.83 0.81 0.81 0.78
25 0.82 0.82 0.81 0.80
26 0.82 0.82 0.81 0.80
27 0.82 0.82 0.81 0.80
28 0.82 0.82 0.81 0.80
29 0.82 0.82 0.81 0.80
30 0.82 0.82 0.81 0.80
31 0.82 0.82 0.81 0.80
32 0.82 0.82 0.81 0.80
33 0.82 0.82 0.81 0.80
34 0.82 0.82 0.81 0.80
Diagnostics 3: k, waer300, waer400, waer600, waer999
1 1.00 1.00 1.00 1.00
2 1.00 1.00 1.00 1.00
3 1.00 1.00 1.00 1.00
4 1.00 1.00 1.00 1.00
5 1.00 1.00 1.00 1.00
6 1.00 1.00 1.00 1.00
7 1.00 1.00 1.00 1.00
8 1.00 1.00 1.00 1.00
9 1.00 1.00 1.00 1.00
10 1.00 1.00 1.00 1.00
11 1.00 1.00 1.00 1.00
12 1.00 1.00 1.00 1.00
13 1.00 1.00 1.00 1.00
14 1.00 1.00 1.00 1.00
15 1.00 1.00 1.00 1.00
16 1.00 1.00 1.00 1.00
17 1.00 1.00 1.00 1.00
18 1.00 1.00 1.00 1.00
19 1.00 1.00 1.00 1.00
20 1.00 1.00 1.00 1.00
21 1.00 1.00 1.00 1.00
22 1.00 1.00 1.00 1.00
23 1.00 1.00 1.00 1.00
24 1.00 1.00 1.00 1.00
25 1.00 1.00 1.00 1.00
26 1.00 1.00 1.00 1.00
27 1.00 1.00 1.00 1.00
28 1.00 1.00 1.00 1.00
29 1.00 1.00 1.00 1.00
30 1.00 1.00 1.00 1.00
31 1.00 1.00 1.00 1.00
32 1.00 1.00 1.00 1.00
33 1.00 1.00 1.00 1.00
34 1.00 1.00 1.00 1.00
Diagnostics 4: k, ssaal, asyal, taual
0 1.00 0.75 0.00
1 1.00 0.75 0.00
2 1.00 0.75 0.02
3 1.00 0.75 2.99
4 1.00 0.75 2.90
5 1.00 0.75 0.08
6 1.00 0.75 0.00
7 1.00 0.75 0.00
8 1.00 0.75 0.00
9 1.00 0.00 0.00
10 1.00 0.77********
11 1.00 0.77********
12 1.00 0.77********
13 1.00 0.77********
14 1.00 0.77********
15 1.00 0.77********
16 1.00 0.77********
17 1.00 0.77********
18 1.00 0.77********
19 1.00 0.77********
20 1.00 0.77********
21 1.00 0.77********
22 1.00 0.77********
23 1.00 0.77********
24 1.00 0.77********
25 1.00 0.74********
26 1.00 0.74********
27 1.00 0.75********
28 1.00 0.71********
29 1.00 0.80********
30 1.00 0.72********
31 1.00 0.78********
32 1.00 0.77********
33 1.00 0.80********
34 1.00 0.77********
-------------------------
-------------------------
In my case, too, wrf.exe is still running. So, why is it happening? Can these warnings cause any trouble with the output? If so, then how can I solve this issue? Did anyone solve it? I checked my biogenic emission input files, and they are generated from the same wrfiput files. So, I think this is not the problem that Elham had. So, what could be the reason? I suspect this could be due to the radiation options. I set the following radiation options in my namelist.input file:
25****************************************
26****************************************
27****************************************
28****************************************
29****************************************
30****************************************
31****************************************
32****************************************
33****************************************
34****************************************
Diagnostics 2: k, gaer300, gaer400, gaer600, gaer999
1 0.82 0.82 0.81 0.80
2 0.82 0.82 0.81 0.80
3 0.82 0.82 0.81 0.80
4 0.82 0.82 0.81 0.80
5 0.82 0.82 0.81 0.80
6 0.82 0.82 0.81 0.80
7 0.82 0.82 0.81 0.80
8 0.82 0.82 0.81 0.80
9 0.82 0.82 0.81 0.80
10 0.82 0.82 0.81 0.80
11 0.82 0.82 0.81 0.80
12 0.82 0.82 0.81 0.80
13 0.82 0.82 0.81 0.80
14 0.82 0.82 0.81 0.80
15 0.82 0.82 0.81 0.80
16 0.81 0.82 0.80 0.77
17 0.81 0.82 0.80 0.77
18 0.81 0.82 0.80 0.77
19 0.81 0.81 0.79 0.75
20 0.83 0.81 0.81 0.78
21 0.82 0.81 0.79 0.76
22 0.82 0.81 0.81 0.77
23 0.82 0.82 0.81 0.80
24 0.83 0.81 0.81 0.78
25 0.82 0.82 0.81 0.80
26 0.82 0.82 0.81 0.80
27 0.82 0.82 0.81 0.80
28 0.82 0.82 0.81 0.80
29 0.82 0.82 0.81 0.80
30 0.82 0.82 0.81 0.80
31 0.82 0.82 0.81 0.80
32 0.82 0.82 0.81 0.80
33 0.82 0.82 0.81 0.80
34 0.82 0.82 0.81 0.80
Diagnostics 3: k, waer300, waer400, waer600, waer999
1 1.00 1.00 1.00 1.00
2 1.00 1.00 1.00 1.00
3 1.00 1.00 1.00 1.00
4 1.00 1.00 1.00 1.00
5 1.00 1.00 1.00 1.00
6 1.00 1.00 1.00 1.00
7 1.00 1.00 1.00 1.00
8 1.00 1.00 1.00 1.00
9 1.00 1.00 1.00 1.00
10 1.00 1.00 1.00 1.00
11 1.00 1.00 1.00 1.00
12 1.00 1.00 1.00 1.00
13 1.00 1.00 1.00 1.00
14 1.00 1.00 1.00 1.00
15 1.00 1.00 1.00 1.00
16 1.00 1.00 1.00 1.00
17 1.00 1.00 1.00 1.00
18 1.00 1.00 1.00 1.00
19 1.00 1.00 1.00 1.00
20 1.00 1.00 1.00 1.00
21 1.00 1.00 1.00 1.00
22 1.00 1.00 1.00 1.00
23 1.00 1.00 1.00 1.00
24 1.00 1.00 1.00 1.00
25 1.00 1.00 1.00 1.00
26 1.00 1.00 1.00 1.00
27 1.00 1.00 1.00 1.00
28 1.00 1.00 1.00 1.00
29 1.00 1.00 1.00 1.00
30 1.00 1.00 1.00 1.00
31 1.00 1.00 1.00 1.00
32 1.00 1.00 1.00 1.00
33 1.00 1.00 1.00 1.00
34 1.00 1.00 1.00 1.00
Diagnostics 4: k, ssaal, asyal, taual
0 1.00 0.75 0.00
1 1.00 0.75 0.00
2 1.00 0.75 0.02
3 1.00 0.75 2.99
4 1.00 0.75 2.90
5 1.00 0.75 0.08
6 1.00 0.75 0.00
7 1.00 0.75 0.00
8 1.00 0.75 0.00
9 1.00 0.00 0.00
10 1.00 0.77********
11 1.00 0.77********
12 1.00 0.77********
13 1.00 0.77********
14 1.00 0.77********
15 1.00 0.77********
16 1.00 0.77********
17 1.00 0.77********
18 1.00 0.77********
19 1.00 0.77********
20 1.00 0.77********
21 1.00 0.77********
22 1.00 0.77********
23 1.00 0.77********
24 1.00 0.77********
25 1.00 0.74********
26 1.00 0.74********
27 1.00 0.75********
28 1.00 0.71********
29 1.00 0.80********
30 1.00 0.72********
31 1.00 0.78********
32 1.00 0.77********
33 1.00 0.80********
34 1.00 0.77********
-------------------------
-------------------------
In my case, too, wrf.exe is still running. So, why is it happening? Can these warnings cause any trouble with the output? If so, then how can I solve this issue? Did anyone solve it? I checked my biogenic emission input files, and they are generated from the same wrfiput files. So, I think this is not the problem that Elham had. So, what could be the reason? I suspect this could be due to the radiation options. I set the following radiation options in my namelist.input file:
&physics
ra_lw_physics = 4, 4,
ra_sw_physics = 4, 4,
radt = 30, 30,
I learned that these warnings are coming every 30 minutes, and then again, the printing in the output files has been going on. So is it because of 'radt' option? I have gone through namelist.input: Best Practices from WRF USERS PAGE (the https://www2.mmm.ucar.edu/wrf/users/namelist_best_prac_wrf.html), and found that it is recommended to use 1 minute per km of dx (e.g., set to 10 for a 10km parent domain). However, I set here 30, which is 2*dx (dx = 15 km for the parent domain, in my case). So, is it because of this? Or is anything else wrong? Can anyone please guide me in this regard? That will be really helpful for me. Thank you in advance.
ra_lw_physics = 4, 4,
ra_sw_physics = 4, 4,
radt = 30, 30,
I learned that these warnings are coming every 30 minutes, and then again, the printing in the output files has been going on. So is it because of 'radt' option? I have gone through namelist.input: Best Practices from WRF USERS PAGE (the https://www2.mmm.ucar.edu/wrf/users/namelist_best_prac_wrf.html), and found that it is recommended to use 1 minute per km of dx (e.g., set to 10 for a 10km parent domain). However, I set here 30, which is 2*dx (dx = 15 km for the parent domain, in my case). So, is it because of this? Or is anything else wrong? Can anyone please guide me in this regard? That will be really helpful for me. Thank you in advance.
With regards,
Ankan
Mary Barth
Feb 6, 2023, 10:07:47 PM2/6/23
to Ankan Sarkar, wrf-chem-bio_emiss, st...@ucar.edu, elham.ba...@gmail.com, Gabriele Pfister
Ankan,
As noted by the earlier email, these are printouts from diagnostic calculations, specifically from the optical driver. The first part is trying to print the aerosol optical depth (tauaer) but the values are so large it doesn't fit the format of the print statement. There is something wrong -- way too much aerosol occurring somewhere. (i=27, j=42)
I suggest reviewing your wrfinput file and/or the last wrfout
file that is generated before these warnings start. Look at all
the aerosol concentrations to see if there is something
anomalously wrong.
Also look at bug reports. I think some things have been cleaned up in recent WRF-chem versions, but your version may not have the correction. One example that we came across is with biomass burning: See https://forum.mmm.ucar.edu/threads/plumerise-bug.9561/
Good luck,
Mary
-- ^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^ Mary Barth Phone: 303-497-8186 Senior Scientist email: bar...@ucar.edu National Center for Atmospheric Research P.O. Box 3000 Boulder, CO 80307 https://staff.ucar.edu/users/barthm SIMA Web Site: http://sima.ucar.edu MUSICA Web Site: http://www2.acom.ucar.edu/sections/multi-scale-infrastructure-chemistry-modeling-musica DC3 Web Site: http://www2.acom.ucar.edu/dc3/ +--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+
Reply all
Reply to author
Forward
Message has been deleted
0 new messages