[RegCNET] بØ(R)صوص: CBMZ chemistry crash

[James Ciarlo`]

Yes, it's just what I feared. I'll try doing just that.

On 1 August 2015 at 15:35, Solmon Fabien <fsolmon at ictp.it> wrote:

> Hello James
>
> if you do not save chemall all the fast (not transported) species are
> reset to zero , the solver will recalculate a new equilibrium at every
> time step .. in fact with this fix you might still generate unrealistic
> concentrations for some untransported species but since you reset to zero
> there are  no further consequence and the code is stable , however it does
> not mean that it performs a good chemistry ... I would try to screen the
> to the solver inputs just before it starts exploding, and try to play with
> chemistry step dtchem.. try to reduce it from 900 s to dt maybe ..
>
>  Fabien
> > 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
> >>
> >>
> > _______________________________________________
> > RegCNET mailing list
> > RegCNET at lists.ictp.it
> > https://lists.ictp.it/mailman/listinfo/regcnet
>
>
> --
> Fabien Solmon
> Research Scientist
> +39 040 2240206 (office) | Earth System Physics
> +39 040 2240449 (fax)    | ICTP (UNESCO-IAEA)
> esp.ictp.it              | Strada Costiera 11
>                          | 34014 Trieste, Italy
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150802/58a06b68/attachment-0001.html>