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 |
---|