Changeset 1571 for trunk/LMDZ.MARS/libf/aeronomars/photochemistry_asis.F90

Timestamp:

Jul 7, 2016, 9:12:29 AM (8 years ago)

Author:

flefevre

Message:

Optimisation des reglages du solveur chimique ASIS.

File:

• trunk/LMDZ.MARS/libf/aeronomars/photochemistry_asis.F90 (modified) (6 diffs)

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

@@ -201 +201 @@
 
    call define_dt(nesp, dt_corrected, dt_guess, ctimestep, cold(:), c(ilev,:),  &
-                  mat1, prod, loss)
+                  mat1, prod, loss, dens(ilev))
 
    if (time + dt_corrected > ptimestep) then
@@ -263 +263 @@
 !================================================================
 
- subroutine define_dt(nesp, dtnew, dtold, ctimestep, cold, ccur, mat1, prod, loss)
+ subroutine define_dt(nesp, dtnew, dtold, ctimestep, cold, ccur, mat1, &
+                      prod, loss, dens)
 
 !================================================================
@@ -282 +283 @@
 real (kind = 8), dimension(nesp,nesp) :: mat1
 real (kind = 8), dimension(nesp)      :: prod, loss
+real (kind = 8)                       :: dens
 
 ! output
@@ -291 +293 @@
 real (kind = 8), dimension(nesp)      :: cnew
 real (kind = 8), dimension(nesp,nesp) :: mat
-real (kind = 8) :: ratio, e, es, coef
+real (kind = 8) :: atol, ratio, e, es, coef
 
 integer                  :: code, iesp, iter
@@ -301 +303 @@
 
 real (kind = 8), parameter :: dtmin   = 10.      ! minimum time step (s)
-real (kind = 8), parameter :: atol    = 1.e4     ! recommended value : 1.e4
+real (kind = 8), parameter :: vmrtol  = 1.e-11   ! absolute tolerance on vmr
 real (kind = 8), parameter :: rtol    = 1./0.1   ! 1/rtol recommended value : 0.1-0.02
 integer,         parameter :: niter   = 3        ! number of iterations
@@ -310 +312 @@
 dttest = dtold   ! dttest = dtold = dt_guess
 
+atol = vmrtol*dens ! absolute tolerance in molecule.cm-3
+
 do iter = 1,niter
 
 if (fast_guess) then
 
-! first guess : fast implicit method
+! first guess : fast semi-implicit method
 
    do iesp = 1, nesp