Changeset 1708

Timestamp:

Apr 28, 2017, 9:15:01 AM (8 years ago)

Author:

flefevre

Message:

Photochimie : le solveur chimique ASIS est desormais solveur par defaut.
Pour plus de details sur la methode voir Cariolle et al, Geosci. Model Dev., 10, 1467-1485, 2017.
(avec des vrais bouts de Mars dedans!)

Location:

trunk/LMDZ.MARS/libf

Files:

• aeronomars/photochemistry_asis.F90 (modified) (2 diffs)
• phymars/physiq_mod.F (modified) (1 diff)

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

@@ -6 +6 @@
 !     ------
 !
-!     Version: 10/11/2015
+!     Version: 27/04/2017
+!
+!     ASIS scheme : for details on the method see
+!     Cariolle et al., Geosci. Model Dev., 10, 1467-1485, 2017.
 !
 !*****************************************************************
@@ -304 +307 @@
 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.1   ! 1/rtol recommended value : 0.1-0.02
+real (kind = 8), parameter :: rtol    = 1./0.05   ! 1/rtol recommended value : 0.1-0.02
 integer,         parameter :: niter   = 3        ! number of iterations
 real (kind = 8), parameter :: coefmax = 2.

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -1330 +1330 @@
      $                       surfdust, surfice)
 !           call photochemistry
-
-            asis = .false.
+!           asis = .false.  euler-backward scheme 
+!           asis = .true.   asis scheme 
+
+            asis = .true.
 
             if (asis) then