Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
SUBROUTINE ICE4_TENDENCIES(KSIZE, KIB, KIE, KIT, KJB, KJE, KJT, KKB, KKE, KKT, KKL, &
&KRR, LDSOFT, PCOMPUTE, &
&OWARM, HSUBG_RC_RR_ACCR, HSUBG_RR_EVAP, HSUBG_AUCV_RC, HSUBG_PR_PDF, &
&PEXN, PRHODREF, PLVFACT, PLSFACT, LDMICRO, K1, K2, K3, &
&PPRES, PCF, PSIGMA_RC, &
&PCIT, &
&PT, PTHT, &
&PRVT, PRCT, PRRT, PRIT, PRST, PRGT, PRHT, PRRT3D, &
&PRVHENI_MR, PRRHONG_MR, PRIMLTC_MR, PRSRIMCG_MR, &
&PRCHONI, PRVDEPS, PRIAGGS, PRIAUTS, PRVDEPG, &
&PRCAUTR, PRCACCR, PRREVAV, &
&PRCRIMSS, PRCRIMSG, PRSRIMCG, PRRACCSS, PRRACCSG, PRSACCRG, PRSMLTG, PRCMLTSR, &
&PRICFRRG, PRRCFRIG, PRICFRR, PRCWETG, PRIWETG, PRRWETG, PRSWETG, &
&PRCDRYG, PRIDRYG, PRRDRYG, PRSDRYG, PRWETGH, PRWETGH_MR, PRGMLTR, &
&PRCWETH, PRIWETH, PRSWETH, PRGWETH, PRRWETH, &
&PRCDRYH, PRIDRYH, PRSDRYH, PRRDRYH, PRGDRYH, PRDRYHG, PRHMLTR, &
&PRCBERI, &
&PRS_TEND, PRG_TEND, PRH_TEND, PSSI, &
&PA_TH, PA_RV, PA_RC, PA_RR, PA_RI, PA_RS, PA_RG, PA_RH, &
&PB_TH, PB_RV, PB_RC, PB_RR, PB_RI, PB_RS, PB_RG, PB_RH, &
&PHLC_HCF, PHLC_LCF, PHLC_HRC, PHLC_LRC, PRAINFR)
!!
!!** PURPOSE
!! -------
!! Computes the tendencies
!!
!! AUTHOR
!! ------
!! S. Riette from the splitting of rain_ice source code (nov. 2014)
!!
!! MODIFICATIONS
!! -------------
!!
!
!
!* 0. DECLARATIONS
! ------------
!
USE MODD_CST
USE MODD_RAIN_ICE_PARAM
USE MODD_RAIN_ICE_DESCR
USE MODD_PARAM_ICE, ONLY : CSNOWRIMING
USE MODI_ICE4_NUCLEATION
USE MODE_ICE4_RRHONG, ONLY: ICE4_RRHONG
USE MODE_ICE4_RIMLTC, ONLY: ICE4_RIMLTC
USE MODE_ICE4_RSRIMCG_OLD, ONLY: ICE4_RSRIMCG_OLD
USE MODE_ICE4_COMPUTE_PDF, ONLY: ICE4_COMPUTE_PDF
USE MODE_ICE4_RAINFR_VERT, ONLY: ICE4_RAINFR_VERT
USE MODE_ICE4_SLOW, ONLY: ICE4_SLOW
USE MODE_ICE4_WARM, ONLY: ICE4_WARM
USE MODE_ICE4_FAST_RS, ONLY: ICE4_FAST_RS
USE MODE_ICE4_FAST_RG, ONLY: ICE4_FAST_RG
USE MODE_ICE4_FAST_RH, ONLY: ICE4_FAST_RH
USE MODE_ICE4_FAST_RI, ONLY: ICE4_FAST_RI
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
USE PARKIND1, ONLY : JPRB
USE YOMHOOK , ONLY : LHOOK, DR_HOOK
!
IMPLICIT NONE
!
!* 0.1 Declarations of dummy arguments :
!
INTEGER, INTENT(IN) :: KSIZE, KIB, KIE, KIT, KJB, KJE, KJT, KKB, KKE, KKT, KKL
INTEGER, INTENT(IN) :: KRR
LOGICAL, INTENT(IN) :: LDSOFT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PCOMPUTE
LOGICAL, INTENT(IN) :: OWARM
CHARACTER*80, INTENT(IN) :: HSUBG_RC_RR_ACCR
CHARACTER*80, INTENT(IN) :: HSUBG_RR_EVAP
CHARACTER(len=4), INTENT(IN) :: HSUBG_AUCV_RC
CHARACTER*80, INTENT(IN) :: HSUBG_PR_PDF ! pdf for subgrid precipitation
REAL, DIMENSION(KSIZE), INTENT(IN) :: PEXN
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRHODREF
REAL, DIMENSION(KSIZE), INTENT(IN) :: PLVFACT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PLSFACT
LOGICAL, DIMENSION(KIT,KJT,KKT), INTENT(IN) :: LDMICRO
INTEGER, DIMENSION(KSIZE), INTENT(IN) :: K1
INTEGER, DIMENSION(KSIZE), INTENT(IN) :: K2
INTEGER, DIMENSION(KSIZE), INTENT(IN) :: K3
REAL, DIMENSION(KSIZE), INTENT(IN) :: PPRES
REAL, DIMENSION(KSIZE), INTENT(IN) :: PCF
REAL, DIMENSION(KSIZE), INTENT(IN) :: PSIGMA_RC
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PCIT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PTHT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRVT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRCT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRRT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRIT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRST
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRGT
REAL, DIMENSION(KSIZE), INTENT(IN) :: PRHT
REAL, DIMENSION(KIT,KJT,KKT), INTENT(IN) :: PRRT3D
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRVHENI_MR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRHONG_MR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIMLTC_MR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSRIMCG_MR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCHONI
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRVDEPS
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIAGGS
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIAUTS
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRVDEPG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCAUTR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCACCR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRREVAV
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCRIMSS
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCRIMSG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSRIMCG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRACCSS
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRACCSG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSACCRG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSMLTG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCMLTSR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRICFRRG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRCFRIG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRICFRR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCWETG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIWETG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRWETG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSWETG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCDRYG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIDRYG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRDRYG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSDRYG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRWETGH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRWETGH_MR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRGMLTR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCWETH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIWETH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSWETH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRGWETH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRWETH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCDRYH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRIDRYH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRSDRYH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRRDRYH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRGDRYH
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRDRYHG
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRHMLTR
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PRCBERI
REAL, DIMENSION(KSIZE, 8), INTENT(INOUT) :: PRS_TEND
REAL, DIMENSION(KSIZE, 8), INTENT(INOUT) :: PRG_TEND
REAL, DIMENSION(KSIZE, 10), INTENT(INOUT) :: PRH_TEND
REAL, DIMENSION(KSIZE), INTENT(INOUT) :: PSSI
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_TH
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RV
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RC
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RR
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RI
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RS
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RG
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PA_RH
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_TH
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RV
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RC
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RR
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RI
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RS
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RG
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PB_RH
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PHLC_HCF
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PHLC_LCF
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PHLC_HRC
REAL, DIMENSION(KSIZE), INTENT(OUT) :: PHLC_LRC
REAL, DIMENSION(KIT,KJT,KKT), INTENT(OUT) :: PRAINFR ! Rain fraction
!
!* 0.2 declaration of local variables
!
REAL, DIMENSION(KSIZE) :: ZRVT, ZRCT, ZRRT, ZRIT, ZRST, ZRGT, &
& ZT, ZTHT, &
& ZZW, &
& ZKA, ZDV, ZAI, ZCJ, &
& ZRF, &
& ZLBDAR, ZLBDAS, ZLBDAG, ZLBDAH, ZLBDAR_RF, &
& ZRGSI, ZRGSI_MR
REAL, DIMENSION(KSIZE) :: PHLI_HCF, PHLI_LCF, PHLI_HRI, PHLI_LRI
CHARACTER(len=80) :: HSUBG_AUCV_RI='NONE'
REAL, DIMENSION(KIT,KJT,KKT) :: ZRRT3D, ZRST3D, ZRGT3D, ZRHT3D
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
INTEGER :: JL
REAL, DIMENSION(KSIZE) :: ZWETG ! 1. if graupel growths in wet mode, 0. otherwise
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (LHOOK) CALL DR_HOOK('ICE4_TENDENCIES', 0, ZHOOK_HANDLE)
PA_TH(:)=0.
PA_RV(:)=0.
PA_RC(:)=0.
PA_RR(:)=0.
PA_RI(:)=0.
PA_RS(:)=0.
PA_RG(:)=0.
PA_RH(:)=0.
PB_TH(:)=0.
PB_RV(:)=0.
PB_RC(:)=0.
PB_RR(:)=0.
PB_RI(:)=0.
PB_RS(:)=0.
PB_RG(:)=0.
PB_RH(:)=0.
!
DO JL=1, KSIZE
ZRVT(JL)=PRVT(JL)
ZRCT(JL)=PRCT(JL)
ZRRT(JL)=PRRT(JL)
ZRIT(JL)=PRIT(JL)
ZRST(JL)=PRST(JL)
ZRGT(JL)=PRGT(JL)
ZTHT(JL)=PTHT(JL)
ZT(JL)=PT(JL)
ENDDO
IF(LDSOFT) THEN
PRVHENI_MR(:)=0.
PRRHONG_MR(:)=0.
PRIMLTC_MR(:)=0.
PRSRIMCG_MR(:)=0.
ELSE
!
!* 2. COMPUTES THE SLOW COLD PROCESS SOURCES
! --------------------------------------
CALL ICE4_NUCLEATION(KSIZE, LDSOFT, PCOMPUTE==1., &
ZTHT, PPRES, PRHODREF, PEXN, PLSFACT, ZT, &
ZRVT, &
PCIT, PRVHENI_MR, PB_TH, PB_RV, PB_RI)
DO JL=1, KSIZE
ZRIT(JL)=ZRIT(JL) + PRVHENI_MR(JL)
ZRVT(JL)=ZRVT(JL) - PRVHENI_MR(JL)
ZTHT(JL)=ZTHT(JL) + PRVHENI_MR(JL)*PLSFACT(JL)
ZT(JL) = ZTHT(JL) * PEXN(JL)
ENDDO
!
!* 3.3 compute the spontaneous freezing source: RRHONG
!
CALL ICE4_RRHONG(KSIZE, LDSOFT, PCOMPUTE, &
&PEXN, PLVFACT, PLSFACT, &
&ZT, ZRRT, &
&ZTHT, &
&PRRHONG_MR, PB_TH, PB_RR, PB_RG)
DO JL=1, KSIZE
PB_RG(JL) = PB_RG(JL) + PRRHONG_MR(JL)
PB_RR(JL) = PB_RR(JL) - PRRHONG_MR(JL)
PB_TH(JL) = PB_TH(JL) + PRRHONG_MR(JL)*(PLSFACT(JL)-PLVFACT(JL))
ENDDO
DO JL=1, KSIZE
ZRGT(JL) = ZRGT(JL) + PRRHONG_MR(JL)
ZRRT(JL) = ZRRT(JL) - PRRHONG_MR(JL)
ZTHT(JL) = ZTHT(JL) + PRRHONG_MR(JL)*(PLSFACT(JL)-PLVFACT(JL)) ! f(L_f*(RRHONG))
ZT(JL) = ZTHT(JL) * PEXN(JL)
ENDDO
!
!* 7.1 cloud ice melting
!
CALL ICE4_RIMLTC(KSIZE, LDSOFT, PCOMPUTE, &
&PEXN, PLVFACT, PLSFACT, &
&ZT, &
&ZTHT, ZRIT, &
&PRIMLTC_MR)
DO JL=1, KSIZE
PB_RC(JL) = PB_RC(JL) + PRIMLTC_MR(JL)
PB_RI(JL) = PB_RI(JL) - PRIMLTC_MR(JL)
PB_TH(JL) = PB_TH(JL) - PRIMLTC_MR(JL)*(PLSFACT(JL)-PLVFACT(JL))
ENDDO
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
DO JL=1, KSIZE
ZRCT(JL) = ZRCT(JL) + PRIMLTC_MR(JL)
ZRIT(JL) = ZRIT(JL) - PRIMLTC_MR(JL)
ZTHT(JL) = ZTHT(JL) - PRIMLTC_MR(JL)*(PLSFACT(JL)-PLVFACT(JL)) ! f(L_f*(-RIMLTC))
ZT(JL) = ZTHT(JL) * PEXN(JL)
ENDDO
!
! 5.1.6 riming-conversion of the large sized aggregates into graupel (old parametrisation)
!
IF(CSNOWRIMING=='OLD ') THEN
ZLBDAS(:)=0.
WHERE(ZRST(:)>0.)
ZLBDAS(:) = MIN(XLBDAS_MAX, XLBS*(PRHODREF(:)*MAX(ZRST(:), XRTMIN(5)))**XLBEXS)
END WHERE
CALL ICE4_RSRIMCG_OLD(KSIZE, LDSOFT, PCOMPUTE==1., &
&PRHODREF, &
&ZLBDAS, &
&ZT, ZRCT, ZRST, &
&PRSRIMCG_MR, PB_RS, PB_RG)
DO JL=1, KSIZE
ZRST(JL) = ZRST(JL) - PRSRIMCG_MR(JL)
ZRGT(JL) = ZRGT(JL) + PRSRIMCG_MR(JL)
ENDDO
ELSE
PRSRIMCG_MR(:) = 0.
ENDIF
ENDIF
!
!* Derived fields
!
IF(KSIZE>0) THEN
IF(.NOT. LDSOFT) THEN
ZZW(:) = EXP(XALPI-XBETAI/ZT(:)-XGAMI*ALOG(ZT(:)))
DO JL=1, KSIZE
PSSI(JL) = ZRVT(JL)*( PPRES(JL)-ZZW(JL) ) / ( XEPSILO * ZZW(JL) ) - 1.0
! Supersaturation over ice
ZKA(JL) = 2.38E-2 + 0.0071E-2*(ZT(JL)-XTT) ! k_a
ZDV(JL) = 0.211E-4*(ZT(JL)/XTT)**1.94 * (XP00/PPRES(JL)) ! D_v
ZAI(JL) = (XLSTT+(XCPV-XCI)*(ZT(JL)-XTT))**2 / (ZKA(JL)*XRV*ZT(JL)**2) &
+ ( XRV*ZT(JL) ) / (ZDV(JL)*ZZW(JL))
ZCJ(JL) = XSCFAC*PRHODREF(JL)**0.3 / SQRT(1.718E-5+0.0049E-5*(ZT(JL)-XTT))
ENDDO
ENDIF
!
!Cloud water split between high and low content part is done here
CALL ICE4_COMPUTE_PDF(KSIZE, HSUBG_AUCV_RC, HSUBG_AUCV_RI, HSUBG_PR_PDF,&
PRHODREF, ZRCT, ZRIT, PCF, ZT, PSIGMA_RC,&
PHLC_HCF, PHLC_LCF, PHLC_HRC, PHLC_LRC,&
PHLI_HCF, PHLI_LCF, PHLI_HRI, PHLI_LRI, ZRF)
IF(HSUBG_RC_RR_ACCR=='PRFR') THEN
!Diagnostic of precipitation fraction
PRAINFR(:,:,:)=UNPACK(ZRF(:), MASK=LDMICRO(:,:,:), FIELD=PRAINFR(:,:,:))
ZRRT3D(:,:,:)=PRRT3D(:,:,:)-UNPACK(PRRHONG_MR(:), MASK=LDMICRO(:,:,:), FIELD=0.)
ZRST3D(:,:,:)=0. !need phasing
ZRGT3D(:,:,:)=0. !need phasing
IF (KRR==7) ZRHT3D(:,:,:)=0. !need phasing
CALL ICE4_RAINFR_VERT(KIB, KIE, KIT, KJB, KJE, KJT, KKB, KKE, KKT, KKL, PRAINFR(:,:,:),
&ZRRT3D(:,:,:), ZRST3D(:,:,:), ZRGT3D(:,:,:), ZRHT3D(:,:,:))
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
DO JL=1,KSIZE
ZRF(JL)=PRAINFR(K1(JL), K2(JL), K3(JL))
END DO
ELSE
PRAINFR(:,:,:)=1.
ZRF(:)=1.
ENDIF
!
!* compute the slope parameters
!
ZLBDAS(:)=0.
WHERE(ZRST(:)>0.)
ZLBDAS(:) = MIN(XLBDAS_MAX, XLBS*(PRHODREF(:)*MAX(ZRST(:), XRTMIN(5)))**XLBEXS)
END WHERE
ZLBDAG(:)=0.
WHERE(ZRGT(:)>0.)
ZLBDAG(:) = XLBG*(PRHODREF(:)*MAX(ZRGT(:), XRTMIN(6)))**XLBEXG
END WHERE
!ZLBDAR will be used when we consider rain diluted over the grid box
ZLBDAR(:)=0.
WHERE(ZRRT(:)>0.)
ZLBDAR(:) = XLBR*( PRHODREF(:)*MAX( ZRRT(:), XRTMIN(3)))**XLBEXR
END WHERE
!ZLBDAR_RF is used when we consider rain concentrated in its fraction
IF (HSUBG_RC_RR_ACCR=='PRFR') THEN
ZLBDAR_RF(:)=0.
WHERE(ZRRT(:)>0. .AND. ZRF(:)>0.)
ZLBDAR_RF(:) = XLBR*( PRHODREF(:) *MAX( ZRRT(:)/ZRF(:) , XRTMIN(3)))**XLBEXR
END WHERE
ELSE
ZLBDAR_RF(:) = ZLBDAR(:)
ENDIF
IF(KRR==7) THEN
ZLBDAH(:)=0.
WHERE(PRHT(:)>0.)
ZLBDAH(:) = XLBH*(PRHODREF(:)*MAX(PRHT(:), XRTMIN(7)))**XLBEXH
END WHERE
ENDIF
ENDIF
!
!
CALL ICE4_SLOW(KSIZE, LDSOFT, PCOMPUTE, PRHODREF, ZT, &
&PSSI, PLVFACT, PLSFACT, &
&ZRVT, ZRCT, ZRIT, ZRST, ZRGT, &
&ZLBDAS, ZLBDAG, &
&ZAI, ZCJ, PHLI_HCF, ZRIT, & !!!!!!! en attendant phasage plus complet
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
&PRCHONI, PRVDEPS, PRIAGGS, PRIAUTS, PRVDEPG, &
&PA_TH, PA_RV, PA_RC, PA_RI, PA_RS, PA_RG)
!
!-------------------------------------------------------------------------------
!
!
!* 3. COMPUTES THE SLOW WARM PROCESS SOURCES
! --------------------------------------
!
!
IF(OWARM) THEN ! Check if the formation of the raindrops by the slow
! warm processes is allowed
CALL ICE4_WARM(KSIZE, LDSOFT, PCOMPUTE, HSUBG_RC_RR_ACCR, HSUBG_RR_EVAP, &
&PRHODREF, PLVFACT, ZT, PPRES, ZTHT,&
&ZLBDAR, ZLBDAR_RF, ZKA, ZDV, ZCJ, &
&PHLC_LCF, PHLC_HCF, PHLC_LRC, PHLC_HRC, &
&PCF, ZRF, &
&ZRVT, ZRCT, ZRRT, &
&PRCAUTR, PRCACCR, PRREVAV, &
&PA_TH, PA_RV, PA_RC, PA_RR)
ELSE
PRCAUTR(:)=0.
PRCACCR(:)=0.
PRREVAV(:)=0.
END IF
!
!-------------------------------------------------------------------------------
!
!
!* 4. COMPUTES THE FAST COLD PROCESS SOURCES FOR r_s
! ----------------------------------------------
!
CALL ICE4_FAST_RS(KSIZE, KSIZE, LDSOFT, PCOMPUTE, &
&PRHODREF, PLVFACT, PLSFACT, PPRES, &
&ZDV, ZKA, ZCJ, &
&ZLBDAR, ZLBDAS, &
&ZT, ZRVT, ZRCT, ZRRT, ZRST, &
&PRIAGGS, &
&PRCRIMSS, PRCRIMSG, PRSRIMCG, &
&PRRACCSS, PRRACCSG, PRSACCRG, PRSMLTG, &
&PRCMLTSR, &
&PRS_TEND, &
&PA_TH, PA_RC, PA_RR, PA_RS, PA_RG)
!
!-------------------------------------------------------------------------------
!
!
!* 5. COMPUTES THE FAST COLD PROCESS SOURCES FOR r_g
! ------------------------------------------------------
!
DO JL=1, KSIZE
ZRGSI(JL) = PRVDEPG(JL) + PRSMLTG(JL) + PRRACCSG(JL) + &
& PRSACCRG(JL) + PRCRIMSG(JL) + PRSRIMCG(JL)
ZRGSI_MR(JL) = PRRHONG_MR(JL) + PRSRIMCG_MR(JL)
ENDDO
CALL ICE4_FAST_RG(KSIZE, KSIZE, LDSOFT, PCOMPUTE, KRR, &
&PRHODREF, PLVFACT, PLSFACT, PPRES, &
&ZDV, ZKA, ZCJ, PCIT, &
&ZLBDAR, ZLBDAS, ZLBDAG, &
&ZT, ZRVT, ZRCT, ZRRT, ZRIT, ZRST, ZRGT, &
&ZRGSI, ZRGSI_MR(:), &
&ZWETG, &
&PRICFRRG, PRRCFRIG, PRICFRR, PRCWETG, PRIWETG, PRRWETG, PRSWETG, &
&PRCDRYG, PRIDRYG, PRRDRYG, PRSDRYG, PRWETGH, PRWETGH_MR, PRGMLTR, &
&PRG_TEND, &
&PA_TH, PA_RC, PA_RR, PA_RI, PA_RS, PA_RG, PA_RH, PB_RG, PB_RH)
!
!-------------------------------------------------------------------------------
!
!
!* 6. COMPUTES THE FAST COLD PROCESS SOURCES FOR r_h
! ----------------------------------------------
!
IF (KRR==7) THEN
CALL ICE4_FAST_RH(KSIZE, KSIZE, LDSOFT, PCOMPUTE, ZWETG, &
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
&PRHODREF, PLVFACT, PLSFACT, PPRES, &
&ZDV, ZKA, ZCJ, &
&ZLBDAS, ZLBDAG, ZLBDAR, ZLBDAH, &
&ZT, ZRVT, ZRCT, ZRRT, ZRIT, ZRST, ZRGT, PRHT, &
&PRCWETH, PRIWETH, PRSWETH, PRGWETH, PRRWETH, &
&PRCDRYH, PRIDRYH, PRSDRYH, PRRDRYH, PRGDRYH, PRDRYHG, PRHMLTR, &
&PRH_TEND, &
&PA_TH, PA_RC, PA_RR, PA_RI, PA_RS, PA_RG, PA_RH)
ELSE
PRCWETH(:)=0.
PRIWETH(:)=0.
PRSWETH(:)=0.
PRGWETH(:)=0.
PRRWETH(:)=0.
PRCDRYH(:)=0.
PRIDRYH(:)=0.
PRSDRYH(:)=0.
PRRDRYH(:)=0.
PRGDRYH(:)=0.
PRDRYHG(:)=0.
PRHMLTR(:)=0.
END IF
!
!-------------------------------------------------------------------------------
!
!
!* 7. COMPUTES SPECIFIC SOURCES OF THE WARM AND COLD CLOUDY SPECIES
! -------------------------------------------------------------
!
CALL ICE4_FAST_RI(KSIZE, LDSOFT, PCOMPUTE, &
&PRHODREF, PLVFACT, PLSFACT, &
&ZAI, ZCJ, PCIT, &
&PSSI, &
&ZRCT, ZRIT, &
&PRCBERI, PA_TH, PA_RC, PA_RI)
!
IF (LHOOK) CALL DR_HOOK('ICE4_TENDENCIES', 1, ZHOOK_HANDLE)
!
END SUBROUTINE ICE4_TENDENCIES