Changeset 1398 for LMDZ4/trunk/libf/phylmd/concvl.F

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/concvl.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/concvl.F

@@ -5 +5 @@
      .             d_t,d_q,d_u,d_v,d_tra,
      .             rain, snow, kbas, ktop, sigd,
-     .             upwd,dnwd,dnwdbis,Ma,mip,Vprecip,
+     .             cbmf,upwd,dnwd,dnwdbis,Ma,mip,Vprecip,
      .             cape,cin,tvp,Tconv,iflag,
      .             pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,
@@ -120 +120 @@
        INTEGER i,k,itra
        REAL qs(klon,klev),qs_wake(klon,klev)
-cLF       REAL cbmf(klon)
+       REAL cbmf(klon)
 cLF       SAVE cbmf
-       REAL, SAVE, ALLOCATABLE :: cbmf(:)
-c$OMP THREADPRIVATE(cbmf)!       
+cIM/JYG      REAL, SAVE, ALLOCATABLE :: cbmf(:)
+ccc$OMP THREADPRIVATE(cbmf)!       
        REAL cbmflast(klon)
        INTEGER ifrst
@@ -185 +185 @@
 c Allocate some variables LF 04/2008
 c
-        allocate(cbmf(klon))
+cIM/JYG allocate(cbmf(klon))
         allocate(ql(klon,klev))
         allocate(t1(klon,klev))
@@ -341 +341 @@
      $              dtime,Ma,upwd,dnwd,dnwdbis,qcondc,wd,cape,
      $              da,phi,mp)
-      
+c
+      DO i = 1,klon
+        cbmf(i) = Ma(i,kbas(i))
+      ENDDO
+c
       else