Context Navigation

cosp_output_mod.F90 @ 5202

Last change on this file since 5202 was 4619, checked in by yann meurdesoif, 16 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.

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

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

Download in other formats:

Original Format