1 | MODULE lmdz_lscp_ini |
---|
2 | |
---|
3 | IMPLICIT NONE |
---|
4 | |
---|
5 | ! PARAMETERS for lscp: |
---|
6 | !-------------------- |
---|
7 | |
---|
8 | REAL RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RV, RG, RPI, EPS_W |
---|
9 | !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RV, RG, RPI, EPS_W) |
---|
10 | |
---|
11 | REAL, SAVE, PROTECTED :: seuil_neb=0.001 ! cloud fraction threshold: a cloud can precipitate when exceeded |
---|
12 | !$OMP THREADPRIVATE(seuil_neb) |
---|
13 | |
---|
14 | REAL, SAVE, PROTECTED :: min_neb_th=1e-10 ! a cloud produced by bi-gaussian really exists when exceeded |
---|
15 | !$OMP THREADPRIVATE(min_neb_th) |
---|
16 | |
---|
17 | REAL, SAVE, PROTECTED :: min_frac_thermals=1.e-10 ! minimum thermals fraction for use of bigaussian distribution |
---|
18 | !$OMP THREADPRIVATE(min_frac_thermals) |
---|
19 | |
---|
20 | INTEGER, SAVE :: lunout, prt_level ! Logical unit number and level for standard output |
---|
21 | !$OMP THREADPRIVATE(lunout,prt_level) |
---|
22 | |
---|
23 | INTEGER, SAVE, PROTECTED :: niter_lscp=5 ! number of iterations to calculate autoconversion to precipitation |
---|
24 | !$OMP THREADPRIVATE(niter_lscp) |
---|
25 | |
---|
26 | INTEGER, SAVE, PROTECTED :: iflag_evap_prec=1 ! precipitation evaporation flag. 0: nothing, 1: "old way", |
---|
27 | ! 2: Max cloud fraction above to calculate the max of reevaporation |
---|
28 | ! >=4: LTP'method i.e. evaporation in the clear-sky fraction of the mesh only |
---|
29 | ! pay attention that iflag_evap_prec=4 may lead to unrealistic and overestimated |
---|
30 | ! evaporation. Use 5 instead |
---|
31 | !$OMP THREADPRIVATE(iflag_evap_prec) |
---|
32 | |
---|
33 | REAL, SAVE, PROTECTED :: t_coup=234.0 ! temperature threshold which determines the phase |
---|
34 | ! for which the saturation vapor pressure is calculated |
---|
35 | !$OMP THREADPRIVATE(t_coup) |
---|
36 | REAL, SAVE, PROTECTED :: DDT0=0.01 ! iteration parameter |
---|
37 | !$OMP THREADPRIVATE(DDT0) |
---|
38 | |
---|
39 | REAL, SAVE, PROTECTED :: ztfondue=278.15 ! parameter to calculate melting fraction of precipitation |
---|
40 | !$OMP THREADPRIVATE(ztfondue) |
---|
41 | |
---|
42 | REAL, SAVE, PROTECTED :: temp_nowater=235.15 ! temperature below which liquid water no longer exists |
---|
43 | !$OMP THREADPRIVATE(temp_nowater) |
---|
44 | |
---|
45 | REAL, SAVE, PROTECTED :: a_tr_sca(4) ! Variables for tracers temporary: alpha parameter for scavenging, 4 possible scavenging processes |
---|
46 | !$OMP THREADPRIVATE(a_tr_sca) |
---|
47 | |
---|
48 | REAL, SAVE, PROTECTED :: min_frac_th_cld=1.e-10 ! minimum thermal fraction to compute a thermal cloud fraction |
---|
49 | !$OMP THREADPRIVATE(min_frac_th_cld) |
---|
50 | |
---|
51 | LOGICAL, SAVE, PROTECTED :: ok_radocond_snow=.false. ! take into account the mass of ice precip in the cloud ice content seen by radiation |
---|
52 | !$OMP THREADPRIVATE(ok_radocond_snow) |
---|
53 | |
---|
54 | REAL, SAVE, PROTECTED :: t_glace_min=258.0 ! lower-bound temperature parameter for cloud phase determination |
---|
55 | !$OMP THREADPRIVATE(t_glace_min) |
---|
56 | |
---|
57 | REAL, SAVE, PROTECTED :: t_glace_max=273.15 ! upper-bound temperature parameter for cloud phase determination |
---|
58 | !$OMP THREADPRIVATE(t_glace_max) |
---|
59 | |
---|
60 | REAL, SAVE, PROTECTED :: exposant_glace=1.0 ! parameter for cloud phase determination |
---|
61 | !$OMP THREADPRIVATE(exposant_glace) |
---|
62 | |
---|
63 | INTEGER, SAVE, PROTECTED :: iflag_vice=0 ! which expression for ice crystall fall velocity |
---|
64 | !$OMP THREADPRIVATE(iflag_vice) |
---|
65 | |
---|
66 | INTEGER, SAVE, PROTECTED :: iflag_t_glace=0 ! which expression for cloud phase partitioning |
---|
67 | !$OMP THREADPRIVATE(iflag_t_glace) |
---|
68 | |
---|
69 | INTEGER, SAVE, PROTECTED :: iflag_cloudth_vert=0 ! option for determining cloud fraction and content in convective boundary layers |
---|
70 | !$OMP THREADPRIVATE(iflag_cloudth_vert) |
---|
71 | |
---|
72 | INTEGER, SAVE, PROTECTED :: iflag_gammasat=0 ! which threshold for homogeneous nucleation below -40oC |
---|
73 | !$OMP THREADPRIVATE(iflag_gammasat) |
---|
74 | |
---|
75 | INTEGER, SAVE, PROTECTED :: iflag_rain_incloud_vol=0 ! use of volume cloud fraction for rain autoconversion |
---|
76 | !$OMP THREADPRIVATE(iflag_rain_incloud_vol) |
---|
77 | |
---|
78 | INTEGER, SAVE, PROTECTED :: iflag_bergeron=0 ! bergeron effect for liquid precipitation treatment |
---|
79 | !$OMP THREADPRIVATE(iflag_bergeron) |
---|
80 | |
---|
81 | INTEGER, SAVE, PROTECTED :: iflag_fisrtilp_qsat=0 ! qsat adjustment (iterative) during autoconversion |
---|
82 | !$OMP THREADPRIVATE(iflag_fisrtilp_qsat) |
---|
83 | |
---|
84 | INTEGER, SAVE, PROTECTED :: iflag_pdf=0 ! type of subgrid scale qtot pdf |
---|
85 | !$OMP THREADPRIVATE(iflag_pdf) |
---|
86 | |
---|
87 | INTEGER, SAVE, PROTECTED :: iflag_icefrac=0 ! which phase partitioning function to use |
---|
88 | !$OMP THREADPRIVATE(iflag_icefrac) |
---|
89 | |
---|
90 | INTEGER, SAVE, PROTECTED :: iflag_autoconversion=0 ! autoconversion option |
---|
91 | !$OMP THREADPRIVATE(iflag_autoconversion) |
---|
92 | |
---|
93 | LOGICAL, SAVE, PROTECTED :: ok_bug_phase_lscp=.true. ! bug on phase partitioning after precipitation processes |
---|
94 | !$OMP THREADPRIVATE(ok_bug_phase_lscp) |
---|
95 | |
---|
96 | LOGICAL, SAVE, PROTECTED :: reevap_ice=.false. ! no liquid precip for T< threshold |
---|
97 | !$OMP THREADPRIVATE(reevap_ice) |
---|
98 | |
---|
99 | REAL, SAVE, PROTECTED :: cld_lc_lsc=2.6e-4 ! liquid autoconversion coefficient, stratiform rain |
---|
100 | !$OMP THREADPRIVATE(cld_lc_lsc) |
---|
101 | |
---|
102 | REAL, SAVE, PROTECTED :: cld_lc_con=2.6e-4 ! liquid autoconversion coefficient, convective rain |
---|
103 | !$OMP THREADPRIVATE(cld_lc_con) |
---|
104 | |
---|
105 | REAL, SAVE, PROTECTED :: cld_tau_lsc=3600. ! liquid autoconversion timescale, stratiform rain |
---|
106 | !$OMP THREADPRIVATE(cld_tau_lsc) |
---|
107 | |
---|
108 | REAL, SAVE, PROTECTED :: cld_tau_con=3600. ! liquid autoconversion timescale, convective rain |
---|
109 | !$OMP THREADPRIVATE(cld_tau_con) |
---|
110 | |
---|
111 | REAL, SAVE, PROTECTED :: cld_expo_lsc=2. ! liquid autoconversion threshold exponent, stratiform rain |
---|
112 | !$OMP THREADPRIVATE(cld_expo_lsc) |
---|
113 | |
---|
114 | REAL, SAVE, PROTECTED :: cld_expo_con=2. ! liquid autoconversion threshold exponent, convective rain |
---|
115 | !$OMP THREADPRIVATE(cld_expo_con) |
---|
116 | |
---|
117 | REAL, SAVE, PROTECTED :: ffallv_lsc=1. ! tuning coefficient crystal fall velocity, stratiform |
---|
118 | !$OMP THREADPRIVATE(ffallv_lsc) |
---|
119 | |
---|
120 | REAL, SAVE, PROTECTED :: ffallv_con=1. ! tuning coefficient crystal fall velocity, convective |
---|
121 | !$OMP THREADPRIVATE(ffallv_con) |
---|
122 | |
---|
123 | REAL, SAVE, PROTECTED :: coef_eva=2e-5 ! tuning coefficient liquid precip evaporation |
---|
124 | !$OMP THREADPRIVATE(coef_eva) |
---|
125 | |
---|
126 | REAL, SAVE, PROTECTED :: coef_sub ! tuning coefficient ice precip sublimation |
---|
127 | !$OMP THREADPRIVATE(coef_sub) |
---|
128 | |
---|
129 | REAL, SAVE, PROTECTED :: expo_eva=0.5 ! tuning coefficient liquid precip evaporation |
---|
130 | !$OMP THREADPRIVATE(expo_eva) |
---|
131 | |
---|
132 | REAL, SAVE, PROTECTED :: expo_sub ! tuning coefficient ice precip sublimation |
---|
133 | !$OMP THREADPRIVATE(expo_sub) |
---|
134 | |
---|
135 | REAL, SAVE, PROTECTED :: cice_velo=1.645 ! factor in the ice fall velocity formulation |
---|
136 | !$OMP THREADPRIVATE(cice_velo) |
---|
137 | |
---|
138 | REAL, SAVE, PROTECTED :: dice_velo=0.16 ! exponent in the ice fall velocity formulation |
---|
139 | !$OMP THREADPRIVATE(dice_velo) |
---|
140 | |
---|
141 | REAL, SAVE, PROTECTED :: dist_liq=300. ! typical deph of cloud-top liquid layer in mpcs |
---|
142 | !$OMP THREADPRIVATE(dist_liq) |
---|
143 | |
---|
144 | REAL, SAVE, PROTECTED :: tresh_cl=0.0 ! cloud fraction threshold for cloud top search |
---|
145 | !$OMP THREADPRIVATE(tresh_cl) |
---|
146 | |
---|
147 | !--Parameters for condensation and ice supersaturation |
---|
148 | |
---|
149 | LOGICAL, SAVE, PROTECTED :: ok_ice_supersat=.FALSE. ! activates the condensation scheme that allows for ice supersaturation |
---|
150 | !$OMP THREADPRIVATE(ok_ice_supersat) |
---|
151 | |
---|
152 | LOGICAL, SAVE, PROTECTED :: ok_unadjusted_clouds=.FALSE. ! if True, relax the saturation adjustment assumption for ice clouds |
---|
153 | !$OMP THREADPRIVATE(ok_unadjusted_clouds) |
---|
154 | |
---|
155 | LOGICAL, SAVE, PROTECTED :: ok_weibull_warm_clouds=.FALSE. ! if True, the weibull condensation scheme replaces the lognormal condensation scheme at positive temperatures |
---|
156 | !$OMP THREADPRIVATE(ok_weibull_warm_clouds) |
---|
157 | |
---|
158 | INTEGER, SAVE, PROTECTED :: iflag_cloud_sublim_pdf=4 ! iflag for the distribution of water inside ice clouds |
---|
159 | !$OMP THREADPRIVATE(iflag_cloud_sublim_pdf) |
---|
160 | |
---|
161 | REAL, SAVE, PROTECTED :: depo_coef_cirrus=.7 ! [-] deposition coefficient for growth of ice crystals in cirrus clouds |
---|
162 | !$OMP THREADPRIVATE(depo_coef_cirrus) |
---|
163 | |
---|
164 | REAL, SAVE, PROTECTED :: capa_cond_cirrus=.5 ! [-] capacitance factor for growth/sublimation of ice crystals in cirrus clouds |
---|
165 | !$OMP THREADPRIVATE(capa_cond_cirrus) |
---|
166 | |
---|
167 | REAL, SAVE, PROTECTED :: std_subl_pdf_lscp=2. ! [%] standard deviation of the gaussian distribution of water inside ice clouds |
---|
168 | !$OMP THREADPRIVATE(std_subl_pdf_lscp) |
---|
169 | |
---|
170 | REAL, SAVE, PROTECTED :: mu_subl_pdf_lscp=1./3. ! [-] shape factor of the gamma distribution of water inside ice clouds |
---|
171 | !$OMP THREADPRIVATE(mu_subl_pdf_lscp) |
---|
172 | |
---|
173 | REAL, SAVE, PROTECTED :: beta_pdf_lscp=1.E-3 ! [SI] tuning coefficient for the standard deviation of the PDF of water vapor in the clear sky region |
---|
174 | !$OMP THREADPRIVATE(beta_pdf_lscp) |
---|
175 | |
---|
176 | REAL, SAVE, PROTECTED :: temp_thresh_pdf_lscp=189. ! [K] factor for the PDF fit of water vapor in UTLS - below this temperature, water vapor is homogeneously distributed in the clear sky region |
---|
177 | !$OMP THREADPRIVATE(temp_thresh_pdf_lscp) |
---|
178 | |
---|
179 | REAL, SAVE, PROTECTED :: k0_pdf_lscp=3.01 ! [-] factor for the PDF fit of water vapor in UTLS |
---|
180 | !$OMP THREADPRIVATE(k0_pdf_lscp) |
---|
181 | |
---|
182 | REAL, SAVE, PROTECTED :: kappa_pdf_lscp=0.0192 ! [] factor for the PDF fit of water vapor in UTLS |
---|
183 | !$OMP THREADPRIVATE(kappa_pdf_lscp) |
---|
184 | |
---|
185 | REAL, SAVE, PROTECTED :: std100_pdf_lscp=4.08 ! [%] standard deviation at RHliq=100% of the PDF fit of water vapor in UTLS |
---|
186 | !$OMP THREADPRIVATE(std100_pdf_lscp) |
---|
187 | |
---|
188 | REAL, SAVE, PROTECTED :: a_homofreez=2.349 ! [-] factor for the Koop homogeneous freezing fit |
---|
189 | !$OMP THREADPRIVATE(a_homofreez) |
---|
190 | |
---|
191 | REAL, SAVE, PROTECTED :: b_homofreez=259. ! [K] factor for the Koop homogeneous freezing fit |
---|
192 | !$OMP THREADPRIVATE(b_homofreez) |
---|
193 | |
---|
194 | REAL, SAVE, PROTECTED :: delta_hetfreez=1. ! [-] value between 0 and 1 to simulate for heterogeneous freezing. |
---|
195 | !$OMP THREADPRIVATE(delta_hetfreez) |
---|
196 | |
---|
197 | REAL, SAVE, PROTECTED :: coef_mixing_lscp=1.E-3 ! [-] tuning coefficient for the mixing process |
---|
198 | !$OMP THREADPRIVATE(coef_mixing_lscp) |
---|
199 | |
---|
200 | REAL, SAVE, PROTECTED :: coef_shear_lscp=0.72 ! [-] additional coefficient for the shearing process (subprocess of the mixing process) |
---|
201 | !$OMP THREADPRIVATE(coef_shear_lscp) |
---|
202 | |
---|
203 | REAL, SAVE, PROTECTED :: chi_mixing_lscp=1. ! [-] factor for the macro distribution of ISSRs wrt clouds in a gridbox |
---|
204 | !$OMP THREADPRIVATE(chi_mixing_lscp) |
---|
205 | !--End of the parameters for condensation and ice supersaturation |
---|
206 | |
---|
207 | !--Parameters for aviation |
---|
208 | |
---|
209 | LOGICAL, SAVE, PROTECTED :: ok_plane_contrail ! activates the contrails parameterisation |
---|
210 | !$OMP THREADPRIVATE(ok_plane_contrail) |
---|
211 | |
---|
212 | REAL, SAVE, PROTECTED :: aspect_ratio_contrails=.1 ! [-] aspect ratio of the contrails clouds |
---|
213 | !$OMP THREADPRIVATE(aspect_ratio_contrails) |
---|
214 | |
---|
215 | REAL, SAVE, PROTECTED :: coef_mixing_contrails ! [-] tuning coefficient for the contrails mixing process |
---|
216 | !$OMP THREADPRIVATE(coef_mixing_contrails) |
---|
217 | |
---|
218 | REAL, SAVE, PROTECTED :: coef_shear_contrails ! [-] additional coefficient for the contrails shearing process (subprocess of the contrails mixing process) |
---|
219 | !$OMP THREADPRIVATE(coef_shear_contrails) |
---|
220 | |
---|
221 | REAL, SAVE, PROTECTED :: chi_mixing_contrails ! [-] factor for the macro distribution of ISSRs wrt clouds in a gridbox |
---|
222 | !$OMP THREADPRIVATE(chi_mixing_contrails) |
---|
223 | |
---|
224 | REAL, SAVE, PROTECTED :: rm_ice_crystals_contrails=7.5E-6 ! [m] geometric radius of ice crystals in contrails |
---|
225 | !$OMP THREADPRIVATE(rm_ice_crystals_contrails) |
---|
226 | |
---|
227 | REAL, SAVE, PROTECTED :: EI_H2O_aviation=1.25 ! [kgH2O/kg] emission index of water vapor for a given fuel type |
---|
228 | !$OMP THREADPRIVATE(EI_H2O_aviation) |
---|
229 | |
---|
230 | REAL, SAVE, PROTECTED :: qheat_fuel_aviation=43.2E6 ! [J/kg] specific combustion heat for a given fuel type |
---|
231 | !$OMP THREADPRIVATE(qheat_fuel_aviation) |
---|
232 | |
---|
233 | REAL, SAVE, PROTECTED :: prop_efficiency_aviation=.3 ! [-] average propulsion efficiency of aircraft |
---|
234 | !$OMP THREADPRIVATE(prop_efficiency_aviation) |
---|
235 | |
---|
236 | REAL, SAVE, PROTECTED :: linear_contrails_lifetime=10800. ! [s] timescale of the lifetime of linear contrails |
---|
237 | !$OMP THREADPRIVATE(linear_contrails_lifetime) |
---|
238 | |
---|
239 | REAL, SAVE, PROTECTED :: initial_width_contrails=200. ! [m] initial width of the linear contrails formed |
---|
240 | !$OMP THREADPRIVATE(initial_width_contrails) |
---|
241 | !--End of the parameters for aviation |
---|
242 | |
---|
243 | !--Parameters for poprecip |
---|
244 | LOGICAL, SAVE, PROTECTED :: ok_poprecip=.FALSE. ! use the processes-oriented formulation of precipitations |
---|
245 | !$OMP THREADPRIVATE(ok_poprecip) |
---|
246 | |
---|
247 | LOGICAL, SAVE, PROTECTED :: ok_corr_vap_evasub=.FALSE. ! use the corrected version of clear-sky water vapor for the evap / subl processes |
---|
248 | !$OMP THREADPRIVATE(ok_corr_vap_evasub) |
---|
249 | |
---|
250 | REAL, SAVE, PROTECTED :: cld_lc_lsc_snow=2.e-5 ! snow autoconversion coefficient, stratiform. default from Chaboureau and PInty 2006 |
---|
251 | !$OMP THREADPRIVATE(cld_lc_lsc_snow) |
---|
252 | |
---|
253 | REAL, SAVE, PROTECTED :: cld_lc_con_snow=2.e-5 ! snow autoconversion coefficient, convective |
---|
254 | !$OMP THREADPRIVATE(cld_lc_con_snow) |
---|
255 | |
---|
256 | REAL, SAVE, PROTECTED :: rain_int_min=1.e-5 ! Minimum local rain intensity [mm/s] before the decrease in associated precipitation fraction |
---|
257 | !$OMP THREADPRIVATE(rain_int_min) |
---|
258 | |
---|
259 | REAL, SAVE, PROTECTED :: thresh_precip_frac=1.E-6 ! precipitation fraction threshold [-] |
---|
260 | !$OMP THREADPRIVATE(thresh_precip_frac) |
---|
261 | |
---|
262 | REAL, SAVE, PROTECTED :: capa_crystal=1. ! Crystal capacitance (shape factor) for lscp_icefrac_turb [-] |
---|
263 | !$OMP THREADPRIVATE(capa_crystal) |
---|
264 | |
---|
265 | REAL, SAVE, PROTECTED :: naero5=0.5 ! Number concentration of aerosol larger than 0.5 microns [scm-3] |
---|
266 | !$OMP THREADPRIVATE(naero5) |
---|
267 | |
---|
268 | REAL, SAVE, PROTECTED :: gamma_snwretro = 0. ! Proportion of snow taken into account in ice retroaction in lscp_icefrac_turb [-] |
---|
269 | !$OMP THREADPRIVATE(gamma_snwretro) |
---|
270 | |
---|
271 | REAL, SAVE, PROTECTED :: gamma_taud = 1. ! Tuning coeff for Lagrangian decorrelation timescale in lscp_icefrac_turb [-] |
---|
272 | !$OMP THREADPRIVATE(gamma_taud) |
---|
273 | |
---|
274 | REAL, SAVE, PROTECTED :: gamma_col=1. ! Tuning coefficient for rain collection efficiency (poprecip) [-] |
---|
275 | !$OMP THREADPRIVATE(gamma_col) |
---|
276 | |
---|
277 | REAL, SAVE, PROTECTED :: gamma_agg=1. ! Tuning coefficient for snow aggregation efficiency (poprecip) [-] |
---|
278 | !$OMP THREADPRIVATE(gamma_agg) |
---|
279 | |
---|
280 | REAL, SAVE, PROTECTED :: gamma_rim=1. ! Tuning coefficient for riming efficiency (poprecip) [-] |
---|
281 | !$OMP THREADPRIVATE(gamma_rim) |
---|
282 | |
---|
283 | REAL, SAVE, PROTECTED :: gamma_melt=1. ! Tuning coefficient for snow melting efficiency (poprecip) [-] |
---|
284 | !$OMP THREADPRIVATE(gamma_melt) |
---|
285 | |
---|
286 | REAL, SAVE, PROTECTED :: gamma_freez=1. ! Tuning coefficient for rain collision freezing efficiency (poprecip) [-] |
---|
287 | !$OMP THREADPRIVATE(gamma_freez) |
---|
288 | |
---|
289 | REAL, SAVE, PROTECTED :: rho_rain=1000. ! Rain density [kg/m3] |
---|
290 | !$OMP THREADPRIVATE(rho_rain) |
---|
291 | |
---|
292 | REAL, SAVE, PROTECTED :: rho_ice=920. ! Ice density [kg/m3] |
---|
293 | !$OMP THREADPRIVATE(rho_ice) |
---|
294 | |
---|
295 | REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! Rain droplets radius (poprecip) [m] |
---|
296 | !$OMP THREADPRIVATE(r_rain) |
---|
297 | |
---|
298 | REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! Ice crystals radius (poprecip) [m] |
---|
299 | !$OMP THREADPRIVATE(r_snow) |
---|
300 | |
---|
301 | REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! Snow autoconversion minimal timescale (when liquid) [s] |
---|
302 | !$OMP THREADPRIVATE(tau_auto_snow_min) |
---|
303 | |
---|
304 | REAL, SAVE, PROTECTED :: tau_auto_snow_max=1000. ! Snow autoconversion minimal timescale (when only ice) [s] |
---|
305 | !$OMP THREADPRIVATE(tau_auto_snow_max) |
---|
306 | |
---|
307 | REAL, SAVE, PROTECTED :: eps=1.E-10 ! Treshold 0 [-] |
---|
308 | !$OMP THREADPRIVATE(eps) |
---|
309 | |
---|
310 | REAL, SAVE, PROTECTED :: alpha_freez=4. ! Slope of exponential for immersion freezing timescale [-] |
---|
311 | !$OMP THREADPRIVATE(alpha_freez) |
---|
312 | |
---|
313 | REAL, SAVE, PROTECTED :: beta_freez=0.1 ! Inv.time immersion freezing [s-1] |
---|
314 | !$OMP THREADPRIVATE(beta_freez) |
---|
315 | |
---|
316 | REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! Rain fall velocity [m/s] |
---|
317 | !$OMP THREADPRIVATE(rain_fallspeed) |
---|
318 | |
---|
319 | REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! Rain fall velocity in clear sky [m/s] |
---|
320 | !$OMP THREADPRIVATE(rain_fallspeed_clr) |
---|
321 | |
---|
322 | REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! Rain fall velocity in cloudy sky [m/s] |
---|
323 | !$OMP THREADPRIVATE(rain_fallspeed_cld) |
---|
324 | |
---|
325 | REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! Snow fall velocity [m/s] |
---|
326 | !$OMP THREADPRIVATE(snow_fallspeed) |
---|
327 | |
---|
328 | REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! Snow fall velocity in clear sky [m/s] |
---|
329 | !$OMP THREADPRIVATE(snow_fallspeed_clr) |
---|
330 | |
---|
331 | REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! Snow fall velocity in cloudy sky [m/s] |
---|
332 | !$OMP THREADPRIVATE(snow_fallspeed_cld) |
---|
333 | !--End of the parameters for poprecip |
---|
334 | |
---|
335 | ! Two parameters used for lmdz_lscp_old only |
---|
336 | INTEGER, SAVE, PROTECTED :: iflag_oldbug_fisrtilp=0, fl_cor_ebil |
---|
337 | !$OMP THREADPRIVATE(iflag_oldbug_fisrtilp,fl_cor_ebil) |
---|
338 | |
---|
339 | CONTAINS |
---|
340 | |
---|
341 | SUBROUTINE lscp_ini(dtime, lunout_in, prt_level_in, ok_ice_supersat_in, ok_plane_contrail_in, & |
---|
342 | iflag_ratqs, fl_cor_ebil_in, & |
---|
343 | RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, RVTMP2_in, & |
---|
344 | RTT_in, RD_in, RV_in, RG_in, RPI_in, EPS_W_in) |
---|
345 | |
---|
346 | |
---|
347 | USE ioipsl_getin_p_mod, ONLY : getin_p |
---|
348 | USE lmdz_cloudth_ini, ONLY : cloudth_ini |
---|
349 | |
---|
350 | REAL, INTENT(IN) :: dtime |
---|
351 | INTEGER, INTENT(IN) :: lunout_in,prt_level_in,iflag_ratqs,fl_cor_ebil_in |
---|
352 | LOGICAL, INTENT(IN) :: ok_ice_supersat_in, ok_plane_contrail_in |
---|
353 | |
---|
354 | REAL, INTENT(IN) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in |
---|
355 | REAL, INTENT(IN) :: RVTMP2_in, RTT_in, RD_in, RV_in, RG_in, RPI_in, EPS_W_in |
---|
356 | character (len=20) :: modname='lscp_ini_mod' |
---|
357 | character (len=80) :: abort_message |
---|
358 | |
---|
359 | |
---|
360 | lunout=lunout_in |
---|
361 | prt_level=prt_level_in |
---|
362 | fl_cor_ebil=fl_cor_ebil_in |
---|
363 | |
---|
364 | ok_ice_supersat=ok_ice_supersat_in |
---|
365 | ok_plane_contrail=ok_plane_contrail_in |
---|
366 | |
---|
367 | RG=RG_in |
---|
368 | RD=RD_in |
---|
369 | RV=RV_in |
---|
370 | RCPD=RCPD_in |
---|
371 | RLVTT=RLVTT_in |
---|
372 | RLSTT=RLSTT_in |
---|
373 | RLMLT=RLMLT_in |
---|
374 | RTT=RTT_in |
---|
375 | RV=RV_in |
---|
376 | RVTMP2=RVTMP2_in |
---|
377 | RPI=RPI_in |
---|
378 | EPS_W=EPS_W_in |
---|
379 | |
---|
380 | |
---|
381 | |
---|
382 | CALL getin_p('niter_lscp',niter_lscp) |
---|
383 | CALL getin_p('iflag_evap_prec',iflag_evap_prec) |
---|
384 | CALL getin_p('seuil_neb',seuil_neb) |
---|
385 | CALL getin_p('ok_radocond_snow',ok_radocond_snow) |
---|
386 | CALL getin_p('t_glace_max',t_glace_max) |
---|
387 | CALL getin_p('t_glace_min',t_glace_min) |
---|
388 | CALL getin_p('exposant_glace',exposant_glace) |
---|
389 | CALL getin_p('iflag_vice',iflag_vice) |
---|
390 | CALL getin_p('iflag_t_glace',iflag_t_glace) |
---|
391 | CALL getin_p('iflag_cloudth_vert',iflag_cloudth_vert) |
---|
392 | CALL getin_p('iflag_gammasat',iflag_gammasat) |
---|
393 | CALL getin_p('iflag_rain_incloud_vol',iflag_rain_incloud_vol) |
---|
394 | CALL getin_p('iflag_bergeron',iflag_bergeron) |
---|
395 | CALL getin_p('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat) |
---|
396 | CALL getin_p('iflag_pdf',iflag_pdf) |
---|
397 | CALL getin_p('iflag_icefrac',iflag_icefrac) |
---|
398 | CALL getin_p('reevap_ice',reevap_ice) |
---|
399 | CALL getin_p('cld_lc_lsc',cld_lc_lsc) |
---|
400 | CALL getin_p('cld_lc_con',cld_lc_con) |
---|
401 | CALL getin_p('cld_lc_lsc_snow',cld_lc_lsc_snow) |
---|
402 | CALL getin_p('cld_lc_con_snow',cld_lc_con_snow) |
---|
403 | CALL getin_p('cld_tau_lsc',cld_tau_lsc) |
---|
404 | CALL getin_p('cld_tau_con',cld_tau_con) |
---|
405 | CALL getin_p('cld_expo_lsc',cld_expo_lsc) |
---|
406 | CALL getin_p('cld_expo_con',cld_expo_con) |
---|
407 | CALL getin_p('ffallv_lsc',ffallv_lsc) |
---|
408 | CALL getin_p('ffallv_lsc',ffallv_con) |
---|
409 | CALL getin_p('coef_eva',coef_eva) |
---|
410 | coef_sub=coef_eva |
---|
411 | CALL getin_p('coef_eva_i',coef_sub) |
---|
412 | CALL getin_p('coef_sub',coef_sub) |
---|
413 | CALL getin_p('expo_eva',expo_eva) |
---|
414 | expo_sub=expo_eva |
---|
415 | CALL getin_p('expo_sub',expo_sub) |
---|
416 | CALL getin_p('iflag_autoconversion',iflag_autoconversion) |
---|
417 | CALL getin_p('dist_liq',dist_liq) |
---|
418 | CALL getin_p('tresh_cl',tresh_cl) |
---|
419 | CALL getin_p('capa_crystal',capa_crystal) |
---|
420 | CALL getin_p('naero5',naero5) |
---|
421 | CALL getin_p('gamma_snwretro',gamma_snwretro) |
---|
422 | CALL getin_p('gamma_taud',gamma_taud) |
---|
423 | CALL getin_p('iflag_oldbug_fisrtilp',iflag_oldbug_fisrtilp) |
---|
424 | CALL getin_p('temp_nowater',temp_nowater) |
---|
425 | CALL getin_p('ok_bug_phase_lscp',ok_bug_phase_lscp) |
---|
426 | ! for poprecip |
---|
427 | CALL getin_p('ok_poprecip',ok_poprecip) |
---|
428 | CALL getin_p('ok_corr_vap_evasub',ok_corr_vap_evasub) |
---|
429 | CALL getin_p('rain_int_min',rain_int_min) |
---|
430 | CALL getin_p('gamma_agg',gamma_agg) |
---|
431 | CALL getin_p('gamma_col',gamma_col) |
---|
432 | CALL getin_p('gamma_rim',gamma_rim) |
---|
433 | CALL getin_p('gamma_freez',gamma_freez) |
---|
434 | CALL getin_p('gamma_melt',gamma_melt) |
---|
435 | CALL getin_p('tau_auto_snow_max',tau_auto_snow_max) |
---|
436 | CALL getin_p('tau_auto_snow_min',tau_auto_snow_min) |
---|
437 | CALL getin_p('r_snow',r_snow) |
---|
438 | CALL getin_p('rain_fallspeed',rain_fallspeed) |
---|
439 | rain_fallspeed_clr=rain_fallspeed |
---|
440 | rain_fallspeed_cld=rain_fallspeed |
---|
441 | CALL getin_p('rain_fallspeed_clr',rain_fallspeed_clr) |
---|
442 | CALL getin_p('rain_fallspeed_cld',rain_fallspeed_cld) |
---|
443 | CALL getin_p('snow_fallspeed',snow_fallspeed) |
---|
444 | snow_fallspeed_clr=snow_fallspeed |
---|
445 | snow_fallspeed_cld=snow_fallspeed |
---|
446 | CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr) |
---|
447 | CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld) |
---|
448 | ! for condensation and ice supersaturation |
---|
449 | CALL getin_p('ok_unadjusted_clouds',ok_unadjusted_clouds) |
---|
450 | CALL getin_p('ok_weibull_warm_clouds',ok_weibull_warm_clouds) |
---|
451 | CALL getin_p('iflag_cloud_sublim_pdf',iflag_cloud_sublim_pdf) |
---|
452 | CALL getin_p('depo_coef_cirrus',depo_coef_cirrus) |
---|
453 | CALL getin_p('capa_cond_cirrus',capa_cond_cirrus) |
---|
454 | CALL getin_p('std_subl_pdf_lscp',std_subl_pdf_lscp) |
---|
455 | CALL getin_p('mu_subl_pdf_lscp',mu_subl_pdf_lscp) |
---|
456 | CALL getin_p('beta_pdf_lscp',beta_pdf_lscp) |
---|
457 | CALL getin_p('temp_thresh_pdf_lscp',temp_thresh_pdf_lscp) |
---|
458 | CALL getin_p('k0_pdf_lscp',k0_pdf_lscp) |
---|
459 | CALL getin_p('kappa_pdf_lscp',kappa_pdf_lscp) |
---|
460 | CALL getin_p('std100_pdf_lscp',std100_pdf_lscp) |
---|
461 | CALL getin_p('a_homofreez',a_homofreez) |
---|
462 | CALL getin_p('b_homofreez',b_homofreez) |
---|
463 | CALL getin_p('delta_hetfreez',delta_hetfreez) |
---|
464 | CALL getin_p('coef_mixing_lscp',coef_mixing_lscp) |
---|
465 | CALL getin_p('coef_shear_lscp',coef_shear_lscp) |
---|
466 | CALL getin_p('chi_mixing_lscp',chi_mixing_lscp) |
---|
467 | ! for aviation |
---|
468 | CALL getin_p('aspect_ratio_contrails',aspect_ratio_contrails) |
---|
469 | coef_mixing_contrails=coef_mixing_lscp |
---|
470 | CALL getin_p('coef_mixing_contrails',coef_mixing_contrails) |
---|
471 | coef_shear_contrails=coef_shear_lscp |
---|
472 | CALL getin_p('coef_shear_contrails',coef_shear_contrails) |
---|
473 | chi_mixing_contrails=chi_mixing_lscp |
---|
474 | CALL getin_p('chi_mixing_contrails',chi_mixing_contrails) |
---|
475 | CALL getin_p('rm_ice_crystals_contrails',rm_ice_crystals_contrails) |
---|
476 | CALL getin_p('EI_H2O_aviation',EI_H2O_aviation) |
---|
477 | CALL getin_p('qheat_fuel_aviation',qheat_fuel_aviation) |
---|
478 | CALL getin_p('prop_efficiency_aviation',prop_efficiency_aviation) |
---|
479 | CALL getin_p('linear_contrails_lifetime',linear_contrails_lifetime) |
---|
480 | CALL getin_p('initial_width_contrails',initial_width_contrails) |
---|
481 | |
---|
482 | |
---|
483 | |
---|
484 | WRITE(lunout,*) 'lscp_ini, niter_lscp:', niter_lscp |
---|
485 | WRITE(lunout,*) 'lscp_ini, iflag_evap_prec:', iflag_evap_prec |
---|
486 | WRITE(lunout,*) 'lscp_ini, seuil_neb:', seuil_neb |
---|
487 | WRITE(lunout,*) 'lscp_ini, ok_radocond_snow:', ok_radocond_snow |
---|
488 | WRITE(lunout,*) 'lscp_ini, t_glace_max:', t_glace_max |
---|
489 | WRITE(lunout,*) 'lscp_ini, t_glace_min:', t_glace_min |
---|
490 | WRITE(lunout,*) 'lscp_ini, exposant_glace:', exposant_glace |
---|
491 | WRITE(lunout,*) 'lscp_ini, iflag_vice:', iflag_vice |
---|
492 | WRITE(lunout,*) 'lscp_ini, iflag_t_glace:', iflag_t_glace |
---|
493 | WRITE(lunout,*) 'lscp_ini, iflag_cloudth_vert:', iflag_cloudth_vert |
---|
494 | WRITE(lunout,*) 'lscp_ini, iflag_gammasat:', iflag_gammasat |
---|
495 | WRITE(lunout,*) 'lscp_ini, iflag_rain_incloud_vol:', iflag_rain_incloud_vol |
---|
496 | WRITE(lunout,*) 'lscp_ini, iflag_bergeron:', iflag_bergeron |
---|
497 | WRITE(lunout,*) 'lscp_ini, iflag_fisrtilp_qsat:', iflag_fisrtilp_qsat |
---|
498 | WRITE(lunout,*) 'lscp_ini, iflag_pdf', iflag_pdf |
---|
499 | WRITE(lunout,*) 'lscp_ini, iflag_icefrac', iflag_icefrac |
---|
500 | WRITE(lunout,*) 'lscp_ini, reevap_ice', reevap_ice |
---|
501 | WRITE(lunout,*) 'lscp_ini, cld_lc_lsc', cld_lc_lsc |
---|
502 | WRITE(lunout,*) 'lscp_ini, cld_lc_con', cld_lc_con |
---|
503 | WRITE(lunout,*) 'lscp_ini, cld_lc_lsc_snow', cld_lc_lsc_snow |
---|
504 | WRITE(lunout,*) 'lscp_ini, cld_lc_con_snow', cld_lc_con_snow |
---|
505 | WRITE(lunout,*) 'lscp_ini, cld_tau_lsc', cld_tau_lsc |
---|
506 | WRITE(lunout,*) 'lscp_ini, cld_tau_con', cld_tau_con |
---|
507 | WRITE(lunout,*) 'lscp_ini, cld_expo_lsc', cld_expo_lsc |
---|
508 | WRITE(lunout,*) 'lscp_ini, cld_expo_con', cld_expo_con |
---|
509 | WRITE(lunout,*) 'lscp_ini, ffallv_lsc', ffallv_lsc |
---|
510 | WRITE(lunout,*) 'lscp_ini, ffallv_con', ffallv_con |
---|
511 | WRITE(lunout,*) 'lscp_ini, coef_eva', coef_eva |
---|
512 | WRITE(lunout,*) 'lscp_ini, coef_sub', coef_sub |
---|
513 | WRITE(lunout,*) 'lscp_ini, expo_eva', expo_eva |
---|
514 | WRITE(lunout,*) 'lscp_ini, expo_sub', expo_sub |
---|
515 | WRITE(lunout,*) 'lscp_ini, iflag_autoconversion', iflag_autoconversion |
---|
516 | WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq |
---|
517 | WRITE(lunout,*) 'lscp_ini, tresh_cl', tresh_cl |
---|
518 | WRITE(lunout,*) 'lscp_ini, capa_crystal', capa_crystal |
---|
519 | WRITE(lunout,*) 'lscp_ini, naero5', naero5 |
---|
520 | WRITE(lunout,*) 'lscp_ini, gamma_snwretro', gamma_snwretro |
---|
521 | WRITE(lunout,*) 'lscp_ini, gamma_taud', gamma_taud |
---|
522 | WRITE(lunout,*) 'lscp_ini, iflag_oldbug_fisrtilp', iflag_oldbug_fisrtilp |
---|
523 | WRITE(lunout,*) 'lscp_ini, fl_cor_ebil', fl_cor_ebil |
---|
524 | WRITE(lunout,*) 'lscp_ini, temp_nowater', temp_nowater |
---|
525 | WRITE(lunout,*) 'lscp_ini, ok_bug_phase_lscp', ok_bug_phase_lscp |
---|
526 | ! for poprecip |
---|
527 | WRITE(lunout,*) 'lscp_ini, ok_poprecip', ok_poprecip |
---|
528 | WRITE(lunout,*) 'lscp_ini, ok_corr_vap_evasub', ok_corr_vap_evasub |
---|
529 | WRITE(lunout,*) 'lscp_ini, rain_int_min:', rain_int_min |
---|
530 | WRITE(lunout,*) 'lscp_ini, gamma_agg:', gamma_agg |
---|
531 | WRITE(lunout,*) 'lscp_ini, gamma_col:', gamma_col |
---|
532 | WRITE(lunout,*) 'lscp_ini, gamma_rim:', gamma_rim |
---|
533 | WRITE(lunout,*) 'lscp_ini, gamma_freez:', gamma_freez |
---|
534 | WRITE(lunout,*) 'lscp_ini, gamma_melt:', gamma_melt |
---|
535 | WRITE(lunout,*) 'lscp_ini, tau_auto_snow_max:',tau_auto_snow_max |
---|
536 | WRITE(lunout,*) 'lscp_ini, tau_auto_snow_min:',tau_auto_snow_min |
---|
537 | WRITE(lunout,*) 'lscp_ini, r_snow:', r_snow |
---|
538 | WRITE(lunout,*) 'lscp_ini, rain_fallspeed_clr:', rain_fallspeed_clr |
---|
539 | WRITE(lunout,*) 'lscp_ini, rain_fallspeed_cld:', rain_fallspeed_cld |
---|
540 | WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr |
---|
541 | WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld |
---|
542 | ! for condensation and ice supersaturation |
---|
543 | WRITE(lunout,*) 'lscp_ini, ok_ice_supersat:', ok_ice_supersat |
---|
544 | WRITE(lunout,*) 'lscp_ini, ok_unadjusted_clouds:', ok_unadjusted_clouds |
---|
545 | WRITE(lunout,*) 'lscp_ini, ok_weibull_warm_clouds:', ok_weibull_warm_clouds |
---|
546 | WRITE(lunout,*) 'lscp_ini, iflag_cloud_sublim_pdf:', iflag_cloud_sublim_pdf |
---|
547 | WRITE(lunout,*) 'lscp_ini, depo_coef_cirrus:', depo_coef_cirrus |
---|
548 | WRITE(lunout,*) 'lscp_ini, capa_cond_cirrus:', capa_cond_cirrus |
---|
549 | WRITE(lunout,*) 'lscp_ini, std_subl_pdf_lscp:', std_subl_pdf_lscp |
---|
550 | WRITE(lunout,*) 'lscp_ini, mu_subl_pdf_lscp:', mu_subl_pdf_lscp |
---|
551 | WRITE(lunout,*) 'lscp_ini, beta_pdf_lscp:', beta_pdf_lscp |
---|
552 | WRITE(lunout,*) 'lscp_ini, temp_thresh_pdf_lscp:', temp_thresh_pdf_lscp |
---|
553 | WRITE(lunout,*) 'lscp_ini, k0_pdf_lscp:', k0_pdf_lscp |
---|
554 | WRITE(lunout,*) 'lscp_ini, kappa_pdf_lscp:', kappa_pdf_lscp |
---|
555 | WRITE(lunout,*) 'lscp_ini, std100_pdf_lscp:', std100_pdf_lscp |
---|
556 | WRITE(lunout,*) 'lscp_ini, a_homofreez:', a_homofreez |
---|
557 | WRITE(lunout,*) 'lscp_ini, b_homofreez:', b_homofreez |
---|
558 | WRITE(lunout,*) 'lscp_ini, delta_hetfreez', delta_hetfreez |
---|
559 | WRITE(lunout,*) 'lscp_ini, coef_mixing_lscp:', coef_mixing_lscp |
---|
560 | WRITE(lunout,*) 'lscp_ini, coef_shear_lscp:', coef_shear_lscp |
---|
561 | WRITE(lunout,*) 'lscp_ini, chi_mixing_lscp:', chi_mixing_lscp |
---|
562 | ! for aviation |
---|
563 | WRITE(lunout,*) 'lscp_ini, aspect_ratio_contrails:', aspect_ratio_contrails |
---|
564 | WRITE(lunout,*) 'lscp_ini, coef_mixing_contrails:', coef_mixing_contrails |
---|
565 | WRITE(lunout,*) 'lscp_ini, coef_shear_contrails:', coef_shear_contrails |
---|
566 | WRITE(lunout,*) 'lscp_ini, chi_mixing_contrails:', chi_mixing_contrails |
---|
567 | WRITE(lunout,*) 'lscp_ini, rm_ice_crystals_contrails:', rm_ice_crystals_contrails |
---|
568 | WRITE(lunout,*) 'lscp_ini, EI_H2O_aviation:', EI_H2O_aviation |
---|
569 | WRITE(lunout,*) 'lscp_ini, qheat_fuel_aviation:', qheat_fuel_aviation |
---|
570 | WRITE(lunout,*) 'lscp_ini, prop_efficiency_aviation:', prop_efficiency_aviation |
---|
571 | WRITE(lunout,*) 'lscp_ini, linear_contrails_lifetime:', linear_contrails_lifetime |
---|
572 | WRITE(lunout,*) 'lscp_ini, initial_width_contrails:', initial_width_contrails |
---|
573 | |
---|
574 | |
---|
575 | |
---|
576 | ! check for precipitation sub-time steps |
---|
577 | IF (ABS(dtime/REAL(niter_lscp)-360.0).GT.0.001) THEN |
---|
578 | WRITE(lunout,*) 'lscp: it is not expected, see Z.X.Li', dtime |
---|
579 | WRITE(lunout,*) 'I would prefer a 6 min sub-timestep' |
---|
580 | ENDIF |
---|
581 | |
---|
582 | ! check consistency between numerical resolution of autoconversion |
---|
583 | ! and other options |
---|
584 | |
---|
585 | IF (iflag_autoconversion .EQ. 2) THEN |
---|
586 | IF ((iflag_vice .NE. 0) .OR. (niter_lscp .GT. 1)) THEN |
---|
587 | abort_message = 'in lscp, iflag_autoconversion=2 requires iflag_vice=0 and niter_lscp=1' |
---|
588 | CALL abort_physic (modname,abort_message,1) |
---|
589 | ENDIF |
---|
590 | ENDIF |
---|
591 | |
---|
592 | |
---|
593 | IF ( ok_weibull_warm_clouds .AND. .NOT. ok_ice_supersat ) THEN |
---|
594 | abort_message = 'in lscp, ok_weibull_warm_clouds=y needs ok_ice_supersat=y' |
---|
595 | CALL abort_physic (modname,abort_message,1) |
---|
596 | ENDIF |
---|
597 | |
---|
598 | IF ( ok_unadjusted_clouds .AND. .NOT. ok_ice_supersat ) THEN |
---|
599 | abort_message = 'in lscp, ok_unadjusted_clouds=y needs ok_ice_supersat=y' |
---|
600 | CALL abort_physic (modname,abort_message,1) |
---|
601 | ENDIF |
---|
602 | |
---|
603 | |
---|
604 | !AA Temporary initialisation |
---|
605 | a_tr_sca(1) = -0.5 |
---|
606 | a_tr_sca(2) = -0.5 |
---|
607 | a_tr_sca(3) = -0.5 |
---|
608 | a_tr_sca(4) = -0.5 |
---|
609 | |
---|
610 | CALL cloudth_ini(iflag_cloudth_vert,iflag_ratqs) |
---|
611 | |
---|
612 | RETURN |
---|
613 | |
---|
614 | END SUBROUTINE lscp_ini |
---|
615 | |
---|
616 | END MODULE lmdz_lscp_ini |
---|