[RegCNET] Problem in RegCM3 aerosol simulation.

Abdourahamane Konare konarea at yahoo.com
Wed Apr 4 13:46:47 CEST 2007


Hi Dilip,
cd PreProc/DATA/AERGLOB/orig

check the factors of conversion in conv_bin.f as I suggested to you
change them based on the attached document of the previous mail
run conv_bin.x, this will create the new AEROSOL.dat file which will be used in ICBC
copy in AERGLOB the AEROSOL.dat or ( check the link to it from the AEROSOL.f )
so when running icbc.x, you will have the new aerosol dataset that will be used for your simulations.

Konare

PS if you don't change any thing it will used the AEROSOL.dat that already exist when downloading the model and data.







Dilip Ganguly <dganguly at prl.res.in> wrote: Hi Konare,
Many thanks for your reply. However, I am not able to understand where 
should I find the conv_bin.f routine. I can understand that there is a 
mistake of conversion factor but how do I correct my aerosol data. Even 
Solmon was saying the same thing to me. I am not generating the aerosol 
data from any inventory database but I am using whatever is provided in 
the ICTP website.

Long back Bi had sent an email mentioning some improvements in 
AEROSOL.dat. I have downloaded the aerosol data (AEROSOL.dat and 
AEROSOL.ctl) from the ICTP website after this announcement and I expected 
the aerosol data to be corrected for the same bug. However, it may so 
happen that the aerosol data available in the ICTP website was not updated 
for the bug related to conversion factor.

I would like to know if the aerosol data currently available in the ICTP 
website is corrected for conversion factor problem. If not then from where 
I can get the corrected aerosol data.

Thanks in advance.

Dilip.


On Wed, 4 Apr 2007, Abdourahamane Konare wrote:

> Hi Dilip,
> The problem you mentioned at your second point is coming from the conv_bin.f  routine in AEROSOL/orig. There is a mistake in the  conversion factors : ! kg/yr/1x1 cell ==> kg/m2/s, have a look in the file attached ( from the previous version) and replace accordingly the coefficients  a(i,j)=a(i,j)/(0.5*(ax1+ax2)*by)/86400/365   in conv_bin.f  by the corresponding one in the file attached, this will deal with the second point.  The problem you mentioned in point one which has been dealt before however if it remains let me know then i will send to you my Main which works.
> 
> Best regards
> 
> 
> Konare
> Dilip Ganguly  wrote: Hi everyone,
> 
> As pointed out by Solmon, I have made the correction in my
> RegCM/Main/aermix.f code. However, even after making this change, I am
> getting exactly the same result as my previous model output obtained
> without incorporating this correction.
> 
> I am surprised to find that the values of aerosol extinction (aext8) and
> all other tracer mixing ratios are exactly same at all atmospheric levels.
> There is absolutely no difference in their magnitudes at different heights
> (sigma pressure levels). This is not expected and it can not be true.
> 
> Another surprising thing which I have noticed is that aerosol extinction
> values are controlled by the dust loading alone. None of the other tracers
> such as BC or OC seem to have any affect on the aerosol extinction.
> Aerosol extinction values (aext8) obtained from RegCM simulation are
> unbelievably low over the regions where there is high column burden of BC
> and OC aerosols. This is a contradiction. I am not able to understand if I
> am making any mistakes while setting various parameters in my simulations
> or BC/OC aerosols simulated by the model does not radiatively interact in
> the RegCM3 model and hence it is not possible to study the effect of these
> aerosols on climate using the present version of the model.
> 
> I have attached my domain.param and regcm.in files along with this mail so
> that it may be checked if I have set all the parameters correctly for
> including dust/aerosols in my model simulation.
> 
> I made following changes in the RegCM/Main/parame file to include all the
> ten tracers
> 
> integer ntr                      ! tracer number
>       parameter( ntr = 10)
>       integer nbin                     ! bins number ( dust)
>       parameter(nbin = 4)
> --------------------------------------------------------------------------------------------
> 
> Can anybody tell, what is going wrong in my RegCM3 aerosol simulation.
> 
> Thanks in advance.
> 
> Dilip.
> 
> 
> -- 
> Dr. Dilip Ganguly
> Postdoctoral Fellow
> Space and Atmospheric Sciences Division
> Physical Research Laboratory
> Navrangpura
> Ahmedabad-380009
> Gujarat, India
> Phone:+91 79 26314558, Fax: +91 79 26314659
> Alternate Email: ganguly.dilip at gmail.com
> visit me at: http://www.prl.res.in/~dganguly_______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
> 
> 
>  
> ---------------------------------
> Be a PS3 game guru.
> Get your game face on with the latest PS3 news and previews at Yahoo! Games.

Dilip. 
  
     ---------------------------------------------------------------------
        Dr. DILIP GANGULY                               
        Space And Atmospheric Sciences Division      
        Physical Research Laboratory                
        Ahmedabad 380009
 Gujarat, India
        Phone: +91 79 26314558    
 Fax:   +91 79 26314659
 Alternate Email: ganguly.dilip at gmail.com          
        URL: http://www.prl.res.in/~dganguly
     ---------------------------------------------------------------------




 
---------------------------------
Food fight? Enjoy some healthy debate
in the Yahoo! Answers Food & Drink Q&A.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20070404/413b569e/attachment-0002.html>


More information about the RegCNET mailing list