Changeset 1398 for LMDZ4/trunk/libf/phylmd/physiq.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/physiq.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -1410 +1410 @@
            ema_pcb(i)  = 0.
            ema_pct(i)  = 0.
-           ema_workcbmf(i) = 0.
+c          ema_workcbmf(i) = 0.
           ENDDO
 cIM15/11/02 rajout initialisation ibas_con,itop_con cf. SB =>BEG
@@ -1435 +1435 @@
 c================================================================================
 
-         ENDIF
+         ENDIF !debut
 
            DO i=1,klon
@@ -2122 +2122 @@
      .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
      .        rain_con, snow_con, ibas_con, itop_con, sigd,
-     .        upwd,dnwd,dnwd0,
+     .        ema_cbmf,upwd,dnwd,dnwd0,
      .        Ma,mip,Vprecip,cape,cin,tvp,Tconv,iflagctrl,
      .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd,
@@ -2194 +2194 @@
               endif
             endif
-          ENDDO
-          DO i = 1, klon
-            ema_cbmf(i) = ema_workcbmf(i)
-          ENDDO      
+          ENDDO     
       ELSE IF (iflag_con.eq.0) THEN
           write(lunout,*) 'On n appelle pas la convection'
@@ -3573 +3570 @@
       DO k = 1, klev
       DO i = 1, klon
-        theta(i,k)=t(i,k)*(100000./pplay(i,k))**(RD/RCPD)
+cJYG/IM theta en debut du pas de temps
+cJYG/IM       theta(i,k)=t(i,k)*(100000./pplay(i,k))**(RD/RCPD)
+cJYG/IM theta en fin de pas de temps de physique
+        theta(i,k)=t_seri(i,k)*(100000./pplay(i,k))**(RD/RCPD)
       ENDDO
       ENDDO