source: LMDZ5/trunk/libf/phymar/PHY_Atm_AT_ALLOC.f90 @ 2222

Last change on this file since 2222 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics
File size: 5.0 KB
Line 
1
2      subroutine PHY_Atm_AT_ALLOC
3
4!------------------------------------------------------------------------------+
5!                                                         Sun 16-Jun-2013  MAR |
6!                                                                              |
7!     subroutine PHY_Atm_AT_ALLOC  allocates prognostic variables of           |
8!                Atmospheric Turbulence Scheme used by MAR                     |
9!                                                                              |
10!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
11!           Last Modification by H. Gallee,               Sun 16-Jun-2013      |
12!                                                                              |
13!------------------------------------------------------------------------------+
14
15
16      use Mod_Real
17      use Mod_PHY____grd
18      use Mod_PHY_AT_kkl
19
20
21      IMPLICIT NONE
22
23
24
25! =================================
26! ALLOCATION Mod_PHY_AT_kkl - BEGIN
27! =================================
28
29      allocate  ( var_AT(kcolp,mzpp) )                       ! Dummy            to Diffuse                                [x]
30      allocate  ( Ac0_AT(      mzp ) )                       ! Tridiagonal Matrix Coefficient A: Common Factor        [m2/s3]
31      allocate  ( Cc0_AT(      mzp ) )                       ! Tridiagonal Matrix Coefficient C: Common Factor        [m2/s3]
32      allocate  ( Kz0_AT(      mzp ) )                       ! Vertical  Turbulent Diffusion Coefficient (MINIMUM)     [m2/s]
33      allocate  ( Ac__AT(kcolp,mzp ) )                       ! Tridiagonal Matrix Coefficient A: Common Factor (t)     [s/m2]
34      allocate  ( Cc__AT(kcolp,mzp ) )                       ! Tridiagonal Matrix Coefficient C: Common Factor (t)     [s/m2]
35      allocate  ( Kz__AT(kcolp,mzp ) )                       ! Vertical  Turbulent Diffusion Coefficient               [m2/s]
36      allocate  ( Kzm_AT(kcolp,mzp ) )                       ! Vertical  Turbulent Diffusion Coefficient (Momentum)    [m2/s]
37      allocate  ( Kzh_AT(kcolp,mzp ) )                       ! Vertical  Turbulent Diffusion Coefficient (Scalars)     [m2/s]
38      allocate  ( Kzh0AT(kcolp,mzp ) )                       ! Vertical  Turbulent Diffusion Coefficient (Scalars)     [m2/s]
39      allocate  ( A___AT(kcolp,mzp ) )                       ! Tridiagonal Matrix Coefficient A                           [-]
40      allocate  ( B___AT(kcolp,mzp ) )                       ! Tridiagonal Matrix Coefficient B                           [-]
41      allocate  ( C___AT(kcolp,mzp ) )                       ! Tridiagonal Matrix Coefficient C                           [-]
42      allocate  ( D___AT(kcolp,mzp ) )                       ! Independant Term               D                           [x]
43      allocate  ( P___AT(      mzp ) )                       ! Auxiliary   Term               P                           [-]
44      allocate  ( Q___AT(      mzp ) )                       ! Auxiliary   Term               Q                           [-]
45      allocate  ( X___AT(      mzp ) )                       ! Auxiliary   Unknown            X                           [x]
46
47      allocate  ( LMO_AT(kcolp     ) )                       ! Monin-Obukhov     Length         (Grid Cell Average)       [m]
48      allocate  ( zi__AT(kcolp     ) )                       ! Inversion         Height         (Grid Cell Average)[m a.g.l.]
49      allocate  ( TKE_AT(kcolp,mzp ) )                       ! Turbulent Kinetic Energy                               [m2/s2]
50      allocate  ( eps_AT(kcolp,mzp ) )                       ! Turbulent Kinetic Energy Dissipation                   [m2/s3]
51      allocate  ( TrT_AT(kcolp,mzp ) )                       ! Turbulent Kinetic Energy Transport                     [m2/s3]
52
53      allocate  ( dua_AT(kcolp,mzp ) )                       ! Wind Speed (x-direc.) Tendency                          [m/s2]
54      allocate  ( dva_AT(kcolp,mzp ) )                       ! Wind Speed (y-direc.) Tendency                          [m/s2]
55      allocate  ( dpktAT(kcolp,mzp ) )                       ! Potential Temperature Tendency, divided by p0**(R/Cp)   [x/s]
56      allocate  ( dqv_AT(kcolp,mzp ) )                       ! Specific  Humidity    Tendency                      [kg/kg/s]
57      allocate  ( dqw_AT(kcolp,mzp ) )                       ! Cloud Droplets Concen.Tendency                      [kg/kg/s]
58      allocate  ( dqi_AT(kcolp,mzp ) )                       ! Cloud Crystals Concen.Tendency                      [kg/kg/s]
59      allocate  ( dCi_AT(kcolp,mzp ) )                       ! CCNi           Concen.Tendency                          [1/s]
60      allocate  ( dqs_AT(kcolp,mzp ) )                       ! Snow Particles Concen.Tendency                      [kg/kg/s]
61      allocate  ( dqr_AT(kcolp,mzp ) )                       ! Rain Drops     Concen.Tendency                      [kg/kg/s]
62
63! =================================
64! ALLOCATION Mod_PHY_AT_kkl -   END
65! =================================
66
67
68      end subroutine PHY_Atm_AT_ALLOC
Note: See TracBrowser for help on using the repository browser.