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

Timestamp:

May 14, 2020, 4:11:59 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Follow-up of the last commit for the transport of the isotopic ratio: forgot the initialisation of zq0() in callsedim_mod.F !
MV

File:

• trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F (modified) (1 diff)

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

@@ -338 +338 @@
 c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       zqi(1:ngrid,1:nlay,1:nq)=pq(1:ngrid,1:nlay,1:nq)
+     &                         +pdqfi(1:ngrid,1:nlay,1:nq)*ptimestep
+      zq0(1:ngrid,1:nlay,1:nq)=pq(1:ngrid,1:nlay,1:nq) !MVals: keep the input value
      &                         +pdqfi(1:ngrid,1:nlay,1:nq)*ptimestep
       zt(1:ngrid,1:nlay)=pt(1:ngrid,1:nlay)