| 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 |
|---|
