! $Header$

                                                                        ! RomP

    evap1, ep1, epmlmmm1, eplamm1, & ! RomP

    wdtraina1, wdtrainm1, & ! RomP

    epmax_diag1) ! epmax_cape

  USE dimphy

  IMPLICIT NONE

  ! .............................START PROLOGUE............................

  ! All argument names (except len,nd,ntra,nloc,delt and the flags) have a

  ! "1" appended.

  ! The "1" is removed for the corresponding compressed (local) variables.

  ! PARAMETERS:

  ! Name            Type         Usage            Description

  ! ----------      ----------     -------  ----------------------------

  ! len           Integer        Input        first (i) dimension

  ! nd            Integer        Input        vertical (k) dimension

  ! ndp1          Integer        Input        nd + 1

  ! ntra          Integer        Input        number of tracors

  ! iflag_con     Integer        Input        version of convect (3/4)

  ! t1            Real           Input        temperature

  ! q1            Real           Input        specific hum

  ! qs1           Real           Input        sat specific hum

  ! u1            Real           Input        u-wind

  ! v1            Real           Input        v-wind

  ! tra1          Real           Input        tracors

  ! p1            Real           Input        full level pressure

  ! ph1           Real           Input        half level pressure

  ! iflag1        Integer        Output       flag for Emanuel conditions

  ! ft1           Real           Output       temp tend

  ! fq1           Real           Output       spec hum tend

  ! fu1           Real           Output       u-wind tend

  ! fv1           Real           Output       v-wind tend

  ! ftra1         Real           Output       tracor tend

  ! precip1       Real           Output       precipitation

  ! VPrecip1      Real           Output       vertical profile of

  ! precipitations

  ! cbmf1         Real           Output       cloud base mass flux

  ! sig1          Real           In/Out       section adiabatic updraft

  ! w01           Real           In/Out       vertical velocity within adiab

  ! updraft

  ! delt          Real           Input        time step

  ! Ma1           Real           Output       mass flux adiabatic updraft

  ! upwd1         Real           Output       total upward mass flux

  ! (adiab+mixed)

  ! dnwd1         Real           Output       saturated downward mass flux

  ! (mixed)

  ! dnwd01        Real           Output       unsaturated downward mass flux

  ! qcondc1       Real           Output       in-cld mixing ratio of

  ! condensed water

  ! wd1           Real           Output       downdraft velocity scale for

  ! sfc fluxes

  ! cape1         Real           Output       CAPE

  ! wdtrainA1     Real           Output   precipitation detrained from

  ! adiabatic draught;

  ! used in tracer transport (cvltr)

  ! wdtrainM1     Real           Output   precipitation detrained from mixed

  ! draughts;

  ! used in tracer transport (cvltr)

  ! da1           Real           Output   used in tracer transport (cvltr)

  ! phi1          Real           Output   used in tracer transport (cvltr)

  ! mp1           Real           Output   used in tracer transport (cvltr)

  ! phi21         Real           Output   used in tracer transport (cvltr)

  ! d1a1          Real           Output   used in tracer transport (cvltr)

  ! dam1          Real           Output   used in tracer transport (cvltr)

  ! evap1         Real           Output

  ! ep1           Real           Output

  ! sij1        Real           Output

  ! elij1         Real           Output

  ! S. Bony, Mar 2002:

  ! * Several modules corresponding to different physical processes

  ! * Several versions of convect may be used:

  ! - iflag_con=3: version lmd  (previously named convect3)

  ! - iflag_con=4: version 4.3b (vect. version, previously convect1/2)

  ! + tard: 	- iflag_con=5: version lmd with ice (previously named convectg)

  ! S. Bony, Oct 2002:

  ! * Vectorization of convect3 (ie version lmd)

  ! ..............................END PROLOGUE.............................

  ! Input

  INTEGER len

  INTEGER nd

  INTEGER ndp1

  INTEGER noff

  INTEGER iflag_con

  INTEGER ntra

  REAL delt

  REAL t1(len, nd)

  REAL q1(len, nd)

  REAL qs1(len, nd)

  REAL u1(len, nd)

  REAL v1(len, nd)

  REAL tra1(len, nd, ntra)

  REAL p1(len, nd)

  REAL ph1(len, ndp1)

  ! Output

  INTEGER iflag1(len)

  REAL ft1(len, nd)

  REAL fq1(len, nd)

  REAL fu1(len, nd)

  REAL fv1(len, nd)

  REAL ftra1(len, nd, ntra)

  REAL precip1(len)

  REAL cbmf1(len)

  REAL sig1(klon, klev)

  REAL w01(klon, klev)

  REAL vprecip1(len, nd+1)

  REAL evap1(len, nd) !RomP

  REAL ep1(len, nd) !RomP

  REAL ma1(len, nd)

  REAL upwd1(len, nd)

  REAL dnwd1(len, nd)

  REAL dnwd01(len, nd)

  REAL qcondc1(len, nd) ! cld

  REAL wd1(len) ! gust

  REAL cape1(len)

  ! RomP >>>

  REAL wdtraina1(len, nd), wdtrainm1(len, nd)

  REAL sij1(len, nd, nd), elij1(len, nd, nd)

  REAL da1(len, nd), phi1(len, nd, nd), mp1(len, nd)

  REAL phi21(len, nd, nd)

  REAL d1a1(len, nd), dam1(len, nd)

  REAL epmlmmm1(len, nd, nd), eplamm1(len, nd)

  ! RomP <<<

  REAL epmax_diag1 (len) ! epmax_cape

  ! -------------------------------------------------------------------

  ! Original Prologue by Kerry Emanuel.

  ! -------------------------------------------------------------------

  ! --- ARGUMENTS

  ! -------------------------------------------------------------------

  ! --- On input:

  ! t:   Array of absolute temperature (K) of dimension ND, with first

  ! index corresponding to lowest model level. Note that this array

  ! will be altered by the subroutine if dry convective adjustment

  ! occurs and if IPBL is not equal to 0.

  ! q:   Array of specific humidity (gm/gm) of dimension ND, with first

  ! index corresponding to lowest model level. Must be defined

  ! at same grid levels as T. Note that this array will be altered

  ! if dry convective adjustment occurs and if IPBL is not equal to 0.

  ! qs:  Array of saturation specific humidity of dimension ND, with first

  ! index corresponding to lowest model level. Must be defined

  ! at same grid levels as T. Note that this array will be altered

  ! if dry convective adjustment occurs and if IPBL is not equal to 0.

  ! u:   Array of zonal wind velocity (m/s) of dimension ND, witth first

  ! index corresponding with the lowest model level. Defined at

  ! same levels as T. Note that this array will be altered if

  ! dry convective adjustment occurs and if IPBL is not equal to 0.

  ! v:   Same as u but for meridional velocity.

  ! tra: Array of passive tracer mixing ratio, of dimensions (ND,NTRA),

  ! where NTRA is the number of different tracers. If no

  ! convective tracer transport is needed, define a dummy

  ! input array of dimension (ND,1). Tracers are defined at

  ! same vertical levels as T. Note that this array will be altered

  ! if dry convective adjustment occurs and if IPBL is not equal to 0.

  ! p:   Array of pressure (mb) of dimension ND, with first

  ! index corresponding to lowest model level. Must be defined

  ! at same grid levels as T.

  ! ph:  Array of pressure (mb) of dimension ND+1, with first index

  ! corresponding to lowest level. These pressures are defined at

  ! levels intermediate between those of P, T, Q and QS. The first

  ! value of PH should be greater than (i.e. at a lower level than)

  ! the first value of the array P.

  ! nl:  The maximum number of levels to which convection can penetrate, plus

  ! 1.

  ! NL MUST be less than or equal to ND-1.

  ! delt: The model time step (sec) between calls to CONVECT

  ! ----------------------------------------------------------------------------

  ! ---   On Output:

  ! iflag: An output integer whose value denotes the following:

  ! VALUE   INTERPRETATION

  ! -----   --------------

  ! 0     Moist convection occurs.

  ! 1     Moist convection occurs, but a CFL condition

  ! on the subsidence warming is violated. This

  ! does not cause the scheme to terminate.

  ! 2     Moist convection, but no precip because ep(inb) lt 0.0001

  ! 3     No moist convection because new cbmf is 0 and old cbmf is 0.

  ! 4     No moist convection; atmosphere is not

  ! unstable

  ! 6     No moist convection because ihmin le minorig.

  ! 7     No moist convection because unreasonable

  ! parcel level temperature or specific humidity.

  ! 8     No moist convection: lifted condensation

  ! level is above the 200 mb level.

  ! 9     No moist convection: cloud base is higher

  ! then the level NL-1.

  ! ft:   Array of temperature tendency (K/s) of dimension ND, defined at

  ! same

  ! grid levels as T, Q, QS and P.

  ! fq:   Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,

  ! defined at same grid levels as T, Q, QS and P.

  ! fu:   Array of forcing of zonal velocity (m/s^2) of dimension ND,

  ! defined at same grid levels as T.

  ! fv:   Same as FU, but for forcing of meridional velocity.

  ! ftra: Array of forcing of tracer content, in tracer mixing ratio per

  ! second, defined at same levels as T. Dimensioned (ND,NTRA).

  ! precip: Scalar convective precipitation rate (mm/day).

  ! VPrecip: Vertical profile of convective precipitation (kg/m2/s).

  ! wd:   A convective downdraft velocity scale. For use in surface

  ! flux parameterizations. See convect.ps file for details.

  ! tprime: A convective downdraft temperature perturbation scale (K).

  ! For use in surface flux parameterizations. See convect.ps

  ! file for details.

  ! qprime: A convective downdraft specific humidity

  ! perturbation scale (gm/gm).

  ! For use in surface flux parameterizations. See convect.ps

  ! file for details.

  ! cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST

  ! BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT

  ! ITS NEXT CALL. That is, the value of CBMF must be "remembered"

  ! by the calling program between calls to CONVECT.

  ! det:   Array of detrainment mass flux of dimension ND.

  ! -------------------------------------------------------------------

  ! Local arrays

  INTEGER i, k, n, il, j

  INTEGER icbmax

  INTEGER icb1(klon)

  INTEGER inb1(klon)

  REAL pnk1(klon)

  REAL qsnk1(klon)

  REAL clw1(klon, klev)

  INTEGER ncum

  ! (local) compressed fields:

  ! ym      integer nloc

  ! ym      parameter (nloc=klon) ! pour l'instant

