Changeset 2213 for trunk/LMDZ.MARS/libf/aeronomars/calchim_mod.F90

Timestamp:

Jan 8, 2020, 12:01:49 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in calchim (only had an impact when computations include the
ionosphere), parameters should not be reset at every call.
EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/calchim_mod.F90 (modified) (2 diffs)

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

@@ -191 +191 @@
 !     for output:
 
-      logical :: output             ! to issue calls to writediagfi and stats
+      logical,save :: output        ! to issue calls to writediagfi and stats
       real :: jo3_3d(ngrid,nlayer)  ! Photodissociation rate O3->O1D (s-1)
       real :: jh2o_3d(ngrid,nlayer)  ! Photodissociation rate H2O->H+OH (s-1)
@@ -200 +200 @@
 
 !=======================================================================
+!     initialization of the chemistry (first call only)
+!=======================================================================
+
+      if (firstcall) then
+
+!=======================================================================
 !     main dashboard for the chemistry
 !=======================================================================
 
-      unichim = .true.     ! true : unified chemistry ! false : separate models in lower and upper atmosphere
-      jonline = .true.     ! true : on-line calculation of photodissociation rates ! false : lookup table
-      ionchem = .false.    ! switch for ionospheric chemistry
-      depos   = .false.    ! switch for dry deposition
-      output  = .false.    ! true : triggers writing of specific outputs (photolysis rates, emission rates, etc)
-
-!=======================================================================
-!     initialization of the chemistry (first call only)
-!=======================================================================
-
-      if (firstcall) then
+        unichim = .true.     ! true : unified chemistry ! false : separate models in lower and upper atmosphere
+        jonline = .true.     ! true : on-line calculation of photodissociation rates ! false : lookup table
+        ionchem = .false.    ! switch for ionospheric chemistry
+        depos   = .false.    ! switch for dry deposition
+        output  = .false.    ! true : triggers writing of specific outputs (photolysis rates, emission rates, etc)
 
          if (photochem) then