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

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/tracer_mod.F90 (modified) (1 diff)

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

@@ -88 +88 @@
 !-----------------------------------------------------------------------
 
+  contains
+  
+    subroutine ini_tracer_mod(nq,tname)
+      implicit none
+      
+      integer,intent(in) :: nq ! number of tracers
+      character(len=*),intent(in) :: tname(nq) ! tracer names
+      
+      integer :: iq, count
+      character(len=20) :: txt ! to store some text
+      
+      ! set dimension and tracer names
+      nqmx=nq
+#ifndef MESOSCALE
+      allocate(noms(nq))
+      do iq=1,nq
+        noms(iq)=tname(iq)
+        write(*,*) "tracer_mod names : ", trim(noms(iq))
+      enddo
+      
+      ! 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(*,*) "ini_tracer_mod: tracers seem to follow old naming ", &
+                   "convention (q01,q02,...)"
+        write(*,*) "you should run newstart to rename them"
+        stop
+      endif
+#endif
+            
+      ! allocate module arrays:
+      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))
+      
+    end subroutine ini_tracer_mod
+
 end module tracer_mod