Changeset 2377 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Oct 20, 2015, 3:48:18 PM (9 years ago)

Author:

musat

Message:

Correct zfull and zhalf CMIP6' diagnostics :
zhalf correspond to the interface levels (and paprs) and include

both the top of the model atmosphere and surface levels

zfull is the actual height above mean sea level and corresponds to the

geopotential height (and pplay levels)

ASima/AI/IM

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -979 +979 @@
        CALL histwrite_phy(o_pres, pplay)
        CALL histwrite_phy(o_paprs, paprs(:,1:klev))
-       IF (vars_defined) THEN
-          DO i=1, klon
-             zx_tmp_fi3d1(i,1)= pphis(i)/RG
-             !020611   zx_tmp_fi3d(i,1)= pphis(i)/RG
-          ENDDO
-          DO k=1, klev
-             !020611        DO k=1, klev-1
-             DO i=1, klon
-                !020611         zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (t_seri(i,k) *RD * 
-                zx_tmp_fi3d1(i,k+1)= zx_tmp_fi3d1(i,k) - (t_seri(i,k) *RD *  &
-                     (paprs(i,k+1) - paprs(i,k))) / ( pplay(i,k) * RG ) 
-             ENDDO
-          ENDDO
-       ENDIF
-       CALL histwrite_phy(o_zfull,zx_tmp_fi3d1(:,2:klevp1))
-       !020611    $o_zfull%name,itau_w,zx_tmp_fi3d)
+       CALL histwrite_phy(o_zfull,zphi/RG)
 
        IF (vars_defined)  THEN
-          DO i=1, klon
-             zx_tmp_fi3d(i,1)= pphis(i)/RG - ( &
-                  (t_seri(i,1)+zxtsol(i))/2. *RD * &
-                  (pplay(i,1) - paprs(i,1)))/( (paprs(i,1)+pplay(i,1))/2.* RG)
-          ENDDO
-          DO k=1, klev-1
-             DO i=1, klon
-                zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - ( &
-                     (t_seri(i,k)+t_seri(i,k+1))/2. *RD *  &
-                     (pplay(i,k+1) - pplay(i,k))) / ( paprs(i,k) * RG ) 
-             ENDDO
-          ENDDO
+        zx_tmp_fi3d(:,1)= pphis(:)/RG
+        DO k = 2, klev
+         DO i = 1, klon
+            zx_tmp_fi3d(i,k) = zphi(i,k)/RG + &
+                          (zphi(i,k)-zphi(i,k-1))/RG * &
+                          (paprs(i,k)-pplay(i,k))/(pplay(i,k)-pplay(i,k-1))
+         ENDDO
+        ENDDO
        ENDIF
        CALL histwrite_phy(o_zhalf, zx_tmp_fi3d)