[1607] | 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 |
---|