source: LMDZ6/trunk/libf/phylmd/cospv2/cosp_config.F90 @ 3958

Last change on this file since 3958 was 3723, checked in by idelkadi, 4 years ago

Debugging COSP v2 for simulators Calipso, Parasol, Cloudsat

File size: 25.4 KB
Line 
1! %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2! Copyright (c) 2015, Regents of the University of Colorado
3! All rights reserved.
4!
5! Redistribution and use in source and binary forms, with or without modification, are
6! permitted provided that the following conditions are met:
7!
8! 1. Redistributions of source code must retain the above copyright notice, this list of
9!    conditions and the following disclaimer.
10!
11! 2. Redistributions in binary form must reproduce the above copyright notice, this list
12!    of conditions and the following disclaimer in the documentation and/or other
13!    materials provided with the distribution.
14!
15! 3. Neither the name of the copyright holder nor the names of its contributors may be
16!    used to endorse or promote products derived from this software without specific prior
17!    written permission.
18!
19! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
20! EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
21! MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL
22! THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
23! SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
24! OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
25! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
26! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
27! OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
28!
29! History:
30! Jul 2007 - A. Bodas-Salcedo - Initial version
31! Jul 2008 - A. Bodas-Salcedo - Added definitions of ISCCP axes
32! Oct 2008 - H. Chepfer       - Added PARASOL_NREFL
33! Jun 2010 - R. Marchand      - Modified to support quickbeam V3, added ifdef for 
34!                               hydrometeor definitions
35! May 2015 - D. Swales        - Tidied up. Set up appropriate fields during initialization.
36! June 2015- D. Swales        - Moved hydrometeor class variables to hydro_class_init in
37!                               the module quickbeam_optics.
38! Mar 2016 - D. Swales        - Added scops_ccfrac. Was previously hardcoded in prec_scops.f90. 
39! Mar 2018 - R. Guzman        - Added LIDAR_NTYPE for the OPAQ diagnostics
40! Apr 2018 - R. Guzman        - Added parameters for GROUND LIDAR and ATLID simulators
41!
42! %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
43
44MODULE MOD_COSP_CONFIG
45    USE COSP_KINDS, ONLY: wp,dp
46    IMPLICIT NONE
47
48   ! #####################################################################################
49   ! Common COSP information
50   ! #####################################################################################
51    character(len=32) ::   &
52         COSP_VERSION              ! COSP Version ID (set in cosp_interface_init)
53    real(wp),parameter ::  &
54         R_UNDEF      = -1.0E30, & ! Missing value
55         R_GROUND     = -1.0E20, & ! Flag for below ground results
56         scops_ccfrac = 0.05       ! Fraction of column (or subcolumn) covered with convective
57                                   ! precipitation (default is 5%). *NOTE* This quantity may vary
58                                   ! between modeling centers.
59    logical :: &
60         use_vgrid                 ! True=Use new grid for L3 CLOUDAT and CALIPSO
61    integer,parameter ::   &
62         SR_BINS = 15,           & ! Number of bins (backscattering coefficient) in CALOPSO LIDAR simulator.
63         N_HYDRO = 9               ! Number of hydrometeor classes used by quickbeam radar simulator.
64
65    ! Mode debug ou prod (AI 0302018)
66    logical :: ok_debug_cosp = .true.
67    ! #################################################################################### 
68    ! Joint histogram bin-boundaries
69    ! tau is used by ISCCP and MISR
70    ! pres is used by ISCCP
71    ! hgt is used by MISR
72    ! ReffLiq is used by MODIS
73    ! ReffIce is used by MODIS
74    ! *NOTE* ALL JOINT-HISTOGRAM BIN BOUNDARIES ARE DECLARED AND DEFINED HERE IN
75    !        COSP_CONFIG, WITH THE EXCEPTION OF THE TAU AXIS USED BY THE MODIS SIMULATOR,
76    !        WHICH IS SET DURING INITIALIZATION IN COSP_INTERFACE_INIT.
77    ! ####################################################################################
78    ! Optical depth bin axis
79    integer,parameter :: &
80         ntau=7 
81    real(wp),parameter,dimension(ntau+1) :: &
82       tau_binBounds = (/0.0, 0.3, 1.3, 3.6, 9.4, 23., 60., 10000./)
83    real(wp),parameter,dimension(ntau) :: &
84         tau_binCenters = (/0.15, 0.80, 2.45, 6.5, 16.2, 41.5, 100.0/)
85    real(wp),parameter,dimension(2,ntau) :: &
86         tau_binEdges = reshape(source=(/0.0, 0.3,  0.3,  1.3,  1.3,  3.6,      3.6,     &
87                                         9.4, 9.4, 23.0, 23.0, 60.0, 60.0, 100000.0/),   &
88                                         shape=(/2,ntau/))
89
90    ! Optical depth bin axes (ONLY USED BY MODIS SIMULATOR IN v1.4)
91    integer :: l,k
92    integer,parameter :: &
93         ntauV1p4 = 6
94    real(wp),parameter,dimension(ntauV1p4+1) :: &
95         tau_binBoundsV1p4 = (/0.3, 1.3, 3.6, 9.4, 23., 60., 10000./)
96    real(wp),parameter,dimension(2,ntauV1p4) :: &
97         tau_binEdgesV1p4 = reshape(source =(/tau_binBoundsV1p4(1),((tau_binBoundsV1p4(k),l=1,2),   &
98                                             k=2,ntauV1p4),100000._wp/),shape = (/2,ntauV1p4/))
99    real(wp),parameter,dimension(ntauV1p4) :: &
100         tau_binCentersV1p4 = (tau_binEdgesV1p4(1,:)+tau_binEdgesV1p4(2,:))/2._wp 
101   
102    ! Cloud-top height pressure bin axis
103    integer,parameter :: &
104         npres = 7     
105    real(wp),parameter,dimension(npres+1) :: &
106         pres_binBounds = (/0., 180., 310., 440., 560., 680., 800., 10000./)
107    real(wp),parameter,dimension(npres) :: &
108         pres_binCenters = (/90000., 74000., 62000., 50000., 37500., 24500., 9000./)   
109    real(wp),parameter,dimension(2,npres) :: &
110         pres_binEdges = reshape(source=(/100000.0, 80000.0, 80000.0, 68000.0, 68000.0,    &
111                                           56000.0, 56000.0, 44000.0, 44000.0, 31000.0,    &
112                                           31000.0, 18000.0, 18000.0,     0.0/),           &
113                                           shape=(/2,npres/))
114
115    ! Cloud-top height bin axis #1
116    integer,parameter :: &
117         nhgt = 16
118    real(wp),parameter,dimension(nhgt+1) :: &
119         hgt_binBounds = (/-.99,0.,0.5,1.,1.5,2.,2.5,3.,4.,5.,7.,9.,11.,13.,15.,17.,99./)
120    real(wp),parameter,dimension(nhgt) :: &
121         hgt_binCenters = 1000*(/0.,0.25,0.75,1.25,1.75,2.25,2.75,3.5,4.5,6.,8.,10.,12.,   &
122         14.5,16.,18./) 
123    real(wp),parameter,dimension(2,nhgt) :: &
124         hgt_binEdges = 1000.0*reshape(source=(/-99.0, 0.0, 0.0, 0.5, 0.5, 1.0, 1.0, 1.5,  &
125                                                  1.5, 2.0, 2.0, 2.5, 2.5, 3.0, 3.0, 4.0,  &
126                                                  4.0, 5.0, 5.0, 7.0, 7.0, 9.0, 9.0,11.0,  &
127                                                  11.0,13.0,13.0,15.0,15.0,17.0,17.0,99.0/),&
128                                                  shape=(/2,nhgt/))   
129
130    ! Liquid and Ice particle bins for MODIS joint histogram of optical-depth and particle
131    ! size
132    integer :: i,j
133    integer,parameter :: &
134         nReffLiq = 6, & ! Number of bins for tau/ReffLiq joint-histogram
135         nReffIce = 6    ! Number of bins for tau/ReffICE joint-histogram
136    real(wp),parameter,dimension(nReffLiq+1) :: &
137         reffLIQ_binBounds = (/0., 8e-6, 1.0e-5, 1.3e-5, 1.5e-5, 2.0e-5, 3.0e-5/)
138    real(wp),parameter,dimension(nReffIce+1) :: &
139         reffICE_binBounds = (/0., 1.0e-5, 2.0e-5, 3.0e-5, 4.0e-5, 6.0e-5, 9.0e-5/)
140    real(wp),parameter,dimension(2,nReffICE) :: &
141         reffICE_binEdges = reshape(source=(/reffICE_binBounds(1),((reffICE_binBounds(k),  &
142                                    l=1,2),k=2,nReffICE),reffICE_binBounds(nReffICE+1)/),  &
143                                    shape = (/2,nReffICE/))
144    real(wp),parameter,dimension(2,nReffLIQ) :: &
145         reffLIQ_binEdges = reshape(source=(/reffLIQ_binBounds(1),((reffLIQ_binBounds(k),  &
146                                    l=1,2),k=2,nReffLIQ),reffLIQ_binBounds(nReffICE+1)/),  &
147                                    shape = (/2,nReffLIQ/))             
148    real(wp),parameter,dimension(nReffICE) :: &
149         reffICE_binCenters = (reffICE_binEdges(1,:)+reffICE_binEdges(2,:))/2._wp
150    real(wp),parameter,dimension(nReffLIQ) :: &
151         reffLIQ_binCenters = (reffLIQ_binEdges(1,:)+reffLIQ_binEdges(2,:))/2._wp
152
153    ! #################################################################################### 
154    ! Constants used by RTTOV.
155    ! #################################################################################### 
156    integer,parameter :: &
157         RTTOV_MAX_CHANNELS = 20
158    character(len=256),parameter :: &
159         rttovDir = '/homedata/rguzman/CALIPSO/RTTOV/rttov_11.3/'
160    ! #################################################################################### 
161    ! Constants used by the PARASOL simulator.
162    ! #################################################################################### 
163    integer,parameter :: &
164         PARASOL_NREFL = 5,  & ! Number of angles in LUT
165         PARASOL_NTAU  = 7     ! Number of optical depths in LUT
166    real(wp),parameter,dimension(PARASOL_NREFL) :: &
167         PARASOL_SZA = (/0.0, 20.0, 40.0, 60.0, 80.0/)
168    REAL(WP),parameter,dimension(PARASOL_NTAU) :: &
169         PARASOL_TAU = (/0., 1., 5., 10., 20., 50., 100./)
170   
171    ! LUTs
172    REAL(WP),parameter,dimension(PARASOL_NREFL,PARASOL_NTAU) :: &
173         ! LUT for liquid particles
174         rlumA = reshape(source=(/ 0.03,     0.03,     0.03,     0.03,     0.03,         &
175                                   0.090886, 0.072185, 0.058410, 0.052498, 0.034730,     &
176                                   0.283965, 0.252596, 0.224707, 0.175844, 0.064488,     &
177                                   0.480587, 0.436401, 0.367451, 0.252916, 0.081667,     &
178                                   0.695235, 0.631352, 0.509180, 0.326551, 0.098215,     &
179                                   0.908229, 0.823924, 0.648152, 0.398581, 0.114411,     &
180                                   1.0,      0.909013, 0.709554, 0.430405, 0.121567/),   &
181                                   shape=(/PARASOL_NREFL,PARASOL_NTAU/)),                &
182         ! LUT for ice particles                                     
183         rlumB = reshape(source=(/ 0.03,     0.03,     0.03,     0.03,     0.03,         &
184                                   0.092170, 0.087082, 0.083325, 0.084935, 0.054157,     &
185                                   0.311941, 0.304293, 0.285193, 0.233450, 0.089911,     &
186                                   0.511298, 0.490879, 0.430266, 0.312280, 0.107854,     &
187                                   0.712079, 0.673565, 0.563747, 0.382376, 0.124127,     &
188                                   0.898243, 0.842026, 0.685773, 0.446371, 0.139004,     &
189                                   0.976646, 0.912966, 0.737154, 0.473317, 0.145269/),   &
190                                   shape=(/PARASOL_NREFL,PARASOL_NTAU/)) 
191
192    ! ####################################################################################
193    ! ISCCP simulator tau/CTP joint histogram information
194    ! ####################################################################################
195    integer,parameter :: &
196         numISCCPTauBins  = ntau, &              ! Number of optical depth bins
197         numISCCPPresBins = npres                ! Number of pressure bins     
198    real(wp),parameter,dimension(ntau+1) :: &
199         isccp_histTau = tau_binBounds           ! Joint-histogram boundaries (optical depth)
200    real(wp),parameter,dimension(npres+1) :: &
201         isccp_histPres = pres_binBounds         ! Joint-histogram boundaries (cloud pressure)
202    real(wp),parameter,dimension(ntau) :: &
203         isccp_histTauCenters = tau_binCenters   ! Joint histogram bin centers (optical depth)
204    real(wp),parameter,dimension(npres) :: &   
205         isccp_histPresCenters = pres_binCenters ! Joint histogram bin centers (cloud pressure)
206    real(wp),parameter,dimension(2,ntau) :: &
207         isccp_histTauEdges = tau_binEdges       ! Joint histogram bin edges (optical depth)
208    real(wp),parameter,dimension(2,npres) :: &   
209         isccp_histPresEdges = pres_binEdges     ! Joint histogram bin edges (cloud pressure)   
210   
211    ! ####################################################################################
212    ! MISR simulator tau/CTH joint histogram information
213    ! ####################################################################################
214    integer,parameter ::  &
215         numMISRHgtBins = nhgt, &             ! Number of cloud-top height bins
216         numMISRTauBins = ntau                ! Number of optical depth bins
217    ! Joint histogram boundaries
218    real(wp),parameter,dimension(numMISRHgtBins+1) :: &
219         misr_histHgt = hgt_binBounds         ! Joint-histogram boundaries (cloud height)
220    real(wp),parameter,dimension(numMISRTauBins+1) :: &
221         misr_histTau = tau_binBounds         ! Joint-histogram boundaries (optical-depth)
222    real(wp),parameter,dimension(numMISRHgtBins) :: &
223         misr_histHgtCenters = hgt_binCenters ! Joint-histogram bin centers (cloud height)
224    real(wp),parameter,dimension(2,numMISRHgtBins) :: &
225         misr_histHgtEdges = hgt_BinEdges     ! Joint-histogram bin edges (cloud height)
226 
227    ! ####################################################################################
228    ! MODIS simulator tau/CTP joint histogram information
229    ! ####################################################################################
230    integer,parameter :: &
231         numMODISPresBins = npres                    ! Number of pressure bins for joint-histogram   
232    real(wp),parameter,dimension(numMODISPresBins + 1) :: &
233         modis_histPres = 100*pres_binBounds         ! Joint-histogram boundaries (cloud pressure)
234    real(wp),parameter,dimension(2, numMODISPresBins) :: &
235         modis_histPresEdges = 100*pres_binEdges     ! Joint-histogram bin edges (cloud pressure)
236    real(wp),parameter,dimension(numMODISPresBins) :: &
237         modis_histPresCenters = 100*pres_binCenters ! Joint-histogram bin centers (cloud pressure)
238
239    ! For the MODIS simulator we want to preserve the ability for cospV1.4.0 to use the
240    ! old histogram bin boundaries for optical depth, so these are set up in initialization.
241    integer :: &
242         numMODISTauBins          ! Number of tau bins for joint-histogram
243    real(wp),save,allocatable,dimension(:) :: &
244         modis_histTau            ! Joint-histogram boundaries (optical depth)
245    !$OMP THREADPRIVATE(modis_histTau)
246    real(wp),save,allocatable,dimension(:,:) :: &
247         modis_histTauEdges       ! Joint-histogram bin edges (optical depth)
248    !$OMP THREADPRIVATE(modis_histTauEdges)
249    real(wp),save,allocatable,dimension(:) :: &
250         modis_histTauCenters     ! Joint-histogram bin centers (optical depth)
251    !$OMP THREADPRIVATE(modis_histTauCenters)
252   
253    ! ####################################################################################
254    ! MODIS simulator tau/ReffICE and tau/ReffLIQ joint-histogram information
255    ! ####################################################################################
256    ! Ice
257    integer,parameter :: &
258         numMODISReffIceBins = nReffIce                ! Number of bins for joint-histogram
259    real(wp),parameter,dimension(nReffIce+1) :: &
260         modis_histReffIce = reffICE_binBounds         ! Effective radius bin boundaries
261    real(wp),parameter,dimension(nReffIce) :: &
262         modis_histReffIceCenters = reffICE_binCenters ! Effective radius bin centers
263    real(wp),parameter,dimension(2,nReffICE) :: &
264         modis_histReffIceEdges = reffICE_binEdges     ! Effective radius bin edges
265       
266    ! Liquid
267    integer,parameter :: &
268         numMODISReffLiqBins = nReffLiq                ! Number of bins for joint-histogram
269    real(wp),parameter,dimension(nReffLiq+1) :: &
270         modis_histReffLiq = reffLIQ_binBounds         ! Effective radius bin boundaries
271    real(wp),parameter,dimension(nReffLiq) :: &
272         modis_histReffLiqCenters = reffICE_binCenters ! Effective radius bin centers
273    real(wp),parameter,dimension(2,nReffICE) :: &
274         modis_histReffLiqEdges = reffLIQ_binEdges     ! Effective radius bin edges
275
276    ! ####################################################################################
277    ! CLOUDSAT reflectivity histogram information
278    ! ####################################################################################
279    integer,parameter :: &
280       CLOUDSAT_DBZE_BINS     =   15, & ! Number of dBZe bins in histogram (cfad)
281       CLOUDSAT_DBZE_MIN      = -100, & ! Minimum value for radar reflectivity
282       CLOUDSAT_DBZE_MAX      =   80, & ! Maximum value for radar reflectivity
283       CLOUDSAT_CFAD_ZE_MIN   =  -50, & ! Lower value of the first CFAD Ze bin
284       CLOUDSAT_CFAD_ZE_WIDTH =    5    ! Bin width (dBZe)
285
286    real(wp),parameter,dimension(CLOUDSAT_DBZE_BINS+1) :: &
287         cloudsat_histRef = (/CLOUDSAT_DBZE_MIN,(/(i, i=int(CLOUDSAT_CFAD_ZE_MIN+CLOUDSAT_CFAD_ZE_WIDTH),&
288                             int(CLOUDSAT_CFAD_ZE_MIN+(CLOUDSAT_DBZE_BINS-1)*CLOUDSAT_CFAD_ZE_WIDTH),    &
289                             int(CLOUDSAT_CFAD_ZE_WIDTH))/),CLOUDSAT_DBZE_MAX/)
290    real(wp),parameter,dimension(2,CLOUDSAT_DBZE_BINS) :: &
291         cloudsat_binEdges = reshape(source=(/cloudsat_histRef(1),((cloudsat_histRef(k), &
292                                   l=1,2),k=2,CLOUDSAT_DBZE_BINS),cloudsat_histRef(CLOUDSAT_DBZE_BINS+1)/),&
293                                   shape = (/2,CLOUDSAT_DBZE_BINS/))     
294    real(wp),parameter,dimension(CLOUDSAT_DBZE_BINS) :: &
295         cloudsat_binCenters = (cloudsat_binEdges(1,:)+cloudsat_binEdges(2,:))/2._wp
296   
297    ! Parameters for Cloudsat near-surface precipitation diagnostics.
298    ! Precipitation classes.
299    integer, parameter :: &
300         nCloudsatPrecipClass = 10
301    integer, parameter :: &
302         pClass_noPrecip      = 0, & ! No precipitation
303         pClass_Rain1         = 1, & ! Rain possible
304         pClass_Rain2         = 2, & ! Rain probable
305         pClass_Rain3         = 3, & ! Rain certain
306         pClass_Snow1         = 4, & ! Snow possible
307         pClass_Snow2         = 5, & ! Snow certain
308         pClass_Mixed1        = 6, & ! Mixed-precipitation possible
309         pClass_Mixed2        = 7, & ! Mixed-precipitation certain
310         pClass_Rain4         = 8, & ! Heavy rain
311         pClass_default       = 9    ! Default
312    ! Reflectivity bin boundaries, used by decision tree to classify precipitation type.
313    real(wp), dimension(4),parameter :: &
314         Zenonbinval =(/0._wp, -5._wp, -7.5_wp, -15._wp/)
315    real(wp), dimension(6),parameter :: &
316         Zbinvallnd = (/10._wp, 5._wp, 2.5_wp, -2.5_wp, -5._wp, -15._wp/)
317    ! Vertical level index(Nlvgrid) for Cloudsat precipitation occurence/frequency diagnostics.
318    ! Level 39 of Nlvgrid(40) is 480-960m.
319!    integer, parameter :: &    !PREC_BUG
320!         cloudsat_preclvl = 39 !PREC_BUG
321   
322    ! ####################################################################################
323    ! Parameters used by the CALIPSO LIDAR simulator
324    ! ####################################################################################
325    ! CALISPO backscatter histogram bins
326    real(wp),parameter ::     &
327       S_cld       = 5.0,     & ! Threshold for cloud detection
328       S_att       = 0.01,    & !
329       S_cld_att   = 30.        ! Threshold for undefined cloud phase detection
330    real(wp),parameter,dimension(SR_BINS+1) :: &
331         calipso_histBsct = (/-1.,0.01,1.2,3.0,5.0,7.0,10.0,15.0,20.0,25.0,30.0,40.0,50.0,   &
332                              60.0,80.0,999./)         ! Backscatter histogram bins
333    real(wp),parameter,dimension(2,SR_BINS) :: &
334         calipso_binEdges = reshape(source=(/calipso_histBsct(1),((calipso_histBsct(k),  &
335                                    l=1,2),k=2,SR_BINS),calipso_histBsct(SR_BINS+1)/),   &
336                                    shape = (/2,SR_BINS/))     
337    real(wp),parameter,dimension(SR_BINS) :: &
338         calipso_binCenters = (calipso_binEdges(1,:)+calipso_binEdges(2,:))/2._wp 
339
340    integer,parameter  ::     &
341       LIDAR_NTEMP = 40, &
342       LIDAR_NCAT  = 4,  & ! Number of categories for cloudtop heights (high/mid/low/tot)
343       LIDAR_NTYPE = 3     ! Number of categories for OPAQ (opaque/thin cloud + z_opaque)
344    real(wp),parameter,dimension(LIDAR_NTEMP) :: &
345       LIDAR_PHASE_TEMP=                                                                 &
346       (/-91.5,-88.5,-85.5,-82.5,-79.5,-76.5,-73.5,-70.5,-67.5,-64.5,                    &
347         -61.5,-58.5,-55.5,-52.5,-49.5,-46.5,-43.5,-40.5,-37.5,-34.5,                    &
348         -31.5,-28.5,-25.5,-22.5,-19.5,-16.5,-13.5,-10.5, -7.5, -4.5,                    &
349          -1.5,  1.5,  4.5,  7.5, 10.5, 13.5, 16.5, 19.5, 22.5, 25.5/)
350    real(wp),parameter,dimension(2,LIDAR_NTEMP) :: &
351       LIDAR_PHASE_TEMP_BNDS=reshape(source=                                             &
352          (/-273.15, -90., -90., -87., -87., -84., -84., -81., -81., -78.,               &
353             -78.,   -75., -75., -72., -72., -69., -69., -66., -66., -63.,               &
354             -63.,   -60., -60., -57., -57., -54., -54., -51., -51., -48.,               &
355             -48.,   -45., -45., -42., -42., -39., -39., -36., -36., -33.,               &
356             -33.,   -30., -30., -27., -27., -24., -24., -21., -21., -18.,               &
357             -18.,   -15., -15., -12., -12.,  -9.,  -9.,  -6.,  -6.,  -3.,               &
358              -3.,     0.,   0.,   3.,   3.,   6.,   6.,   9.,   9.,  12.,               &
359              12.,    15.,  15.,  18.,  18.,  21.,  21.,  24.,  24., 100. /),            &
360              shape=(/2,40/))       
361
362    ! ####################################################################################
363    ! Parameters used by the GROUND LIDAR simulator
364    ! ####################################################################################
365    ! GROUND LIDAR backscatter histogram bins
366!    real(wp),parameter ::     &
367!       S_cld       = 5.0,     & ! Threshold for cloud detection
368!       S_att       = 0.01,    & !
369!       S_cld_att   = 30.        ! Threshold for undefined cloud phase detection
370    real(wp),parameter,dimension(SR_BINS+1) :: &
371         grLidar532_histBsct = (/-1.,0.01,1.2,3.0,5.0,7.0,10.0,15.0,20.0,25.0,30.0,40.0,50.0,  &
372                                 60.0,80.0,999./)         ! Backscatter histogram bins
373    real(wp),parameter,dimension(2,SR_BINS) :: &
374         grLidar532_binEdges = reshape(source=(/grLidar532_histBsct(1),((grLidar532_histBsct(k),  &
375                                    l=1,2),k=2,SR_BINS),grLidar532_histBsct(SR_BINS+1)/),   &
376                                    shape = (/2,SR_BINS/))     
377    real(wp),parameter,dimension(SR_BINS) :: &
378         grLidar532_binCenters = (grLidar532_binEdges(1,:)+grLidar532_binEdges(2,:))/2._wp 
379
380!    integer,parameter  ::     &
381!       LIDAR_NCAT  = 4       ! Number of categories for cloudtop heights (high/mid/low/tot)
382
383    ! ####################################################################################
384    ! Parameters used by the ATLID LIDAR simulator
385    ! ####################################################################################
386    ! ATLID LIDAR backscatter histogram bins
387    real(wp),parameter ::     &
388       S_cld_atlid       = 1.74,    & ! Threshold for cloud detection
389       S_att_atlid       = 0.01,    & !
390       S_cld_att_atlid   = 6.67        ! Threshold for undefined cloud phase detection
391    real(wp),parameter,dimension(SR_BINS+1) :: &
392         atlid_histBsct = (/-1.,0.01,1.03,1.38,1.74,2.07,2.62,3.65,4.63,5.63,6.67,8.8,11.25,  &
393                                 13.2,17.2,999./)         ! Backscatter histogram bins
394    real(wp),parameter,dimension(2,SR_BINS) :: &
395         atlid_binEdges = reshape(source=(/atlid_histBsct(1),((atlid_histBsct(k),  &
396                                    l=1,2),k=2,SR_BINS),atlid_histBsct(SR_BINS+1)/),   &
397                                    shape = (/2,SR_BINS/))     
398    real(wp),parameter,dimension(SR_BINS) :: &
399         atlid_binCenters = (atlid_binEdges(1,:)+atlid_binEdges(2,:))/2._wp 
400
401!    integer,parameter  ::     &
402!       LIDAR_NCAT  = 4       ! Number of categories for cloudtop heights (high/mid/low/tot)
403
404    ! ####################################################################################
405    ! New vertical grid used by CALIPSO and CLOUDSAT L3 (set up during initialization)
406    ! ####################################################################################
407    integer :: &
408         Nlvgrid      ! Number of levels in New grid
409    real(wp),dimension(:), save, allocatable :: &
410       vgrid_zl,  & ! New grid bottoms
411       vgrid_zu,  & ! New grid tops
412       vgrid_z      ! New grid center
413    !$OMP THREADPRIVATE(vgrid_zl,vgrid_zu,vgrid_z)
414
415    ! ####################################################################################
416    ! Vertical grid used by CALIPSO and CLOUDSAT L3 (LMDZ parameter)
417    ! ####################################################################################
418
419    integer, parameter :: &
420         niv_sorties = 40
421    real(wp),parameter,dimension(niv_sorties) :: &
422         vgrid_z_in = (/240.0, 720.0, 1200.0, 1680.0, 2160.0, 2640.0, 3120.0, 3600.0, &
423                        4080.0, 4560.0, 5040.0, 5520.0, 6000.0, 6480.0, 6960.0, 7440.0, &
424                        7920.0, 8400.0, 8880.0, 9360.0, 9840.0, 10320.0, 10800.0, &
425                        11280.0, 11760.0, 12240.0, 12720.0, 13200.0, 13680.0, 14160.0, &
426                        14640.0, 15120.0, 15600.0, 16080.0, 16560.0, 17040.0, 17520.0, &
427                        18000.0, 18480.0, 18960.0/)
428
429END MODULE MOD_COSP_CONFIG
Note: See TracBrowser for help on using the repository browser.