source: LMDZ5/trunk/libf/phymar/PHY_Atm_CM_ALLOC.f90 @ 5450

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

Inclusion de la physique de MAR


Integration of MAR physics

File size: 6.8 KB
RevLine 
[2089]1
2      subroutine PHY_Atm_CM_ALLOC
3
4!------------------------------------------------------------------------------+
5!                                                         Sun  9-Jun-2013  MAR |
6!                                                                              |
7!     subroutine PHY_Atm_CM_ALLOC  allocates prognostic variables of           |
8!                Cloud Microphysical    Scheme used by MAR                     |
9!                                                                              |
10!     version 3.p.4.1 created by H. Gallee,               Tue 19-Mar-2013      |
11!           Last Modification by H. Gallee,               Sun  9-Jun-2013      |
12!                                                                              |
13!------------------------------------------------------------------------------+
14
15
16      use Mod_Real
17      use Mod_PHY____grd
18      use Mod_PHY_CM_kkl
19
20
21      IMPLICIT NONE
22
23
24
25! =================================
26! ALLOCATION Mod_PHY_CM_kkl - BEGIN
27! =================================
28
29      allocate  ( Ta__CM(kcolp,mzp ) )                       !  Air   Temperature                                          [K]
30      allocate  ( qvswCM(kcolp,mzpp) )                       !  Saturation Specific Humidity     (over liquid water)   [kg/kg]
31      allocate  ( qvsiCM(kcolp,mzpp) )                       !  Saturation Specific Humidity     (over ice)            [kg/kg]
32      allocate  ( qw__CM(kcolp,mzpp) )                       !  Cloud Droplets      Concentration                      [kg/kg]
33      allocate  ( qwd_CM(kcolp,mzp ) )                       !  Cloud Droplets      Concentration Variation            [kg/kg]
34      allocate  ( CCNwCM(kcolp,mzp ) )                       !  Cloud Droplets      Number                              [-/m3]
35      allocate  ( qi__CM(kcolp,mzpp) )                       !  Cloud Ice Particles Concentration                      [kg/kg]
36      allocate  ( qid_CM(kcolp,mzp ) )                       !  Cloud Ice Particles Concentration Variation            [kg/kg]
37      allocate  ( CCNiCM(kcolp,mzp ) )                       !  Cloud Ice Particles Number                              [-/m3]
38      allocate  ( CFraCM(kcolp,mzp ) )                       !  Cloud               Fraction                            [-]   
39      allocate  ( qs__CM(kcolp,mzpp) )                       !  Snow      Particles Concentration                      [kg/kg]
40! #qg allocate  ( qg__CM(kcolp,mzpp) )                       !  Graupels            Concentration                      [kg/kg]
41      allocate  ( qr__CM(kcolp,mzpp) )                       !  Rain  Drops         Concentration                      [kg/kg]
42      allocate  ( HLatCM(kcolp,mzp ) )                       !  Latent Heat Release                                     [W/m2]
43
44      allocate  ( uss_CM(kcolp) )                            !  Snow  Particles     Turbulent Surface Flux            [kg m/s]
45      allocate  ( Ice0CM(kcolp) )                            !  Ice C.Accumulation (time t-dt)                        [m w.e.]
46      allocate  ( ICE_CM(kcolp) )                            !  Ice C.Accumulation (time t   )                        [m w.e.]
47      allocate  ( Sno0CM(kcolp) )                            !  Snow  Accumulation (time t-dt, before snow erosion)   [m w.e.]
48      allocate  ( SnobCM(kcolp) )                            !  Snow  Accumulation (time t-dt, after  snow erosion)   [m w.e.]
49      allocate  ( SnowCM(kcolp) )                            !  Snow  Accumulation (time t   )                        [m w.e.]
50      allocate  ( Rai0CM(kcolp) )                            !  Rain  Accumulation (time t-dt)                        [m w.e.]
51      allocate  ( RainCM(kcolp) )                            !  Rain  Accumulation (time t   )                        [m w.e.]
52
53      allocate  ( dpktCM(kcolp,mzp ) )                       !  Reduced Potential Temperature TENDENCY                  [KX/s]
54      allocate  ( dqv_CM(kcolp,mzp ) )                       !  Specific          Humidity    TENDENCY               [kg/kg/s]
55      allocate  ( dqw_CM(kcolp,mzp ) )                       !  Cloud Droplets Concentration  TENDENCY               [kg/kg/s]
56      allocate  ( dqi_CM(kcolp,mzp ) )                       !  Cloud Crystals Concentration  TENDENCY               [kg/kg/s]
57      allocate  ( dqs_CM(kcolp,mzp ) )                       !  Snow Particles Concentration  TENDENCY               [kg/kg/s]
58      allocate  ( dqr_CM(kcolp,mzp ) )                       !  Rain Drops     Concentration  TENDENCY               [kg/kg/s]
59! #cw allocate  ( dCw_CM(kcolp,mzp ) )                       !  CCN            Concentration  TENDENCY                   [1/s]
60      allocate  ( dCi_CM(kcolp,mzp ) )                       !  CIN            Concentration  TENDENCY                   [1/s]
61      allocate  ( dCF_CM(kcolp,mzp ) )                       !  Cloud Fraction                TENDENCY                   [1/s]
62
63      allocate  ( wat0EW(kcolp) )                            !  Total Precipitable  Water  in the  Air Column         [m w.e.]
64      allocate  ( wat1EW(kcolp) )                            !  Total Precipitable  Water  in the  Air Column         [m w.e.]
65      allocate  ( wat2EW(kcolp) )                            !  Total Precipitable  Water  in the  Air Column         [m w.e.]
66      allocate  ( watfEW(kcolp) )                            !  Water Flux (Atm. --> Srf.) during 1 Time Step         [m w.e.]
67      allocate  ( enr0EW(kcolp) )                            !  Total Energy (Sens. +Lat.) in the  Air Column         [m w.e.]
68      allocate  ( enr1EW(kcolp) )                            !  Total Energy (Sens. +Lat.) in the  Air Column         [m w.e.]
69      allocate  ( enr2EW(kcolp) )                            !  Total Energy (Sens. +Lat.) in the  Air Column         [m w.e.]
70      allocate  ( mphyEW(kcolp) )                            !                                                        [m w.e.]
71
72
73!  Isotopes Proxies
74!  ~~~~~~~~~~~~~~~~
75      allocate  ( Hcd_CM(kcolp) )                            ! latent heat release                                   [mm w.e.]
76      allocate  ( Tcd_CM(kcolp) )                            ! latent heat release weighted Air Temperature                [K]
77      allocate  ( Zcd_CM(kcolp) )                            ! latent heat release weighted Altitude                       [m]
78      allocate  ( Hsb_CM(kcolp) )                            ! latent heat absorb.                                   [mm w.e.]
79      allocate  ( Tsb_CM(kcolp) )                            ! latent heat absorb. weighted Air Temperature                [K]
80      allocate  ( Zsb_CM(kcolp) )                            ! latent heat absorb. weighted Altitude                       [m]
81
82
83
84! =================================
85! ALLOCATION Mod_PHY_CM_kkl -   END
86! =================================
87
88
89
90      end subroutine PHY_Atm_CM_ALLOC
Note: See TracBrowser for help on using the repository browser.