Changeset 1495 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 10, 2015, 4:45:56 AM (9 years ago)

Author:

flefevre

Message:

Implementation du coeur chimique Adaptative Semi-Implicit Symmetric (ASIS)
Par defaut non actif pour l'instant (variable logique asis = .false.)

Location:

trunk/LMDZ.MARS/libf

Files:

• aeronomars/calchim_asis.F90 (added)
• aeronomars/photochemistry_asis.F90 (added)
• aeronomars/photolysis_asis.F90 (added)
• aeronomars/types_asis.F90 (added)
• phymars/physiq.F (modified) (3 diffs)

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

@@ -202 +202 @@
       real rsedcloud(ngrid,nlayer) ! Cloud sedimentation radius
       real rhocloud(ngrid,nlayer)  ! Cloud density (kg.m-3)
-      REAL surfdust(ngrid,nlayer) ! dust surface area (m2/m3, if photochemistry)
-      REAL surfice(ngrid,nlayer)  !  ice surface area (m2/m3, if photochemistry)
       REAL inertiesoil(ngrid,nsoilmx) ! Time varying subsurface
                                       ! thermal inertia (J.s-1/2.m-2.K-1)
                                       ! (used only when tifeedback=.true.)
-
+c     Variables used by the photochemistry
+      logical :: asis             ! true  : adaptative semi-implicit symmetric (asis) chemical solver
+                                  ! false : euler backward chemical solver
+      REAL surfdust(ngrid,nlayer) ! dust surface area (m2/m3, if photochemistry)
+      REAL surfice(ngrid,nlayer)  !  ice surface area (m2/m3, if photochemistry)
 c     Variables used by the slope model
       REAL sl_ls, sl_lct, sl_lat
@@ -1220 +1222 @@
      $                       surfdust, surfice)
 !           call photochemistry
+
+            asis = .false.
+
+            if (asis) then
+            call calchim_asis(ngrid,nlayer,nq,
+     &                   ptimestep,zplay,zplev,pt,pdt,dist_sol,mu0,
+     $                   zzlev,zzlay,zday,pq,pdq,zdqchim,zdqschim,
+     $                   zdqcloud,zdqscloud,tauref,co2ice,
+     $                   pu,pdu,pv,pdv,surfdust,surfice)
+            else
             call calchim(ngrid,nlayer,nq,
      &                   ptimestep,zplay,zplev,pt,pdt,dist_sol,mu0,
@@ -1225 +1237 @@
      $                   zdqcloud,zdqscloud,tauref,co2ice,
      $                   pu,pdu,pv,pdv,surfdust,surfice)
+            end if
 
            ! increment values of tracers: