Changeset 2832 for LMDZ5

Timestamp:

Mar 28, 2017, 3:39:31 PM (7 years ago)

Author:

acozic

Message:

transfer aps and bps to INCA / in phys_output_mod replace Ahyb and Bhyb by aps and bps as decide by email (poihl / 28/03/2017)

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phys_output_mod.F90 (modified) (5 diffs)
• physiq_mod.F90 (modified) (3 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
+    USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps
     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)                 :: Ahyb, Bhyb, Alt
+    REAL, DIMENSION(klev)                 :: Alt
     CHARACTER(LEN=4), DIMENSION(nlevSTD)  :: clevSTD
     REAL, DIMENSION(nlevSTD)              :: rlevSTD
@@ -291 +291 @@
     zdtime_moy = dtime         ! Frequence ou l on moyenne
 
-    ! Calcul des Ahyb, Bhyb et Alt
+    ! Calcul  Alt
     DO k=1,klev
-       Ahyb(k)=(ap(k)+ap(k+1))/2.
-       Bhyb(k)=(bp(k)+bp(k+1))/2.
+!       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,*)'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
@@ -343 +343 @@
             levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
     CALL wxios_add_vaxis("Ahyb", &
-            levmax(iff) - levmin(iff) + 1, Ahyb)
+            levmax(iff) - levmin(iff) + 1, aps)
     CALL wxios_add_vaxis("Bhyb", &
-            levmax(iff) - levmin(iff) + 1, Bhyb)
+            levmax(iff) - levmin(iff) + 1, bps)
     CALL wxios_add_vaxis("Alt", &
             levmax(iff) - levmin(iff) + 1, Alt)
@@ -416 +416 @@
 !!!! Composantes de la coordonnee sigma-hybride 
           CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &
-               levmax(iff) - levmin(iff) + 1,Ahyb,nvertap(iff))
+               levmax(iff) - levmin(iff) + 1,aps,nvertap(iff))
 
           CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &
-               levmax(iff) - levmin(iff) + 1,Bhyb,nvertbp(iff))
+               levmax(iff) - levmin(iff) + 1,bps,nvertbp(iff))
 
           CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &                       

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -249 +249 @@
     use FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
     use ACAMA_GWD_rando_m, only: ACAMA_GWD_rando
+    USE VERTICAL_LAYERS_MOD, ONLY: aps,bps
+
 
     IMPLICIT none
@@ -659 +661 @@
     REAL beta_prec_fisrt(klon,klev) ! taux de conv de l'eau cond (fisrt)
     ! RomP <<<
+    REAL          :: calday
 
     !IM cf FH pour Tiedtke 080604
@@ -4243 +4246 @@
             pphi, &
             pphis, &
-            zx_rh)
+            zx_rh, &
+            aps, bps)
 
        CALL VTe(VTinca)