Context Navigation

cosp_output_mod.F90 @ 2884

Last change on this file since 2884 was 2822, checked in by idelkadi, 8 years ago

Ecriture des champs de sorties Cosp 4D dans les fichiers Netcdf en utilisant les fonctionalites Xios.

ajout d'une routine histwrite_cosp4D appelee si les sorties sont gerees par Xios (dans cosp_output_write.F90)
definitions des differents axes et grilles 4D pour Cosp (dans context_lmdz.xml)

Traitement des valeurs indefinies par Xios et mise a zero si l'on utilise Ioipsl.

Correction d'erreur sur les sorties Modis

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format