Changeset 1621 for trunk/LMDZ.GENERIC/libf/dyn3d/disvert.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.GENERIC/libf/dyn3d/disvert.F (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/dyn3d/disvert.F

@@ -5 +5 @@
       USE callkeys_mod, ONLY: kastprof,pceil
       USE comvert_mod, ONLY: ap,bp,sig,nivsigs,nivsig,pa,preff,
-     .                  aps,bps,presnivs,pseudoalt
+     .                  aps,bps,presnivs,pseudoalt,scaleheight
 
 c    Auteur :  F. Forget Y. Wanherdrick, P. Levan
@@ -31 +31 @@
       INTEGER l,ll
       REAL snorm
-      REAL alpha,beta,gama,delta,deltaz,h,quoi,quand
+      REAL alpha,beta,gama,delta,deltaz,quoi,quand
       REAL zsig(llm)
       INTEGER np,ierr
@@ -108 +108 @@
          PRINT*,'****************************'
 
-         READ(99,*) h
+         READ(99,*) scaleheight
          do l=1,llm
             read(99,*) zsig(l)
@@ -117 +117 @@
          if(autozlevs)then
             open(91,file="z2sig.def",form='formatted')
-            read(91,*) h
+            read(91,*) scaleheight
             DO l=1,llm-2
                read(91,*) Hmax
@@ -126 +126 @@
             print*,'Auto-shifting h in disvert.F to:'
 !            h = Hmax / log(psurf/100.0)
-            h = Hmax / log(psurf/pceil)
-            print*,'h = ',h,' km'
+            scaleheight = Hmax / log(psurf/pceil)
+            print*,'h = ',scaleheight,' km'
         endif
        
         sig(1)=1
         do l=2,llm
-           sig(l) = 0.5 * ( exp(-zsig(l)/h) + exp(-zsig(l-1)/h) )
+           sig(l) = 0.5 * ( exp(-zsig(l)/scaleheight) + 
+     &                      exp(-zsig(l-1)/scaleheight) )
         end do
         sig(llm+1)=0
@@ -223 +224 @@
       DO l = 1, llm
        presnivs(l) = aps(l)+bps(l)*preff
-       pseudoalt(l) = -h*log(presnivs(l)/preff)
+       pseudoalt(l) = -scaleheight*log(presnivs(l)/preff)
       ENDDO