Changeset 2786 for LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

Timestamp:

Jan 30, 2017, 11:24:47 AM (8 years ago)

Author:

Ehouarn Millour

Message:

Further convergence with planetary models:

• move vertical_layers_mod.F90 to phy_common and call init_vertical_layers in inigeomphy_mod
• add aps and bps (mid-layer coordinates) to vertical_layers_mod.F90
• same adaptations for the 1D case

EM

File:

• LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

@@ -40 +40 @@
    USE physiq_mod, ONLY: physiq
    USE comvert_mod, ONLY: presnivs, ap, bp, dpres,nivsig, nivsigs, pa, &
-                          preff
+                          preff, aps, bps, pseudoalt, scaleheight
+   USE vertical_layers_mod, ONLY: init_vertical_layers
    USE temps_mod, ONLY: annee_ref, calend, day_end, day_ini, day_ref, &
                         itau_dyn, itau_phy, start_time
@@ -634 +635 @@
         call disvert0(pa,preff,ap,bp,dpres,presnivs,nivsigs,nivsig)
         print *,'On utilise disvert0'
+        aps(1:llm)=0.5*(ap(1:llm)+ap(2:llm+1))
+        bps(1:llm)=0.5*(bp(1:llm)+bp(2:llm+1))
+        scaleheight=8.
+        pseudoalt(1:llm)=-scaleheight*log(presnivs(1:llm)/preff)
       ELSE
         call disvert()
@@ -641 +646 @@
       ENDIF
       ! initialize ap,bp, etc. in vertical_layers_mod
+      CALL init_vertical_layers(llm,preff,scaleheight, &
+                                ap,bp,aps,bps,presnivs,pseudoalt)
+
       sig_s=presnivs/preff
       plev =ap+bp*psurf