source: LMDZ6/trunk/libf/phylmd/cosp2/cosp_output_mod.F90 @ 4848

Last change on this file since 4848 was 4619, checked in by yann meurdesoif, 17 months ago

Suppress usage of preprocessing key CPP_XIOS.
Wrapper file is used to suppress XIOS symbol when xios is not linked and not used (-io ioipsl)
The CPP_XIOS key is replaced in model by "using_xios" boolean variable to switch between IOIPSL or XIOS output.

YM

File size: 25.6 KB
Line 
1! A.Idelkadi sept 2013 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2!! Module pour declarer et initialiser les parametres de controle des fichiers de sorties et des champs a sortir
3!! La routine cosp_output_open (appelee 1 seule fois dans phy_cosp.F90) permet :
4!! de creer les fichiers avec leurs grilles horizontales et verticales
5!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
6  MODULE cosp_output_mod
7
8!  USE MOD_COSP_CONSTANTS
9!  USE MOD_COSP_TYPES
10
11!  use MOD_COSP_INTERFACE_v1p4
12!  use MOD_COSP_CONFIG
13!  use MOD_COSP_Modis_Simulator, only : cosp_modis
14!  use mod_modis_sim, only : numMODISReffIceBins, reffICE_binCenters, &
15!                            numMODISReffLiqBins, reffLIQ_binCenters
16     USE COSP_KINDS, ONLY: wp,dp
17     IMPLICIT NONE
18! cosp_output_mod
19      INTEGER          :: i
20!!!!!!! Controle des fichier de sorties Cosp !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
21      LOGICAL, DIMENSION(3), SAVE  :: cosp_outfilekeys
22      INTEGER, DIMENSION(3), SAVE  :: cosp_nidfiles
23!$OMP THREADPRIVATE(cosp_outfilekeys, cosp_nidfiles)
24      INTEGER, DIMENSION(3), SAVE  :: nhoricosp,nvert,nvertmcosp,nvertcol,nvertbze, &
25                                      nvertsratio,nvertisccp,nvertp,nverttemp,nvertmisr, &
26                                      nvertReffIce,nvertReffLiq,nverttau
27      REAL, DIMENSION(3), SAVE                :: zoutm_cosp
28!$OMP THREADPRIVATE(nhoricosp, nvert,nvertmcosp,nvertcol,nvertsratio,nvertbze,nvertisccp,nvertp,zoutm_cosp,nverttemp,nvertmisr)
29!$OMP THREADPRIVATE(nvertReffIce,nvertReffLiq,nverttau)
30      REAL, SAVE                   :: zdtimemoy_cosp
31!$OMP THREADPRIVATE(zdtimemoy_cosp)
32      CHARACTER(LEN=20), DIMENSION(3), SAVE  :: cosp_outfiletypes
33      CHARACTER(LEN=20), DIMENSION(3), SAVE  :: cosp_outfilenames
34      REAL, DIMENSION(3), SAVE               :: cosp_ecritfiles
35!$OMP THREADPRIVATE(cosp_outfiletypes, cosp_outfilenames, cosp_ecritfiles)
36
37!!!!  Controle des variables a sortir dans les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
38  TYPE ctrl_outcosp
39     LOGICAL,DIMENSION(3)                 :: cles             !!! Sortir ou non le champs
40     CHARACTER(len=20)                    :: name       
41     CHARACTER(len=150)                   :: description      !!! Nom
42     CHARACTER(len=20)                    :: unit             !!! Unite
43     CHARACTER(len=20),DIMENSION(3)  :: cosp_typeecrit        !!! Operation (ave, inst, ...)
44  END TYPE ctrl_outcosp
45
46! CALIPSO vars
47  TYPE(ctrl_outcosp), SAVE :: o_cllcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
48         "cllcalipso", "Lidar Low-level Cloud Fraction", "1", (/ ('', i=1, 3) /))                                   
49  TYPE(ctrl_outcosp), SAVE :: o_clmcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
50         "clmcalipso", "Lidar Mid-level Cloud Fraction", "1", (/ ('', i=1, 3) /))
51  TYPE(ctrl_outcosp), SAVE :: o_clhcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
52         "clhcalipso", "Lidar High-level Cloud Fraction", "1", (/ ('', i=1, 3) /))
53  TYPE(ctrl_outcosp), SAVE :: o_cltcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
54         "cltcalipso", "Lidar Total Cloud Fraction", "1", (/ ('', i=1, 3) /))
55  TYPE(ctrl_outcosp), SAVE :: o_clcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
56         "clcalipso", "Lidar Cloud Fraction (532 nm)", "1", (/ ('', i=1, 3) /))
57  TYPE(ctrl_outcosp), SAVE :: o_cfad_lidarsr532 = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
58         "cfad_lidarsr532", "Lidar Scattering Ratio CFAD (532 nm)", "1", (/ ('', i=1, 3) /))   
59  TYPE(ctrl_outcosp), SAVE :: o_parasol_refl = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
60         "parasol_refl", "PARASOL-like mono-directional reflectance","1", (/ ('', i=1, 3) /))
61  TYPE(ctrl_outcosp), SAVE :: o_parasol_crefl = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &             
62         "parasol_crefl", "PARASOL-like mono-directional reflectance (integral)","1", (/ ('', i=1, 3) /))                 
63  TYPE(ctrl_outcosp), SAVE :: o_Ncrefl = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
64         "Ncrefl", "Nb PARASOL-like mono-directional reflectance (integral)","1", (/ ('', i=1, 3) /))
65  TYPE(ctrl_outcosp), SAVE :: o_atb532 = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
66         "atb532", "Lidar Attenuated Total Backscatter (532 nm)","1", (/ ('', i=1, 3) /))
67  TYPE(ctrl_outcosp), SAVE :: o_beta_mol532 = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
68         "beta_mol532", "Lidar Molecular Backscatter (532 nm)","m-1 sr-1", (/ ('', i=1, 3) /))
69!! AI  11 2015
70  TYPE(ctrl_outcosp), SAVE :: o_cllcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
71         "cllcalipsoice", "CALIPSO Ice-Phase Low Level Cloud Fraction", "%", (/ ('', i=1, 3) /)) 
72  TYPE(ctrl_outcosp), SAVE :: o_cllcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
73         "cllcalipsoliq", "CALIPSO Liq-Phase Low Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
74  TYPE(ctrl_outcosp), SAVE :: o_clmcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
75         "clmcalipsoice", "CALIPSO Ice-Phase Mid Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
76  TYPE(ctrl_outcosp), SAVE :: o_clmcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
77         "clmcalipsoliq", "CALIPSO Liq-Phase Mid Level Cloud Fraction", "%", (/ ('', i=1, 3) /))                 
78  TYPE(ctrl_outcosp), SAVE :: o_clhcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
79         "clhcalipsoice", "CALIPSO Ice-Phase High Level Cloud Fraction", "%", (/ ('', i=1, 3) /)) 
80  TYPE(ctrl_outcosp), SAVE :: o_clhcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
81         "clhcalipsoliq", "CALIPSO Liq-Phase High Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
82  TYPE(ctrl_outcosp), SAVE :: o_cltcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
83         "cltcalipsoice", "CALIPSO Ice-Phase Tot Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
84  TYPE(ctrl_outcosp), SAVE :: o_cltcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
85         "cltcalipsoliq", "CALIPSO Liq-Phase Tot Level Cloud Fraction", "%", (/ ('', i=1, 3) /))                 
86  TYPE(ctrl_outcosp), SAVE :: o_cllcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
87         "cllcalipsoun", "CALIPSO Undefined-Phase Low Level Cloud Fraction", "%", (/ ('', i=1, 3) /)) 
88  TYPE(ctrl_outcosp), SAVE :: o_clmcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
89         "clmcalipsoun", "CALIPSO Undefined-Phase Mid Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
90  TYPE(ctrl_outcosp), SAVE :: o_clhcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
91         "clhcalipsoun", "CALIPSO Undefined-Phase High Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
92  TYPE(ctrl_outcosp), SAVE :: o_cltcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
93         "cltcalipsoun", "CALIPSO Undefined-Phase Tot Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
94  TYPE(ctrl_outcosp), SAVE :: o_clcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
95         "clcalipsoice", "Lidar Ice-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
96  TYPE(ctrl_outcosp), SAVE :: o_clcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
97         "clcalipsoliq", "Lidar Liq-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
98  TYPE(ctrl_outcosp), SAVE :: o_clcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
99         "clcalipsoun", "Lidar Undef-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))     
100  TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmpice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
101         "clcalipsotmpice", "Lidar Ice-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
102  TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmpliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
103         "clcalipsotmpliq", "Lidar Liq-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
104  TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmpun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
105         "clcalipsotmpun", "Lidar Undef-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
106  TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
107         "clcalipsotmp", "Lidar Cloud Fraction", "%", (/ ('', i=1, 3) /))
108
109  TYPE(ctrl_outcosp), SAVE :: o_clopaquecalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &  !OPAQ
110         "clopaquecalipso", "Lidar Opaque Cloud Fraction", "%", (/ ('', i=1, 3) /))             !OPAQ
111  TYPE(ctrl_outcosp), SAVE :: o_clthincalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &    !OPAQ
112         "clthincalipso", "Lidar Thin Cloud Fraction", "%", (/ ('', i=1, 3) /))                 !OPAQ
113  TYPE(ctrl_outcosp), SAVE :: o_clzopaquecalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
114         "clzopaquecalipso", "Lidar mean opacity altitude", "m", (/ ('', i=1, 3) /))            !OPAQ
115  TYPE(ctrl_outcosp), SAVE :: o_clcalipsoopaque = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &  !OPAQ
116         "clcalipsoopaque", "Lidar Opaque profile Cloud Fraction", "%", (/ ('', i=1, 3) /))     !OPAQ
117  TYPE(ctrl_outcosp), SAVE :: o_clcalipsothin = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &    !OPAQ
118         "clcalipsothin", "Lidar Thin profile Cloud Fraction", "%", (/ ('', i=1, 3) /))         !OPAQ
119  TYPE(ctrl_outcosp), SAVE :: o_clcalipsozopaque = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
120         "clcalipsozopaque", "Lidar z_opaque Fraction", "%", (/ ('', i=1, 3) /))                !OPAQ
121  TYPE(ctrl_outcosp), SAVE :: o_clcalipsoopacity = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
122         "clcalipsoopacity", "Lidar opacity Fraction", "%", (/ ('', i=1, 3) /))                 !OPAQ
123
124  TYPE(ctrl_outcosp), SAVE :: o_proftemp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &         !TIBO
125         "proftemp", "Temperature profiles (40 lev)", "K", (/ ('', i=1, 3) /))                  !TIBO
126  TYPE(ctrl_outcosp), SAVE :: o_profSR = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &           !TIBO
127         "profSR", "Lidar Scattering Ratio profiles (532 nm)", "1", (/ ('', i=1, 3) /))         !TIBO
128
129! Radar Cloudsat
130  TYPE(ctrl_outcosp), SAVE :: o_cfadDbze94 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
131         "cfadDbze94", "CloudSat Radar Reflectivity CFAD", "%", (/ ('', i=1, 3) /))
132  TYPE(ctrl_outcosp), SAVE :: o_dbze94 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
133         "dbze94", "CloudSat Radar Reflectivity", "%", (/ ('', i=1, 3) /))
134
135! Calipso + Cloudsat
136  TYPE(ctrl_outcosp), SAVE :: o_clcalipso2 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
137         "clcalipso2", "CALIPSO Cloud Fraction Undetected by CloudSat", "1", (/ ('', i=1, 3) /))
138  TYPE(ctrl_outcosp), SAVE :: o_cltlidarradar = ctrl_outcosp((/ .TRUE., .TRUE.,.TRUE. /), &         
139         "cltlidarradar", "Lidar and Radar Total Cloud Fraction", "%", (/ ('', i=1, 3) /))
140     
141! ISCCP vars
142  TYPE(ctrl_outcosp), SAVE :: o_sunlit = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
143         "sunlit", "1 for day points, 0 for nightime","1",(/ ('', i=1, 3) /))                   
144  TYPE(ctrl_outcosp), SAVE :: o_clisccp2 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
145         "clisccp2", "Cloud Fraction as Calculated by the ISCCP Simulator","%", (/ ('', i=1, 3) /))
146  TYPE(ctrl_outcosp), SAVE :: o_boxtauisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
147         "boxtauisccp", "Optical Depth in Each Column as Calculated by the ISCCP Simulator","1", (/ ('', i=1, 3) /))
148  TYPE(ctrl_outcosp), SAVE :: o_boxptopisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
149         "boxptopisccp", "Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator","Pa", (/ ('', i=1, 3) /))
150  TYPE(ctrl_outcosp), SAVE :: o_tclisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
151          "tclisccp", "Total Cloud Fraction as Calculated by the ISCCP Simulator", "%", (/ ('', i=1, 3) /))
152  TYPE(ctrl_outcosp), SAVE :: o_ctpisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
153          "ctpisccp", "Mean Cloud Top Pressure as Calculated by the ISCCP Simulator", "Pa", (/ ('', i=1, 3) /))
154  TYPE(ctrl_outcosp), SAVE :: o_tauisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
155          "tauisccp", "Optical Depth as Calculated by the ISCCP Simulator", "1", (/ ('', i=1, 3) /))
156  TYPE(ctrl_outcosp), SAVE :: o_albisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
157          "albisccp", "Mean Cloud Albedo as Calculated by the ISCCP Simulator", "1", (/ ('', i=1, 3) /))
158  TYPE(ctrl_outcosp), SAVE :: o_meantbisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
159          "meantbisccp", " Mean all-sky 10.5 micron brightness temperature as calculated &
160           by the ISCCP Simulator","K", (/ ('', i=1, 3) /))
161  TYPE(ctrl_outcosp), SAVE :: o_meantbclrisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
162          "meantbclrisccp", "Mean clear-sky 10.5 micron brightness temperature as calculated &
163           by the ISCCP Simulator","K", (/ ('', i=1, 3) /))
164
165! MISR simulator
166  TYPE(ctrl_outcosp), SAVE :: o_clMISR = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
167         "clMISR", "Cloud Fraction as Calculated by the MISR Simulator","%", (/ ('', i=1, 3) /))
168
169! MODIS simulator
170  TYPE(ctrl_outcosp), SAVE :: o_cllmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
171         "cllmodis", "MODIS Low-level Cloud Fraction", "%", (/ ('', i=1, 3) /))                                   
172  TYPE(ctrl_outcosp), SAVE :: o_clmmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
173         "clmmodis", "MODIS Mid-level Cloud Fraction", "%", (/ ('', i=1, 3) /))
174  TYPE(ctrl_outcosp), SAVE :: o_clhmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
175         "clhmodis", "MODIS High-level Cloud Fraction", "%", (/ ('', i=1, 3) /))
176  TYPE(ctrl_outcosp), SAVE :: o_cltmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
177         "cltmodis", "MODIS Total Cloud Fraction", "%", (/ ('', i=1, 3) /))
178  TYPE(ctrl_outcosp), SAVE :: o_clwmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
179         "clwmodis", "MODIS Cloud Fraction water mean", "%", (/ ('', i=1, 3) /))
180  TYPE(ctrl_outcosp), SAVE :: o_climodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
181         "climodis", "MODIS Cloud Fraction ice mean", "%", (/ ('', i=1, 3) /))
182  TYPE(ctrl_outcosp), SAVE :: o_tautmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
183         "tautmodis", "MODIS Optical_Thickness_Total_Mean", "1", (/ ('', i=1, 3) /))                                   
184  TYPE(ctrl_outcosp), SAVE :: o_tauwmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
185         "tauwmodis", "MODIS Optical_Thickness_Water_Mean", "1", (/ ('', i=1, 3) /))
186  TYPE(ctrl_outcosp), SAVE :: o_tauimodis= ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
187         "tauimodis", "MODIS Optical_Thickness_Ice_Mean", "1", (/ ('', i=1, 3) /))
188  TYPE(ctrl_outcosp), SAVE :: o_tautlogmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
189         "tautlogmodis", "MODIS Optical_Thickness_Total_logMean", "1", (/ ('', i=1, 3) /))                                   
190  TYPE(ctrl_outcosp), SAVE :: o_tauwlogmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
191         "tauwlogmodis", "MODIS Optical_Thickness_Water_logMean", "1", (/ ('', i=1, 3) /))
192  TYPE(ctrl_outcosp), SAVE :: o_tauilogmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
193         "tauilogmodis", "MODIS Optical_Thickness_Ice_logMean", "1", (/ ('', i=1, 3) /))
194  TYPE(ctrl_outcosp), SAVE :: o_reffclwmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
195         "reffclwmodis", "Modis Cloud_Particle_Size_Water_Mean", "m", (/ ('', i=1, 3) /))
196  TYPE(ctrl_outcosp), SAVE :: o_reffclimodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
197         "reffclimodis", "Modis Cloud_Particle_Size_Ice_Mean", "m", (/ ('', i=1, 3) /))
198  TYPE(ctrl_outcosp), SAVE :: o_pctmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
199         "pctmodis", "Modis Cloud_Top_Pressure_Total_Mean", "Pa", (/ ('', i=1, 3) /))
200  TYPE(ctrl_outcosp), SAVE :: o_lwpmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
201         "lwpmodis", "Modis Liquid_Water_Path_Mean", "kg m-2", (/ ('', i=1, 3) /))
202  TYPE(ctrl_outcosp), SAVE :: o_iwpmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
203         "iwpmodis", "Modis Ice_Water_Path_Mean", "kg m-2", (/ ('', i=1, 3) /))
204  TYPE(ctrl_outcosp), SAVE :: o_clmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
205         "clmodis", "MODIS Cloud Area Fraction", "%", (/ ('', i=1, 3) /))
206  TYPE(ctrl_outcosp), SAVE :: o_crimodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
207         "crimodis", "Optical_Thickness_vs_ReffIce from Modis", "%", (/ ('',i=1, 3) /))         
208  TYPE(ctrl_outcosp), SAVE :: o_crlmodis = ctrl_outcosp((/ .TRUE., .TRUE.,.TRUE. /), &
209         "crlmodis", "Optical_Thickness_vs_ReffLiq from Modis", "%", (/ ('',i=1, 3) /))         
210
211! Rttovs simulator
212  TYPE(ctrl_outcosp), SAVE :: o_tbrttov = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
213         "tbrttov", "Rttovs Cloud Area Fraction", "%", (/ ('', i=1, 3) /))
214
215! Scops and others
216  TYPE(ctrl_outcosp), SAVE :: o_fracout = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
217         "fracout", "Subcolumn output from SCOPS", "%", (/ ('', i=1, 3) /))
218
219  LOGICAL, SAVE :: cosp_varsdefined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
220  REAL, SAVE  :: Cosp_fill_value
221!$OMP THREADPRIVATE(Cosp_fill_value)
222 
223
224CONTAINS
225
226!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
227!!!!!!!!! Ouverture des fichier et definition des  axes!!!!!!!!
228  !! histbeg, histvert
229!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
230
231  SUBROUTINE cosp_output_open(Nlevlmdz, Ncolumns, presnivs, dtime, freq_cosp, &
232                              ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml,  &
233                              ecrit_mth, ecrit_day, ecrit_hf, use_vgrid, vgrid)
234  use MOD_COSP_INTERFACE_v1p4, only :  cosp_vgrid
235!  use MOD_COSP
236  use mod_cosp_config, only : DBZE_BINS, SR_BINS, CFAD_ZE_MIN, PARASOL_NREFL, &
237                              CFAD_ZE_WIDTH,vgrid_zl,vgrid_zu,vgrid_z,PARASOL_SZA, &
238                              isccp_histPresCenters,tau_binCenters, LIDAR_NTEMP, &
239                              LIDAR_PHASE_TEMP,misr_histHgtCenters,numMISRHgtBins, &
240                              numMODISReffIceBins,reffICE_binCenters, &
241                              numMODISReffLiqBins, reffLIQ_binCenters, pres_binCenters
242
243  USE iophy
244  USE ioipsl
245  USE phys_cal_mod
246  USE time_phylmdz_mod, ONLY: day_ref, annee_ref, day_ini, start_time, itau_phy
247  USE print_control_mod, ONLY: lunout
248    ! ug Pour les sorties XIOS
249  USE wxios
250
251  IMPLICIT NONE
252
253!!! Variables d'entree
254  integer                  :: Nlevlmdz, Ncolumns      ! Number of levels
255  real,dimension(Nlevlmdz) :: presnivs
256  real                     :: dtime, freq_cosp, ecrit_day, ecrit_hf, ecrit_mth
257  logical                  :: use_vgrid
258  logical                  :: ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml
259  type(cosp_vgrid)   :: vgrid   ! Information on vertical grid of stats
260!  type(cosp_lidarstats) :: stlidar ! Summary statistics from lidar simulator
261
262!!! Variables locales
263  integer                   :: idayref, iff, ii
264  real                      :: zjulian,zjulian_start
265  real(wp),dimension(Ncolumns)  :: column_ax
266  real(wp),dimension(DBZE_BINS) ::  dbze_ax
267  CHARACTER(LEN=20), DIMENSION(3)  :: chfreq = (/ '1day', '1d  ', '3h  ' /)           
268  real(wp),parameter,dimension(SR_BINS) :: sratio_ax = (/0.005, &
269                                                  0.605,2.09,4.,6., &
270                                          8.5,12.5,17.5,22.5,27.5,35.,45.,55.,70.,50040./)
271
272!!! Variables d'entree
273
274! ug Variables utilisées pour récupérer le calendrier pour xios
275  INTEGER :: x_an, x_mois, x_jour
276  REAL :: x_heure
277  INTEGER :: ini_an, ini_mois, ini_jour
278  REAL :: ini_heure
279
280    WRITE(lunout,*) 'Debut cosp_output_mod.F90'
281    print*,'cosp_varsdefined',cosp_varsdefined
282    ! Initialisations (Valeurs par defaut)
283
284!! Definition valeurs axes
285    do ii=1,Ncolumns
286      column_ax(ii) = real(ii)
287    enddo
288
289    do i=1,DBZE_BINS
290     dbze_ax(i) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(i - 0.5)
291    enddo
292 
293    cosp_outfilenames(1) = 'histmthCOSP'
294    cosp_outfilenames(2) = 'histdayCOSP'
295    cosp_outfilenames(3) = 'histhfCOSP'
296
297    cosp_outfiletypes(1) = 'ave(X)'
298    cosp_outfiletypes(2) = 'ave(X)'
299    cosp_outfiletypes(3) = 'ave(X)'
300
301    cosp_outfilekeys(1) = ok_mensuelCOSP
302    cosp_outfilekeys(2) = ok_journeCOSP
303    cosp_outfilekeys(3) = ok_hfCOSP
304
305    cosp_ecritfiles(1) = mth_len*86400.
306    cosp_ecritfiles(2) = 1.*86400.
307    cosp_ecritfiles(3) = 0.125*86400.
308
309! Lecture des parametres dans output.def ou config.def
310
311    CALL getin('cosp_outfilenames',cosp_outfilenames)
312    CALL getin('cosp_outfilekeys',cosp_outfilekeys)
313    CALL getin('cosp_ecritfiles',cosp_ecritfiles)
314    CALL getin('cosp_outfiletypes',cosp_outfiletypes)
315
316    WRITE(lunout,*)'cosp_outfilenames=',cosp_outfilenames
317    WRITE(lunout,*)'cosp_outfilekeys=',cosp_outfilekeys
318    WRITE(lunout,*)'cosp_ecritfiles=',cosp_ecritfiles
319    WRITE(lunout,*)'cosp_outfiletypes=',cosp_outfiletypes
320   
321    idayref = day_ref
322    CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
323    CALL ymds2ju(annee_ref, 1, day_ini, start_time, zjulian_start)
324
325    IF (using_xios) THEN
326
327  ! recuperer la valeur indefine Xios
328  !    CALL xios_get_field_attr("clcalipso",default_value=Cosp_fill_value)
329  !         Cosp_fill_value=missing_val
330      Cosp_fill_value=0.
331      print*,'Cosp_fill_value=',Cosp_fill_value
332
333      CALL wxios_add_vaxis("height", vgrid%Nlvgrid, vgrid%z)
334      print*,'wxios_add_vaxis vgrid%Nlvgrid, vgrid%z',vgrid%Nlvgrid,vgrid%z
335
336      WRITE(lunout,*) 'wxios_add_vaxis height_mlev, Nlevlmdz vgrid%mz ', &
337                     Nlevlmdz,vgrid%mz
338      CALL wxios_add_vaxis("height_mlev", Nlevlmdz, vgrid%mz)
339
340      WRITE(lunout,*) 'wxios_add_vaxis sza, PARASOL_NREFL ', &
341                       PARASOL_NREFL, PARASOL_SZA
342      CALL wxios_add_vaxis("sza", PARASOL_NREFL, PARASOL_SZA)
343
344      WRITE(lunout,*) 'wxios_add_vaxis pressure2 ',7,pres_binCenters
345      CALL wxios_add_vaxis("pressure2", 7, pres_binCenters)
346
347      WRITE(lunout,*) 'wxios_add_vaxis column ',Ncolumns,column_ax
348      CALL wxios_add_vaxis("column", Ncolumns, column_ax)
349
350      WRITE(lunout,*) 'wxios_add_vaxis temp LIDAR_NTEMP, LIDAR_PHASE_TEMP ', &
351                    LIDAR_NTEMP, LIDAR_PHASE_TEMP
352      CALL wxios_add_vaxis("temp", LIDAR_NTEMP, LIDAR_PHASE_TEMP)
353
354      WRITE(lunout,*) 'wxios_add_vaxis cth16 numMISRHgtBins, misr_histHgtCenters ', &
355                      numMISRHgtBins, misr_histHgtCenters
356      CALL wxios_add_vaxis("cth16", numMISRHgtBins, misr_histHgtCenters)
357
358      WRITE(lunout,*) 'wxios_add_vaxis dbze DBZE_BINS, dbze_ax ', &
359                       DBZE_BINS, dbze_ax
360      CALL wxios_add_vaxis("dbze", DBZE_BINS, dbze_ax)
361
362      WRITE(lunout,*) 'wxios_add_vaxis scatratio SR_BINS, sratio_ax', &
363                     SR_BINS, sratio_ax
364      CALL wxios_add_vaxis("scatratio", SR_BINS, sratio_ax)
365
366      WRITE(lunout,*) 'wxios_add_vaxis ReffIce numMODISReffIceBins, &
367                     reffICE_binCenters',numMODISReffIceBins, reffICE_binCenters
368     CALL wxios_add_vaxis("ReffIce", numMODISReffIceBins, reffICE_binCenters)
369
370     WRITE(lunout,*) 'wxios_add_vaxis ReffLiq numMODISReffLiqBins, &
371                    reffLIQ_binCenters', numMODISReffLiqBins, reffLIQ_binCenters
372     CALL wxios_add_vaxis("ReffLiq", numMODISReffLiqBins, reffLIQ_binCenters)
373
374     WRITE(lunout,*) 'wxios_add_vaxis 7, tau_binCenters', &
375                      7, tau_binCenters
376     CALL wxios_add_vaxis("tau", 7, tau_binCenters)
377
378  ENDIF
379
380    zdtimemoy_cosp = freq_COSP         ! Frequence ou l on moyenne
381
382    DO iff=1,3
383       zoutm_cosp(iff) = cosp_ecritfiles(iff) ! Frequence ou l on ecrit en seconde
384
385       IF (cosp_outfilekeys(iff)) THEN
386           CALL histbeg_phy_all(cosp_outfilenames(iff),itau_phy,zjulian,&
387             dtime,nhoricosp(iff),cosp_nidfiles(iff))
388!           print*,'histbeg_phy nhoricosp(iff),cosp_nidfiles(iff)', &
389!                    nhoricosp(iff),cosp_nidfiles(iff)
390
391        IF(using_xios) THEN
392          IF (.not. ok_all_xml) then
393           WRITE(lunout,*) 'wxios_add_file ',cosp_outfilenames(iff)
394           CALL wxios_add_file(cosp_outfilenames(iff),chfreq(iff),10)
395          ENDIF
396        ENDIF
397
398#ifndef CPP_IOIPSL_NO_OUTPUT
399! Definition de l'axe vertical
400       if (use_vgrid) then
401! Axe vertical Cosp 40 niveaux (en m)
402      CALL histvert(cosp_nidfiles(iff),"height","height","m",vgrid%Nlvgrid,vgrid%z,nvert(iff))
403       else
404! Axe vertical modele LMDZ presnivs
405      CALL histvert(cosp_nidfiles(iff),"presnivs","Vertical levels","Pa",vgrid%Nlvgrid,presnivs,nvert(iff),"down")
406       endif
407! Axe vertical niveaux modele (en m)
408      CALL histvert(cosp_nidfiles(iff),"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertmcosp(iff))
409
410      CALL histvert(cosp_nidfiles(iff),"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp(iff))
411
412      CALL histvert(cosp_nidfiles(iff),"pressure2","pressure","mb",7,pres_binCenters,nvertisccp(iff),"down")
413
414      CALL histvert(cosp_nidfiles(iff),"column","column","count",Ncolumns,column_ax,nvertcol(iff)) !DBUG
415
416      CALL histvert(cosp_nidfiles(iff),"temp","temperature","C",LIDAR_NTEMP,LIDAR_PHASE_TEMP,nverttemp(iff))
417
418      CALL histvert(cosp_nidfiles(iff),"cth","altitude","m",numMISRHgtBins,misr_histHgtCenters,nvertmisr(iff))
419 
420      CALL histvert(cosp_nidfiles(iff),"ReffIce","Effective_particle_size_Ice","microns",numMODISReffIceBins, reffICE_binCenters, &
421                    nvertReffIce(iff))                                         
422     
423      CALL histvert(cosp_nidfiles(iff),"ReffLiq","Effective_particle_size_Liq","microns",numMODISReffLiqBins, reffLIQ_binCenters, &                                 
424                    nvertReffLiq(iff))
425
426      CALL histvert(cosp_nidfiles(iff),"dbze","equivalent_reflectivity_factor","dBZ",DBZE_BINS,dbze_ax,nvertbze(iff))
427     
428      CALL histvert(cosp_nidfiles(iff),"scatratio","backscattering_ratio","1",SR_BINS,sratio_ax,nvertsratio(iff))
429
430      CALL histvert(cosp_nidfiles(iff),"tau","cloud optical depth","1",7,tau_binCenters,nverttau(iff))
431     
432!!! Valeur indefinie en cas IOIPSL
433     Cosp_fill_value=0.
434
435#endif
436
437      ENDIF
438  ENDDO
439
440    end SUBROUTINE cosp_output_open
441
442 END MODULE cosp_output_mod
Note: See TracBrowser for help on using the repository browser.