MODULE lmdz_wake_ini

IMPLICIT NONE

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

  ! But : Decrire le comportement des poches froides apparaissant dans les

  ! grands systemes convectifs, et fournir l'energie disponible pour

  ! le declenchement de nouvelles colonnes convectives.

  ! State variables :

  ! deltatw    : temperature difference between wake and off-wake regions

  ! deltaqw    : specific humidity difference between wake and off-wake regions

  ! sigmaw     : fractional area covered by wakes.

  ! wdens      : number of wakes per unit area

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

  ! Declaration de variables

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

  ! Variables a fixer

!jyg<

!!  REAL, SAVE                                            :: stark, wdens_ref, coefgw, alpk

  INTEGER, SAVE, PROTECTED                                         :: prt_level

  REAL, SAVE, PROTECTED, DIMENSION(2)                              :: wdens_ref

  REAL, SAVE, PROTECTED                                            :: stark, coefgw, alpk, wk_pupper

!>jyg

  REAL, SAVE, PROTECTED                                            :: crep_upper, crep_sol

  !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, wk_pupper, crep_upper, crep_sol)

  REAL, SAVE, PROTECTED                                            :: tau_cv

  !$OMP THREADPRIVATE(tau_cv)

  REAL, SAVE, PROTECTED                                            :: rzero, aa0 ! minimal wake radius and area

  !$OMP THREADPRIVATE(rzero, aa0)

  LOGICAL, SAVE, PROTECTED                                         :: flag_wk_check_trgl

  !$OMP THREADPRIVATE(flag_wk_check_trgl)

  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_act

  !$OMP THREADPRIVATE(iflag_wk_act)

  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_check_trgl

  !$OMP THREADPRIVATE(iflag_wk_check_trgl)

  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_pop_dyn

  !$OMP THREADPRIVATE(iflag_wk_pop_dyn)

  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_profile

  !$OMP THREADPRIVATE(iflag_wk_profile)

  REAL, SAVE, PROTECTED                                            :: wdensmin

  !$OMP THREADPRIVATE(wdensmin)

  REAL, SAVE, PROTECTED                                            :: sigmad, hwmin, wapecut, cstart

  !$OMP THREADPRIVATE(sigmad, hwmin, wapecut, cstart)

  REAL, SAVE, PROTECTED                                            :: sigmaw_max

  !$OMP THREADPRIVATE(sigmaw_max)

  REAL, SAVE, PROTECTED                                            :: dens_rate

  !$OMP THREADPRIVATE(dens_rate)

  REAL, SAVE, PROTECTED                                            :: epsilon_loc

  !$OMP THREADPRIVATE(epsilon_loc)

  REAL, SAVE, PROTECTED                                            :: epsim1,RG,RD

  !$OMP THREADPRIVATE(epsim1,RG,RD)

CONTAINS

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

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

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

  ! *

  ! WAKE                                                        *

  ! retour a un Pupper fixe                                *

  ! *

  ! written by   :  GRANDPEIX Jean-Yves   09/03/2000            *

  ! modified by :   ROEHRIG Romain        01/29/2007            *

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

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

  ! Initialisations

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

  USE ioipsl_getin_p_mod, ONLY : getin_p

  real eps

  integer, intent(in) :: prt_lev

  real, intent(in) :: rg_in,rd_in,rv_in

  ! Essais d'initialisation avec sigmaw = 0.02 et hw = 10.

!!  DATA wdensmin/1.e-12/

  ! cc nrlmd

  ! cc

  ! Longueur de maille (en m)

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

  ! ALON = 3.e5

  ! alon = 1.E6

  ! Configuration de coefgw,stark,wdens (22/02/06 by YU Jingmei)

  ! coefgw : Coefficient pour les ondes de gravite

  ! stark : Coefficient k dans Cstar=k*sqrt(2*WAPE)

  ! wdens : Densite surfacique de poche froide

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

  ! cc nrlmd      coefgw=10

  ! coefgw=1

  ! wdens0 = 1.0/(alon**2)

  ! cc nrlmd      wdens = 1.0/(alon**2)

  ! cc nrlmd      stark = 0.50

  ! CRtest

  ! cc nrlmd      alpk=0.1

  ! alpk = 1.0

  ! alpk = 0.5

  ! alpk = 0.05

  ! Get wapecut from parameter file

  wapecut = 1.

  ! cc nrlmd Lecture du fichier wake_param.data

  ! cc nrlmd Lecture du fichier wake_param.data

!jyg<

!!  wdens_ref=8.E-12

!!  CALL getin_p('wdens_ref',wdens_ref)

!>jyg

!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!  Population dynamics parameters    !!!!!!!!!!!!!!!!!!!!!!!!!!!!

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

                                                    ! and wdens prognostic

                                            ! 1: act(:)=1.

                                            ! 2: act(:)=f(Wape)

                                                    ! and wdens prognostic

!

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

!jyg<

!!  WRITE(*,*) 'wdens_ref=', wdens_ref

!>jyg

END SUBROUTINE wake_ini

END MODULE lmdz_wake_ini