Changeset 1398 for LMDZ4/trunk/libf/phylmd/write_histmthNMC.h

Timestamp:

Jun 4, 2010, 6:56:18 PM (14 years ago)

Author:

musat

Message:

Last corrections for CMIP5:

• Add O3 at standard level files histmthNMC.nc
• Add positive attribute "down" for vertical axes for all output files
• Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
• Correct units for vertical axes : mb instead of hPa
• Add mass flux at the bottom of clouds
• Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
• Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
• Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
• Fields albs, albslw become alb1, alb2 in output files
• Correct title for rugs_* fields
• Correct units for pbase and ptop are Pa (not mb)
• Correct ndayrain field

FH/JLD/JYG/MAF/IM

File:

• LMDZ4/trunk/libf/phylmd/write_histmthNMC.h (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/write_histmthNMC.h

@@ -33 +33 @@
         DO i=1, klon
          IF(tnondef(i,k,1).NE.missing_val) THEN
-          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,1))/ecrit_hf2mth
+          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,1))/freq_moyNMC(1)
          ELSE
           zx_tmp_fiNC(i,k) = missing_val
@@ -60 +60 @@
        CALL histwrite_phy(nid_mthnmc,"TxT",itau_w,T2sumSTD(:,:,1))
 c
+       DO k=1, nlevSTD
+        DO i=1, klon
+         IF(O3sumSTD(i,k,1).NE.missing_val) THEN
+          zx_tmp_fiNC(i,k) = O3sumSTD(i,k,1) * 1.e+9
+         ELSE
+          zx_tmp_fiNC(i,k) = missing_val
+         ENDIF
+        ENDDO
+       ENDDO !k=1, nlevSTD
+       CALL histwrite_phy(nid_mthnmc,"tro3",itau_w,
+     $ zx_tmp_fiNC)
+c
+       if (read_climoz == 2) THEN
+       DO k=1, nlevSTD
+        DO i=1, klon
+         IF(O3daysumSTD(i,k,1).NE.missing_val) THEN
+          zx_tmp_fiNC(i,k) = O3daysumSTD(i,k,1) * 1.e+9
+         ELSE
+          zx_tmp_fiNC(i,k) = missing_val
+         ENDIF
+        ENDDO
+       ENDDO !k=1, nlevSTD
+c
+        CALL histwrite_phy(nid_mthnmc,"tro3_daylight",itau_w,
+     $  zx_tmp_fiNC)
+       endif 
+c
       if (ok_sync) then
 c$OMP MASTER