Context Navigation

cosp_output_mod.F90 @ 3543

Last change on this file since 3543 was 3435, checked in by Laurent Fairhead, 6 years ago

"Historic" :-) commit merging the physics branch used for DYNAMICO with the LMDZ trunk.
The same physics branch can now be used seamlessly with the traditional lon-lat LMDZ
dynamical core and DYNAMICO.
Testing consisted in running a lon-lat LMDZ bucket simulation with the NPv6.1 physics package
with the original trunk sources and the merged sources. Tests were succesful in the sense that
numeric continuity was preserved in the restart files from both simulation. Further tests
included running both versions of the physics codes for one year in a LMDZOR setting in which
the restart files also came out identical.

Caution:

as the physics package now manages unstructured grids, grid information needs to be transmitted

to the surface scheme ORCHIDEE. This means that the interface defined in surf_land_orchidee_mod.F90
is only compatible with ORCHIDEE version orchidee2.1 and later versions. If previous versions of
ORCHIDEE need to be used, the CPP key ORCHIDEE_NOUNSTRUCT needs to be set at compilation time.
This is done automatically if makelmdz/makelmdz_fcm are called with the veget orchidee2.0 switch

due to a limitation in XIOS, the time at which limit conditions will be read in by DYNAMICO will be

delayed by one physic timestep with respect to the time it is read in by the lon-lat model. This is caused
by the line

IF (MOD(itime-1, lmt_pas) == 0 .OR. (jour_lu /= jour .AND. grid_type /= unstructured)) THEN ! time to read

in limit_read_mod.F90

Work still needed on COSP integration and XML files for DYNAMICO

