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

Timestamp:

Feb 8, 2017, 4:36:25 PM (8 years ago)

Author:

emillour

Message:

Mars GCM:

• Added possibility to run with an Helium "he" tracer (to be initialized at constant value of 3.6e-7 kg/kg_air, i.e. the 4ppm of Krasnopolsky 1996 EUVE satellite, using newstart).
• corrected case for CH4 in aeronomars/photochemistry.F (was using undefined indexes when there is no CH4).
• updated aki/cpi coefficients for Argon used to compute mean atmospheric Cp and thermal conductivity in aeronomars/concentrations.F

JYC+EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• dyn1d/meteoriticflux_mod.F90 (deleted)
• initracer.F (modified) (2 diffs)
• physiq_mod.F (modified) (2 diffs)
• tracer_mod.F90 (modified) (1 diff)

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

@@ -88 +88 @@
       igcm_no2=0
       igcm_n2d=0
+      igcm_he=0
       igcm_co2plus=0
       igcm_oplus=0
@@ -238 +239 @@
           count=count+1
         endif
+        if (noms(iq).eq."he") then
+          igcm_he=iq
+          mmol(igcm_he)=4.
+          count=count+1
+        endif
         if (noms(iq).eq."co2plus") then
           igcm_co2plus=iq

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

@@ -20 +20 @@
      &                      rho_ice_co2,nuiceco2_sed,nuiceco2_ref,
      &                      igcm_dust_mass, igcm_dust_number, igcm_h2o2,
-     &                      nuice_ref, rho_ice, rho_dust, ref_r0
+     &                      nuice_ref, rho_ice, rho_dust, ref_r0,
+     &                      igcm_he
       use comsoil_h, only: inertiedat, ! soil thermal inertia
      &                     tsoil, nsoilmx ! number of subsurface layers
@@ -2336 +2337 @@
          endif ! of if (tracer.and.(igcm_co2.ne.0))
 
+
+        ! Output He tracer, if there is one
+        if (tracer.and.(igcm_he.ne.0)) then
+          call WRITEDIAGFI(ngrid,"he","helium mass mixing ratio",
+     &                     "kg/kg",3,zq(1,1,igcm_he))
+          vmr=zq(1:ngrid,1:nlayer,igcm_he)
+     &       *mmean(1:ngrid,1:nlayer)/mmol(igcm_he)
+          call WRITEDIAGFI(ngrid,'vmr_he','helium vol. mixing ratio',
+     &                       'mol/mol',3,vmr)
+        endif
+
 c        ----------------------------------------------------------
 c        Outputs of the water cycle

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

@@ -73 +73 @@
       integer,save :: igcm_no2
       integer,save :: igcm_n2d
+      integer,save :: igcm_he
       integer,save :: igcm_ch4
       ! Ions