#define nloc klon

  ! RomP >>>

  REAL sigd(nloc)

  ! RomP <<<

  ! -------------------------------------------------------------------

  ! --- SET CONSTANTS AND PARAMETERS

  ! -------------------------------------------------------------------

  ! print *, '-> cv_driver'      !jyg

  ! -- set simulation flags:

  ! (common cvflag)

  CALL cv_flag(0)

  ! -- set thermodynamical constants:

  ! (common cvthermo)

  ! -- set convect parameters

  ! includes microphysical parameters and parameters that

  ! control the rate of approach to quasi-equilibrium)

  ! (common cvparam)

  END IF

  END IF

  ! ---------------------------------------------------------------------

  ! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS

  ! ---------------------------------------------------------------------

  ! ym

    END DO

  END IF

  ! RomP >>>

  ! RomP <<<

  ! --------------------------------------------------------------------

  ! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY

  ! --------------------------------------------------------------------

    ! print*,'Emanuel version 30 '

    CALL cv30_prelim(len, nd, ndp1, t1, q1, p1, ph1 & ! nd->na

  END IF

    CALL cv_prelim(len, nd, ndp1, t1, q1, p1, ph1, lv1, cpn1, tv1, gz1, h1, &

  END IF

  ! --------------------------------------------------------------------

  ! --- CONVECTIVE FEED

  ! --------------------------------------------------------------------

    CALL cv30_feed(len, nd, t1, q1, qs1, p1, ph1, hm1, gz1 & !

                                                             ! nd->na

  END IF

    CALL cv_feed(len, nd, t1, q1, qs1, p1, hm1, gz1, nk1, icb1, icbmax, &

  END IF

  ! --------------------------------------------------------------------

  ! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part

  ! (up through ICB for convect4, up through ICB+1 for convect3)

  ! Calculates the lifted parcel virtual temperature at nk, the

  ! actual temperature, and the adiabatic liquid water content.

  ! --------------------------------------------------------------------

    CALL cv30_undilute1(len, nd, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1 & ! nd->na

  END IF

    CALL cv_undilute1(len, nd, t1, q1, qs1, gz1, p1, nk1, icb1, icbmax, tp1, &

  END IF

  ! -------------------------------------------------------------------

  ! --- TRIGGERING

  ! -------------------------------------------------------------------

    CALL cv30_trigger(len, nd, icb1, plcl1, p1, th1, tv1, tvp1 & !

                                                                 ! nd->na

  END IF

  END IF

  ! =====================================================================

  ! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY

  ! =====================================================================

    END IF

  END DO

  ! print*,'cv_driver : klon, ncum = ',len,ncum

    ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

    ! --- COMPRESS THE FIELDS

    ! (-> vectorization over convective gridpoints)

    ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

      CALL cv30_compress(len, nloc, ncum, nd, ntra, iflag1, nk1, icb1, icbs1, &

        plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, v1, &

        gz1, th1, tra1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, sig1, &

        w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, buoybase, t, &

        q, qs, u, v, gz, th, tra, h, lv, cpn, p, ph, tv, tp, tvp, clw, sig, &

    END IF

      CALL cv_compress(len, nloc, ncum, nd, iflag1, nk1, icb1, cbmf1, plcl1, &

        tnk1, qnk1, gznk1, t1, q1, qs1, u1, v1, gz1, h1, lv1, cpn1, p1, ph1, &

        tv1, tp1, tvp1, clw1, iflag, nk, icb, cbmf, plcl, tnk, qnk, gznk, t, &

    END IF

    ! -------------------------------------------------------------------

    ! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :

    ! ---   FIND THE REST OF THE LIFTED PARCEL TEMPERATURES

    ! ---   &

    ! ---   COMPUTE THE PRECIPITATION EFFICIENCIES AND THE

    ! ---   FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD

    ! ---   &

    ! ---   FIND THE LEVEL OF NEUTRAL BUOYANCY

    ! -------------------------------------------------------------------

      CALL cv30_undilute2(nloc, ncum, nd, icb, icbs, nk & !na->nd

        , tnk, qnk, gznk, t, q, qs, gz, p, h, tv, lv, pbase, buoybase, plcl, &

    END IF

      CALL cv_undilute2(nloc, ncum, nd, icb, nk, tnk, qnk, gznk, t, q, qs, &

    END IF

    ! -------------------------------------------------------------------

    ! --- CLOSURE

    ! -------------------------------------------------------------------

      CALL cv30_closure(nloc, ncum, nd, icb, inb & ! na->nd

      ! epmax_cape

      call cv30_epmax_fn_cape(nloc,ncum,nd &

                ,cape,ep,hp,icb,inb,clw,nk,t,h,lv &

        ! on écrase ep et recalcule hp

    END IF

      CALL cv_closure(nloc, ncum, nd, nk, icb, tv, tvp, p, ph, dph, plcl, &

    END IF

    ! -------------------------------------------------------------------

    ! --- MIXING

    ! -------------------------------------------------------------------

      CALL cv30_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb & !

                                                                ! na->nd

        , ph, t, q, qs, u, v, tra, h, lv, qnk, hp, tv, tvp, ep, clw, m, sig, &

    END IF

      CALL cv_mixing(nloc, ncum, nd, icb, nk, inb, inbis, ph, t, q, qs, u, v, &

        h, lv, qnk, hp, tv, tvp, ep, clw, cbmf, m, ment, qent, uent, vent, &

    END IF

    ! -------------------------------------------------------------------

    ! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS

    ! -------------------------------------------------------------------

      ! RomP >>>

      CALL cv30_unsat(nloc, ncum, nd, nd, ntra, icb, inb & ! na->nd

        , t, q, qs, gz, u, v, tra, p, ph, th, tv, lv, cpn, ep, sigp, clw, m, &

        ment, elij, delt, plcl, mp, qp, up, vp, trap, wt, water, evap, b, &

      ! RomP <<<

    END IF

      CALL cv_unsat(nloc, ncum, nd, inb, t, q, qs, gz, u, v, p, ph, h, lv, &

    END IF

    ! -------------------------------------------------------------------

    ! --- YIELD

    ! (tendencies, precipitation, variables of interface with other

    ! processes, etc)

    ! -------------------------------------------------------------------

      CALL cv30_yield(nloc, ncum, nd, nd, ntra & ! na->nd

        , icb, inb, delt, t, q, u, v, tra, gz, p, ph, h, hp, lv, cpn, th, ep, &

        clw, m, tp, mp, qp, up, vp, trap, wt, water, evap, b, ment, qent, &

        uent, vent, nent, elij, traent, sig, tv, tvp, iflag, precip, vprecip, &

        ft, fq, fu, fv, ftra, upwd, dnwd, dnwd0, ma, mike, tls, tps, qcondc, &

    END IF

      CALL cv_yield(nloc, ncum, nd, nk, icb, inb, delt, t, q, u, v, gz, p, &

        ph, h, hp, lv, cpn, ep, clw, frac, m, mp, qp, up, vp, wt, water, &

        evap, ment, qent, uent, vent, nent, elij, tv, tvp, iflag, wd, qprime, &

    END IF

    ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

    ! --- passive tracers

    ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

      ! RomP >>>

      CALL cv30_tracer(nloc, len, ncum, nd, nd, ment, sij, da, phi, phi2, &

      ! RomP <<<

    END IF

    ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

    ! --- UNCOMPRESS THE FIELDS

    ! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

    ! set iflag1 =42 for non convective points

    END DO

      CALL cv30_uncompress(nloc, len, ncum, nd, ntra, idcum, iflag, precip, &

        vprecip, evap, ep, sig, w0 & !RomP

        , ft, fq, fu, fv, ftra, inb, ma, upwd, dnwd, dnwd0, qcondc, wd, cape, &

        da, phi, mp, phi2, d1a, dam, sij & !RomP

        , elij, clw, epmlmmm, eplamm & !RomP

        , wdtraina, wdtrainm,epmax_diag &     !RomP

        , iflag1, precip1, vprecip1, evap1, ep1, sig1, w01 & !RomP

        , ft1, fq1, fu1, fv1, ftra1, inb1, ma1, upwd1, dnwd1, dnwd01, &

        qcondc1, wd1, cape1, da1, phi1, mp1, phi21, d1a1, dam1, sij1 & !RomP

        , elij1, clw1, epmlmmm1, eplamm1 & !RomP

    END IF

      CALL cv_uncompress(nloc, len, ncum, nd, idcum, iflag, precip, cbmf, ft, &

        fq, fu, fv, ma, qcondc, iflag1, precip1, cbmf1, ft1, fq1, fu1, fv1, &

    END IF

  END IF ! ncum>0

  ! print *, 'fin cv_driver ->'      !jyg

END SUBROUTINE cv_driver

! ==================================================================

  USE ioipsl_getin_p_mod, ONLY : getin_p

  IMPLICIT NONE

  ! Argument : iflag_ice_thermo : ice thermodynamics is taken into account if

  ! iflag_ice_thermo >=1

  INTEGER iflag_ice_thermo

  include "cvflag.h"

  ! -- si .TRUE., on rend la gravite plus explicite et eventuellement

  ! differente de 10.0 dans convect3:

  !

! si icvflag_Tpa=0, alors la fraction de glace dans l'ascendance adiabatique est

  ! fonction de la temperature de l'environnement et la temperature de l'ascendance est

  ! calculee en deux it�rations, une en supposant qu'il n'y a pas de glace et l'autre

  ! en ajoutant la glace (ancien sch�ma d'Arnaud Jam).

! si icvflag_Tpa=1, alors la fraction de glace dans l'ascendance adiabatique est

  ! fonction de la temperature de l'environnement et la temperature de l'ascendance est

  ! calculee en une seule iteration.

! si icvflag_Tpa=2, alors la fraction de glace dans l'ascendance adiabatique est

  ! fonction de la temperature de l'ascendance et la temperature de l'ascendance est

  ! calculee en une seule iteration.

END SUBROUTINE cv_flag

! ==================================================================

  IMPLICIT NONE

  ! -------------------------------------------------------------

  ! Set thermodynamical constants for convectL

  ! -------------------------------------------------------------

  include "YOMCST.h"

  include "cvthermo.h"

  INTEGER iflag_con

  ! original set from convect:

  ELSE

    ! constants consistent with LMDZ:

    ! ori      t0  = RTT

    ! maf       grav= 10.    ! implicitely or explicitely used in convect3

  END IF

  ! ginv=1.0/g

END SUBROUTINE cv_thermo