source: LMDZ5/trunk/libf/phylmd/add_pbl_tend.F90 @ 2079

Last change on this file since 2079 was 1998, checked in by fhourdin, 11 years ago

Reactivation du calcul d'un zmax continu pour les thermiques

(thermcell_height.F90, thermcell_plume.F90, thermcell_dry.F90)
ouvre la voie à la réactivation d'une fermeture humide des thermiques
iflag_thermals_closure=2
(conf_phys_m.F90, thermcell.h, thermcell_main.F90)

Modification liée à la conservation de l'eau

(add_phys_tend.F90, add_pbl_tend.F90, physiq.F90)

Modifications liées au déclenchement stochastique

  1. possibilité de revenir à la Ale déterministe pour le criter ALE>|CIN| iflag_trig_bl=2, 1 par défaut)
  2. possibilité d'activer une fermeture statistique où ALP est divisé par la probabilité de déclenchement iflag_clos_bl=1 (0 par défaut, ancienne option 1 passée en =2)

Modification de l'entrainemement dans la version "stratocumulus" du

modèle du thermique (quand iflag_thermals_ed=8).
(modifie thermcell_plume.F90)

Catherine, Jean-Yves et Frédéric
File size: 1.6 KB
Line 
1SUBROUTINE add_pbl_tend(zdu, zdv, zdt, zdq, zdql, paprs, text)
2  ! ======================================================================
3  ! Ajoute les tendances de couche limite, soit determinees par la
4  ! parametrisation
5  ! physique, soit forcees,  aux variables d etat de la dynamique t_seri,
6  ! q_seri ...
7  ! ======================================================================
8
9
10  ! ======================================================================
11  ! Declarations
12  ! ======================================================================
13
14  USE dimphy
15  USE phys_local_var_mod
16  USE phys_state_var_mod
17  IMPLICIT NONE
18  include "dimensions.h"
19  REAL hthturb_gcssold(llm)
20  REAL hqturb_gcssold(llm)
21  REAL dtime_frcg
22  LOGICAL turb_fcg_gcssold
23  COMMON /turb_forcing/dtime_frcg, hthturb_gcssold, hqturb_gcssold, &
24    turb_fcg_gcssold
25
26  ! Arguments :
27  ! ------------
28  REAL zdu(klon, klev), zdv(klon, klev)
29  REAL zdt(klon, klev), zdq(klon, klev), zdql(klon, klev)
30  CHARACTER *(*) text
31  REAL paprs(klon,klev+1)
32
33  ! Local :
34  ! --------
35  REAL zzdt(klon, klev), zzdq(klon, klev)
36  INTEGER i, k
37
38  IF (turb_fcg_gcssold) THEN
39    DO k = 1, klev
40      DO i = 1, klon
41        zzdt(i, k) = hthturb_gcssold(k)*dtime_frcg
42        zzdq(i, k) = hqturb_gcssold(k)*dtime_frcg
43      END DO
44    END DO
45    PRINT *, ' add_pbl_tend, dtime_frcg ', dtime_frcg
46    PRINT *, ' add_pbl_tend, zzdt ', zzdt
47    PRINT *, ' add_pbl_tend, zzdq ', zzdq
48    CALL add_phys_tend(zdu, zdv, zzdt, zzdq, zdql, paprs, text)
49  ELSE
50    CALL add_phys_tend(zdu, zdv, zdt, zdq, zdql, paprs, text)
51  END IF
52
53
54  RETURN
55END SUBROUTINE add_pbl_tend
Note: See TracBrowser for help on using the repository browser.