Changeset 2332 for trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.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/dyn1d/testphys1d.F (modified) (6 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, nqperes,nqdesc,iqfils,nqfils,
-     &                    iqpere, nqdesc_tot
+      use infotrac, only: nqtot, tname, nqperes,nqfils
       use comsoil_h, only: volcapa, layer, mlayer, inertiedat, nsoilmx
       use comgeomfi_h, only: sinlat, ini_fillgeom
@@ -268 +267 @@
 
        ! 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))
+       ALLOCATE(nqfils(nqtot))
        nqperes=0
-       nqfils(:)=0
-       nqdesc(:)=0
-       iqfils(:,:)=0
-       iqpere(:)=0
-       nqdesc_tot=0    
+       nqfils(:)=0  
        DO iq=1,nqtot
        if (tnom_transp(iq) == 'air') then
@@ -283 +276 @@
      &          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?
@@ -296 +288 @@
      &   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         
+             nqfils(ipere)=nqfils(ipere)+1         
              continu=.false.
            else !if (tnom_transp(iq) == tnom_0(ipere)) then
@@ -314 +304 @@
        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:
@@ -673 +638 @@
       call phys_state_var_init(1,llm,nq,tname,
      .          day0,time,daysec,dtphys,rad,g,r,cpp,
-     .          nqdesc,iqfils,nqperes,nqfils)! MVals: variables isotopes
+     .          nqperes,nqfils)! MVals: variables isotopes
       call ini_fillgeom(1,latitude,longitude,(/1.0/))
       call conf_phys(1,llm,nq)