[1998] | 1 | SUBROUTINE add_pbl_tend(zdu, zdv, zdt, zdq, zdql, paprs, text) |
---|
[1992] | 2 | ! ====================================================================== |
---|
| 3 | ! Ajoute les tendances de couche limite, soit determinees par la |
---|
| 4 | ! parametrisation |
---|
| 5 | ! physique, soit forcees, aux variables d etat de la dynamique t_seri, |
---|
| 6 | ! q_seri ... |
---|
| 7 | ! ====================================================================== |
---|
[1973] | 8 | |
---|
| 9 | |
---|
[1992] | 10 | ! ====================================================================== |
---|
| 11 | ! Declarations |
---|
| 12 | ! ====================================================================== |
---|
[1973] | 13 | |
---|
[1992] | 14 | USE dimphy |
---|
| 15 | USE phys_local_var_mod |
---|
| 16 | USE phys_state_var_mod |
---|
| 17 | IMPLICIT NONE |
---|
| 18 | include "dimensions.h" |
---|
| 19 | REAL hthturb_gcssold(llm) |
---|
| 20 | REAL hqturb_gcssold(llm) |
---|
| 21 | REAL dtime_frcg |
---|
| 22 | LOGICAL turb_fcg_gcssold |
---|
| 23 | COMMON /turb_forcing/dtime_frcg, hthturb_gcssold, hqturb_gcssold, & |
---|
| 24 | turb_fcg_gcssold |
---|
[1973] | 25 | |
---|
[1992] | 26 | ! Arguments : |
---|
| 27 | ! ------------ |
---|
| 28 | REAL zdu(klon, klev), zdv(klon, klev) |
---|
| 29 | REAL zdt(klon, klev), zdq(klon, klev), zdql(klon, klev) |
---|
| 30 | CHARACTER *(*) text |
---|
[1998] | 31 | REAL paprs(klon,klev+1) |
---|
[1973] | 32 | |
---|
[1992] | 33 | ! Local : |
---|
| 34 | ! -------- |
---|
| 35 | REAL zzdt(klon, klev), zzdq(klon, klev) |
---|
| 36 | INTEGER i, k |
---|
[1973] | 37 | |
---|
[1992] | 38 | IF (turb_fcg_gcssold) THEN |
---|
| 39 | DO k = 1, klev |
---|
| 40 | DO i = 1, klon |
---|
| 41 | zzdt(i, k) = hthturb_gcssold(k)*dtime_frcg |
---|
| 42 | zzdq(i, k) = hqturb_gcssold(k)*dtime_frcg |
---|
| 43 | END DO |
---|
| 44 | END DO |
---|
| 45 | PRINT *, ' add_pbl_tend, dtime_frcg ', dtime_frcg |
---|
| 46 | PRINT *, ' add_pbl_tend, zzdt ', zzdt |
---|
| 47 | PRINT *, ' add_pbl_tend, zzdq ', zzdq |
---|
[1998] | 48 | CALL add_phys_tend(zdu, zdv, zzdt, zzdq, zdql, paprs, text) |
---|
[1992] | 49 | ELSE |
---|
[1998] | 50 | CALL add_phys_tend(zdu, zdv, zdt, zdq, zdql, paprs, text) |
---|
[1992] | 51 | END IF |
---|
| 52 | |
---|
| 53 | |
---|
| 54 | RETURN |
---|
| 55 | END SUBROUTINE add_pbl_tend |
---|