Changeset 2164

Timestamp:

Sep 25, 2019, 10:05:39 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Add the possibility to super-cycle chemistry computations. Chemistry will
be computed every "ichemistry" physics step (but chemistry tendencies are
added every physics step). Default is ichemistry=1.
AB+EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/calchim_mod.F90 (modified) (1 diff)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -2752 +2752 @@
 - Externalize in "physiq" the computation of solar zenithal angle
   (it should be computed at every physics timestep, regardless of iradia)
+
+== 25/09/2019 == AB+EM
+Add the possibility to super-cycle chemistry computations. Chemistry will
+be computed every "ichemistry" physics step (but chemistry tendencies are
+added every physics step). Default is ichemistry=1.

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

@@ -3 +3 @@
   IMPLICIT NONE
 
+  INTEGER, SAVE :: ichemistry ! compute chemistry every ichemistry physics step
   REAL,SAVE,ALLOCATABLE :: zdqchim(:,:,:) ! Tendancy on pq due to photochemistry
   REAL,SAVE,ALLOCATABLE :: zdqschim(:,:) ! Tendancy on qsurf due to photochemistry

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

@@ -44 +44 @@
      &                      ini_scatterers,tauvis
       use datafile_mod, only: datadir
+      use calchim_mod, only: ichemistry
 
       IMPLICIT NONE
@@ -688 +689 @@
          write(*,*) "tituscap",tituscap
                      
-
+! Chemistry:
          write(*,*) "photochemistry: include chemical species"
          photochem=.false. ! default value
          call getin("photochem",photochem)
          write(*,*) " photochem = ",photochem
+         
+         write(*,*) "Compute chemistry (if photochem is .true.)",
+     &   "every ichemistry physics step (default: ichemistry=1)"
+         ichemistry=1
+         call getin("ichemistry",ichemistry)
+         write(*,*) " ichemistry = ",ichemistry
 
 

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

@@ -15 +15 @@
 
       use watercloud_mod, only: watercloud, zdqcloud, zdqscloud
-      use calchim_mod, only: calchim, zdqchim, zdqschim
+      use calchim_mod, only: calchim, ichemistry, zdqchim, zdqschim
       use watersat_mod, only: watersat
       use co2condens_mod, only: co2condens
@@ -1615 +1615 @@
          IF (photochem) then
 
+           if (modulo(icount-1,ichemistry).eq.0) then
+           ! compute chemistry every ichemistry physics step
+
 !           dust and ice surface area
             call surfacearea(ngrid, nlayer, naerkind,
@@ -1627 +1630 @@
      $                   zdqcloud,zdqscloud,tau(:,1),co2ice,
      $                   pu,pdu,pv,pdv,surfdust,surfice)
+
+            endif ! of if (modulo(icount-1,ichemistry).eq.0)
 
            ! increment values of tracers: