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
RevLine 
[1335]1!
2! phys_local_var_mod.F90 1327 2010-03-17 15:33:56Z idelkadi $
3
[1938]4MODULE phys_output_var_mod
[1335]5
[1938]6  use dimphy
7  ! Variables outputs pour les ecritures des sorties
8  !======================================================================
9  !
10  !
11  !======================================================================
12  ! Declaration des variables
[1335]13
[1938]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
[2042]19  REAL, ALLOCATABLE :: bils_ech(:) ! Contribution of energy conservation
[1938]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
[2042]25  !$OMP THREADPRIVATE(bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
[1753]26
[1865]27
[1938]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
[1761]37
[1938]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)
[1791]47
[1938]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)
[1791]51
[1938]52  INTEGER, SAVE:: levmin(nfiles) = 1
53  INTEGER, SAVE:: levmax(nfiles)
54  !$OMP THREADPRIVATE(levmin, levmax)
[1791]55
[1938]56  REAL, SAVE                :: zdtime_moy
57  !$OMP THREADPRIVATE(zdtime_moy)
[1797]58
[1938]59  LOGICAL, SAVE :: vars_defined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
[1807]60
[2333]61  REAL, allocatable:: zustr_gwd_hines(:), zvstr_gwd_hines(:) ! (klon)
62  REAL, allocatable:: zustr_gwd_front(:), zvstr_gwd_front(:) ! (klon)
[1938]63  REAL, allocatable:: zustr_gwd_rando(:), zvstr_gwd_rando(:) ! (klon)
64  !$OMP THREADPRIVATE(zustr_gwd_rando, zvstr_gwd_rando)
65
[1791]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
[1938]73
[1335]74CONTAINS
75
[1938]76  !======================================================================
77  SUBROUTINE phys_output_var_init
78    use dimphy
[1335]79
[1938]80    IMPLICIT NONE
[1335]81
[1938]82    include "clesphys.h"
[1335]83
[1938]84    !------------------------------------------------
[1335]85
[1938]86    allocate(snow_o(klon), zfra_o(klon))
87    allocate(itau_con(klon))
[2042]88    allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
[1335]89
[2333]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))
[1938]93    IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
[1335]94
[1938]95  END SUBROUTINE phys_output_var_init
[1335]96
[1938]97  !======================================================================
98  SUBROUTINE phys_output_var_end
99    use dimphy
100    IMPLICIT NONE
101
102    deallocate(snow_o,zfra_o,itau_con)
[2042]103    deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
[1938]104
105  END SUBROUTINE phys_output_var_end
106
[1335]107END MODULE phys_output_var_mod
Note: See TracBrowser for help on using the repository browser.