Changeset 2322 for trunk/LMDZ.MARS/libf/phymars/tracer_mod.F90

Timestamp:

May 13, 2020, 5:14:02 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Extent of the transport of the isotopic ratio implemented in the dynamics to all the Van Leer transport schemes used in the physics (for now it only
concerns the tracer HDO).

• libf/dynphy_lonlat/phymars/: iniphysiq_mod.F90: transmission of the content of variables describing the isotopes defined in the dynamics (precisely by dyn3d_common/infotrac.F90,

which reads traceur.def) to the physics

• libf/phymars/: phys_state_var_init_mod.F90, tracer_mod.F : initialisation of the variables describing the isotopes in the physics callsedim_mod.F: implementation of the transport of the isotopic ratio in the Van Leer scheme used for sedimentation (applies to hdo ice) co2condens_mod.F: implementation of the transport of the isotopic ratio in the Van Leer scheme used for condensation of CO2 (applies to hdo ice and

vapour)
MV

File:

• trunk/LMDZ.MARS/libf/phymars/tracer_mod.F90 (modified) (3 diffs)

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

@@ -95 +95 @@
       ! other tracers
       integer,save :: igcm_ar_n2 ! for simulations using co2 +neutral gas
-
+      ! MVals: isotopes
+      integer, save                 :: nqperes ! numbers of tracers defined as "peres"
+      integer, allocatable, save    :: nqfils(:) ! numbers of sons ("fils") of the considered tracer
+      integer, allocatable, save    :: iqfils(:,:) ! indice of a son, ex: iqfils(nqfils(ipere),ipere)
+      real, parameter               :: qperemin=1.e-16 ! threschold for the "pere" mixing ratio qpere to calculate Ratio=qfils/qpere
+      real, parameter               :: masseqmin=1.e-16 ! threschold for the "pere" transporting masse
+      !integer, allocatable, save    :: nqdesc(:) ! number of sons + all gran-sons over all generations: not useful for now in the martian case as there are no gran-sons
 
 !-----------------------------------------------------------------------
@@ -101 +107 @@
   contains
   
-    subroutine ini_tracer_mod(nq,tname)
+    subroutine ini_tracer_mod(nq,tname,dyn_nqdesc,dyn_iqfils,dyn_nqperes,dyn_nqfils)! MVals: variables isotopes
       implicit none
       
       integer,intent(in) :: nq ! number of tracers
       character(len=*),intent(in) :: tname(nq) ! tracer names
+      !MVals: variables isotopes
+      integer,intent(in) :: dyn_nqperes
+      integer,intent(in) :: dyn_nqfils(nq)
+      integer,intent(in) :: dyn_nqdesc(nq)
+      integer,intent(in) :: dyn_iqfils(nq,nq)
       
       integer :: iq, count
@@ -117 +128 @@
         write(*,*) "tracer_mod names : ", trim(noms(iq))
       enddo
+
+      !MVals: isotopes variables initialisation
+      do iq=1,nq
+        if (dyn_nqfils(iq).ne.dyn_nqdesc(iq)) then
+          write(*,*) ' for now all descendants must be sons: check the', &
+                     '  relations between tracers in traceur.def !'
+          call abort_physic("ini_tracer_mod","relatives pattern between tracers not accepted",1)
+        endif
+      enddo
+      allocate(nqfils(nq))!,nqdesc(nq))    
+      allocate(iqfils(nq,nq))
+      nqperes=dyn_nqperes   
+      nqfils(:)=dyn_nqfils(:)
+      iqfils(:,:)=dyn_iqfils(:,:)
      
 #ifndef MESOSCALE