source: LMDZ5/trunk/libf/phylmd/add_pbl_tend.F @ 1985

Last change on this file since 1985 was 1973, checked in by fhourdin, 11 years ago

Le add_pbl_tend est passé de phy1d à phylmd et appelé dans phylmd/physiq.F90
si klon_glo=1
add_pbl_tend is moved from phy1d to phylmd and called in phylmd/physiq.F90
if klon_glo=1

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/
22     s  dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold
23
24! Arguments :
25!------------
26      REAL zdu(klon,klev),zdv(klon,klev)
27      REAL zdt(klon,klev),zdq(klon,klev),zdql(klon,klev)
28      CHARACTER*(*) text
29
30! Local :
31!--------
32      REAL zzdt(klon,klev),zzdq(klon,klev)
33      INTEGER i,k
34!
35      IF (Turb_fcg_gcssold) THEN
36       DO k = 1,klev
37        DO i = 1,klon
38         zzdt(i,k) = hthturb_gcssold(k)*dtime_frcg
39         zzdq(i,k) = hqturb_gcssold(k)*dtime_frcg
40        ENDDO
41       ENDDO
42       print *,' add_pbl_tend, dtime_frcg ', dtime_frcg
43       print *,' add_pbl_tend, zzdt ', zzdt
44       print *,' add_pbl_tend, zzdq ', zzdq
45       CALL add_phys_tend (zdu,zdv,zzdt,zzdq,zdql,text)
46      ELSE
47       CALL add_phys_tend (zdu,zdv, zdt, zdq,zdql,text)
48      ENDIF
49!
50
51      RETURN
52      END
Note: See TracBrowser for help on using the repository browser.