source: LMDZ5/trunk/libf/phylmd/phys_output_var_mod.F90 @ 2339

Last change on this file since 2339 was 2333, checked in by lguez, 9 years ago

New parameterization of gravity wave drag due to front/jet systems, by

  1. de la Camara and F. Lott. The new Camara-Lott parameterization

replaces the Hines parameterization so it is activated if not ok_hines
and ok_gwd_rando.

Also changed distribution of phase speeds in FLOTT_GWD_rando, from
uniform to Gaussian. Bug fix in sugwd_strato. Bug fix in the arguments
of the call to add_phys_tend for methane oxydation.

For the new Camara-Lott parameterization, we need to compute relative
vorticity in calfis and pass it as a new argument "rot" to
physiq. Interpolation of relative vorticity to the physics grid is not
optimal for now: it is not weighted by cell areas.

Alvaro de la Camara, Fran\c{}cois Lott
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 4.2 KB
Line 
1!
2! phys_local_var_mod.F90 1327 2010-03-17 15:33:56Z idelkadi $
3
4MODULE phys_output_var_mod
5
6  use dimphy
7  ! Variables outputs pour les ecritures des sorties
8  !======================================================================
9  !
10  !
11  !======================================================================
12  ! Declaration des variables
13
14  REAL, SAVE, ALLOCATABLE :: snow_o(:), zfra_o(:)
15  !$OMP THREADPRIVATE(snow_o, zfra_o)
16  INTEGER, SAVE, ALLOCATABLE ::  itau_con(:)       ! Nombre de pas ou rflag <= 1
17  !$OMP THREADPRIVATE(itau_con)
18  REAL, ALLOCATABLE :: bils_ec(:) ! Contribution of energy conservation
19  REAL, ALLOCATABLE :: bils_ech(:) ! Contribution of energy conservation
20  REAL, ALLOCATABLE :: bils_tke(:) ! Contribution of energy conservation
21  REAL, ALLOCATABLE :: bils_diss(:) ! Contribution of energy conservation
22  REAL, ALLOCATABLE :: bils_kinetic(:) ! bilan de chaleur au sol, kinetic
23  REAL, ALLOCATABLE :: bils_enthalp(:) ! bilan de chaleur au sol
24  REAL, ALLOCATABLE :: bils_latent(:) ! bilan de chaleur au sol
25  !$OMP THREADPRIVATE(bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
26
27
28  ! ug Plein de variables venues de phys_output_mod
29  INTEGER, PARAMETER                           :: nfiles = 9
30  LOGICAL, DIMENSION(nfiles), SAVE             :: clef_files
31  LOGICAL, DIMENSION(nfiles), SAVE             :: clef_stations
32  INTEGER, DIMENSION(nfiles), SAVE             :: lev_files
33  INTEGER, DIMENSION(nfiles), SAVE             :: nid_files
34  INTEGER, DIMENSION(nfiles), SAVE  :: nnid_files
35  !$OMP THREADPRIVATE(clef_files, clef_stations, lev_files,nid_files,nnid_files)
36  INTEGER, DIMENSION(nfiles), SAVE :: nnhorim
37
38  INTEGER, DIMENSION(nfiles), SAVE :: nhorim, nvertm
39  INTEGER, DIMENSION(nfiles), SAVE :: nvertap, nvertbp, nvertAlt
40  REAL, DIMENSION(nfiles), SAVE                :: zoutm
41  CHARACTER(LEN=20), DIMENSION(nfiles), SAVE   :: type_ecri
42  !$OMP THREADPRIVATE(nnhorim, nhorim, nvertm, zoutm,type_ecri)
43  CHARACTER(LEN=20), DIMENSION(nfiles), SAVE  :: type_ecri_files, phys_out_filetypes
44  !$OMP THREADPRIVATE(type_ecri_files, phys_out_filetypes)
45  CHARACTER(LEN=20), DIMENSION(nfiles), SAVE  :: phys_out_filenames
46  !$OMP THREADPRIVATE(phys_out_filenames)
47
48  ! swaero_diag : flag indicates if it is necessary to do calculation for some aerosol diagnostics
49  LOGICAL, SAVE                                :: swaero_diag=.FALSE.
50  !$OMP THREADPRIVATE(swaero_diag)
51
52  INTEGER, SAVE:: levmin(nfiles) = 1
53  INTEGER, SAVE:: levmax(nfiles)
54  !$OMP THREADPRIVATE(levmin, levmax)
55
56  REAL, SAVE                :: zdtime_moy
57  !$OMP THREADPRIVATE(zdtime_moy)
58
59  LOGICAL, SAVE :: vars_defined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
60
61  REAL, allocatable:: zustr_gwd_hines(:), zvstr_gwd_hines(:) ! (klon)
62  REAL, allocatable:: zustr_gwd_front(:), zvstr_gwd_front(:) ! (klon)
63  REAL, allocatable:: zustr_gwd_rando(:), zvstr_gwd_rando(:) ! (klon)
64  !$OMP THREADPRIVATE(zustr_gwd_rando, zvstr_gwd_rando)
65
66  TYPE ctrl_out
67     INTEGER,DIMENSION(nfiles)            :: flag
68     CHARACTER(len=20)                    :: name
69     CHARACTER(len=150)                   :: description
70     CHARACTER(len=20)                    :: unit
71     CHARACTER(len=20),DIMENSION(nfiles)  :: type_ecrit
72  END TYPE ctrl_out
73
74CONTAINS
75
76  !======================================================================
77  SUBROUTINE phys_output_var_init
78    use dimphy
79
80    IMPLICIT NONE
81
82    include "clesphys.h"
83
84    !------------------------------------------------
85
86    allocate(snow_o(klon), zfra_o(klon))
87    allocate(itau_con(klon))
88    allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
89
90    IF (ok_hines) allocate(zustr_gwd_hines(klon), zvstr_gwd_hines(klon))
91    IF (.not.ok_hines.and.ok_gwd_rando) &
92                  allocate(zustr_gwd_front(klon), zvstr_gwd_front(klon))
93    IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
94
95  END SUBROUTINE phys_output_var_init
96
97  !======================================================================
98  SUBROUTINE phys_output_var_end
99    use dimphy
100    IMPLICIT NONE
101
102    deallocate(snow_o,zfra_o,itau_con)
103    deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
104
105  END SUBROUTINE phys_output_var_end
106
107END MODULE phys_output_var_mod
Note: See TracBrowser for help on using the repository browser.