Changeset 2467 for trunk/LMDZ.MARS/libf/phymars

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/phymars/physiq_mod.F (modified) (3 diffs)

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

@@ -342 +342 @@
       REAL zdvmolvis(ngrid,nlayer)
       real zdqmoldiff(ngrid,nlayer,nq)
+      real*8 PhiEscH,PhiEscH2,PhiEscD
 
       REAL dwatercap(ngrid), dwatercap_dif(ngrid)     ! (kg/m-2)
@@ -1887 +1888 @@
      $     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)
 
         DO l=1,nlayer
@@ -3304 +3306 @@
      $           3,zdtnirco2)
 
+            !H, H2 and D escape fluxes
+
+            call WRITEDIAGFI(ngrid,"PhiH","H escape flux","s-1",
+     $           0,PhiEscH)
+            call WRITEDIAGFI(ngrid,"PhiH2","H2 escape flux","s-1",
+     $           0,PhiEscH2)
+            call WRITEDIAGFI(ngrid,"PhiD","D escape flux","s-1",
+     $           0,PhiEscD)
+
+!            call wstats(ngrid,"PhiH","H escape flux","s-1",
+!     $           0,PhiEscH)
+!            call wstats(ngrid,"PhiH2","H2 escape flux","s-1",
+!     $           0,PhiEscH2)
+!            call wstats(ngrid,"PhiD","D escape flux","s-1",
+!     $           0,PhiEscD)
+            
 !            call wstats(ngrid,"q15um","15 um cooling","K/s",
 !     $           3,zdtnlte)