source: LMDZ5/branches/LF-private/libf/phy1d/add_pbl_tend.F @ 5440

Last change on this file since 5440 was 1607, checked in by lguez, 13 years ago

Import 1D files. Files added to directory "phy1d" were in directories:

lmdz1d_source_20111207/phy1d_source
lmdz1d_source_20111207/phy1d_source_upd

extracted from:

http://www.lmd.jussieu.fr/~jyg/lmdz1d_source_20111207.tar.gz

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.