source: LMDZ5/trunk/libf/phylmd/add_pbl_tend.F90 @ 1992

Last change on this file since 1992 was 1992, checked in by lguez, 10 years ago

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.
File size: 1.5 KB
Line 
1SUBROUTINE add_pbl_tend(zdu, zdv, zdt, zdq, zdql, 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
32  ! Local :
33  ! --------
34  REAL zzdt(klon, klev), zzdq(klon, klev)
35  INTEGER i, k
36
37  IF (turb_fcg_gcssold) THEN
38    DO k = 1, klev
39      DO i = 1, klon
40        zzdt(i, k) = hthturb_gcssold(k)*dtime_frcg
41        zzdq(i, k) = hqturb_gcssold(k)*dtime_frcg
42      END DO
43    END DO
44    PRINT *, ' add_pbl_tend, dtime_frcg ', dtime_frcg
45    PRINT *, ' add_pbl_tend, zzdt ', zzdt
46    PRINT *, ' add_pbl_tend, zzdq ', zzdq
47    CALL add_phys_tend(zdu, zdv, zzdt, zzdq, zdql, text)
48  ELSE
49    CALL add_phys_tend(zdu, zdv, zdt, zdq, zdql, text)
50  END IF
51
52
53  RETURN
54END SUBROUTINE add_pbl_tend
Note: See TracBrowser for help on using the repository browser.