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

Last change on this file since 2529 was 2529, checked in by oboucher, 9 years ago

Reviewing the way swaero_diag is computed to correct an old bug that
makes topswad and co diagnostics wrong for the first timestep...
And introducing swaero_diag in rrtm/recmwf_aero so that double and
triple radiation calls when diagnostics are not needed to save time

  • 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.5 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
[2529]49  !--OB: this needs to be set to TRUE by default and changed back to FALSE after first radiation call
50  !--    and corrected back to TRUE based on output requests
51  LOGICAL, SAVE                                :: swaero_diag=.TRUE.
[1938]52  !$OMP THREADPRIVATE(swaero_diag)
[1791]53
[1938]54  INTEGER, SAVE:: levmin(nfiles) = 1
55  INTEGER, SAVE:: levmax(nfiles)
56  !$OMP THREADPRIVATE(levmin, levmax)
[1791]57
[1938]58  REAL, SAVE                :: zdtime_moy
59  !$OMP THREADPRIVATE(zdtime_moy)
[1797]60
[1938]61  LOGICAL, SAVE :: vars_defined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
[1807]62
[2333]63  REAL, allocatable:: zustr_gwd_hines(:), zvstr_gwd_hines(:) ! (klon)
64  REAL, allocatable:: zustr_gwd_front(:), zvstr_gwd_front(:) ! (klon)
[1938]65  REAL, allocatable:: zustr_gwd_rando(:), zvstr_gwd_rando(:) ! (klon)
[2341]66  !$OMP THREADPRIVATE(zustr_gwd_hines, zvstr_gwd_hines)
67  !$OMP THREADPRIVATE(zustr_gwd_front, zvstr_gwd_front)
[1938]68  !$OMP THREADPRIVATE(zustr_gwd_rando, zvstr_gwd_rando)
69
[1791]70  TYPE ctrl_out
71     INTEGER,DIMENSION(nfiles)            :: flag
72     CHARACTER(len=20)                    :: name
73     CHARACTER(len=150)                   :: description
74     CHARACTER(len=20)                    :: unit
75     CHARACTER(len=20),DIMENSION(nfiles)  :: type_ecrit
76  END TYPE ctrl_out
[1938]77
[1335]78CONTAINS
79
[1938]80  !======================================================================
81  SUBROUTINE phys_output_var_init
82    use dimphy
[1335]83
[1938]84    IMPLICIT NONE
[1335]85
[1938]86    include "clesphys.h"
[1335]87
[1938]88    !------------------------------------------------
[1335]89
[1938]90    allocate(snow_o(klon), zfra_o(klon))
91    allocate(itau_con(klon))
[2042]92    allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
[1335]93
[2333]94    IF (ok_hines) allocate(zustr_gwd_hines(klon), zvstr_gwd_hines(klon))
95    IF (.not.ok_hines.and.ok_gwd_rando) &
96                  allocate(zustr_gwd_front(klon), zvstr_gwd_front(klon))
[1938]97    IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
[1335]98
[1938]99  END SUBROUTINE phys_output_var_init
[1335]100
[1938]101  !======================================================================
102  SUBROUTINE phys_output_var_end
103    use dimphy
104    IMPLICIT NONE
105
106    deallocate(snow_o,zfra_o,itau_con)
[2042]107    deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
[1938]108
109  END SUBROUTINE phys_output_var_end
110
[1335]111END MODULE phys_output_var_mod
Note: See TracBrowser for help on using the repository browser.