[RegCNET] Ø¨Ø(R)ØµÙˆØµ: CBMZ chemistry crash

Sat Aug 15 10:50:51 CEST 2015

Hi again everyone,

I have tried reducing the dt stepwise (since I could not find dtchem) and
have so far brought it down to 30. I am afraid that the model is still
crashing in the same manner.

I am currently scouring the code to try to understand what may be causing
this error, as it has now also popped up in my CB6 simulation (not just
CBMZ). I am attaching the error again for reference.

Thanks and regards,
James

On 2 August 2015 at 14:18, James Ciarlo` <james.ciarlo at physics.org> wrote:

> Hi again,
>
> I have tried to use the dtchem but I could not find it anywhere in the
> User Guide or the model code, even after doing an svn update.
>
> On 2 August 2015 at 12:11, James Ciarlo` <james.ciarlo at physics.org> wrote:
>
>> 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/20150815/0d083e56/attachment-0001.html>
-------------- next part --------------
 $$$ Jvalue min/max NO2surf :  0.10299E-03 0.10321E-03
 $$$ 1999-02-07 03:00:00 UTC, ktau   =           106920
 $$$ 1st, 2nd time deriv of ps   =  0.10338E-04 0.21996E-06
 $$$  no. of points w/convection =    3201
 SRF variables written at 1999-02-07 03:00:00 UTC
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.1957239150519D-13
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.9786195752595D-14
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.1957239150519D-13
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.1957239150519D-13
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.2357280358500D-13
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.5893200896249D-14
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.1118394380616D-13
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.2795985951539D-14
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.5591971903078D-14
 DLSODE-  Warning..internal T (=R1) and H (=R2) are
       such that in the machine, T + H = T on the next step
       (H = step size). Solver will continue anyway
      In above,  R1 =  0.2941501835651D+03   R2 =  0.1397992975770D-14
 DLSODE-  Above warning has been issued I1 times.
       It will not be issued again for this problem
      In above message,  I1 =        10
 DLSODE-  At T (=R1) and step size H (=R2), the
       corrector convergence failed repeatedly
       or with ABS(H) = HMIN
      In above,  R1 =  0.2941501835651D+03   R2 =  0.2490238543513-155
 Abort called by computing node          153 at 2015-08-14 15:38:12.137 +0200