Changeset 1621 for trunk/LMDZ.MARS/libf/phymars/initracer.F

Timestamp:

Oct 25, 2016, 9:23:21 AM (8 years ago)

Author:

emillour

Message:

Further work on full dynamics/physics separation.

LMDZ.COMMON:

• added phy_common/vertical_layers_mod.F90 to store information on vertical grid. This is where routines in the physics should get the information.
• The contents of vertical_layers_mod intialized via dynphy_lonlat/inigeomphy_mod.F90.

LMDZ.MARS:

• physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
• created an "ini_tracer_mod" routine in module "tracer_mod" for a cleaner initialization of the later.
• removed some purely dynamics-related outputs (etot0, zoom parameters, etc.) from diagfi.nc and stats.nc outputs as these informations are not available in the physics.

LMDZ.GENERIC:

• physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
• added nqtot to tracer_h.F90.
• removed some purely dynamics-related outputs (etot0, zoom parameters, etc.) from diagfi.nc and stats.nc outputs as these informations are not available in the physics.

LMDZ.VENUS:

• physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
• added infotrac_phy.F90 to store information on tracers in the physics. Initialized via iniphysiq.
• added cpdet_phy_mod.F90 to store t2tpot etc. functions to be used in the physics. Initialized via iniphysiq. IMPORTANT: there are some hard-coded constants! These should match what is in cpdet_mod.F90 in the dynamics.
• got rid of references to moyzon_mod module within the physics. The required variables (tmoy, plevmoy) are passed to the physics as arguments to physiq.

LMDZ.TITAN:

• added infotrac_phy.F90 to store information on tracers in the physics. Initialized via iniphysiq.
• added cpdet_phy_mod.F90 to store t2tpot etc. functions to be used in the physics.
• Extra work required to completely decouple physics and dynamics: moyzon_mod should be cleaned up and information passed from dynamics to physics as as arguments. Likewise moyzon_ch and moyzon_mu should not be queried from logic_mod (which is in the dynamics).

EM

File:

• trunk/LMDZ.MARS/libf/phymars/initracer.F (modified) (2 diffs)

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

@@ -1 +1 @@
       SUBROUTINE initracer(ngrid,nq,qsurf)
 
-#ifndef MESOSCALE
-       use infotrac, only: tname
-#endif
        use tracer_mod
        USE comcstfi_h
@@ -54 +51 @@
 c-----------------------------------------------------------------------
 
-! Initialization: allocate arrays in tracer_mod
-      allocate(mmol(nq))
-      allocate(radius(nq))
-      allocate(rho_q(nq))
-      allocate(alpha_lift(nq))
-      allocate(alpha_devil(nq))
-      allocate(igcm_dustbin(nq))
-      allocate(nqdust(nq))
-
-#ifndef MESOSCALE
-      allocate(noms(nq))
-! Initialization: get tracer names from the dynamics and check if we are
-!                 using 'old' tracer convention ('q01',q02',...)
-!                 or new convention (full tracer names)
-      ! check if tracers have 'old' names
-
-      count=0
-      do iq=1,nq
-        txt=" "
-        write(txt,'(a1,i2.2)') 'q',iq
-        if (txt.eq.tname(iq)) then
-          count=count+1
-        endif
-      enddo ! of do iq=1,nq
-      
-      if (count.eq.nq) then
-        write(*,*) "initracer: tracers seem to follow old naming ",
-     &             "convention (q01,q02,...)"
-        write(*,*) "you should run newstart to rename them"
-        stop
-      endif
-
-      ! copy tracer names from dynamics
-      do iq=1,nq
-        noms(iq)=tname(iq)
-        write(*,*) "initracer names : ", noms(iq)
-      enddo
-#endif
 
 c------------------------------------------------------------