SUBROUTINE alpale ( debut, itap, dtime, paprs, omega, t_seri,   &
                    alp_offset, it_wape_prescr,  wape_prescr, fip_prescr, &
                    ale_bl_prescr, alp_bl_prescr, &
                    wake_pe, wake_fip,  &
                    Ale_bl, Ale_bl_trig, Alp_bl, &
                    Ale, Alp, Ale_wake, Alp_wake )

! **************************************************************
! *
! ALPALE                                                       *
! *
! *
! written by   : Jean-Yves Grandpeix, 12/05/2016              *
! modified by :                                               *
! **************************************************************

  USE dimphy
  USE ioipsl_getin_p_mod, ONLY: getin_p
  USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level
  USE phys_local_var_mod, ONLY: zw2       ! Variables internes non sauvegardees de la physique

  IMPLICIT NONE

!================================================================
! Auteur(s)   : Jean-Yves Grandpeix, 12/05/2016
! Objet : Sums up all contributions to Ale and Alp
!================================================================

! Input arguments
!----------------
  LOGICAL, INTENT(IN)                                        :: debut
  INTEGER, INTENT(IN)                                        :: itap
  REAL, INTENT(IN)                                           :: dtime
  INTEGER, INTENT(IN)                                        :: it_wape_prescr
  REAL, INTENT(IN)                                           :: wape_prescr, fip_prescr
  REAL, INTENT(IN)                                           :: Ale_bl_prescr, Alp_bl_prescr
  REAL, INTENT(IN)                                           :: alp_offset
  REAL, DIMENSION(klon,klev+1), INTENT(IN)                   :: paprs
  REAL, DIMENSION(klon,klev), INTENT(IN)                     :: t_seri
  REAL, DIMENSION(klon,klev), INTENT(IN)                     :: omega
  REAL, DIMENSION(klon), INTENT(IN)                          :: wake_pe, wake_fip
  REAL, DIMENSION(klon), INTENT(IN)                          :: Ale_bl, Ale_bl_trig, Alp_bl


! Output arguments
!----------------
  REAL, DIMENSION(klon), INTENT(OUT)                         :: Ale, Alp
  REAL, DIMENSION(klon), INTENT(OUT)                         :: Ale_wake, Alp_wake

  include "alpale.h"
  include "YOMCST.h"
  include "YOETHF.h"

