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

Timestamp:

Feb 16, 2021, 12:31:33 PM (4 years ago)

Author:

emillour

Message:

Mars GCM:

• Adding the deuterium chemistry now that the HDO cycle is included.
• Chemistry still works as before if deuterium tracers are not present.
• Added handling of hdo in molecular diffusion (moldiff_red).

FGG+JYC+EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• initracer.F (modified) (2 diffs)
• tracer_mod.F90 (modified) (1 diff)

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

@@ -84 +84 @@
       igcm_o3=0
       igcm_h=0
+      igcm_d=0
+      igcm_hd=0
       igcm_h2=0
+      igcm_od=0
+      igcm_do2=0
+      igcm_hdo2=0
       igcm_oh=0
       igcm_ho2=0
@@ -365 +370 @@
           count=count+1
         endif
+        if (noms(iq).eq."od") then
+          igcm_od=iq
+          mmol(igcm_od)=18.
+          count=count+1
+        endif
+        if (noms(iq).eq."d") then
+           igcm_d=iq
+           mmol(igcm_d)=2.
+           count=count+1
+        endif
+        if (noms(iq).eq."hd") then
+           igcm_hd=iq
+           mmol(igcm_hd)=3.
+           count=count+1
+        endif
+        if (noms(iq).eq."do2") then
+           igcm_do2=iq
+           mmol(igcm_do2)=34.
+           count=count+1
+        endif
+        if (noms(iq).eq."hdo2") then
+           igcm_hdo2=iq
+           mmol(igcm_hdo2)=35.
+           count=count+1
+        endif
         if (noms(iq).eq."co2_ice") then
           igcm_co2_ice=iq

trunk/LMDZ.MARS/libf/phymars/tracer_mod.F90

@@ -77 +77 @@
       integer,save :: igcm_he
       integer,save :: igcm_ch4
+      !Deuterated species derived from HDO
+      integer,save :: igcm_od
+      integer,save :: igcm_d
+      integer,save :: igcm_hd
+      integer,save :: igcm_do2
+      integer,save :: igcm_hdo2
       ! Ions
       integer,save :: igcm_co2plus