Changeset 1412 for LMDZ4

Timestamp:

Jul 9, 2010, 3:06:40 PM (14 years ago)

Author:

idelkadi

Message:

Nettoyage dans physiq.F des anciennes versions des options iflag_cldcon=5 ou 6.

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (21 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -1 +1 @@
 ! $Id$
-!
 c#define IO_DEBUG
 
@@ -262 +261 @@
       CHARACTER*4 bb2
       CHARACTER*2 bb3
-c
 
       real twriteSTD(klon,nlevSTD,nfiles)
@@ -865 +863 @@
       REAL rflag(klon)          ! flag fonctionnement de convect
       INTEGER iflagctrl(klon)          ! flag fonctionnement de convect
+
 c -- convect43:
       INTEGER ntra              ! nb traceurs pour convect4.3
@@ -1233 +1232 @@
       call phys_state_var_init(read_climoz)
       call phys_output_var_init
-
       print*, '================================================='
 cIM for NMC files
@@ -1253 +1251 @@
 c         pmflxr=0.
 c         pmflxs=0.
+
         itau_con=0
         first=.false.
@@ -1424 +1423 @@
            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
@@ -1449 +1448 @@
 c================================================================================
 
-         ENDIF
+         ENDIF !debut
 
            DO i=1,klon
@@ -2171 +2170 @@
      .        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,
@@ -2190 +2189 @@
 
           ELSE ! ok_cvl
+
 c MAF conema3 ne contient pas les traceurs
           CALL conema3 (dtime,
@@ -2229 +2229 @@
 
           DO i = 1, klon
-            ema_pcb(i)  = pbase(i)
+            ema_pcb(i)  = paprs(i,ibas_con(i))
           ENDDO
           DO i = 1, klon
-
 ! L'idicage de itop_con peut cacher un pb potentiel
 ! FH sous la dictee de JYG, CR
@@ -2243 +2242 @@
               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'
@@ -2463 +2459 @@
 c  ==============
 
-! Dans le cas où on active les thermiques, on fait partir l'ajustement
+! Dans le cas où on active les thermiques, on fait partir l'ajustement
 ! a partir du sommet des thermiques.
 ! Dans le cas contraire, on demarre au niveau 1.
@@ -2547 +2543 @@
      s   (klon,klev,q_seri,zqsat,clwcon0th,ptconvth,ratqsc,rnebcon0th)
          if(prt_level.ge.9) print*,' CLOUDS_GNO OK'
-
-      endif
+       
+       endif
 
 c   ratqs stables
@@ -2611 +2607 @@
 
       if (iflag_cldcon.eq.1 .or.iflag_cldcon.eq.2
-     s    .or.iflag_cldcon.eq.4) then
+     s    .or.iflag_cldcon.ge.4) then
 
 ! On ajoute une constante au ratqsc*2 pour tenir compte de 
@@ -2646 +2642 @@
       endif
 
-      print*,'PHSYIQ NUAGES4'
 
 c
@@ -2842 +2837 @@
       else
 ! Ancienne version
-         cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)
-         cldliq(:,:)=cldliq(:,:)+rnebcon(:,:)*clwcon(:,:)
+      cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)
+      cldliq(:,:)=cldliq(:,:)+rnebcon(:,:)*clwcon(:,:)
       endif
 
@@ -3312 +3307 @@
 !     s        ref_liq,ref_ice
           call phys_cosp(itap,dtime,freq_cosp,
-     $                 ecrit_mth,ecrit_day,ecrit_hf,overlap,
-     $                   klon,klev,rlon,rlat,presnivs,
+     $                   ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP,
+     $                   ecrit_mth,ecrit_day,ecrit_hf,
+     $                   klon,klev,rlon,rlat,presnivs,overlap,
      $                   ref_liq,ref_ice,
      $                   pctsrf(:,is_ter)+pctsrf(:,is_lic),
@@ -3322 +3318 @@
      $                   pmflxr(:,1:klev),pmflxs(:,1:klev),
      $                   mr_ozone,cldtau, cldemi)
+
 !     L          calipso2D,calipso3D,cfadlidar,parasolrefl,atb,betamol,
 !     L          cfaddbze,clcalipso2,dbze,cltlidarradar,
@@ -3469 +3466 @@
       wwriteSTD(:,:,4)=wlevSTD(:,:)
 c
+cIM initialisation 5eme fichier de sortie 
 cIM ajoute 5eme niveau 170310 BEG
       twriteSTD(:,:,5)=tlevSTD(:,:)
@@ -3581 +3579 @@
 cIM global posePB#include "write_bilKP_ave.h"
 c
+
 c Sauvegarder les valeurs de t et q a la fin de la physique:
 c
@@ -3639 +3638 @@
       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