[RegCNET] Problem in RegCM3 aerosol simulation.
Prodromos Zanis
pzanis at geol.uoa.gr
Thu Apr 5 11:15:29 CEST 2007
Dear Dilip and Konare
I was following your exchange of emails and I could not find either the
directory RegCM/PreProc/DATA/AERGLOB/orig. I would be grateful if you
could also send me the updated Aerosol.data and ctl files.
With kind regards
Prodromos Zanis
> Hi Dilip,
>
> I'll send to you later today the AEROSOL.dat AEROSOL.ctl obtained after
> the changes and we will see.
>
> Best regards,
>
> Konare
>
> Dilip Ganguly <dganguly at prl.res.in> wrote: Hi,
> Thanks Konare for your reply. But, I do not have any directory like
> RegCM/PreProc/DATA/AERGLOB/orig
> and hence I am not able to find the routines conv_bin.f and conv_bin.x
>
> When I run the datalinker.x, RegCM/PreProc/DATA/AERGLOB/ directory is
> created which contains link to the files AEROSOL.dat and AEROSOL.ctl.
>
> So, I do not know how to incorporate the changes you are suggesting to
> correct the conversion factor problem in aerosol data.
>
> If anyone has the corrected aerosol data then I request him/her to put it
> on the website so that I may download it. Best thing would be if Bi does
> it on the ICTP RegCM webpage.
>
> Konare, can you send me the corrected aerosol data.
>
> Thanks.
>
> Dilip.
>
>
> On Wed, 4 Apr 2007, Abdourahamane Konare wrote:
>
>>
>> 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 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
> 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
> ---------------------------------------------------------------------
>
>
>
>
>
> ---------------------------------
> The fish are biting.
> Get more visitors on your site using Yahoo! Search
> Marketing._______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>
--
****************************************************
Dr. Prodromos Zanis
Centre for Atmospheric Physics and Climatology
Academy of Athens
3rd September 131, Athens 11251, Greece
Tel. +30 210 8832048
Fax: +30 210 8832048
e-mail: pzanis at geol.uoa.gr
Web address: http://users.auth.gr/~zanis/
*****************************************************
More information about the RegCNET
mailing list