Ignore:
Timestamp:
Dec 14, 2001, 3:39:31 PM (23 years ago)
Author:
lmdzadmin
Message:

Correction bug pour KE vectorise FH, JYG
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

    r295 r301  
    689689cKE43
    690690c Initialisation pour la convection de K.E. (sb):
    691          IF (iflag_con.EQ.4) THEN
     691         IF (iflag_con.GE.3) THEN
    692692
    693693         PRINT*, "*** Convection de Kerry Emanuel 4.3  "
     
    12281228     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    12291229     .                "ave(X)", zsto,zout)
    1230 c
    12311230cKE43
    1232       IF (iflag_con .EQ. 4) THEN ! sb
     1231      IF (iflag_con .GE. 3) THEN ! sb
    12331232c
    12341233         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg",
     
    13351334
    13361335         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
    1337      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     1336     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    13381337     .                "ave(X)", zsto,zout)
    13391338
    13401339         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
    1341      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     1340     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    13421341     .                "ave(X)", zsto,zout)
    13431342
     
    14151414c
    14161415cKE43
    1417       IF (iflag_con.EQ.4) THEN ! (sb)
     1416      IF (iflag_con.GE.3) THEN ! (sb)
    14181417c
    14191418         CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s",
     
    19601959         itop_con(i) = klev+1 - kctop(i)
    19611960      ENDDO
    1962       ELSE IF (iflag_con.EQ.3) THEN
    1963           stop'reactiver le call conlmd dans physiq.F'
    1964 c     CALL conccm (dtime,paprs,pplay,t_seri,q_seri,conv_q,
    1965 c    s             d_t_con, d_q_con,
    1966 c    s             rain_con, snow_con, ibas_con, itop_con)
    1967 cKE43
    1968       ELSE IF (iflag_con.EQ.4) THEN
     1961      ELSE IF (iflag_con.GE.3) THEN
    19691962c nb of tracers for the KE convection:
    19701963          if (nqmax .GE. 4) then
     
    19731966              ntra = 1
    19741967          endif
    1975 cke43 (arguments inutiles enleves => des SAVE dans conema43?)
    1976 c$$$          CALL conema43(dtime,paprs,pplay,t_seri,q_seri,
    1977 c$$$     $        u_seri,v_seri,tr_seri,nbtr,
    1978 c$$$     .        ema_workcbmf,
    1979 c$$$     .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
    1980 c$$$     .        wdn, tdn, qdn,
    1981 c$$$     .        rain_con, snow_con, ibas_con, itop_con,
    1982 c$$$     .        upwd,dnwd,dnwd0,bas,top,Ma,cape,tvp,rflag,
    1983 c$$$     .        pbase
    1984 c$$$     .        ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,
    1985 c$$$     .        pori_con,plcl_con,dtma_con,dtlcl_con)
    1986           if (0.eq.1) then ! vectorise
     1968          if (iflag_con.eq.4) then ! vectorise
    19871969          CALL conemav (dtime,paprs,pplay,t_seri,q_seri,
    19881970     .        u_seri,v_seri,tr_seri,nbtr,
     
    19931975c    .        Ma,cape,tvp,(/(nint(rflag(i)),i=1,size(rflag))/),
    19941976     .        Ma,cape,tvp,iflagctrl,
     1977     .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr)
     1978
     1979          else
     1980
     1981          CALL conema (dtime,paprs,pplay,t_seri,q_seri,
     1982     .        u_seri,v_seri,tr_seri,nbtr,
     1983     .        ema_work1,ema_work2,
     1984     .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
     1985     .        rain_con, snow_con, ibas_con, itop_con,
     1986     .        upwd,dnwd,dnwd0,bas,top,
     1987     .        Ma,cape,tvp,rflag,
    19951988     .       pbase
    19961989     .        ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr)
    1997 
    1998           else
    1999 
    2000           CALL conema (dtime,paprs,pplay,t_seri,q_seri,
    2001      $        u_seri,v_seri,tr_seri,nbtr,
    2002      .        ema_work1,ema_work2,
    2003      .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
    2004 c$$$     .        wdn, tdn, qdn,
    2005      .        rain_con, snow_con, ibas_con, itop_con,
    2006      .        upwd,dnwd,dnwd0,bas,top,Ma,cape,tvp,rflag,
    2007      .        pbase
    2008      .        ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr)
    2009 c$$$     .        pori_con,plcl_con,dtma_con,dtlcl_con)
    20101990          endif
    20111991          DO i = 1, klon
     
    20722052      IF (nqmax.GT.2) THEN !--melange convectif de traceurs
    20732053c
    2074           IF (iflag_con .NE. 2 .AND.  iflag_con .NE. 4 ) THEN
     2054          IF (iflag_con .LT. 2 .AND.  iflag_con .GT. 4 ) THEN
    20752055              PRINT*, 'Pour l instant, seul conflx fonctionne ',
    20762056     $            'avec traceurs', iflag_con
    20772057              PRINT*,' Mettre iflag_con',
    2078      $            ' = 2  ou 4 dans run.def et repasser'
     2058     $            ' = 2, 3 ou 4 dans run.def et repasser'
    20792059              CALL abort
    20802060              ENDIF
     
    23482328C Pascale : il faut quand meme apeller phytrac car il gere les sorties
    23492329cKE43       des traceurs => il faut donc mettre des flags a .false.
    2350       IF (iflag_con.EQ.4) THEN
     2330      IF (iflag_con.GE.3) THEN
    23512331c           on ajoute les tendances calculees par KE43
    23522332c$$$ OM on onhibe la convection sur les traceurs
     
    28932873      CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    28942874cKE43
    2895       IF (iflag_con .EQ. 4) THEN ! sb
     2875      IF (iflag_con .GE. 3) THEN ! sb
    28962876c
    28972877      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
     
    29952975      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
    29962976      CALL histwrite(nid_mth,"ptconv",itap,zx_tmp_3d,
    2997      .                                   iim*(jjm+1)*klev,ndex3d)
     2977     .                                   iim*(jjmp1)*klev,ndex3d)
    29982978c
    29992979      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
    30002980      CALL histwrite(nid_mth,"ratqs",itap,zx_tmp_3d,
    3001      .                                   iim*(jjm+1)*klev,ndex3d)
     2981     .                                   iim*(jjmp1)*klev,ndex3d)
    30022982c
    30032983      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
     
    30723052      ENDIF
    30733053cKE43
    3074       IF (iflag_con.EQ.4) THEN ! (sb)
     3054      IF (iflag_con.GE.3) THEN ! (sb)
    30753055c
    30763056      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
Note: See TracChangeset for help on using the changeset viewer.