| 1 | SUBROUTINE add_pbl_tend(zdu, zdv, zdt, zdq, zdql, paprs, text) |
|---|
| 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 | ! ====================================================================== |
|---|
| 8 | |
|---|
| 9 | |
|---|
| 10 | ! ====================================================================== |
|---|
| 11 | ! Declarations |
|---|
| 12 | ! ====================================================================== |
|---|
| 13 | |
|---|
| 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 |
|---|
| 25 | |
|---|
| 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 |
|---|
| 31 | REAL paprs(klon,klev+1) |
|---|
| 32 | |
|---|
| 33 | ! Local : |
|---|
| 34 | ! -------- |
|---|
| 35 | REAL zzdt(klon, klev), zzdq(klon, klev) |
|---|
| 36 | INTEGER i, k |
|---|
| 37 | |
|---|
| 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 |
|---|
| 48 | CALL add_phys_tend(zdu, zdv, zzdt, zzdq, zdql, paprs, text) |
|---|
| 49 | ELSE |
|---|
| 50 | CALL add_phys_tend(zdu, zdv, zdt, zdq, zdql, paprs, text) |
|---|
| 51 | END IF |
|---|
| 52 | |
|---|
| 53 | |
|---|
| 54 | RETURN |
|---|
| 55 | END SUBROUTINE add_pbl_tend |
|---|
