Changeset 3296 for LMDZ6/trunk/libf/phylmd/macv2sp.F90

Timestamp:

Mar 23, 2018, 10:13:38 PM (6 years ago)

Author:

oboucher

Message:

a USE statement was missing in mo_simple_plumes.F90
z array has been changed to reflect Stephanie's recommendation
that z is the altitude of the middle of the layers

File:

• LMDZ6/trunk/libf/phylmd/macv2sp.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/macv2sp.F90

@@ -36 +36 @@
   REAL,DIMENSION(klon,klev) :: aod_prof, ssa_prof, asy_prof
   REAL,DIMENSION(klon,klev) :: z, dz
-  REAL,DIMENSION(klon)      :: oro, zrho
+  REAL,DIMENSION(klon)      :: oro, zrho, zt
   !
   INTEGER, PARAMETER :: nmon = 12
@@ -63 +63 @@
   ! define the height and dheight arrays
   !
-  z(:,1)  = 0.                                      ! altitude of surface taken as 0
-  DO k = 1, klev-1
-    oro(:) = pphis(:)/RG                            ! surface height in m
-    zrho(:) = pplay(:, k)/t_seri(:, k)/RD           ! air density in kg/m3
-    dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG  ! layer thickness in m
-    z(:,k+1) = z(:,k) + dz(:,k)                     ! height of interfaces in m, starting from 0 at surface
+  oro(:)  = pphis(:)/RG                             ! surface height in m
+  !
+  DO k = 1, klev
+    zrho(:) = pplay(:,k)/t_seri(:,k)/RD                         ! air density in kg/m3
+    dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG              ! layer thickness in m
+    IF (k==1) THEN
+       z(:,1) = oro(:) + (paprs(:,1)-pplay(:,1))/zrho(:)/RG     ! altitude middle of first layer in m
+       zt(:)  = oro(:) + dz(:,1)                                ! altitude top of first layer in m
+    ELSE
+      z(:,k) = zt(:) + (paprs(:,k)-pplay(:,k))/zrho(:)/RG       ! altitude middle of layer k in m
+      zt(:)  = zt(:) + dz(:,k)                                  ! altitude top of layer k in m
+    ENDIF
   ENDDO
   !