[RegCNET] CBMZ chemistry crash

[Shalaby Kamal Ahmed]

Hi everyone
My thought is that, most long lived species (the transported species) in
the scheme is saved, what is left is extremely short lived species like
OH, HO2, O1D and so one, those species generally are created and destroyed
immediately and no real need to save them, this problem happen I think
when there is a small numerical error that could lead to the instability
of the solver. So by suppressing chemall we may not lose the efficiency of
the solver. But I need to check this throughly.
best regards
> Hi everyone,
>
> I have recently tried commenting the chemall sections in the
> mod_che_chemistry.F90, as suggested by Ahmed.
>
>
>           ! 1 : initialise xrin with the concentrations from
>           !     previous chemsolv step
> !          do ic = 1 , totsp
> !            xrin(ic) = chemall(j,i,k,ic)
> !          end do
>
> and
>
>           ! save the concentrations of all species for next chemistry step
> !          do ic = 1 , totsp
> !            chemall(j,i,k,ic) = xrout(ic)
> !          end do
>           ! Store photolysis rates for diagnostic
>
> Dear Ahmed,
>
> The simulation has just run to completion this way. however, I am
> concerned
> that this way, there is no chemistry continuity between time steps. It
> this
> the case? if so, how can we resolve this?
>
> Best regards,
> James
>
>
>
> On 1 August 2015 at 12:49, <zeinab.salah at yahoo.com> wrote:
>
>> Hi James,
>>
>> We try to run the recent version also with CBMZ or DCCB and we get this
>> error also.
>> The file of regcm.in and the error message are in the attachments.
>>
>> Best regards,
>> Zeinab
>>
>>
>>
>> ------------------------------
>> *Ù
Ù†:* James Ciarlo` <james.ciarlo at physics.org>
>> *إلى:* RegCNET <regcnet at lists.ictp.it>
>> *نسخة كربونية:* "ashalaby at ictp.it" <ashalaby at ictp.it>
>> *تاريخ الإرسال:* الثلاثاء 21 يوليو، 2015‏
>> 2:01 Ù

>> *الÙ
وضوع:* [RegCNET] CBMZ chemistry crash
>>
>> Dear RegCNETers,
>>
>> I have installed a recent version of the RegCM and tried to run a CBMZ
>> simulation for 1 year. However, the chemistry crashed and the simulation
>> kept on running producing a massive error log.
>>
>> After including a few lines in the mod_che_chemistry.F90 I got the
>> following information on when the error is happening:
>>
>> xrout, xrin, trac%indcbmz
>>   7.864835401599373E+153   38524086126.5782               10
>>
>> 10 is H2O2.
>>
>> I have attached the regcm.in file and the error obtained. I cannot
>> figure
>> out why this is happening. I was wondering if anyone has any idea on the
>> matter.
>>
>> Thanks and regards,
>> James
>>
>> _______________________________________________
>> RegCNET mailing list
>> RegCNET at lists.ictp.it
>> https://lists.ictp.it/mailman/listinfo/regcnet
>>
>>
>