Changeset 1398 for LMDZ4/trunk/libf/phylmd/calcul_STDlev.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/calcul_STDlev.h (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/calcul_STDlev.h

@@ -107 +107 @@
      &             zx_tmp_fi3d,T2STD)
 
-
+c
+      zx_tmp_fi3d(:,:)=wo(:,:,1) * dobson_u * 1e3 / zmasse / rmo3 * rmd
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
+     &             zx_tmp_fi3d,O3STD)
+c
+      if (read_climoz == 2) THEN
+      zx_tmp_fi3d(:,:)=wo(:,:,2) * dobson_u * 1e3 / zmasse / rmo3 * rmd
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
+     &             zx_tmp_fi3d,O3daySTD)
+      endif
+c
         DO l=1, klev
         DO i=1, klon
@@ -115 +125 @@
         CALL plevel_new(klon,klev,nlevSTD,.true.,zx_tmp_fi3d,rlevSTD,
      &             omega,wlevSTD)
-
 c
 cIM on somme les valeurs toutes les freq_calNMC secondes