EM, YM, LF

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
224	CONTAINS
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
249	#ifdef CPP_XIOS
250	! ug Pour les sorties XIOS
251	USE wxios
252	#endif
253
254	IMPLICIT NONE
255
256	!!! Variables d'entree
257	integer :: Nlevlmdz, Ncolumns ! Number of levels
258	real,dimension(Nlevlmdz) :: presnivs
259	real :: dtime, freq_cosp, ecrit_day, ecrit_hf, ecrit_mth
260	logical :: use_vgrid
261	logical :: ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml
262	type(cosp_vgrid) :: vgrid ! Information on vertical grid of stats
263	! type(cosp_lidarstats) :: stlidar ! Summary statistics from lidar simulator
264
265	!!! Variables locales
266	integer :: idayref, iff, ii
267	real :: zjulian,zjulian_start
268	real(wp),dimension(Ncolumns) :: column_ax
269	real(wp),dimension(DBZE_BINS) :: dbze_ax
270	CHARACTER(LEN=20), DIMENSION(3) :: chfreq = (/ '1day', '1d ', '3h ' /)
271	real(wp),parameter,dimension(SR_BINS) :: sratio_ax = (/0.005, &
272	0.605,2.09,4.,6., &
273	8.5,12.5,17.5,22.5,27.5,35.,45.,55.,70.,50040./)
274
275	!!! Variables d'entree
276
277	#ifdef CPP_XIOS
278	! ug Variables utilisées pour récupérer le calendrier pour xios
279	INTEGER :: x_an, x_mois, x_jour
280	REAL :: x_heure
281	INTEGER :: ini_an, ini_mois, ini_jour
282	REAL :: ini_heure
283	#endif
284
285	WRITE(lunout,*) 'Debut cosp_output_mod.F90'
286	print*,'cosp_varsdefined',cosp_varsdefined
287	! Initialisations (Valeurs par defaut)
288
289	!! Definition valeurs axes
290	do ii=1,Ncolumns
291	column_ax(ii) = real(ii)
292	enddo
293
294	do i=1,DBZE_BINS
295	dbze_ax(i) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(i - 0.5)
296	enddo
297
298	cosp_outfilenames(1) = 'histmthCOSP'
299	cosp_outfilenames(2) = 'histdayCOSP'
300	cosp_outfilenames(3) = 'histhfCOSP'
301
302	cosp_outfiletypes(1) = 'ave(X)'
303	cosp_outfiletypes(2) = 'ave(X)'
304	cosp_outfiletypes(3) = 'ave(X)'
305
306	cosp_outfilekeys(1) = ok_mensuelCOSP
307	cosp_outfilekeys(2) = ok_journeCOSP
308	cosp_outfilekeys(3) = ok_hfCOSP
309
310	cosp_ecritfiles(1) = mth_len*86400.
311	cosp_ecritfiles(2) = 1.*86400.
312	cosp_ecritfiles(3) = 0.125*86400.
313
314	! Lecture des parametres dans output.def ou config.def
315
316	CALL getin('cosp_outfilenames',cosp_outfilenames)
317	CALL getin('cosp_outfilekeys',cosp_outfilekeys)
318	CALL getin('cosp_ecritfiles',cosp_ecritfiles)
319	CALL getin('cosp_outfiletypes',cosp_outfiletypes)
320
321	WRITE(lunout,*)'cosp_outfilenames=',cosp_outfilenames
322	WRITE(lunout,*)'cosp_outfilekeys=',cosp_outfilekeys
323	WRITE(lunout,*)'cosp_ecritfiles=',cosp_ecritfiles
324	WRITE(lunout,*)'cosp_outfiletypes=',cosp_outfiletypes
325
326	idayref = day_ref
327	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
328	CALL ymds2ju(annee_ref, 1, day_ini, start_time, zjulian_start)
329
330	#ifdef CPP_XIOS
331
332	! recuperer la valeur indefine Xios
333	! CALL xios_get_field_attr("clcalipso",default_value=Cosp_fill_value)
334	! Cosp_fill_value=missing_val
335	Cosp_fill_value=0.
336	print*,'Cosp_fill_value=',Cosp_fill_value
337
338	CALL wxios_add_vaxis("height", vgrid%Nlvgrid, vgrid%z)
339	print*,'wxios_add_vaxis vgrid%Nlvgrid, vgrid%z',vgrid%Nlvgrid,vgrid%z
340
341	WRITE(lunout,*) 'wxios_add_vaxis height_mlev, Nlevlmdz vgrid%mz ', &
342	Nlevlmdz,vgrid%mz
343	CALL wxios_add_vaxis("height_mlev", Nlevlmdz, vgrid%mz)
344
345	WRITE(lunout,*) 'wxios_add_vaxis sza, PARASOL_NREFL ', &
346	PARASOL_NREFL, PARASOL_SZA
347	CALL wxios_add_vaxis("sza", PARASOL_NREFL, PARASOL_SZA)
348
349	WRITE(lunout,*) 'wxios_add_vaxis pressure2 ',7,pres_binCenters
350	CALL wxios_add_vaxis("pressure2", 7, pres_binCenters)
351
352	WRITE(lunout,*) 'wxios_add_vaxis column ',Ncolumns,column_ax
353	CALL wxios_add_vaxis("column", Ncolumns, column_ax)
354
355	WRITE(lunout,*) 'wxios_add_vaxis temp LIDAR_NTEMP, LIDAR_PHASE_TEMP ', &
356	LIDAR_NTEMP, LIDAR_PHASE_TEMP
357	CALL wxios_add_vaxis("temp", LIDAR_NTEMP, LIDAR_PHASE_TEMP)
358
359	WRITE(lunout,*) 'wxios_add_vaxis cth16 numMISRHgtBins, misr_histHgtCenters ', &
360	numMISRHgtBins, misr_histHgtCenters
361	CALL wxios_add_vaxis("cth16", numMISRHgtBins, misr_histHgtCenters)
362
363	WRITE(lunout,*) 'wxios_add_vaxis dbze DBZE_BINS, dbze_ax ', &
364	DBZE_BINS, dbze_ax
365	CALL wxios_add_vaxis("dbze", DBZE_BINS, dbze_ax)
366
367	WRITE(lunout,*) 'wxios_add_vaxis scatratio SR_BINS, sratio_ax', &
368	SR_BINS, sratio_ax
369	CALL wxios_add_vaxis("scatratio", SR_BINS, sratio_ax)
370
371	WRITE(lunout,*) 'wxios_add_vaxis ReffIce numMODISReffIceBins, &
372	reffICE_binCenters',numMODISReffIceBins, reffICE_binCenters
373	CALL wxios_add_vaxis("ReffIce", numMODISReffIceBins, reffICE_binCenters)
374
375	WRITE(lunout,*) 'wxios_add_vaxis ReffLiq numMODISReffLiqBins, &
376	reffLIQ_binCenters', numMODISReffLiqBins, reffLIQ_binCenters
377	CALL wxios_add_vaxis("ReffLiq", numMODISReffLiqBins, reffLIQ_binCenters)
378
379	WRITE(lunout,*) 'wxios_add_vaxis 7, tau_binCenters', &
380	7, tau_binCenters
381	CALL wxios_add_vaxis("tau", 7, tau_binCenters)
382
383	#endif
384
385	zdtimemoy_cosp = freq_COSP ! Frequence ou l on moyenne
386
387	DO iff=1,3
388	zoutm_cosp(iff) = cosp_ecritfiles(iff) ! Frequence ou l on ecrit en seconde
389
390	IF (cosp_outfilekeys(iff)) THEN
391	CALL histbeg_phy_all(cosp_outfilenames(iff),itau_phy,zjulian,&
392	dtime,nhoricosp(iff),cosp_nidfiles(iff))
393	! print*,'histbeg_phy nhoricosp(iff),cosp_nidfiles(iff)', &
394	! nhoricosp(iff),cosp_nidfiles(iff)
395
396	#ifdef CPP_XIOS
397	IF (.not. ok_all_xml) then
398	WRITE(lunout,*) 'wxios_add_file ',cosp_outfilenames(iff)
399	CALL wxios_add_file(cosp_outfilenames(iff),chfreq(iff),10)
400	ENDIF
401	#endif
402
403	#ifndef CPP_IOIPSL_NO_OUTPUT
404	! Definition de l'axe vertical
405	if (use_vgrid) then
406	! Axe vertical Cosp 40 niveaux (en m)
407	CALL histvert(cosp_nidfiles(iff),"height","height","m",vgrid%Nlvgrid,vgrid%z,nvert(iff))
408	else
409	! Axe vertical modele LMDZ presnivs
410	CALL histvert(cosp_nidfiles(iff),"presnivs","Vertical levels","Pa",vgrid%Nlvgrid,presnivs,nvert(iff),"down")
411	endif
412	! Axe vertical niveaux modele (en m)
413	CALL histvert(cosp_nidfiles(iff),"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertmcosp(iff))
414
415	CALL histvert(cosp_nidfiles(iff),"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp(iff))
416
417	CALL histvert(cosp_nidfiles(iff),"pressure2","pressure","mb",7,pres_binCenters,nvertisccp(iff),"down")
418
419	CALL histvert(cosp_nidfiles(iff),"column","column","count",Ncolumns,column_ax,nvertcol(iff)) !DBUG
420
421	CALL histvert(cosp_nidfiles(iff),"temp","temperature","C",LIDAR_NTEMP,LIDAR_PHASE_TEMP,nverttemp(iff))
422
423	CALL histvert(cosp_nidfiles(iff),"cth","altitude","m",numMISRHgtBins,misr_histHgtCenters,nvertmisr(iff))
424
425	CALL histvert(cosp_nidfiles(iff),"ReffIce","Effective_particle_size_Ice","microns",numMODISReffIceBins, reffICE_binCenters, &
426	nvertReffIce(iff))
427
428	CALL histvert(cosp_nidfiles(iff),"ReffLiq","Effective_particle_size_Liq","microns",numMODISReffLiqBins, reffLIQ_binCenters, &
429	nvertReffLiq(iff))
430
431	CALL histvert(cosp_nidfiles(iff),"dbze","equivalent_reflectivity_factor","dBZ",DBZE_BINS,dbze_ax,nvertbze(iff))
432
433	CALL histvert(cosp_nidfiles(iff),"scatratio","backscattering_ratio","1",SR_BINS,sratio_ax,nvertsratio(iff))
434
435	CALL histvert(cosp_nidfiles(iff),"tau","cloud optical depth","1",7,tau_binCenters,nverttau(iff))
436
437	!!! Valeur indefinie en cas IOIPSL
438	Cosp_fill_value=0.
439
440	#endif
441
442	ENDIF
443	ENDDO
444
445	end SUBROUTINE cosp_output_open
446
447	END MODULE cosp_output_mod

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

Download in other formats:

Original Format