! $Id: cv3_cine.F90 2416 2015-12-24 11:58:33Z jyg $

  ! **************************************************************

  ! *

  ! CV3_CINE                                                    *

  ! *

  ! *

  ! written by   :   Frederique Cheruy                          *

  ! vectorization:   Jean-Yves Grandpeix, 19/06/2003, 11.54.43  *

  ! modified by :                                               *

  ! **************************************************************

  IMPLICIT NONE

  include "YOMCST.h"

  include "cvthermo.h"

  include "cv3param.h"

  ! input:

  INTEGER ncum, nd, nloc

  INTEGER icb(nloc), inb(nloc)

  REAL pbase(nloc), plcl(nloc)

  REAL p(nloc, nd), ph(nloc, nd+1)

  REAL tv(nloc, nd), tvp(nloc, nd)

  ! output

  REAL cina(nloc), cinb(nloc), plfc(nloc)

  ! local variables

  INTEGER il, i, j, k

  REAL dpmax

  REAL deltap, dcin

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

  ! Initialization

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

  END DO

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

  ! Recompute buoyancies

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

      ! print*,'tvp tv=',tvp(il,k),tv(il,k)

    END DO

  END DO

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

  ! calcul de la flottabilite a LCL (Buoylcl)

  ! ifst = first P-level above lcl

  ! isublcl = highest P-level below lcl.

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

  END DO

    ELSE

    END IF

  END DO

    tvlcl(il) = tv(il, ifst(il)-1) + (tv(il,ifst(il))-tv(il,ifst(il)-1))*( &

  END DO

  END DO

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

  ! premiere couche contenant un  niveau de flotabilite positive

  ! et premiere couche contenant un  niveau de flotabilite negative

  ! au dessus du niveau de condensation

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

  END DO

        ELSE

        END IF

      END IF

    END DO

  END DO

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

  ! When there is no positive buoyancy level, set Plfc, Cina and Cinb

  ! to arbitrary extreme values.

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

    END IF

  END DO

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

  ! -- Two cases : BUOYlcl >= 0 and BUOYlcl < 0.

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

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

  ! -- 1.0 BUOYlcl >=0.

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

  dpmax = 50.

  END DO

  ! 1.1 No inhibition case

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

  ! If buoyancy is positive at LCL and stays positive over a large enough

  ! pressure interval (=DPMAX), inhibition is set to zero,

      END IF

    END IF

  END DO

  ! 1.2 Upper inhibition only case

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

  END DO

  ! 1.2.1 Recompute itop (=1st layer with positive buoyancy above ineg)

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

    END IF

  END DO

        END IF

      END IF

    END DO

  END DO

  ! If there is no layer with positive buoyancy above ineg, set Plfc,

  ! Cina and Cinb to arbitrary extreme values.

    END IF

  END DO

  END DO

      ! 1.2.2  Calcul de la pression du niveau de flot. nulle juste au-dessus

      ! de LCL

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

        ! In order to get P0, one may interpolate linearly buoyancies

        ! between P(ineg) and P(ineg-1).

        p0(il) = (buoy(il,ineg(il))*p(il,ineg(il)-1)-buoy(il,ineg(il)-1)*p(il,ineg(il)))/ &

      ELSE

        ! In order to get P0, one has to interpolate between P(ineg) and

        ! Plcl.

        p0(il) = (buoy(il,ineg(il))*plcl(il)-buoylcl(il)*p(il,ineg(il)))/ &

      END IF

    END IF

  END DO

  ! 1.2.3 Computation of PLFC

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

      plfc(il) = (buoy(il,itop(il))*p(il,itop(il)-1)-buoy(il,itop( &

    END IF

  END DO

  ! 1.2.4 Computation of CINA

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

  ! Upper part of CINA : integral from P(itop-1) to Plfc

    END IF

  END DO

  ! Middle part of CINA : integral from P(ineg) to P(itop-1)

        END IF

      END IF

    END DO

  END DO

  ! Lower part of CINA : integral from P0 to P(ineg)

    END IF

  END DO

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

  ! -- 2.0 BUOYlcl <0.

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

  END DO

  ! 2.0.1 Premiere  couche ou la flotabilite est negative au dessus du sol

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

  ! au cas ou elle existe  sinon ilow=1 (nk apres)

  ! on suppose que la parcelle part de la premiere couche

    END IF

  END DO

        END IF

      END IF

    END DO

  END DO

  ! 2.0.2  Calcul de la pression du niveau de flot. nulle sous le nuage

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

        p0(il) = (buoy(il,ilow(il))*p(il,ilow(il)-1)-buoy(il,ilow( &

      ELSE

      END IF

    END IF

  END DO

  ! 2.1. Computation of CINB

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

    lswitch2(il) = (isublcl(il)==1 .AND. ilow(il)==1) .OR. &

  END DO

  ! c      IF (    (isublcl .EQ. 1 .AND. ilow .EQ. 1)

  ! c     $    .OR.(isublcl .EQ. ilow-1)) THEN

  ! 2.1.1 First case : Plcl just above P0

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

    END IF

  END DO

  END DO

  ! c      ELSE

  ! 2.1.2 Second case : there is at least one P-level between P0 and Plcl

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

  ! Lower part of CINB : integral from P0 to P(ilow)

    END IF

  END DO

  ! Middle part of CINB : integral from P(ilow) to P(isublcl)

  ! c      DO k = ilow,isublcl-1

      END IF

    END DO

  END DO

  ! Upper part of CINB : integral from P(isublcl) to Plcl

      dcin = rd*(buoy(il,isublcl(il))+buoylcl(il))*deltap/ &

    END IF

  END DO

  ! c      ENDIF

  ! 2.2 Computation of CINA

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

  END DO

  ! 2.2.1 FIrst case : Plcl > P(itop-1)

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

  ! In order to get Plfc, one may interpolate linearly buoyancies

  ! between P(itop) and P(itop-1).

      plfc(il) = (buoy(il,itop(il))*p(il,itop(il)-1)-buoy(il,itop( &

    END IF

  END DO

  ! Upper part of CINA : integral from P(itop-1) to Plfc

    END IF

  END DO

  ! Middle part of CINA : integral from P(icb+1) to P(itop-1)

      END IF

    END DO

  END DO

  ! Lower part of CINA : integral from Plcl to P(icb+1)

          dcin = 0.5*rd*(buoy(il,icb(il))+buoy(il,icb(il)+1))*deltap/ &

        END IF

        dcin = rd*(buoylcl(il)+buoy(il,icb(il)))*deltap/ &

      ELSE

        dcin = rd*(buoylcl(il)+buoy(il,icb(il)+1))*deltap/ &

      END IF

    END IF

  END DO

  END DO

  ! c      ELSE

  ! 2.2.2 Second case : Plcl lies between P(itop-1) and P(itop);

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

  ! In order to get Plfc, one has to interpolate between P(itop) and Plcl.

      plfc(il) = (buoy(il,itop(il))*plcl(il)-buoylcl(il)*p(il,itop(il)))/ &

    END IF

  END DO

    END IF

  END DO

  ! c      ENDIF

END SUBROUTINE cv3_cine