source: LMDZ5/trunk/libf/cosp/cosp_constants.F90 @ 2182

Last change on this file since 2182 was 2080, checked in by idelkadi, 10 years ago
  • Rajout des derictives OPENMP pour lire les fichiers namelists
  • Corrections dans la re-ecriture des sorties (traitement des valeurs indefinies)
  • Suppression des parametres non utilises des namelists
  • 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: 6.6 KB
Line 
1! (c) British Crown Copyright 2008, the Met Office.
2! All rights reserved.
3!
4! Redistribution and use in source and binary forms, with or without modification, are permitted
5! provided that the following conditions are met:
6!
7!     * Redistributions of source code must retain the above copyright notice, this list
8!       of conditions and the following disclaimer.
9!     * Redistributions in binary form must reproduce the above copyright notice, this list
10!       of conditions and the following disclaimer in the documentation and/or other materials
11!       provided with the distribution.
12!     * Neither the name of the Met Office nor the names of its contributors may be used
13!       to endorse or promote products derived from this software without specific prior written
14!       permission.
15!
16! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
17! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
18! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
19! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
20! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
21! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER
22! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
23! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
24
25!
26! History:
27! Jul 2007 - A. Bodas-Salcedo - Initial version
28! Jul 2008 - A. Bodas-Salcedo - Added definitions of ISCCP axes
29! Oct 2008 - H. Chepfer       - Added PARASOL_NREFL
30!
31!
32MODULE MOD_COSP_CONSTANTS
33    use netcdf, only: nf90_fill_real
34    IMPLICIT NONE
35   
36    ! Indices to address arrays of LS and CONV hydrometeors
37    integer,parameter :: I_LSCLIQ = 1
38    integer,parameter :: I_LSCICE = 2
39    integer,parameter :: I_LSRAIN = 3
40    integer,parameter :: I_LSSNOW = 4
41    integer,parameter :: I_CVCLIQ = 5
42    integer,parameter :: I_CVCICE = 6
43    integer,parameter :: I_CVRAIN = 7
44    integer,parameter :: I_CVSNOW = 8
45    integer,parameter :: I_LSGRPL = 9
46   
47    ! Missing value
48!!    real,parameter :: R_UNDEF = -1.0E30
49!     real,parameter :: R_UNDEF = 9.96921e+36
50      real,parameter :: R_UNDEF = nf90_fill_real
51
52    ! Number of possible output variables
53    integer,parameter :: N_OUT_LIST = 27
54    ! Value for forward model result from a level that is under the ground
55    real,parameter :: R_GROUND = -1.0E20
56
57    ! Stratiform and convective clouds in frac_out
58    integer, parameter :: I_LSC = 1, & ! Large-scale clouds
59                          I_CVC = 2    ! Convective clouds
60   
61    !--- Radar constants
62    ! CFAD constants
63    integer,parameter :: DBZE_BINS     =   15   ! Number of dBZe bins in histogram (cfad)
64    real,parameter    :: DBZE_MIN      = -100.0 ! Minimum value for radar reflectivity
65    real,parameter    :: DBZE_MAX      =   80.0 ! Maximum value for radar reflectivity
66    real,parameter    :: CFAD_ZE_MIN   =  -50.0 ! Lower value of the first CFAD Ze bin
67    real,parameter    :: CFAD_ZE_WIDTH =    5.0 ! Bin width (dBZe)
68
69   
70    !--- Lidar constants
71    ! CFAD constants
72    integer,parameter :: SR_BINS       =   15
73    integer,parameter :: DPOL_BINS     =   6
74    real,parameter    :: LIDAR_UNDEF   =   999.999
75    ! Other constants
76    integer,parameter :: LIDAR_NCAT    =   4
77    integer,parameter :: PARASOL_NREFL =   5 ! parasol
78    real,parameter,dimension(PARASOL_NREFL) :: PARASOL_SZA = (/0.0, 20.0, 40.0, 6.0, 80.0/)
79!    real,parameter,dimension(PARASOL_NREFL) :: PARASOL_SZA = (/1.0, 2.0, 3.0, 4.0, 5.0/)
80    real,parameter    :: DEFAULT_LIDAR_REFF = 30.0e-6 ! Default lidar effective radius
81   
82    !--- MISR constants
83    integer,parameter :: MISR_N_CTH = 16
84
85    !--- RTTOV constants
86    integer,parameter :: RTTOV_MAX_CHANNELS = 20
87   
88    ! ISCCP tau-Pc axes
89    real,parameter,dimension(7) :: ISCCP_TAU = (/0.15, 0.80, 2.45, 6.5, 16.2, 41.5, 50000.0/)
90    real,parameter,dimension(2,7) :: ISCCP_TAU_BNDS = reshape(source=(/0.0,0.3,0.3,1.30,1.30,3.6,3.6,9.4, &
91                                                      9.4,23.0,23.0,60.0,60.0,100000.0/), shape=(/2,7/))
92   
93!     real,parameter,dimension(7) :: ISCCP_PC = (/9000., 24500., 37500., 50000., 62000., 74000., 90000./)
94!     real,parameter,dimension(2,7) :: ISCCP_PC_BNDS = reshape(source=(/0.0,18000.0,18000.0,31000.0,31000.0, &
95!                                44000.0,44000.0,56000.0,56000.0,68000.0,68000.0,80000.0,80000.0,100000.0/), shape=(/2,7/))
96   
97    real,parameter,dimension(7) :: ISCCP_PC = (/90000., 74000., 62000., 50000., 37500., 24500., 9000./)
98    real,parameter,dimension(2,7) :: ISCCP_PC_BNDS = reshape(source=(/100000.0,80000.0,80000.0,68000.0,68000.0,56000.0 &
99                               ,56000.0,44000.0,44000.0,31000.0,31000.0,18000.0,18000.0,0.0/), shape=(/2,7/))
100   
101    real,parameter,dimension(MISR_N_CTH) :: MISR_CTH = (/ 0., 0.25, 0.75, 1.25, 1.75, 2.25, 2.75, 3.5, &
102                                            4.5, 6., 8., 10., 12., 14.5, 16., 18./)
103    real,parameter,dimension(2,MISR_N_CTH) :: MISR_CTH_BNDS = reshape(source=(/ &
104                                            -99.0,  0.0,       0.0,  0.5,       0.5,  1.0,      1.0,  1.5, &
105                                              1.5,  2.0,       2.0,  2.5,       2.5,  3.0,      3.0,  4.0, &
106                                              4.0,  5.0,       5.0,  7.0,       7.0,  9.0,      9.0, 11.0, &
107                                             11.0, 13.0,      13.0, 15.0,      15.0, 17.0,     17.0, 99.0/), &
108                                             shape=(/2,MISR_N_CTH/))
109           
110    !  Table hclass for quickbeam
111    integer,parameter :: N_HYDRO = 9
112    real :: HCLASS_TYPE(N_HYDRO),HCLASS_COL(N_HYDRO),HCLASS_PHASE(N_HYDRO), &
113            HCLASS_CP(N_HYDRO),HCLASS_DMIN(N_HYDRO),HCLASS_DMAX(N_HYDRO)
114    real :: HCLASS_APM(N_HYDRO),HCLASS_BPM(N_HYDRO),HCLASS_RHO(N_HYDRO), &
115            HCLASS_P1(N_HYDRO),HCLASS_P2(N_HYDRO),HCLASS_P3(N_HYDRO)
116    data HCLASS_TYPE/5,1,2,2,5,1,2,2,2/
117    data HCLASS_COL/1,2,3,4,5,6,7,8,9/
118    data HCLASS_PHASE/0,1,0,1,0,1,0,1,1/
119    data HCLASS_CP/0,0,1,1,0,0,1,1,1/
120    data HCLASS_DMIN/-1,-1,-1,-1,-1,-1,-1,-1,-1/
121    data HCLASS_DMAX/-1,-1,-1,-1,-1,-1,-1,-1,-1/
122    data HCLASS_APM/524,110.8,524,-1,524,110.8,524,-1,-1/
123    data HCLASS_BPM/3,2.91,3,-1,3,2.91,3,-1,-1/
124    data HCLASS_RHO/-1,-1,-1,100,-1,-1,-1,100,400/
125    data HCLASS_P1/-1,-1,8000000.,3000000.,-1,-1,8000000.,3000000.,4000000./
126    data HCLASS_P2/6,40,-1,-1,6,40,-1,-1,-1/
127    data HCLASS_P3/0.3,2,-1,-1,0.3,2,-1,-1,-1/
128
129   
130   
131END MODULE MOD_COSP_CONSTANTS
Note: See TracBrowser for help on using the repository browser.