Error: mieaer /refr/ outside range 1e-3 - 10 refr= -0.10931E+02

[Sreyashi Debnath]

Hello members,

I am getting an error: mieaer /refr/ outside range 1e-3 - 10 refr=    -0.10931E+02, while running WRF-Chem.

Kindly help me out in solving the problem.

Thanking you in anticipation.

Regards,

Sreyashi Debnath

IITM, Pune.

[Sreyashi Debnath]

Hello.!!

I have gone through the link ( https://ruc.noaa.gov/wrf/wrf-chem/FAQ.htm) you have shared and recompiled the model as mentioned over there but still I am getting the error. 

I have noticed one thing, that unlike the example provided in the WRF-Chem FAQ page i.e. 'mieaer /refr/ outside range 1e-3 - 10 refr= 0.10003E+02', the value of 'refr' in my case is always coming negative: e.g.  mieaer /refr/ outside range 1e-3 - 10 refr=    -0.10293E+02 (notice the negative sign before 0.10293E+02).

Can you please suggest what should be done in this case?

Regards,

Sreyashi Debnath

IITM, Pune

On Sat, May 18, 2019 at 4:34 PM Rama Krishna <k.b....@gmail.com> wrote:

HI Sreyashi,

For your problem this may be the solution plz find this

The error is due to the calculated refractive index (10.931) slightly exceeding the preset upper limit (10). The solution is to increase the upper limit a little bit. For example, to revise the upper limit into 15, make the following modification in WRFV3/chem/module_optical_averaging.F, and recompile the model.

Line 4678: if(abs(refr).gt.10.0.or.abs(refr).le.0.001)then

--> if(abs(refr).gt.15.0.or.abs(refr).le.0.001)then

Please also find the link below last answer

https://ruc.noaa.gov/wrf/wrf-chem/FAQ.htm

--
You received this message because you are subscribed to the Google Groups "wrf-chem" group.
To unsubscribe from this group and stop receiving emails from it, send an email to wrf-chem+u...@ucar.edu.

--

 K Rama krishna.

[Rajesh Kumar]

Hi Sreyashi,

Could you please make sure that you don't have any negative PM2.5 emissions in any your domain? One of the students also ran in this problem earlier and we found out it was related to negative PM2.5 emissions in either wrfchemi* or wrffirechemi* files. The negative emissions were resulting from subtraction of BC and OC emissions from total PM2.5 emissions in either anthro_emis or fire_emis preprocessing? This should not happen in theory but I have seen such incidents at isolated grid point before. If there are negative emissions, just set them to zero otherwise they lead to negative concentrations in the model and consequently negative refractive indices. 

Best,

Rajesh

--

Rajesh Kumar, Project Scientist II

Research Applications Laboratory 

National Center for Atmospheric Research 

E-mail: rku...@ucar.edu

Phone: +1-303-497-1413

Web: https://staff.ucar.edu/users/rkumar

[Sreyashi Debnath]

Hello Rajesh Sir,

Thank you for the information. Actually, I did subtract BC and OC emission from total PM2.5 emissions in anthro_emis preprocessing since I came to know that to eliminate over estimation of BC and OC, these two emissions should be subtracted from total PM2.5 emission.

So, could you please suggest what should be done to minimise the over estimation if not subtracting the two species from total emission.??

Regards,

Sreyashi

[Rajesh Kumar]

BC and OC should still be subtracted and in my experience, you will have very small negative values at 1 or 2 grid points. Just set the negative values to zero.

Sent from my iPhone

[Sreyashi Debnath]

Alright.!!

Thanks a lot sir, will try to do accordingly.

Regards,

Sreyashi

[Caterina Mogno]

Dear Rajesh,

regarding this topic, does  the EDGAR-HTAP inputs provided in the WRFchem tools for the Community page already include OC and BC substraction from PM2.5, and PM2.5 substraction from PM10?

Thanks in advance for help!

Kind regards,

Caterina

To unsubscribe from this group and stop receiving emails from it, send an email to wrf-...@ucar.edu.

--

 K Rama krishna.

[Rajesh Kumar]

Dear Caterina,

We did not subtract OC and BC contributions from EDGAR-HTAP inputs provided on the website. However, we do that in the namelist for anthro_emis as follows: 

E_PM2_5 -> PM25 + -1.0*OC + -1.0*BC

Note that in some cases, we have seen that this subtraction leads to very small negative values. It should not happen in theory but occurs mathematically at very few grid boxes over Asia. The negative emissions lead to negative concentrations of the species and blow the model away. So, I'll suggest to check your PM2_5 and PM10 emissions in your wrfchemi* files and make sure that you don't have any negative values. 

Best,

Rajesh

[Caterina Mogno]

Dear Rajesh,

thanks for your quick reply.  In the example namelist for anthro_emis provided for MOZART-MOSAIC ( INP-Examples/ MOZART-MOSAIC.inp)

the emis_map states the following:

'PM25I(a)->0.2*PM2.5',
'PM25J(a)->0.8*PM2.5',

'PM_10(a)->PM10'

Based on what you said, the correct mapping would be:

'PM25I(a)->0.2*PM2.5' - 0.2*BC - 0.2*OC' ,
'PM25J(a)->0.8*PM2.5 - 0.8*BC - 0.8*OC ',

'PM_10(a)->PM10 - PM2.5'

Is that correct both for PM2.5 and PM10?

Thank you very much in advance.

Kind regards,

Caterina

[Rajesh Kumar]

Dear Caterina,

Yes, that is correct. 

Best,

Rajesh

[Gonzalo Ferrada]

Hello everyone,

I am experiencing the same error in version 3.9.1:

calculate MEGAN emissions at ktau, gmtp, tmidh =         1296   5.000000       5.979167

-------------- FATAL CALLED ---------------

FATAL CALLED FROM FILE:  <stdin>  LINE:      76

mieaer /refr/ outside range 1e-3 - 10 refr=    -0.24828E+02

-------------------------------------------

Similar as someone mentioned before, my error is because it generates a negative value. This occurs when the simulation is well advanced (53rd hour).

I have checked all my inputs (chemi, firechemi and biochemi) and none of them have any negative value in their species.

However, I see some aerosols with negative values in my wrfbdy file... is that the error? or does a negative value in wrfbdy means an outflow of aerosols? 

Thanks,

Gonzalo.

On Friday, April 3, 2020 at 11:21:37 AM UTC-5, rkumar wrote:

Dear Caterina,

Yes, that is correct. 

Best,

Rajesh

On Fri, Apr 3, 2020 at 10:19 AM Caterina Mogno <cateri...@gmail.com> wrote:

Dear Rajesh,

thanks for your quick reply.  In the example namelist for anthro_emis provided for MOZART-MOSAIC ( INP-Examples/ MOZART-MOSAIC.inp)

the emis_map states the following:

'PM25I(a)->0.2*PM2.5',
'PM25J(a)->0.8*PM2.5',

'PM_10(a)->PM10'

Based on what you said, the correct mapping would be:

'PM25I(a)->0.2*PM2.5' - 0.2*BC - 0.2*OC' ,
'PM25J(a)->0.8*PM2.5 - 0.8*BC - 0.8*OC ',

'PM_10(a)->PM10 - PM2.5'

Is that correct both for PM2.5 and PM10?

Thank you very much in advance.

Kind regards,

Caterina

On Wed, 1 Apr 2020 at 17:54, Rajesh Kumar <rku...@ucar.edu> wrote:

Dear Caterina,

We did not subtract OC and BC contributions from EDGAR-HTAP inputs provided on the website. However, we do that in the namelist for anthro_emis as follows: 

E_PM2_5 -> PM25 + -1.0*OC + -1.0*BC

Note that in some cases, we have seen that this subtraction leads to very small negative values. It should not happen in theory but occurs mathematically at very few grid boxes over Asia. The negative emissions lead to negative concentrations of the species and blow the model away. So, I'll suggest to check your PM2_5 and PM10 emissions in your wrfchemi* files and make sure that you don't have any negative values. 

Best,

Rajesh

[Mary Barth]

Hello Gonzalo,

I'm not sure if anyone replied. The wrfbdy file contains mixing ratios or concentrations of species (e.g. so4a1_BXS) and tendencies (so4a1_BTXS). The latter are temporal tendencies and could be negative. But the former are concentrations and should be positive.

-- 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  
DC3 Web Site:  http://www2.acom.ucar.edu/dc3/
+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+

[Gonzalo Ferrada]

Thanks Mary!

I resolved the issue by playing with the chemdt. I had chemdt = 2.5 (minutes) and, every time the error happened, I switched to chemdt = 5.0 and finished those simulations with no errors.

Cheers,

--
Gonzalo A. Ferrada, Ph.D. (c)
Center for Global and Regional Environmental Research
Department of Chemical and Biochemical Engineering
The University of Iowa

http://gcarmichgroup.com/people/gonzalo.html 

[Lin Lu]

Hi all,

Hi all,

According to the discussions and answers in this topic, do we need to modify the mozart-mosiac.inp file like this: 

'PM25I(a)->0.2*PM2.5',
'PM25J(a)->0.8*PM2.5',

'PM_10(a)->PM10'

Based on what you said, the correct mapping would be:

'PM25I(a)->0.2*PM2.5' - 0.2*BC - 0.2*OC' ,
'PM25J(a)->0.8*PM2.5 - 0.8*BC - 0.8*OC ',

'PM_10(a)->PM10 - PM2.5'

Since, the mapping in the examples provided on the download page do not include - 0.2*BC - 0.2*OC'  from PM2.5

Is this necessary for all the users? If not, why?

I'd be grateful for the clarification provided.

Thanks a lot in advance.

Regards, Lin 

On Wednesday, April 1, 2020 at 5:54:28 PM UTC+1, rkumar wrote:

Dear Caterina,

We did not subtract OC and BC contributions from EDGAR-HTAP inputs provided on the website. However, we do that in the namelist for anthro_emis as follows: 

E_PM2_5 -> PM25 + -1.0*OC + -1.0*BC

Note that in some cases, we have seen that this subtraction leads to very small negative values. It should not happen in theory but occurs mathematically at very few grid boxes over Asia. The negative emissions lead to negative concentrations of the species and blow the model away. So, I'll suggest to check your PM2_5 and PM10 emissions in your wrfchemi* files and make sure that you don't have any negative values. 

Best,

Rajesh

[Leo Matak]

Dear all,

I am encountering the same issue ( mieaer /refr/ outside range 1e-3 - 10 refr=    -0.10931E+02), and I had seen some of my PM10 or PM2.5 are negative indeed. My question now is related to : "If there are negative emissions, just set them to zero otherwise they lead to negative concentrations in the model and consequently negative refractive indices.", how do you mean? 

In the module 

chem/module_optical_averaging.F

i set the refr=0 if becomes negative? or what?

Thanks,

Leo

[Rajesh Kumar]

Hi Leo,

You should set the negative emissions to zero in wrfchemi* files and don't have to modify anything in the model code. 

Thanks,
Rajesh

Rajesh Kumar, Project Scientist III