Changeset 2150

Timestamp:

Sep 5, 2019, 2:50:33 PM (5 years ago)

Author:

flefevre

Message:

Correction de bug mineur sur calcul du pas de temps photochimique.

File:

• trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90

@@ -352 +352 @@
 real (kind = 8), parameter :: dtmin   = 10.      ! minimum time step (s)
 real (kind = 8), parameter :: vmrtol  = 1.e-11   ! absolute tolerance on vmr
-real (kind = 8), parameter :: rtol    = 1./0.05   ! 1/rtol recommended value : 0.1-0.02
+real (kind = 8), parameter :: rtol    = 0.05     ! rtol recommended value : 0.1-0.02
 integer,         parameter :: niter   = 3        ! number of iterations
 real (kind = 8), parameter :: coefmax = 2.
@@ -406 +406 @@
 do iesp = 1,nesp
    es = 2.*abs((ratio*cnew(iesp) - (1. + ratio)*ccur(iesp) + cold(iesp))   &
-         /(1. + ratio)/max(ccur(iesp),atol))
+         /(1. + ratio)/max(ccur(iesp)*rtol,atol))
 
    if (es > e) then
@@ -412 +412 @@
    end if
 end do
-e = rtol*e
 
 ! timestep correction
@@ -700 +699 @@
 !---  c008: ho2 + ho2 -> h2o2 + o2
 
-!     jpl 2006
-
-!     c008(:) = 3.5e-13*exp(430./t(:))
+!     jpl 2015
+
+!     c008(:) = 3.0e-13*exp(460./t(:))
 
 !     christensen et al., grl, 13, 2002