source: LMDZ5/branches/testing/libf/phy1d/add_pbl_tend.F @ 1910

Last change on this file since 1910 was 1910, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1860:1909 into testing branch

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.6 KB
Line 
1      SUBROUTINE add_pbl_tend (zdu,zdv,zdt,zdq,zdql,text)
2!======================================================================
3! Ajoute les tendances de couche limite, soit determinees par la parametrisation
4! physique, soit forcees,  aux variables d etat de la dynamique t_seri, q_seri ...
5!======================================================================
6
7
8!======================================================================
9! Declarations
10!======================================================================
11
12      use dimphy
13      use phys_local_var_mod
14      use phys_state_var_mod
15      IMPLICIT none
16#include "dimensions.h"
17        real  hthturb_gcssold(llm)
18        real  hqturb_gcssold(llm)
19        real dtime_frcg
20        logical  Turb_fcg_gcssold
21        common /turb_forcing/ dtime_frcg,
22     $           Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold
23! Arguments :
24!------------
25      REAL zdu(klon,klev),zdv(klon,klev)
26      REAL zdt(klon,klev),zdq(klon,klev),zdql(klon,klev)
27      CHARACTER*(*) text
28
29! Local :
30!--------
31      REAL zzdt(klon,klev),zzdq(klon,klev)
32      INTEGER i,k
33!
34      IF (Turb_fcg_gcssold) THEN
35       DO k = 1,klev
36        DO i = 1,klon
37         zzdt(i,k) = hthturb_gcssold(k)*dtime_frcg
38         zzdq(i,k) = hqturb_gcssold(k)*dtime_frcg
39        ENDDO
40       ENDDO
41       print *,' add_pbl_tend, dtime_frcg ', dtime_frcg
42       print *,' add_pbl_tend, zzdt ', zzdt
43       print *,' add_pbl_tend, zzdq ', zzdq
44       CALL add_phys_tend (zdu,zdv,zzdt,zzdq,zdql,text)
45      ELSE
46       CALL add_phys_tend (zdu,zdv, zdt, zdq,zdql,text)
47      ENDIF
48!
49
50      RETURN
51      END
Note: See TracBrowser for help on using the repository browser.