Changeset 2322 for trunk/LMDZ.MARS/libf/phymars/dyn1d

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/dyn1d/testphys1d.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -5 +5 @@
       use dimphy, only : init_dimphy
       use mod_grid_phy_lmdz, only : regular_lonlat
-      use infotrac, only: nqtot, tname
+      use infotrac, only: nqtot, tname, nqperes,nqdesc,iqfils,nqfils,
+     &                    iqpere, nqdesc_tot
       use comsoil_h, only: volcapa, layer, mlayer, inertiedat, nsoilmx
       use comgeomfi_h, only: sinlat, ini_fillgeom
@@ -127 +128 @@
       REAL halfaxe, excentric, Lsperi
       Logical paleomars
+
+c   MVals: isotopes as in the dynamics (CRisi)
+      INTEGER :: ifils,ipere,generation
+      CHARACTER(len=30), ALLOCATABLE, DIMENSION(:) :: tnom_transp ! transporting fluid short name
+      CHARACTER(len=80) :: line ! to store a line of text     
+      INTEGER ierr0
+      LOGICAL :: continu
 
 c=======================================================================
@@ -240 +248 @@
         allocate(dqdyn(nlayer,nq))
         allocate(mqtot(nq))
+        allocate(tnom_transp(nq))
         
         ! read tracer names from file traceur.def
         do iq=1,nq
-          read(90,*,iostat=ierr) tname(iq)
+          read(90,'(80a)',iostat=ierr) line ! store the line from traceur.def
           if (ierr.ne.0) then
             write(*,*) 'testphys1d: error reading tracer names...'
             stop
           endif
+          ! if format is tnom_0, tnom_transp (isotopes)
+          read(line,*,iostat=ierr0) tname(iq),tnom_transp(iq)
+          if (ierr0.ne.0) then
+            read(line,*) tname(iq)
+            tnom_transp(iq)='air'
+          endif
+
         enddo
         close(90)
+
+       ! Isotopes: as in the 3D case we have to determine father/son relations for isotopes and carrying fluid
+       ALLOCATE(nqfils(nqtot),nqdesc(nqtot))    
+       ALLOCATE(iqfils(nqtot,nqtot))   
+       ALLOCATE(iqpere(nqtot))
+       nqperes=0
+       nqfils(:)=0
+       nqdesc(:)=0
+       iqfils(:,:)=0
+       iqpere(:)=0
+       nqdesc_tot=0    
+       DO iq=1,nqtot
+       if (tnom_transp(iq) == 'air') then
+         ! ceci est un traceur père
+         WRITE(*,*) 'Le traceur',iq,', appele ',
+     &          trim(tname(iq)),', est un pere'
+         nqperes=nqperes+1
+         iqpere(iq)=0
+       else !if (tnom_transp(iq) == 'air') then
+         ! ceci est un fils. Qui est son père?
+         WRITE(*,*) 'Le traceur',iq,', appele ',
+     &                trim(tname(iq)),', est un fils'
+         continu=.true.
+         ipere=1
+         do while (continu)           
+           if (tnom_transp(iq) .eq. tname(ipere)) then
+             ! Son père est ipere
+             WRITE(*,*) 'Le traceur',iq,'appele ',
+     &   trim(tname(iq)),' est le fils de ',
+     &   ipere,'appele ',trim(tname(ipere))
+             nqfils(ipere)=nqfils(ipere)+1  
+             iqfils(nqfils(ipere),ipere)=iq
+             iqpere(iq)=ipere         
+             continu=.false.
+           else !if (tnom_transp(iq) == tnom_0(ipere)) then
+             ipere=ipere+1 
+             if (ipere.gt.nqtot) then
+                 WRITE(*,*) 'Le traceur',iq,'appele ',
+     &           trim(tname(iq)),', est orpelin.'
+                 CALL abort_gcm('infotrac_init',
+     &                  'Un traceur est orphelin',1)
+             endif !if (ipere.gt.nqtot) then
+           endif !if (tnom_transp(iq) == tnom_0(ipere)) then
+         enddo !do while (continu)
+       endif !if (tnom_transp(iq) == 'air') then
+       enddo !DO iq=1,nqtot
+       WRITE(*,*) 'nqperes=',nqperes    
+       WRITE(*,*) 'nqfils=',nqfils
+       WRITE(*,*) 'iqpere=',iqpere
+       WRITE(*,*) 'iqfils=',iqfils
+       ! Calculer le nombre de descendants à partir de iqfils et de nbfils
+       DO iq=1,nqtot    
+       generation=0
+       continu=.true.
+       ifils=iq
+       do while (continu)
+          ipere=iqpere(ifils)
+         if (ipere.gt.0) then
+          nqdesc(ipere)=nqdesc(ipere)+1   
+          nqdesc_tot=nqdesc_tot+1      
+          iqfils(nqdesc(ipere),ipere)=iq
+          ifils=ipere
+          generation=generation+1
+         else !if (ipere.gt.0) then
+          continu=.false.
+         endif !if (ipere.gt.0) then
+       enddo !do while (continu)    
+       WRITE(*,*) 'Le traceur ',iq,', appele ',trim(tname(iq)),
+     &               ' est un traceur de generation: ',generation
+       ENDDO !DO iq=1,nqtot
+       WRITE(*,*) 'infotrac: nqdesc=',nqdesc
+       WRITE(*,*) 'iqfils=',iqfils
+       WRITE(*,*) 'nqdesc_tot=',nqdesc_tot
 
         ! initialize tracers here:
@@ -583 +672 @@
       call init_dimphy(1,nlayer) ! Initialize dimphy module
       call phys_state_var_init(1,llm,nq,tname,
-     .          day0,time,daysec,dtphys,rad,g,r,cpp)
+     .          day0,time,daysec,dtphys,rad,g,r,cpp,
+     .          nqdesc,iqfils,nqperes,nqfils)! MVals: variables isotopes
       call ini_fillgeom(1,latitude,longitude,(/1.0/))
       call conf_phys(1,llm,nq)