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