Changeset 2833

Timestamp:

Mar 28, 2017, 4:32:15 PM (7 years ago)

Author:

acozic

Message:

use pseudoalt calculated in the dynamic instead of Alt in phys_output_mod

File:

• LMDZ5/trunk/libf/phylmd/phys_output_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -44 +44 @@
     USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat
     USE print_control_mod, ONLY: prt_level,lunout
-    USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps
+    USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt
     USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
 #ifdef CPP_XIOS
@@ -93 +93 @@
     INTEGER                               :: idayref
     REAL                                  :: zjulian_start, zjulian
-    REAL, DIMENSION(klev)                 :: Alt
     CHARACTER(LEN=4), DIMENSION(nlevSTD)  :: clevSTD
     REAL, DIMENSION(nlevSTD)              :: rlevSTD
@@ -291 +290 @@
     zdtime_moy = dtime         ! Frequence ou l on moyenne
 
-    ! Calcul  Alt
-    DO k=1,klev
-!       Ahyb(k)=(ap(k)+ap(k+1))/2.
-!       Bhyb(k)=(bp(k)+bp(k+1))/2.
-       Alt(k)=log(preff/presnivs(k))*8.
-    ENDDO
-    !          if(prt_level.ge.1) then
-!    WRITE(lunout,*)'Ap Hybrid = ',Ahyb(1:klev)
-!    WRITE(lunout,*)'Bp Hybrid = ',Bhyb(1:klev)
-    WRITE(lunout,*)'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
-    !          ENDIF
 
   ecrit_files(7) = ecrit_files(1)
@@ -347 +335 @@
             levmax(iff) - levmin(iff) + 1, bps)
     CALL wxios_add_vaxis("Alt", &
-            levmax(iff) - levmin(iff) + 1, Alt)
+            levmax(iff) - levmin(iff) + 1, pseudoalt)
    ELSE
     ! NMC files
@@ -422 +410 @@
 
           CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &                       
-               levmax(iff) - levmin(iff) + 1,Alt,nvertAlt(iff))
+               levmax(iff) - levmin(iff) + 1,pseudoalt,nvertAlt(iff))
 
           ELSE