Changeset 2332 for trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F

Timestamp:

May 26, 2020, 12:17:13 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Follow-up of the last commit for the transport of the isotopic ratio: simplification of the transmission of variables from the dynamics to the
physics:

• libf/dynphy_lonlat/phymars/: iniphysiq_mod.F90: transmission of the content of 2 variables describing the isotopes instead of 4 (nqperes: number of tracers "peres", nqfils:

number of tracers "fils")

• libf/phymars/: phys_state_var_init_mod.F90, tracer_mod.F: idem callsedim_mod.F: idem co2condens_mod.F: idem
• libf/phymars/dyn1d: testphys1d.F: idem (the reading interface for traceur.def has been completed to fill the variables nqperes and nqfils).

MV

File:

• trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F (modified) (3 diffs)

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

@@ -23 +23 @@
      &                      igcm_stormdust_number,igcm_topdust_mass, 
      &                      igcm_topdust_number,
-     &                      iqfils,nqfils,qperemin,masseqmin ! MVals: variables isotopes
+     &                      nqfils,qperemin,masseqmin ! MVals: variables isotopes
       USE newsedim_mod, ONLY: newsedim
       USE comcstfi_h, ONLY: g
@@ -536 +536 @@
               !MVals: Loop over the sons ("fils")
               if (nqfils(iq).gt.0) then
-                iq2=iqfils(nqfils(iq),iq) ! for now always nqfils(iq)=1 (special case of HDO only son of H2O)
+                if (iq.eq.igcm_h2o_ice) then
+                 iq2=igcm_hdo_ice
+                else
+                 call abort_physic("callsedim_mod","invalid isotope",1)
+                endif 
                 !MVals: input paramters in vlz_fi for hdo
                 do l=1,nlay
@@ -595 +599 @@
               !MVals: Special case of isotopes: for now only HDO
               if (nqfils(iq).gt.0) then
-               iq2=iqfils(nqfils(iq),iq)
+                if (iq.eq.igcm_h2o_ice) then
+                 iq2=igcm_hdo_ice
+                else
+                 call abort_physic("callsedim_mod","invalid isotope",1)
+                endif
                pdqsed(ig,l,iq2)=(zqi(ig,l,iq2)-
      $            (pq(ig,l,iq2) + pdqfi(ig,l,iq2)*ptimestep))/ptimestep