Changeset 2467 for trunk/LMDZ.MARS/libf/aeronomars/thermosphere.F

Timestamp:

Mar 3, 2021, 1:48:21 PM (4 years ago)

Author:

emillour

Message:

Mars GCM:

• bug fix in the diagnostic of H, H2, and D total escape which was wrong when running in parallel. The escape rates can now also be included as outputs in diagfi.

FGG

File:

• trunk/LMDZ.MARS/libf/aeronomars/thermosphere.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/aeronomars/thermosphere.F

@@ -3 +3 @@
      $     mu0,ptimestep,ptime,zday,tsurf,zzlev,zzlay,
      &     pt,pq,pu,pv,pdt,pdq,
-     $     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff)
+     $     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff,
+     $     PhiEscH,PhiEscH2,PhiEscD)
 
       use conc_mod, only: rnew, cpnew
@@ -34 +35 @@
       REAL,INTENT(out) :: zdvmolvis(ngrid,nlayer)
       REAL,INTENT(out) :: zdqmoldiff(ngrid,nlayer,nq)
+      REAL*8,INTENT(out) :: PhiEscH,PhiEscH2,PhiEscD
 
       INTEGER :: l,ig
@@ -80 +82 @@
         call moldiff_red(ngrid,nlayer,nq,
      &                   pplay,pplev,pt,pdt,pq,pdq,ptimestep,
-     &                   zzlay,zdteuv,zdtconduc,zdqmoldiff)
+     &                   zzlay,zdteuv,zdtconduc,zdqmoldiff,
+     &                   PhiEscH,PhiEscH2,PhiEscD)
       endif