! Local variables
!----------------
  INTEGER                                                    :: i, k
  REAL, DIMENSION(klon)                                      :: www
  REAL, SAVE                                                 :: ale_max=1000.
  REAL, SAVE                                                 :: alp_max=2.
  CHARACTER*20 modname
  CHARACTER*80 abort_message


    !$OMP THREADPRIVATE(ale_max,alp_max)

       ! Calcul de l'energie disponible ALE (J/kg) et de la puissance
       ! disponible ALP (W/m2) pour le soulevement des particules dans
       ! le modele convectif

       do i = 1,klon
          ALE(i) = 0.
          ALP(i) = 0.
       enddo

       !calcul de ale_wake et alp_wake
       if (iflag_wake>=1) then
          if (itap <= it_wape_prescr) then
             do i = 1,klon
                ale_wake(i) = wape_prescr
                alp_wake(i) = fip_prescr
             enddo
          else
             do i = 1,klon
                !jyg  ALE=WAPE au lieu de ALE = 1/2 Cstar**2
                !cc           ale_wake(i) = 0.5*wake_cstar(i)**2
                ale_wake(i) = wake_pe(i)
                alp_wake(i) = wake_fip(i)
             enddo
          endif
       else
          do i = 1,klon
             ale_wake(i) = 0.
             alp_wake(i) = 0.
          enddo
       endif
       !combinaison avec ale et alp de couche limite: constantes si pas
       !de couplage, valeurs calculees dans le thermique sinon
       if (iflag_coupl==0) then
          if (debut.and.prt_level>9) &
               WRITE(lunout,*)'ALE et ALP imposes'
          do i = 1,klon
             !on ne couple que ale
             !           ALE(i) = max(ale_wake(i),Ale_bl(i))
             ALE(i) = max(ale_wake(i),ale_bl_prescr)
             !on ne couple que alp
             !           ALP(i) = alp_wake(i) + Alp_bl(i)
             ALP(i) = alp_wake(i) + alp_bl_prescr
          enddo
       else
          IF(prt_level>9)WRITE(lunout,*)'ALE et ALP couples au thermique'
          !         do i = 1,klon
          !             ALE(i) = max(ale_wake(i),Ale_bl(i))
          ! avant        ALP(i) = alp_wake(i) + Alp_bl(i)
          !             ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
          !         write(20,*)'ALE',ALE(i),Ale_bl(i),ale_wake(i)
          !         write(21,*)'ALP',ALP(i),Alp_bl(i),alp_wake(i)
          !         enddo

          ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
          ! Modif FH 2010/04/27. Sans doute temporaire.
          ! Deux options pour le alp_offset : constant si >?? 0 ou
          ! proportionnel ??a w si <0
          ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
          ! Estimation d'une vitesse verticale effective pour ALP
          if (1==0) THEN
             www(1:klon)=0.
             do k=2,klev-1
                do i=1,klon
                   www(i)=max(www(i),-omega(i,k)*RD*t_seri(i,k) &
                        /(RG*paprs(i,k)) *zw2(i,k)*zw2(i,k))
                   ! if (paprs(i,k)>pbase(i)) then
                   ! calcul approche de la vitesse verticale en m/s
                   !  www(i)=max(www(i),-omega(i,k)*RD*temp(i,k)/(RG*paprs(i,k))
                   !             endif
                   !   Le 0.1 est en gros H / ps = 1e4 / 1e5
                enddo
             enddo
             do i=1,klon
                if (www(i)>0. .and. ale_bl(i)>0. ) www(i)=www(i)/ale_bl(i)
             enddo
          ENDIF


          do i = 1,klon
             ALE(i) = max(ale_wake(i),Ale_bl(i))
             !cc nrlmd le 10/04/2012----------Stochastic triggering------------
             if (iflag_trig_bl>=1) then
                ALE(i) = max(ale_wake(i),Ale_bl_trig(i))
             endif
             !cc fin nrlmd le 10/04/2012
             if (alp_offset>=0.) then
                ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
             else
                abort_message ='Ne pas passer la car www non calcule'
                CALL abort_physic (modname,abort_message,1)

                ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                !                                _                  _
                ! Ajout d'une composante 3 * A * w w'2 a w'3 avec
                ! w=www : w max sous pbase ou A est la fraction
                ! couverte par les ascendances w' on utilise le fait
                ! que A * w'3 = ALP et donc A * w'2 ~ ALP / sqrt(ALE)
                ! (on ajoute 0.1 pour les singularites)
                ALP(i)=alp_wake(i)*(1.+3.*www(i)/( sqrt(ale_wake(i))+0.1) ) &
                     +alp_bl(i)  *(1.+3.*www(i)/( sqrt(ale_bl(i))  +0.1) )
                !    ALP(i)=alp_wake(i)+Alp_bl(i)+alp_offset*min(omega(i,6),0.)
                !             if (alp(i)<0.) then
                !                PRINT*,'ALP ',alp(i),alp_wake(i) &
                !                     ,Alp_bl(i),alp_offset*min(omega(i,6),0.)
                !             endif
             endif
          enddo
          ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

       endif
       do i=1,klon
          if (alp(i)>alp_max) then
             IF(prt_level>9)WRITE(lunout,*)                             &
                  'WARNING SUPER ALP (seuil=',alp_max, &
                  '): i, alp, alp_wake,ale',i,alp(i),alp_wake(i),ale(i)
             alp(i)=alp_max
          endif
          if (ale(i)>ale_max) then
             IF(prt_level>9)WRITE(lunout,*)                             &
                  'WARNING SUPER ALE (seuil=',ale_max, &
                  '): i, alp, alp_wake,ale',i,ale(i),ale_wake(i),alp(i)
             ale(i)=ale_max
          endif
       enddo

       !fin calcul ale et alp
       !=======================================================================


  RETURN
  END

