[RegCNET] بخصوص: CBMZ chemistry crash

[James Ciarlo`]

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150801/7bb8f16e/attachment.html>