source: LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90 @ 4869

Last change on this file since 4869 was 4869, checked in by evignon, 2 months ago

identification d'un bug dans l'ajustement de la temperature
lors de la fonte dans lscp
Ajout d'un flag pour maintenir/enlever le probleme
File size: 18.1 KB
Line 
1MODULE lmdz_lscp_ini
2
3IMPLICIT NONE
4
5  ! PARAMETERS for lscp:
6  !--------------------
7 
8  REAL RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RPI
9  !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RPI)
10 
11  LOGICAL, SAVE, PROTECTED :: ok_bug_fonte_lscp=.true. ! maintains bug on temperature adjustment associated with snow melting
12  !$OMP THREADPRIVATE(ok_bug_fonte_lscp)
13
14  REAL, SAVE, PROTECTED :: seuil_neb=0.001      ! cloud fraction threshold: a cloud really exists when exceeded
15  !$OMP THREADPRIVATE(seuil_neb)
16
17  INTEGER, SAVE :: lunout, prt_level            ! Logical unit number and level for standard output
18  !$OMP THREADPRIVATE(lunout,prt_level)
19
20  INTEGER, SAVE, PROTECTED :: niter_lscp=5      ! number of iterations to calculate autoconversion to precipitation
21  !$OMP THREADPRIVATE(niter_lscp)
22
23  INTEGER, SAVE, PROTECTED :: iflag_evap_prec=1 ! precipitation evaporation flag. 0: nothing, 1: "old way",
24                                                ! 2: Max cloud fraction above to calculate the max of reevaporation
25                                                ! >=4: LTP'method i.e. evaporation in the clear-sky fraction of the mesh only
26                                                ! pay attention that iflag_evap_prec=4 may lead to unrealistic and overestimated
27                                                ! evaporation. Use 5 instead
28  !$OMP THREADPRIVATE(iflag_evap_prec)
29
30  REAL, SAVE, PROTECTED :: t_coup=234.0         ! temperature threshold which determines the phase
31                                                ! for which the saturation vapor pressure is calculated
32  !$OMP THREADPRIVATE(t_coup)
33  REAL, SAVE, PROTECTED :: DDT0=0.01            ! iteration parameter
34  !$OMP THREADPRIVATE(DDT0)
35
36  REAL, SAVE, PROTECTED :: ztfondue=278.15      ! parameter to calculate melting fraction of precipitation
37  !$OMP THREADPRIVATE(ztfondue)
38
39  REAL, SAVE, PROTECTED :: temp_nowater=233.15  ! temperature below which liquid water no longer exists
40  !$OMP THREADPRIVATE(temp_nowater)
41
42  REAL, SAVE, PROTECTED :: a_tr_sca(4)          ! Variables for tracers temporary: alpha parameter for scavenging, 4 possible scavenging processes
43  !$OMP THREADPRIVATE(a_tr_sca)
44 
45  INTEGER, SAVE, PROTECTED ::  iflag_mpc_bl=0   ! flag to activate boundary layer mixed phase cloud param
46  !$OMP THREADPRIVATE(iflag_mpc_bl)
47 
48  LOGICAL, SAVE, PROTECTED :: ok_radocond_snow=.false. ! take into account the mass of ice precip in the cloud ice content seen by radiation
49  !$OMP THREADPRIVATE(ok_radocond_snow)
50
51  REAL, SAVE, PROTECTED :: t_glace_min=258.0    ! lower-bound temperature parameter for cloud phase determination
52  !$OMP THREADPRIVATE(t_glace_min)
53
54  REAL, SAVE, PROTECTED :: t_glace_max=273.15   ! upper-bound temperature parameter for cloud phase determination
55  !$OMP THREADPRIVATE(t_glace_max)
56
57  REAL, SAVE, PROTECTED :: exposant_glace=1.0   ! parameter for cloud phase determination
58  !$OMP THREADPRIVATE(exposant_glace)
59
60  INTEGER, SAVE, PROTECTED :: iflag_vice=0      ! which expression for ice crystall fall velocity
61  !$OMP THREADPRIVATE(iflag_vice)
62
63  INTEGER, SAVE, PROTECTED :: iflag_t_glace=0   ! which expression for cloud phase partitioning
64  !$OMP THREADPRIVATE(iflag_t_glace)
65
66  INTEGER, SAVE, PROTECTED :: iflag_cloudth_vert=0         ! option for determining cloud fraction and content in convective boundary layers
67  !$OMP THREADPRIVATE(iflag_cloudth_vert)
68
69  INTEGER, SAVE, PROTECTED :: iflag_gammasat=0             ! which threshold for homogeneous nucleation below -40oC
70  !$OMP THREADPRIVATE(iflag_gammasat)
71
72  INTEGER, SAVE, PROTECTED :: iflag_rain_incloud_vol=0     ! use of volume cloud fraction for rain autoconversion
73  !$OMP THREADPRIVATE(iflag_rain_incloud_vol)
74
75  INTEGER, SAVE, PROTECTED :: iflag_bergeron=0             ! bergeron effect for liquid precipitation treatment 
76  !$OMP THREADPRIVATE(iflag_bergeron)
77
78  INTEGER, SAVE, PROTECTED :: iflag_fisrtilp_qsat=0        ! qsat adjustment (iterative) during autoconversion
79  !$OMP THREADPRIVATE(iflag_fisrtilp_qsat)
80
81  INTEGER, SAVE, PROTECTED :: iflag_pdf=0                  ! type of subgrid scale qtot pdf
82  !$OMP THREADPRIVATE(iflag_pdf)
83
84  INTEGER, SAVE, PROTECTED :: iflag_autoconversion=0       ! autoconversion option
85  !$OMP THREADPRIVATE(iflag_autoconversion)
86
87  LOGICAL, SAVE, PROTECTED :: reevap_ice=.false.           ! no liquid precip for T< threshold
88  !$OMP THREADPRIVATE(reevap_ice)
89
90  REAL, SAVE, PROTECTED :: cld_lc_lsc=2.6e-4               ! liquid autoconversion coefficient, stratiform rain
91  !$OMP THREADPRIVATE(cld_lc_lsc)
92
93  REAL, SAVE, PROTECTED :: cld_lc_con=2.6e-4                ! liquid autoconversion coefficient, convective rain
94  !$OMP THREADPRIVATE(cld_lc_con)
95
96  REAL, SAVE, PROTECTED :: cld_tau_lsc=3600.                ! liquid autoconversion timescale, stratiform rain
97  !$OMP THREADPRIVATE(cld_tau_lsc)
98
99  REAL, SAVE, PROTECTED :: cld_tau_con=3600.                ! liquid autoconversion timescale, convective rain
100  !$OMP THREADPRIVATE(cld_tau_con)
101
102  REAL, SAVE, PROTECTED :: cld_expo_lsc=2.                  ! liquid autoconversion threshold exponent, stratiform rain
103  !$OMP THREADPRIVATE(cld_expo_lsc)
104
105  REAL, SAVE, PROTECTED :: cld_expo_con=2.                  ! liquid autoconversion threshold exponent, convective rain
106  !$OMP THREADPRIVATE(cld_expo_con)
107
108  REAL, SAVE, PROTECTED :: ffallv_lsc=1.                    ! tuning coefficient crystal fall velocity, stratiform
109  !$OMP THREADPRIVATE(ffallv_lsc)
110
111  REAL, SAVE, PROTECTED :: ffallv_con=1.                    ! tuning coefficient crystal fall velocity, convective
112  !$OMP THREADPRIVATE(ffallv_con)
113
114  REAL, SAVE, PROTECTED :: coef_eva=2e-5                    ! tuning coefficient liquid precip evaporation
115  !$OMP THREADPRIVATE(coef_eva)
116
117  REAL, SAVE, PROTECTED :: coef_sub                        ! tuning coefficient ice precip sublimation
118  !$OMP THREADPRIVATE(coef_sub)
119
120  REAL, SAVE, PROTECTED :: expo_eva=0.5                     ! tuning coefficient liquid precip evaporation
121  !$OMP THREADPRIVATE(expo_eva)
122
123  REAL, SAVE, PROTECTED :: expo_sub                       ! tuning coefficient ice precip sublimation
124  !$OMP THREADPRIVATE(expo_sub)
125
126  REAL, SAVE, PROTECTED :: cice_velo=1.645                  ! factor in the ice fall velocity formulation
127  !$OMP THREADPRIVATE(cice_velo)
128
129  REAL, SAVE, PROTECTED ::  dice_velo=0.16                  ! exponent in the ice fall velocity formulation
130  !$OMP THREADPRIVATE(dice_velo)
131
132  REAL, SAVE, PROTECTED :: dist_liq=300.                    ! typical deph of cloud-top liquid layer in mpcs
133  !$OMP THREADPRIVATE(dist_liq)
134
135  REAL, SAVE, PROTECTED    :: tresh_cl=0.0                  ! cloud fraction threshold for cloud top search
136  !$OMP THREADPRIVATE(tresh_cl)
137
138  !--Parameters for poprecip
139  LOGICAL, SAVE, PROTECTED :: ok_poprecip=.FALSE.           ! use the processes-oriented formulation of precipitations
140  !$OMP THREADPRIVATE(ok_poprecip)
141
142  REAL, SAVE, PROTECTED :: cld_lc_lsc_snow               ! snow autoconversion coefficient, stratiform rain
143  !$OMP THREADPRIVATE(cld_lc_lsc_snow)
144
145  REAL, SAVE, PROTECTED :: cld_lc_con_snow                ! snow autoconversion coefficient, convective rain
146  !$OMP THREADPRIVATE(cld_lc_con_snow)
147
148  REAL, SAVE, PROTECTED :: rain_int_min=0.001               ! Minimum local rain intensity [mm/s] before the decrease in associated precipitation fraction
149  !$OMP THREADPRIVATE(rain_int_min)
150
151  REAL, SAVE, PROTECTED :: thresh_precip_frac=1.E-6         ! precipitation fraction threshold TODO [-]
152  !$OMP THREADPRIVATE(thresh_precip_frac)
153
154  REAL, SAVE, PROTECTED :: gamma_col=1.                     ! A COMMENTER TODO [-]
155  !$OMP THREADPRIVATE(gamma_col)
156
157  REAL, SAVE, PROTECTED :: gamma_agg=1.                     ! A COMMENTER TODO [-]
158  !$OMP THREADPRIVATE(gamma_agg)
159
160  REAL, SAVE, PROTECTED :: gamma_rim=1.                     ! A COMMENTER TODO [-]
161  !$OMP THREADPRIVATE(gamma_rim)
162
163  REAL, SAVE, PROTECTED :: rho_rain=1000.                    ! A COMMENTER TODO [kg/m3]
164  !$OMP THREADPRIVATE(rho_rain)
165
166  REAL, SAVE, PROTECTED :: rho_ice=920.                    ! A COMMENTER TODO [kg/m3]
167  !$OMP THREADPRIVATE(rho_ice)
168
169  REAL, SAVE, PROTECTED :: rho_snow                        ! A COMMENTER TODO [kg/m3]
170  !$OMP THREADPRIVATE(rho_snow)
171
172  REAL, SAVE, PROTECTED :: r_rain=500.E-6                   ! A COMMENTER TODO [m]
173  !$OMP THREADPRIVATE(r_rain)
174
175  REAL, SAVE, PROTECTED :: r_snow=1.E-3                    ! A COMMENTER TODO [m]
176  !$OMP THREADPRIVATE(r_snow)
177
178  REAL, SAVE, PROTECTED :: r_ice=50.E-6                    ! A COMMENTER TODO [m]
179  !$OMP THREADPRIVATE(r_ice)
180
181  REAL, SAVE, PROTECTED :: tau_auto_snow_min=100.          ! A COMMENTER TODO [s]
182  !$OMP THREADPRIVATE(tau_auto_snow_min)
183
184  REAL, SAVE, PROTECTED :: tau_auto_snow_max=1000.          ! A COMMENTER TODO [s]
185  !$OMP THREADPRIVATE(tau_auto_snow_max)
186
187  REAL, SAVE, PROTECTED :: eps=1.E-10                       ! A COMMENTER TODO [-]
188  !$OMP THREADPRIVATE(eps)
189
190  REAL, SAVE, PROTECTED :: air_thermal_conduct=2.4e-2      ! A COMMENTER TODO [-]
191  !$OMP THREADPRIVATE(air_thermal_conduct)
192
193  REAL, SAVE, PROTECTED :: coef_ventil=1.                   ! A COMMENTER TODO [-]
194  !$OMP THREADPRIVATE(coef_ventil)
195
196  REAL, SAVE, PROTECTED :: alpha_freez=4.                 ! A COMMENTER TODO [-]
197  !$OMP THREADPRIVATE(alpha_freez)
198
199  REAL, SAVE, PROTECTED :: beta_freez=0.1                 ! A COMMENTER TODO [m-3.s-1]
200  !$OMP THREADPRIVATE(beta_freez)
201
202  REAL, SAVE, PROTECTED :: gamma_freez=1.                 ! A COMMENTER TODO [-]
203  !$OMP THREADPRIVATE(gamma_freez)
204
205  REAL, SAVE, PROTECTED :: rain_fallspeed=4.              ! A COMMENTER TODO [m/s]
206  !$OMP THREADPRIVATE(rain_fallspeed)
207
208  REAL, SAVE, PROTECTED :: rain_fallspeed_clr              ! A COMMENTER TODO [m/s]
209  !$OMP THREADPRIVATE(rain_fallspeed_clr)
210
211  REAL, SAVE, PROTECTED :: rain_fallspeed_cld             ! A COMMENTER TODO [m/s]
212  !$OMP THREADPRIVATE(rain_fallspeed_cld)
213
214  REAL, SAVE, PROTECTED :: snow_fallspeed=1.             ! A COMMENTER TODO [m/s]
215  !$OMP THREADPRIVATE(snow_fallspeed)
216
217  REAL, SAVE, PROTECTED :: snow_fallspeed_clr             ! A COMMENTER TODO [m/s]
218  !$OMP THREADPRIVATE(snow_fallspeed_clr)
219
220  REAL, SAVE, PROTECTED :: snow_fallspeed_cld             ! A COMMENTER TODO [m/s]
221  !$OMP THREADPRIVATE(snow_fallspeed_cld)
222  !--End of the parameters for poprecip
223
224! Two parameters used for lmdz_lscp_old only
225  INTEGER, SAVE, PROTECTED :: iflag_oldbug_fisrtilp=0, fl_cor_ebil
226  !$OMP THREADPRIVATE(iflag_oldbug_fisrtilp,fl_cor_ebil)
227
228CONTAINS
229
230SUBROUTINE lscp_ini(dtime,lunout_in,prt_level_in,ok_ice_sursat, iflag_ratqs, fl_cor_ebil_in, RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, &
231                    RVTMP2_in, RTT_in,RD_in,RG_in,RPI_in)
232
233
234   USE ioipsl_getin_p_mod, ONLY : getin_p
235   USE ice_sursat_mod, ONLY: ice_sursat_init
236   USE lmdz_cloudth_ini, ONLY : cloudth_ini
237
238   REAL, INTENT(IN)      :: dtime
239   INTEGER, INTENT(IN)   :: lunout_in,prt_level_in,iflag_ratqs,fl_cor_ebil_in
240   LOGICAL, INTENT(IN)   :: ok_ice_sursat
241
242   REAL, INTENT(IN)      :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in
243   REAL, INTENT(IN)      ::  RVTMP2_in, RTT_in, RD_in, RG_in, RPI_in
244   character (len=20) :: modname='lscp_ini_mod'
245   character (len=80) :: abort_message
246
247
248    lunout=lunout_in
249    prt_level=prt_level_in
250    fl_cor_ebil=fl_cor_ebil_in
251
252    RG=RG_in
253    RD=RD_in
254    RCPD=RCPD_in
255    RLVTT=RLVTT_in
256    RLSTT=RLSTT_in
257    RLMLT=RLMLT_in
258    RTT=RTT_in
259    RG=RG_in
260    RVTMP2=RVTMP2_in
261    RPI=RPI_in
262
263
264
265    CALL getin_p('niter_lscp',niter_lscp)
266    CALL getin_p('iflag_evap_prec',iflag_evap_prec)
267    CALL getin_p('seuil_neb',seuil_neb)
268    CALL getin_p('iflag_mpc_bl',iflag_mpc_bl)
269    CALL getin_p('ok_radocond_snow',ok_radocond_snow)
270    CALL getin_p('t_glace_max',t_glace_max)
271    CALL getin_p('t_glace_min',t_glace_min)
272    CALL getin_p('exposant_glace',exposant_glace)
273    CALL getin_p('iflag_vice',iflag_vice)
274    CALL getin_p('iflag_t_glace',iflag_t_glace)
275    CALL getin_p('iflag_cloudth_vert',iflag_cloudth_vert)
276    CALL getin_p('iflag_gammasat',iflag_gammasat)
277    CALL getin_p('iflag_rain_incloud_vol',iflag_rain_incloud_vol)
278    CALL getin_p('iflag_bergeron',iflag_bergeron)
279    CALL getin_p('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat)
280    CALL getin_p('iflag_pdf',iflag_pdf)
281    CALL getin_p('reevap_ice',reevap_ice)
282    CALL getin_p('cld_lc_lsc',cld_lc_lsc)
283    CALL getin_p('cld_lc_con',cld_lc_con)
284    cld_lc_lsc_snow=cld_lc_lsc
285    cld_lc_con_snow=cld_lc_con
286    CALL getin_p('cld_lc_lsc_snow',cld_lc_lsc_snow)
287    CALL getin_p('cld_lc_con_snow',cld_lc_con_snow)
288    CALL getin_p('cld_tau_lsc',cld_tau_lsc)
289    CALL getin_p('cld_tau_con',cld_tau_con)
290    CALL getin_p('cld_expo_lsc',cld_expo_lsc)
291    CALL getin_p('cld_expo_con',cld_expo_con)
292    CALL getin_p('ffallv_lsc',ffallv_lsc)
293    CALL getin_p('ffallv_lsc',ffallv_con)
294    CALL getin_p('coef_eva',coef_eva)
295    coef_sub=coef_eva
296    CALL getin_p('coef_eva_i',coef_sub)
297    CALL getin_p('coef_sub',coef_sub)
298    CALL getin_p('expo_eva',expo_eva)
299    expo_sub=expo_eva
300    CALL getin_p('expo_sub',expo_sub)
301    CALL getin_p('iflag_autoconversion',iflag_autoconversion)
302    CALL getin_p('dist_liq',dist_liq)
303    CALL getin_p('tresh_cl',tresh_cl)
304    CALL getin_p('iflag_oldbug_fisrtilp',iflag_oldbug_fisrtilp)
305    CALL getin_p('ok_poprecip',ok_poprecip)
306    CALL getin_p('rain_int_min',rain_int_min)
307    CALL getin_p('gamma_agg',gamma_agg)
308    CALL getin_p('gamma_col',gamma_col)
309    CALL getin_p('gamma_rim',gamma_rim)
310    CALL getin_p('gamma_freez',gamma_freez)
311    CALL getin_p('r_snow',r_snow)
312    CALL getin_p('rain_fallspeed',rain_fallspeed)
313    rain_fallspeed_clr=rain_fallspeed
314    rain_fallspeed_cld=rain_fallspeed
315    CALL getin_p('rain_fallspeed_clr',rain_fallspeed_clr)
316    CALL getin_p('rain_fallspeed_cld',rain_fallspeed_cld)
317    CALL getin_p('snow_fallspeed',snow_fallspeed)
318    snow_fallspeed_clr=snow_fallspeed
319    snow_fallspeed_cld=snow_fallspeed
320    CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr)
321    CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld)
322    CALL getin_p('ok_bug_fonte_lscp',ok_bug_fonte_lscp)
323
324
325
326    WRITE(lunout,*) 'lscp_ini, niter_lscp:', niter_lscp
327    WRITE(lunout,*) 'lscp_ini, iflag_evap_prec:', iflag_evap_prec
328    WRITE(lunout,*) 'lscp_ini, seuil_neb:', seuil_neb
329    WRITE(lunout,*) 'lscp_ini, iflag_mpc_bl:', iflag_mpc_bl
330    WRITE(lunout,*) 'lscp_ini, ok_radocond_snow:', ok_radocond_snow
331    WRITE(lunout,*) 'lscp_ini, t_glace_max:', t_glace_max
332    WRITE(lunout,*) 'lscp_ini, t_glace_min:', t_glace_min
333    WRITE(lunout,*) 'lscp_ini, exposant_glace:', exposant_glace
334    WRITE(lunout,*) 'lscp_ini, iflag_vice:', iflag_vice
335    WRITE(lunout,*) 'lscp_ini, iflag_t_glace:', iflag_t_glace
336    WRITE(lunout,*) 'lscp_ini, iflag_cloudth_vert:', iflag_cloudth_vert
337    WRITE(lunout,*) 'lscp_ini, iflag_gammasat:', iflag_gammasat
338    WRITE(lunout,*) 'lscp_ini, iflag_rain_incloud_vol:', iflag_rain_incloud_vol
339    WRITE(lunout,*) 'lscp_ini, iflag_bergeron:', iflag_bergeron
340    WRITE(lunout,*) 'lscp_ini, iflag_fisrtilp_qsat:', iflag_fisrtilp_qsat
341    WRITE(lunout,*) 'lscp_ini, iflag_pdf', iflag_pdf
342    WRITE(lunout,*) 'lscp_ini, reevap_ice', reevap_ice
343    WRITE(lunout,*) 'lscp_ini, cld_lc_lsc', cld_lc_lsc
344    WRITE(lunout,*) 'lscp_ini, cld_lc_con', cld_lc_con
345    WRITE(lunout,*) 'lscp_ini, cld_lc_lsc_snow', cld_lc_lsc_snow
346    WRITE(lunout,*) 'lscp_ini, cld_lc_con_snow', cld_lc_con_snow
347    WRITE(lunout,*) 'lscp_ini, cld_tau_lsc', cld_tau_lsc
348    WRITE(lunout,*) 'lscp_ini, cld_tau_con', cld_tau_con
349    WRITE(lunout,*) 'lscp_ini, cld_expo_lsc', cld_expo_lsc
350    WRITE(lunout,*) 'lscp_ini, cld_expo_con', cld_expo_con
351    WRITE(lunout,*) 'lscp_ini, ffallv_lsc', ffallv_lsc
352    WRITE(lunout,*) 'lscp_ini, ffallv_con', ffallv_con
353    WRITE(lunout,*) 'lscp_ini, coef_eva', coef_eva
354    WRITE(lunout,*) 'lscp_ini, coef_sub', coef_sub
355    WRITE(lunout,*) 'lscp_ini, expo_eva', expo_eva
356    WRITE(lunout,*) 'lscp_ini, expo_sub', expo_sub
357    WRITE(lunout,*) 'lscp_ini, iflag_autoconversion', iflag_autoconversion
358    WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq
359    WRITE(lunout,*) 'lscp_ini, tresh_cl', tresh_cl
360    WRITE(lunout,*) 'lscp_ini, iflag_oldbug_fisrtilp', iflag_oldbug_fisrtilp
361    WRITE(lunout,*) 'lscp_ini, fl_cor_ebil', fl_cor_ebil
362    WRITE(lunout,*) 'lscp_ini, ok_poprecip', ok_poprecip
363    WRITE(lunout,*) 'lscp_ini, rain_int_min:', rain_int_min
364    WRITE(lunout,*) 'lscp_ini, gamma_agg:', gamma_agg
365    WRITE(lunout,*) 'lscp_ini, gamma_col:', gamma_col
366    WRITE(lunout,*) 'lscp_ini, gamma_rim:', gamma_rim
367    WRITE(lunout,*) 'lscp_ini, gamma_freez:', gamma_freez
368    WRITE(lunout,*) 'lscp_ini, r_snow:', r_snow
369    WRITE(lunout,*) 'lscp_ini, rain_fallspeed_clr:', rain_fallspeed_clr
370    WRITE(lunout,*) 'lscp_ini, rain_fallspeed_cld:', rain_fallspeed_cld
371    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr
372    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld
373    WRITE(lunout,*) 'lscp_ini, ok_bug_fonte_lscp:', ok_bug_fonte_lscp
374
375
376
377
378
379    ! check for precipitation sub-time steps
380    IF (ABS(dtime/REAL(niter_lscp)-360.0).GT.0.001) THEN
381        WRITE(lunout,*) 'lscp: it is not expected, see Z.X.Li', dtime
382        WRITE(lunout,*) 'I would prefer a 6 min sub-timestep'
383    ENDIF
384
385    ! check consistency between numerical resolution of autoconversion
386    ! and other options
387   
388    IF (iflag_autoconversion .EQ. 2) THEN
389        IF ((iflag_vice .NE. 0) .OR. (niter_lscp .GT. 1)) THEN
390           abort_message = 'in lscp, iflag_autoconversion=2 requires iflag_vice=0 and niter_lscp=1'
391           CALL abort_physic (modname,abort_message,1)
392        ENDIF
393    ENDIF
394
395
396    ! check consistency between the use of the processes-oriented precipitation formulations
397    ! and other options
398   
399    IF (ok_poprecip) THEN
400        IF ((iflag_evap_prec .LT. 4) .OR. (niter_lscp .GT. 1)) THEN
401           abort_message = 'in lscp, ok_poprecip=y requires iflag_evap_prec >= 4 and niter_lscp=1'
402           CALL abort_physic (modname,abort_message,1)
403        ENDIF
404    ENDIF
405
406    !--Initialisations of constants depending on other constants
407    !--rho_snow formula from r_snow (Brandes et al. 2007 - JAMC)
408    rho_snow = 1.e3 * 0.178 * (r_snow * 2. * 1000.)**(-0.922)
409
410
411    !AA Temporary initialisation
412    a_tr_sca(1) = -0.5
413    a_tr_sca(2) = -0.5
414    a_tr_sca(3) = -0.5
415    a_tr_sca(4) = -0.5
416   
417    IF (ok_ice_sursat) CALL ice_sursat_init()
418
419    CALL cloudth_ini(iflag_cloudth_vert,iflag_ratqs)
420
421RETURN
422
423END SUBROUTINE lscp_ini
424
425END MODULE lmdz_lscp_ini
Note: See TracBrowser for help on using the repository browser.