Changeset 619 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Apr 15, 2005, 2:36:17 PM (19 years ago)

Author:

lmdzadmin

Message:

Rajout convection Kerry Emanuel pour traceurs- MAF+JYG

File:

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

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -212 +212 @@
       REAL d_qx(klon,klev,nqmax)
       REAL d_ps(klon)
+      real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
 
       INTEGER klevp1, klevm1
@@ -769 +770 @@
       REAL zx_t, zx_qs, zdelta, zcor, zfra, zlvdcp, zlsdcp
       real zqsat(klon,klev)
-      INTEGER i, k, iq, ig, j, nsrf, ll
+      INTEGER i, k, iq, ig, j, iiq, nsrf, ll
       REAL t_coup
       PARAMETER (t_coup=234.0)
@@ -795 +796 @@
       INTEGER iflagctrl(klon)          ! flag fonctionnement de convect
 c -- convect43:
-      INTEGER ntra              ! nb traceurs pour convect4.3
+      INTEGER ntra              ! nb traceurs 
       REAL pori_con(klon)    ! pressure at the origin level of lifted parcel
       REAL plcl_con(klon),dtma_con(klon),dtlcl_con(klon)
@@ -1153 +1154 @@
 
          WRITE(lunout,*)"*** Convection de Kerry Emanuel 4.3  "
-         WRITE(lunout,*)
-     .      "On va utiliser le melange convectif des traceurs qui"
-         WRITE(lunout,*)"est calcule dans convect4.3"
-         WRITE(lunout,*)" !!! penser aux logical flags de phytrac"
 
           DO i = 1, klon
@@ -1167 +1164 @@
           DO i = 1, klon
            ibas_con(i) = 1
-           itop_con(i) = klev+1
+           itop_con(i) = 1
           ENDDO
 cIM15/11/02 rajout initialisation ibas_con,itop_con cf. SB =>END
@@ -1733 +1730 @@
       ELSE IF (iflag_con.GE.3) THEN
 c nb of tracers for the KE convection:
-          if (nqmax .GE. 4) then
-              ntra = nbtr
-          else
-              ntra = 1 
-          endif
-c
+c MAF la partie traceurs est faite dans phytrac
+c on met ntra=1 pour limiter les appels mais on peut
+c supprimer les calculs / ftra.
+              ntra = 1
 c sb, oct02:
 c Schema de convection modularise et vectorise:
@@ -1747 +1742 @@
           CALL concvl (iflag_con,
      .        dtime,paprs,pplay,t_seri,q_seri,
-     .        u_seri,v_seri,tr_seri,nbtr,
+     .        u_seri,v_seri,tr_seri,ntra,
      .        ema_work1,ema_work2,
      .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
@@ -1753 +1748 @@
      .        upwd,dnwd,dnwd0,
      .        Ma,cape,tvp,iflagctrl,
-     .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd)
+     .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd,
+     .        pmflxr,pmflxs,
+     .        da,phi,mp)
+
 cIM cf. FH
               clwcon0=qcondc
+              pmfu(:,:)=upwd(:,:)+dnwd(:,:)
 
           ELSE ! ok_cvl
-
+c MAF conema3 ne contient pas les traceurs
           CALL conema3 (dtime,
      .        paprs,pplay,t_seri,q_seri,
-     .        u_seri,v_seri,tr_seri,nbtr,
+     .        u_seri,v_seri,tr_seri,ntra,
      .        ema_work1,ema_work2,
      .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
@@ -1885 +1884 @@
       zx_ajustq=.FALSE.
 c
-      IF (nqmax.GT.2) THEN !--melange convectif de traceurs
-c
-          IF (iflag_con .NE. 2 .AND. debut) THEN 
-              WRITE(lunout,*)'Pour l instant, seul conflx fonctionne ',
-     $            'avec traceurs', iflag_con
-              WRITE(lunout,*)' Mettre iflag_con', 
-     $            ' = 2 dans run.def et repasser'
-c              CALL abort
-              ENDIF 
-c
-      ENDIF !--nqmax.GT.2
-c
-c Appeler l'ajustement sec
-c
 c===================================================================
 c Convection seche (thermiques ou ajustement)
@@ -2644 +2629 @@
 c   Calcul  des tendances traceurs
 c====================================================================
-C Pascale : il faut quand meme apeller phytrac car il gere les sorties
-cKE43       des traceurs => il faut donc mettre des flags a .false.
-      IF (iflag_con.GE.3) THEN
-c           on ajoute les tendances calculees par KE43
-cXXX OM on onhibe la convection sur les traceurs
-        DO iq=1, nqmax-2 ! Sandrine a -3 ???
-cXXX OM on inhibe la convection sur les traceur
-cXXX        DO k = 1, nlev
-cXXX        DO i = 1, klon
-cXXX          tr_seri(i,k,iq) = tr_seri(i,k,iq) + d_tr(i,k,iq)
-cXXX        ENDDO
-cXXX        ENDDO
-        WRITE(iqn,'(i2.2)') iq
-        CALL minmaxqfi(tr_seri(1,1,iq),0.,1.e33,'couche lim iq='//iqn)
-        ENDDO
-CMAF modif pour garder info du nombre de traceurs auxquels
-C la physique s'applique
-      ELSE
-CMAF modif pour garder info du nombre de traceurs auxquels
-C la physique s'applique
 C
-      call phytrac (rnpb,
+      call phytrac (iflag_con,rnpb,
      I                   itap, julien, gmtime,
      I                   debut,lafin,
@@ -2680 +2645 @@
      I                   ibas_con,
      I                   pmflxr,pmflxs,prfl,psfl,
-#ifdef INCA
+     I                   da,phi,mp,
+     I                   upwd,dnwd,
+#ifdef INCA_CH4
      I                   flxmass_w,
 #endif
      O                   tr_seri)
-      ENDIF 
 
       IF (offline) THEN