Newer
Older
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
IF (OWARM) THEN ! rain_ice_warm
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'AUTO_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'AUTO_BU_RRR')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'ACCR_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'ACCR_BU_RRR')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'REVA_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRVS(:,:,:)*PRHODJ(:,:,:),6,'REVA_BU_RRV')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'REVA_BU_RRR')
ENDIF
!rain_ice_fast_rs
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'RIM_BU_RTH')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'RIM_BU_RRC')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'RIM_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'RIM_BU_RRG')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'ACC_BU_RTH')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'ACC_BU_RRR')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'ACC_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'ACC_BU_RRG')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'CMEL_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'CMEL_BU_RRG')
!rain_ice_fast_rg
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'CFRZ_BU_RTH')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'CFRZ_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9,'CFRZ_BU_RRI')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'CFRZ_BU_RRG')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'WETG_BU_RTH')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'WETG_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'WETG_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9,'WETG_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'WETG_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'WETG_BU_RRG')
IF (LBUDGET_RH) CALL BUDGET (PRHS(:,:,:)*PRHODJ(:,:,:),12,'WETG_BU_RRH')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'DRYG_BU_RTH')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'DRYG_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'DRYG_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9,'DRYG_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'DRYG_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'DRYG_BU_RRG')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'GMLT_BU_RTH')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'GMLT_BU_RRR')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'GMLT_BU_RRG')
IF(KRR==7) THEN ! rain_ice_fast_rh
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'WETH_BU_RTH')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'WETH_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'WETH_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9,'WETH_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'WETH_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'WETH_BU_RRG')
IF (LBUDGET_RH) CALL BUDGET (PRHS(:,:,:)*PRHODJ(:,:,:),12,'WETH_BU_RRH')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'HMLT_BU_RTH')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8,'HMLT_BU_RRR')
IF (LBUDGET_RH) CALL BUDGET (PRHS(:,:,:)*PRHODJ(:,:,:),12,'HMLT_BU_RRH')
ENDIF
!rain_ice_fast_ri
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'IMLT_BU_RTH')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'IMLT_BU_RRC')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9,'IMLT_BU_RRI')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4,'BERFI_BU_RTH')
IF (LBUDGET_RC) CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7,'BERFI_BU_RRC')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9,'BERFI_BU_RRI')
!
END IF
!
!-------------------------------------------------------------------------------
!
!* 8. COMPUTE THE SEDIMENTATION (RS) SOURCE
! -------------------------------------
!
!* 8.1 time splitting loop initialization
!
ZTSPLITR= PTSTEP / FLOAT(KSPLITR)
!
!
IF (HSEDIM == 'STAT') THEN
CALL RAIN_ICE_SEDIMENTATION_STAT
ELSEIF (HSEDIM == 'SPLI') THEN
CALL RAIN_ICE_SEDIMENTATION_SPLIT
ELSE
WRITE(*,*) ' STOP'
WRITE(*,*) ' NO SEDIMENTATION SCHEME FOR HSEDIM=',HSEDIM
CALL PRINT_MSG(NVERB_FATAL,'GEN','RAIN_ICE','')
!sedimentation of rain fraction
CALL RAINFR_VERT(PRAINFR, PRRS(:,:,:)*PTSTEP)
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
!
!
!-------------------------------------------------------------------------------
!
!-------------------------------------------------------------------------------
!
!
CONTAINS
!
!
!-------------------------------------------------------------------------------
!
!
SUBROUTINE RAIN_ICE_SEDIMENTATION_SPLIT
!
!* 0. DECLARATIONS
! ------------
!
IMPLICIT NONE
!
!* 0.2 declaration of local variables
!
!
INTEGER , DIMENSION(SIZE(GSEDIMC)) :: IC1,IC2,IC3 ! Used to replace the COUNT
INTEGER , DIMENSION(SIZE(GSEDIMR)) :: IR1,IR2,IR3 ! Used to replace the COUNT
INTEGER , DIMENSION(SIZE(GSEDIMS)) :: IS1,IS2,IS3 ! Used to replace the COUNT
INTEGER , DIMENSION(SIZE(GSEDIMI)) :: II1,II2,II3 ! Used to replace the COUNT
INTEGER , DIMENSION(SIZE(GSEDIMG)) :: IG1,IG2,IG3 ! Used to replace the COUNT
INTEGER , DIMENSION(SIZE(GSEDIMH)) :: IH1,IH2,IH3 ! Used to replace the COUNT
INTEGER :: ILENALLOCC,ILENALLOCR,ILENALLOCI,ILENALLOCS,ILENALLOCG,ILENALLOCH
INTEGER :: ILISTLENC,ILISTLENR,ILISTLENI,ILISTLENS,ILISTLENG,ILISTLENH
INTEGER, ALLOCATABLE :: ILISTR(:),ILISTC(:),ILISTI(:),ILISTS(:),ILISTG(:),ILISTH(:)
! Optimization for NEC
!INTEGER, SAVE :: IOLDALLOCC = SIZE(PEXNREF,1)*SIZE(PEXNREF,2)*SIZE(PEXNREF,3)/10
!INTEGER, SAVE :: IOLDALLOCR = SIZE(PEXNREF,1)*SIZE(PEXNREF,2)*SIZE(PEXNREF,3)/10
!INTEGER, SAVE :: IOLDALLOCI = SIZE(PEXNREF,1)*SIZE(PEXNREF,2)*SIZE(PEXNREF,3)/10
!INTEGER, SAVE :: IOLDALLOCS = SIZE(PEXNREF,1)*SIZE(PEXNREF,2)*SIZE(PEXNREF,3)/10
!INTEGER, SAVE :: IOLDALLOCG = SIZE(PEXNREF,1)*SIZE(PEXNREF,2)*SIZE(PEXNREF,3)/10
!INTEGER, SAVE :: IOLDALLOCH = SIZE(PEXNREF,1)*SIZE(PEXNREF,2)*SIZE(PEXNREF,3)/10
INTEGER, SAVE :: IOLDALLOCC = 6000
INTEGER, SAVE :: IOLDALLOCR = 6000
INTEGER, SAVE :: IOLDALLOCI = 6000
INTEGER, SAVE :: IOLDALLOCS = 6000
INTEGER, SAVE :: IOLDALLOCG = 6000
INTEGER, SAVE :: IOLDALLOCH = 6000
!
REAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) :: ZCONC3D ! droplet condensation
!-------------------------------------------------------------------------------
!
!
! O. Initialization of for sedimentation
!
IF (OSEDIC) PINPRC (:,:) = 0.
PINPRR (:,:) = 0.
PINPRR3D (:,:,:) = 0.
PINPRS (:,:) = 0.
PINPRG (:,:) = 0.
IF ( KRR == 7 ) PINPRH (:,:) = 0.
IF (PRESENT(PFPR)) PFPR(:,:,:,:) = 0.
!
!* 1. Parameters for cloud sedimentation
!
IF (OSEDIC) THEN
ZRAY(:,:,:) = 0.
ZLBC(:,:,:) = XLBC(1)
ZFSEDC(:,:,:) = XFSEDC(1)
ZCONC3D(:,:,:)= XCONC_LAND
ZCONC_TMP(:,:)= XCONC_LAND
IF (PRESENT(PSEA)) THEN
ZCONC_TMP(:,:)=PSEA(:,:)*XCONC_SEA+(1.-PSEA(:,:))*XCONC_LAND
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
DO JK=IKTB,IKTE
ZLBC(:,:,JK) = PSEA(:,:)*XLBC(2)+(1.-PSEA(:,:))*XLBC(1)
ZFSEDC(:,:,JK) = (PSEA(:,:)*XFSEDC(2)+(1.-PSEA(:,:))*XFSEDC(1))
ZFSEDC(:,:,JK) = MAX(MIN(XFSEDC(1),XFSEDC(2)),ZFSEDC(:,:,JK))
ZCONC3D(:,:,JK)= (1.-PTOWN(:,:))*ZCONC_TMP(:,:)+PTOWN(:,:)*XCONC_URBAN
ZRAY(:,:,JK) = 0.5*((1.-PSEA(:,:))*GAMMA(XNUC+1.0/XALPHAC)/(GAMMA(XNUC)) + &
PSEA(:,:)*GAMMA(XNUC2+1.0/XALPHAC2)/(GAMMA(XNUC2)))
END DO
ELSE
ZCONC3D(:,:,:) = XCONC_LAND
ZRAY(:,:,:) = 0.5*(GAMMA(XNUC+1.0/XALPHAC)/(GAMMA(XNUC)))
END IF
ZRAY(:,:,:) = MAX(1.,ZRAY(:,:,:))
ZLBC(:,:,:) = MAX(MIN(XLBC(1),XLBC(2)),ZLBC(:,:,:))
ENDIF
!
!* 2. compute the fluxes
!
! optimization by looking for locations where
! the precipitating fields are larger than a minimal value only !!!
! For optimization we consider each variable separately
ZRTMIN(:) = XRTMIN(:) * ZINVTSTEP
IF (OSEDIC) GSEDIMC(:,:,:) = .FALSE.
GSEDIMR(:,:,:) = .FALSE.
GSEDIMI(:,:,:) = .FALSE.
GSEDIMS(:,:,:) = .FALSE.
GSEDIMG(:,:,:) = .FALSE.
IF ( KRR == 7 ) GSEDIMH(:,:,:) = .FALSE.
!
ILENALLOCR = 0
IF (OSEDIC) ILENALLOCC = 0
ILENALLOCI = 0
ILENALLOCS = 0
ILENALLOCG = 0
IF ( KRR == 7 ) ILENALLOCH = 0
!
! ZPiS = Specie i source creating during the current time step
! PRiS = Source of the previous time step
!
IF (OSEDIC) THEN
ZPRCS(:,:,:) = 0.0
ZPRCS(:,:,:) = PRCS(:,:,:)-PRCT(:,:,:)* ZINVTSTEP
PRCS(:,:,:) = PRCT(:,:,:)* ZINVTSTEP
END IF
ZPRRS(:,:,:) = 0.0
ZPRSS(:,:,:) = 0.0
ZPRGS(:,:,:) = 0.0
IF ( KRR == 7 ) ZPRHS(:,:,:) = 0.0
!
ZPRRS(:,:,:) = PRRS(:,:,:)-PRRT(:,:,:)* ZINVTSTEP
ZPRSS(:,:,:) = PRSS(:,:,:)-PRST(:,:,:)* ZINVTSTEP
ZPRGS(:,:,:) = PRGS(:,:,:)-PRGT(:,:,:)* ZINVTSTEP
IF ( KRR == 7 ) ZPRHS(:,:,:) = PRHS(:,:,:)-PRHT(:,:,:)* ZINVTSTEP
PRRS(:,:,:) = PRRT(:,:,:)* ZINVTSTEP
PRSS(:,:,:) = PRST(:,:,:)* ZINVTSTEP
PRGS(:,:,:) = PRGT(:,:,:)* ZINVTSTEP
IF ( KRR == 7 ) PRHS(:,:,:) = PRHT(:,:,:)* ZINVTSTEP
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
1316
1317
!
! PRiS = Source of the previous time step + source created during the subtime
! step
!
DO JN = 1 , KSPLITR
IF( JN==1 ) THEN
IF (OSEDIC) PRCS(:,:,:) = PRCS(:,:,:) + ZPRCS(:,:,:)/KSPLITR
PRRS(:,:,:) = PRRS(:,:,:) + ZPRRS(:,:,:)/KSPLITR
PRSS(:,:,:) = PRSS(:,:,:) + ZPRSS(:,:,:)/KSPLITR
PRGS(:,:,:) = PRGS(:,:,:) + ZPRGS(:,:,:)/KSPLITR
IF ( KRR == 7 ) PRHS(:,:,:) = PRHS(:,:,:) + ZPRHS(:,:,:)/KSPLITR
DO JK = IKTB , IKTE
ZW(:,:,JK) =ZTSPLITR/(PRHODREF(:,:,JK)* PDZZ(:,:,JK))
END DO
ELSE
IF (OSEDIC) PRCS(:,:,:) = PRCS(:,:,:) + ZPRCS(:,:,:)*ZTSPLITR
PRRS(:,:,:) = PRRS(:,:,:) + ZPRRS(:,:,:)*ZTSPLITR
PRSS(:,:,:) = PRSS(:,:,:) + ZPRSS(:,:,:)*ZTSPLITR
PRGS(:,:,:) = PRGS(:,:,:) + ZPRGS(:,:,:)*ZTSPLITR
IF ( KRR == 7 ) PRHS(:,:,:) = PRHS(:,:,:) + ZPRHS(:,:,:)*ZTSPLITR
END IF
!
IF (OSEDIC) GSEDIMC(IIB:IIE,IJB:IJE,IKTB:IKTE) = &
PRCS(IIB:IIE,IJB:IJE,IKTB:IKTE)>ZRTMIN(2)
GSEDIMR(IIB:IIE,IJB:IJE,IKTB:IKTE) = &
PRRS(IIB:IIE,IJB:IJE,IKTB:IKTE)>ZRTMIN(3)
GSEDIMI(IIB:IIE,IJB:IJE,IKTB:IKTE) = &
PRIS(IIB:IIE,IJB:IJE,IKTB:IKTE)>ZRTMIN(4)
GSEDIMS(IIB:IIE,IJB:IJE,IKTB:IKTE) = &
PRSS(IIB:IIE,IJB:IJE,IKTB:IKTE)>ZRTMIN(5)
GSEDIMG(IIB:IIE,IJB:IJE,IKTB:IKTE) = &
PRGS(IIB:IIE,IJB:IJE,IKTB:IKTE)>ZRTMIN(6)
IF ( KRR == 7 ) GSEDIMH(IIB:IIE,IJB:IJE,IKTB:IKTE) = &
PRHS(IIB:IIE,IJB:IJE,IKTB:IKTE)>ZRTMIN(7)
!
IF (OSEDIC) ISEDIMC = COUNTJV( GSEDIMC(:,:,:),IC1(:),IC2(:),IC3(:))
ISEDIMR = COUNTJV( GSEDIMR(:,:,:),IR1(:),IR2(:),IR3(:))
ISEDIMI = COUNTJV( GSEDIMI(:,:,:),II1(:),II2(:),II3(:))
ISEDIMS = COUNTJV( GSEDIMS(:,:,:),IS1(:),IS2(:),IS3(:))
ISEDIMG = COUNTJV( GSEDIMG(:,:,:),IG1(:),IG2(:),IG3(:))
IF ( KRR == 7 ) ISEDIMH = COUNTJV( GSEDIMH(:,:,:),IH1(:),IH2(:),IH3(:))
!
!* 2.1 for cloud
!
IF (OSEDIC) THEN
ZWSED(:,:,:) = 0.
IF( JN==1 ) PRCS(:,:,:) = PRCS(:,:,:) * PTSTEP
IF( ISEDIMC >= 1 ) THEN
IF ( ISEDIMC .GT. ILENALLOCC ) THEN
IF ( ILENALLOCC .GT. 0 ) THEN
DEALLOCATE (ZRCS, ZRHODREFC, ILISTC,ZWLBDC,ZCONC,ZRCT, &
ZZT,ZPRES,ZRAY1D,ZFSEDC1D,ZWLBDA,ZCC )
END IF
ILENALLOCC = MAX (IOLDALLOCC, 2*ISEDIMC )
IOLDALLOCC = ILENALLOCC
ALLOCATE(ZRCS(ILENALLOCC), ZRHODREFC(ILENALLOCC), ILISTC(ILENALLOCC), &
ZWLBDC(ILENALLOCC), ZCONC(ILENALLOCC), ZRCT(ILENALLOCC), ZZT(ILENALLOCC), &
ZPRES(ILENALLOCC), ZRAY1D(ILENALLOCC), ZFSEDC1D(ILENALLOCC), &
ZWLBDA(ILENALLOCC), ZCC(ILENALLOCC) )
END IF
!
DO JL=1,ISEDIMC
ZRCS(JL) = PRCS(IC1(JL),IC2(JL),IC3(JL))
ZRHODREFC(JL) = PRHODREF(IC1(JL),IC2(JL),IC3(JL))
ZWLBDC(JL) = ZLBC(IC1(JL),IC2(JL),IC3(JL))
ZCONC(JL) = ZCONC3D(IC1(JL),IC2(JL),IC3(JL))
ZRCT(JL) = PRCT(IC1(JL),IC2(JL),IC3(JL))
ZZT(JL) = PTHT(IC1(JL),IC2(JL),IC3(JL))
ZPRES(JL) = PPABST(IC1(JL),IC2(JL),IC3(JL))
ZRAY1D(JL) = ZRAY(IC1(JL),IC2(JL),IC3(JL))
ZFSEDC1D(JL) = ZFSEDC(IC1(JL),IC2(JL),IC3(JL))
END DO
!
ILISTLENC = 0
DO JL=1,ISEDIMC
IF( ZRCS(JL) .GT. ZRTMIN(2) ) THEN
ILISTLENC = ILISTLENC + 1
ILISTC(ILISTLENC) = JL
END IF
END DO
DO JJ = 1, ILISTLENC
JL = ILISTC(JJ)
IF (ZRCS(JL) .GT. ZRTMIN(2) .AND. ZRCT(JL) .GT. XRTMIN(2)) THEN
ZWLBDC(JL) = ZWLBDC(JL) * ZCONC(JL) / (ZRHODREFC(JL) * ZRCT(JL))
ZWLBDC(JL) = ZWLBDC(JL)**XLBEXC
ZRAY1D(JL) = ZRAY1D(JL) / ZWLBDC(JL) !! ZRAY : mean diameter=M(1)/2
ZZT(JL) = ZZT(JL) * (ZPRES(JL)/XP00)**(XRD/XCPD)
ZWLBDA(JL) = 6.6E-8*(101325./ZPRES(JL))*(ZZT(JL)/293.15)
ZCC(JL) = XCC*(1.+1.26*ZWLBDA(JL)/ZRAY1D(JL)) !! XCC modified for cloud
ZWSED (IC1(JL),IC2(JL),IC3(JL))= ZRHODREFC(JL)**(-XCEXVT +1 ) * &
ZWLBDC(JL)**(-XDC)*ZCC(JL)*ZFSEDC1D(JL) * ZRCS(JL)
END IF
END DO
END IF
DO JK = IKTB , IKTE
PRCS(:,:,JK) = PRCS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,2)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRC(:,:) = PINPRC(:,:) + ZWSED(:,:,IKB) / XRHOLW / KSPLITR
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
IF( JN==KSPLITR ) THEN
PRCS(:,:,:) = PRCS(:,:,:) * ZINVTSTEP
END IF
END IF
!
!* 2.2 for rain
!
IF( JN==1 ) PRRS(:,:,:) = PRRS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
IF( ISEDIMR >= 1 ) THEN
IF ( ISEDIMR .GT. ILENALLOCR ) THEN
IF ( ILENALLOCR .GT. 0 ) THEN
DEALLOCATE (ZRRS, ZRHODREFR, ILISTR)
END IF
ILENALLOCR = MAX (IOLDALLOCR, 2*ISEDIMR )
IOLDALLOCR = ILENALLOCR
ALLOCATE(ZRRS(ILENALLOCR), ZRHODREFR(ILENALLOCR), ILISTR(ILENALLOCR))
END IF
!
DO JL=1,ISEDIMR
ZRRS(JL) = PRRS(IR1(JL),IR2(JL),IR3(JL))
ZRHODREFR(JL) = PRHODREF(IR1(JL),IR2(JL),IR3(JL))
END DO
!
ILISTLENR = 0
DO JL=1,ISEDIMR
IF( ZRRS(JL) .GT. ZRTMIN(3) ) THEN
ILISTLENR = ILISTLENR + 1
ILISTR(ILISTLENR) = JL
END IF
END DO
DO JJ = 1, ILISTLENR
JL = ILISTR(JJ)
ZWSED (IR1(JL),IR2(JL),IR3(JL))= XFSEDR * ZRRS(JL)**XEXSEDR * &
ZRHODREFR(JL)**(XEXSEDR-XCEXVT)
END DO
END IF
DO JK = IKTB , IKTE
PRRS(:,:,JK) = PRRS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,3)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRR(:,:) = PINPRR(:,:) + ZWSED(:,:,IKB)/XRHOLW/KSPLITR
PINPRR3D(:,:,:) = PINPRR3D(:,:,:) + ZWSED(:,:,1:IKT)/XRHOLW/KSPLITR
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
1411
IF( JN==KSPLITR ) THEN
PRRS(:,:,:) = PRRS(:,:,:) * ZINVTSTEP
END IF
!
!* 2.3 for pristine ice
!
IF( JN==1 ) PRIS(:,:,:) = PRIS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
IF( ISEDIMI >= 1 ) THEN
IF ( ISEDIMI .GT. ILENALLOCI ) THEN
IF ( ILENALLOCI .GT. 0 ) THEN
DEALLOCATE (ZRIS, ZRHODREFI, ILISTI)
END IF
ILENALLOCI = MAX (IOLDALLOCI, 2*ISEDIMI )
IOLDALLOCI = ILENALLOCI
ALLOCATE(ZRIS(ILENALLOCI), ZRHODREFI(ILENALLOCI), ILISTI(ILENALLOCI))
END IF
!
DO JL=1,ISEDIMI
ZRIS(JL) = PRIS(II1(JL),II2(JL),II3(JL))
ZRHODREFI(JL) = PRHODREF(II1(JL),II2(JL),II3(JL))
END DO
!
ILISTLENI = 0
DO JL=1,ISEDIMI
IF( ZRIS(JL) .GT. MAX(ZRTMIN(4),1.0E-7 )) THEN ! limitation of the McF&H formula
ILISTLENI = ILISTLENI + 1
ILISTI(ILISTLENI) = JL
END IF
END DO
DO JJ = 1, ILISTLENI
JL = ILISTI(JJ)
ZWSED (II1(JL),II2(JL),II3(JL))= XFSEDI * ZRIS(JL) * &
ZRHODREFI(JL)**(1.0-XCEXVT) * & ! McF&H
MAX( 0.05E6,-0.15319E6-0.021454E6* &
ALOG(ZRHODREFI(JL)*ZRIS(JL)) )**XEXCSEDI
END DO
END IF
DO JK = IKTB , IKTE
PRIS(:,:,JK) = PRIS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,4)=ZWSED(:,:,JK)
ENDDO
ENDIF
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
1451
1452
1453
1454
1455
IF( JN==KSPLITR ) THEN
PRIS(:,:,:) = PRIS(:,:,:) * ZINVTSTEP
END IF
!
!* 2.4 for aggregates/snow
!
IF( JN==1 ) PRSS(:,:,:) = PRSS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
IF( ISEDIMS >= 1 ) THEN
IF ( ISEDIMS .GT. ILENALLOCS ) THEN
IF ( ILENALLOCS .GT. 0 ) THEN
DEALLOCATE (ZRSS, ZRHODREFS, ILISTS)
END IF
ILENALLOCS = MAX (IOLDALLOCS, 2*ISEDIMS )
IOLDALLOCS = ILENALLOCS
ALLOCATE(ZRSS(ILENALLOCS), ZRHODREFS(ILENALLOCS), ILISTS(ILENALLOCS))
END IF
!
DO JL=1,ISEDIMS
ZRSS(JL) = PRSS(IS1(JL),IS2(JL),IS3(JL))
ZRHODREFS(JL) = PRHODREF(IS1(JL),IS2(JL),IS3(JL))
END DO
!
ILISTLENS = 0
DO JL=1,ISEDIMS
IF( ZRSS(JL) .GT. ZRTMIN(5) ) THEN
ILISTLENS = ILISTLENS + 1
ILISTS(ILISTLENS) = JL
END IF
END DO
DO JJ = 1, ILISTLENS
JL = ILISTS(JJ)
ZWSED (IS1(JL),IS2(JL),IS3(JL))= XFSEDS * ZRSS(JL)**XEXSEDS * &
ZRHODREFS(JL)**(XEXSEDS-XCEXVT)
END DO
END IF
DO JK = IKTB , IKTE
PRSS(:,:,JK) = PRSS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,5)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRS(:,:) = PINPRS(:,:) + ZWSED(:,:,IKB)/XRHOLW/KSPLITR
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
IF( JN==KSPLITR ) THEN
PRSS(:,:,:) = PRSS(:,:,:) * ZINVTSTEP
END IF
!
!* 2.5 for graupeln
!
ZWSED(:,:,:) = 0.
IF( JN==1 ) PRGS(:,:,:) = PRGS(:,:,:) * PTSTEP
IF( ISEDIMG >= 1 ) THEN
IF ( ISEDIMG .GT. ILENALLOCG ) THEN
IF ( ILENALLOCG .GT. 0 ) THEN
DEALLOCATE (ZRGS, ZRHODREFG, ILISTG)
END IF
ILENALLOCG = MAX (IOLDALLOCG, 2*ISEDIMG )
IOLDALLOCG = ILENALLOCG
ALLOCATE(ZRGS(ILENALLOCG), ZRHODREFG(ILENALLOCG), ILISTG(ILENALLOCG))
END IF
!
DO JL=1,ISEDIMG
ZRGS(JL) = PRGS(IG1(JL),IG2(JL),IG3(JL))
ZRHODREFG(JL) = PRHODREF(IG1(JL),IG2(JL),IG3(JL))
END DO
!
ILISTLENG = 0
DO JL=1,ISEDIMG
IF( ZRGS(JL) .GT. ZRTMIN(6) ) THEN
ILISTLENG = ILISTLENG + 1
ILISTG(ILISTLENG) = JL
END IF
END DO
DO JJ = 1, ILISTLENG
JL = ILISTG(JJ)
ZWSED (IG1(JL),IG2(JL),IG3(JL))= XFSEDG * ZRGS(JL)**XEXSEDG * &
ZRHODREFG(JL)**(XEXSEDG-XCEXVT)
END DO
END IF
DO JK = IKTB , IKTE
PRGS(:,:,JK) = PRGS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,6)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRG(:,:) = PINPRG(:,:) + ZWSED(:,:,IKB)/XRHOLW/KSPLITR
IF( JN==KSPLITR ) THEN
1508
1509
1510
1511
1512
1513
1514
1515
1516
1517
1518
1519
1520
1521
1522
1523
1524
1525
1526
1527
1528
1529
1530
1531
1532
1533
1534
1535
1536
1537
1538
1539
1540
1541
1542
1543
1544
1545
1546
PRGS(:,:,:) = PRGS(:,:,:) * ZINVTSTEP
END IF
!
!* 2.6 for hail
!
IF ( KRR == 7 ) THEN
IF( JN==1 ) PRHS(:,:,:) = PRHS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
IF( ISEDIMH >= 1 ) THEN
IF ( ISEDIMH .GT. ILENALLOCH ) THEN
IF ( ILENALLOCH .GT. 0 ) THEN
DEALLOCATE (ZRHS, ZRHODREFH, ILISTH)
END IF
ILENALLOCH = MAX (IOLDALLOCH, 2*ISEDIMH )
IOLDALLOCH = ILENALLOCH
ALLOCATE(ZRHS(ILENALLOCH), ZRHODREFH(ILENALLOCH), ILISTH(ILENALLOCH))
END IF
!
DO JL=1,ISEDIMH
ZRHS(JL) = PRHS(IH1(JL),IH2(JL),IH3(JL))
ZRHODREFH(JL) = PRHODREF(IH1(JL),IH2(JL),IH3(JL))
END DO
!
ILISTLENH = 0
DO JL=1,ISEDIMH
IF( ZRHS(JL) .GT. ZRTMIN(7) ) THEN
ILISTLENH = ILISTLENH + 1
ILISTH(ILISTLENH) = JL
END IF
END DO
DO JJ = 1, ILISTLENH
JL = ILISTH(JJ)
ZWSED (IH1(JL),IH2(JL),IH3(JL))= XFSEDH * ZRHS(JL)**XEXSEDH * &
ZRHODREFH(JL)**(XEXSEDH-XCEXVT)
END DO
END IF
DO JK = IKTB , IKTE
PRHS(:,:,JK) = PRHS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,7)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRH(:,:) = PINPRH(:,:) + ZWSED(:,:,IKB)/XRHOLW/KSPLITR
IF( JN==KSPLITR ) THEN
PRHS(:,:,:) = PRHS(:,:,:) * ZINVTSTEP
END IF
END IF
!
END DO
!
IF (OSEDIC) THEN
IF (ILENALLOCC .GT. 0) DEALLOCATE (ZRCS, ZRHODREFC, &
ILISTC,ZWLBDC,ZCONC,ZRCT, ZZT,ZPRES,ZRAY1D,ZFSEDC1D, ZWLBDA,ZCC)
END IF
IF (ILENALLOCR .GT. 0 ) DEALLOCATE(ZRHODREFR,ZRRS,ILISTR)
IF (ILENALLOCI .GT. 0 ) DEALLOCATE(ZRHODREFI,ZRIS,ILISTI)
IF (ILENALLOCS .GT. 0 ) DEALLOCATE(ZRHODREFS,ZRSS,ILISTS)
IF (ILENALLOCG .GT. 0 ) DEALLOCATE(ZRHODREFG,ZRGS,ILISTG)
IF (KRR == 7 .AND. (ILENALLOCH .GT. 0 )) DEALLOCATE(ZRHODREFH,ZRHS,ILISTH)
!
!* 2.3 budget storage
!
IF (LBUDGET_RC .AND. OSEDIC) &
CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7 ,'SEDI_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8 ,'SEDI_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRIS(:,:,:)*PRHODJ(:,:,:),9 ,'SEDI_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRSS(:,:,:)*PRHODJ(:,:,:),10,'SEDI_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRGS(:,:,:)*PRHODJ(:,:,:),11,'SEDI_BU_RRG')
IF ( KRR == 7 .AND. LBUDGET_RH) &
CALL BUDGET (PRHS(:,:,:)*PRHODJ(:,:,:),12,'SEDI_BU_RRH')
!* 2.4 DROPLET DEPOSITION AT THE 1ST LEVEL ABOVE GROUND
!
IF (LDEPOSC) THEN
GDEP(:,:) = .FALSE.
GDEP(IIB:IIE,IJB:IJE) = PRCS(IIB:IIE,IJB:IJE,IKB) >0
PRCS(:,:,IKB) = PRCS(:,:,IKB) - XVDEPOSC * PRCT(:,:,IKB) / PDZZ(:,:,IKB)
PINPRC(:,:) = PINPRC(:,:) + XVDEPOSC * PRCT(:,:,IKB) * PRHODREF(:,:,IKB) /XRHOLW
PINDEP(:,:) = XVDEPOSC * PRCT(:,:,IKB) * PRHODREF(:,:,IKB) /XRHOLW
END WHERE
END IF
!
!* 2.5 budget storage
!
IF ( LBUDGET_RC .AND. LDEPOSC ) &
CALL BUDGET (PRCS(:,:,:)*PRHODJ(:,:,:),7 ,'DEPO_BU_RRC')
END SUBROUTINE RAIN_ICE_SEDIMENTATION_SPLIT
!
!-------------------------------------------------------------------------------
!
SUBROUTINE RAIN_ICE_SEDIMENTATION_STAT
!
!* 0. DECLARATIONS
! ------------
!
IMPLICIT NONE
!
!* 0.2 declaration of local variables
!
!
REAL :: ZP1,ZP2,ZH,ZZWLBDA,ZZWLBDC,ZZCC
REAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2)) :: ZQP
INTEGER :: JI,JJ,JK
INTEGER :: JCOUNT, JL
INTEGER, DIMENSION(SIZE(PRHODREF,1)*SIZE(PRHODREF,2)) :: I1, I2
REAL, DIMENSION(SIZE(PRHODREF,1),SIZE(PRHODREF,2),SIZE(PRHODREF,3)) :: ZCONC3D ! droplet condensation
!-------------------------------------------------------------------------------
!
!
!* 1. Parameters for cloud sedimentation
!
IF (OSEDIC) THEN
ZRAY(:,:,:) = 0.
ZLBC(:,:,:) = XLBC(1)
ZFSEDC(:,:,:) = XFSEDC(1)
ZCONC3D(:,:,:)= XCONC_LAND
ZCONC_TMP(:,:)= XCONC_LAND
IF (PRESENT(PSEA)) THEN
ZCONC_TMP(:,:)=PSEA(:,:)*XCONC_SEA+(1.-PSEA(:,:))*XCONC_LAND
DO JK=IKTB,IKTE
ZLBC(:,:,JK) = PSEA(:,:)*XLBC(2)+(1.-PSEA(:,:))*XLBC(1)
ZFSEDC(:,:,JK) = (PSEA(:,:)*XFSEDC(2)+(1.-PSEA(:,:))*XFSEDC(1))
ZFSEDC(:,:,JK) = MAX(MIN(XFSEDC(1),XFSEDC(2)),ZFSEDC(:,:,JK))
ZCONC3D(:,:,JK)= (1.-PTOWN(:,:))*ZCONC_TMP(:,:)+PTOWN(:,:)*XCONC_URBAN
ZRAY(:,:,JK) = 0.5*((1.-PSEA(:,:))*GAMMA(XNUC+1.0/XALPHAC)/(GAMMA(XNUC)) + &
PSEA(:,:)*GAMMA(XNUC2+1.0/XALPHAC2)/(GAMMA(XNUC2)))
END DO
ELSE
ZCONC3D(:,:,:) = XCONC_LAND
ZRAY(:,:,:) = 0.5*(GAMMA(XNUC+1.0/XALPHAC)/(GAMMA(XNUC)))
END IF
ZRAY(:,:,:) = MAX(1.,ZRAY(:,:,:))
ZLBC(:,:,:) = MAX(MIN(XLBC(1),XLBC(2)),ZLBC(:,:,:))
ENDIF
IF (LDEPOSC) PINDEP (:,:) = 0.
!
!* 2. compute the fluxes
!
ZRTMIN(:) = XRTMIN(:) * ZINVTSTEP
!
IF (OSEDIC) THEN
ZPRCS(:,:,:) = 0.0
ZPRCS(:,:,:) = PRCS(:,:,:)-PRCT(:,:,:)* ZINVTSTEP
PRCS(:,:,:) = PRCT(:,:,:)* ZINVTSTEP
END IF
ZPRRS(:,:,:) = 0.0
ZPRSS(:,:,:) = 0.0
ZPRGS(:,:,:) = 0.0
IF ( KRR == 7 ) ZPRHS(:,:,:) = 0.0
!
ZPRRS(:,:,:) = PRRS(:,:,:)-PRRT(:,:,:)* ZINVTSTEP
ZPRSS(:,:,:) = PRSS(:,:,:)-PRST(:,:,:)* ZINVTSTEP
ZPRGS(:,:,:) = PRGS(:,:,:)-PRGT(:,:,:)* ZINVTSTEP
IF ( KRR == 7 ) ZPRHS(:,:,:) = PRHS(:,:,:)-PRHT(:,:,:)* ZINVTSTEP
PRRS(:,:,:) = PRRT(:,:,:)* ZINVTSTEP
PRSS(:,:,:) = PRST(:,:,:)* ZINVTSTEP
PRGS(:,:,:) = PRGT(:,:,:)* ZINVTSTEP
IF ( KRR == 7 ) PRHS(:,:,:) = PRHT(:,:,:)* ZINVTSTEP
!
IF (OSEDIC) PRCS(:,:,:) = PRCS(:,:,:) + ZPRCS(:,:,:)
PRRS(:,:,:) = PRRS(:,:,:) + ZPRRS(:,:,:)
PRSS(:,:,:) = PRSS(:,:,:) + ZPRSS(:,:,:)
PRGS(:,:,:) = PRGS(:,:,:) + ZPRGS(:,:,:)
IF ( KRR == 7 ) PRHS(:,:,:) = PRHS(:,:,:) + ZPRHS(:,:,:)
DO JK = IKTB , IKTE
ZW(:,:,JK) =PTSTEP/(PRHODREF(:,:,JK)* PDZZ(:,:,JK) )
END DO
!
!* 2.1 for cloud
!
IF (OSEDIC) THEN
PRCS(:,:,:) = PRCS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
ZWSEDW1(:,:,:) = 0.
ZWSEDW2(:,:,:) = 0.
! calculation of P1, P2 and sedimentation flux
1699
1700
1701
1702
1703
1704
1705
1706
1707
1708
1709
1710
1711
1712
1713
1714
1715
1716
1717
1718
1719
1720
1721
1722
1723
1724
1725
1726
1727
1728
1729
1730
1731
1732
1733
1734
1735
1736
1737
1738
1739
1740
1741
1742
1743
1744
1745
DO JK = IKE , IKB, -1*KKL
!estimation of q' taking into account incomming ZWSED
ZQP(:,:)=ZWSED(:,:,JK+KKL)*ZW(:,:,JK)
JCOUNT=COUNTJV2((PRCS(:,:,JK) > ZRTMIN(2) .AND. PRCT(:,:,JK) > ZRTMIN(2)) .OR. &
(ZQP(:,:) > ZRTMIN(2)),I1(:),I2(:))
DO JL=1, JCOUNT
JI=I1(JL)
JJ=I2(JL)
!calculation of w
! mars 2009 : ajout d'un test
!IF ( PRCS(JI,JJ,JK) > ZRTMIN(2) ) THEN
IF(PRCS(JI,JJ,JK) > ZRTMIN(2) .AND. PRCT(JI,JJ,JK) > ZRTMIN(2)) THEN
ZZWLBDA=6.6E-8*(101325./PPABST(JI,JJ,JK))*(PTHT(JI,JJ,JK)/293.15)
ZZWLBDC=(ZLBC(JI,JJ,JK)*ZCONC3D(JI,JJ,JK) &
&/(PRHODREF(JI,JJ,JK)*PRCT(JI,JJ,JK)))**XLBEXC
ZZCC=XCC*(1.+1.26*ZZWLBDA*ZZWLBDC/ZRAY(JI,JJ,JK)) !! ZCC : Fall speed
ZWSEDW1 (JI,JJ,JK)=PRHODREF(JI,JJ,JK)**(-XCEXVT ) * &
& ZZWLBDC**(-XDC)*ZZCC*ZFSEDC(JI,JJ,JK)
ENDIF
IF ( ZQP(JI,JJ) > ZRTMIN(2) ) THEN
ZZWLBDA=6.6E-8*(101325./PPABST(JI,JJ,JK))*(PTHT(JI,JJ,JK)/293.15)
ZZWLBDC=(ZLBC(JI,JJ,JK)*ZCONC3D(JI,JJ,JK) &
&/(PRHODREF(JI,JJ,JK)*ZQP(JI,JJ)))**XLBEXC
ZZCC=XCC*(1.+1.26*ZZWLBDA*ZZWLBDC/ZRAY(JI,JJ,JK)) !! ZCC : Fall speed
ZWSEDW2 (JI,JJ,JK)=PRHODREF(JI,JJ,JK)**(-XCEXVT ) * &
& ZZWLBDC**(-XDC)*ZZCC*ZFSEDC(JI,JJ,JK)
ENDIF
ENDDO
DO JJ = IJB, IJE
DO JI = IIB, IIE
ZH=PDZZ(JI,JJ,JK)
ZP1 = MIN(1., ZWSEDW1(JI,JJ,JK) * PTSTEP / ZH)
! mars 2009 : correction : ZWSEDW1 => ZWSEDW2
!IF (ZWSEDW1(JI,JJ,JK) /= 0.) THEN
IF (ZWSEDW2(JI,JJ,JK) /= 0.) THEN
ZP2 = MAX(0.,1 - ZH &
& / (PTSTEP*ZWSEDW2(JI,JJ,JK)) )
ELSE
ZP2 = 0.
ENDIF
ZWSED (JI,JJ,JK)=ZP1*PRHODREF(JI,JJ,JK)*&
&ZH*PRCS(JI,JJ,JK)&
&* ZINVTSTEP+ ZP2 * ZWSED (JI,JJ,JK+KKL)
ENDDO
ENDDO
DO JK = IKTB , IKTE
PRCS(:,:,JK) = PRCS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
END DO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,2)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRC(:,:) = ZWSED(:,:,IKB)/XRHOLW ! in m/s
PRCS(:,:,:) = PRCS(:,:,:) * ZINVTSTEP
ENDIF
!
!* 2.2 for rain
!
PRRS(:,:,:) = PRRS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
ZWSEDW1(:,:,:) = 0.
ZWSEDW2(:,:,:) = 0.
! calculation of ZP1, ZP2 and sedimentation flux
1771
1772
1773
1774
1775
1776
1777
1778
1779
1780
1781
1782
1783
1784
1785
1786
1787
1788
1789
1790
1791
1792
1793
1794
1795
1796
1797
1798
1799
1800
1801
1802
1803
1804
DO JK = IKE , IKB, -1*KKL
!estimation of q' taking into account incomming ZWSED
ZQP(:,:)=ZWSED(:,:,JK+KKL)*ZW(:,:,JK)
JCOUNT=COUNTJV2((PRRS(:,:,JK) > ZRTMIN(3)) .OR. &
(ZQP(:,:) > ZRTMIN(3)),I1(:),I2(:))
DO JL=1, JCOUNT
JI=I1(JL)
JJ=I2(JL)
!calculation of w
IF ( PRRS(JI,JJ,JK) > ZRTMIN(3) ) THEN
ZWSEDW1 (JI,JJ,JK)= XFSEDR *PRRS(JI,JJ,JK)**(XEXSEDR-1)* &
PRHODREF(JI,JJ,JK)**(XEXSEDR-XCEXVT-1)
ENDIF
IF ( ZQP(JI,JJ) > ZRTMIN(3) ) THEN
ZWSEDW2 (JI,JJ,JK)= XFSEDR *(ZQP(JI,JJ))**(XEXSEDR-1)* &
PRHODREF(JI,JJ,JK)**(XEXSEDR-XCEXVT-1)
ENDIF
ENDDO
DO JJ = IJB, IJE
DO JI = IIB, IIE
ZH=PDZZ(JI,JJ,JK)
ZP1 = MIN(1., ZWSEDW1(JI,JJ,JK) * PTSTEP / ZH )
IF (ZWSEDW2(JI,JJ,JK) /= 0.) THEN
ZP2 = MAX(0.,1 - ZH &
& / (PTSTEP*ZWSEDW2(JI,JJ,JK)) )
ELSE
ZP2 = 0.
ENDIF
ZWSED (JI,JJ,JK)=ZP1*PRHODREF(JI,JJ,JK)*&
&ZH*PRRS(JI,JJ,JK)&
&* ZINVTSTEP+ ZP2 * ZWSED (JI,JJ,JK+KKL)
ENDDO
ENDDO
DO JK = IKTB , IKTE
PRRS(:,:,JK) = PRRS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
ENDDO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,3)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRR(:,:) = ZWSED(:,:,IKB)/XRHOLW ! in m/s
PINPRR3D(:,:,:) = ZWSED(:,:,1:IKT)/XRHOLW ! in m/s
PRRS(:,:,:) = PRRS(:,:,:) * ZINVTSTEP
!
!* 2.3 for pristine ice
!
PRIS(:,:,:) = PRIS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
ZWSEDW1(:,:,:) = 0.
ZWSEDW2(:,:,:) = 0.
! calculation of ZP1, ZP2 and sedimentation flux
1828
1829
1830
1831
1832
1833
1834
1835
1836
1837
1838
1839
1840
1841
1842
1843
1844
1845
1846
1847
1848
1849
1850
1851
1852
1853
1854
1855
1856
1857
1858
1859
1860
1861
1862
1863
1864
1865
DO JK = IKE , IKB, -1*KKL
!estimation of q' taking into account incomming ZWSED
ZQP(:,:)=ZWSED(:,:,JK+KKL)*ZW(:,:,JK)
JCOUNT=COUNTJV2((PRIS(:,:,JK) > MAX(ZRTMIN(4),1.0E-7 )) .OR. &
(ZQP(:,:) > MAX(ZRTMIN(4),1.0E-7 )),I1(:),I2(:))
DO JL=1, JCOUNT
JI=I1(JL)
JJ=I2(JL)
!calculation of w
IF ( PRIS(JI,JJ,JK) > MAX(ZRTMIN(4),1.0E-7 ) ) THEN
ZWSEDW1 (JI,JJ,JK)= XFSEDI * &
& PRHODREF(JI,JJ,JK)**(XCEXVT) * & ! McF&H
& MAX( 0.05E6,-0.15319E6-0.021454E6* &
& ALOG(PRHODREF(JI,JJ,JK)*PRIS(JI,JJ,JK)) )**XEXCSEDI
ENDIF
IF ( ZQP(JI,JJ) > MAX(ZRTMIN(4),1.0E-7 ) ) THEN
ZWSEDW2 (JI,JJ,JK)= XFSEDI * &
& PRHODREF(JI,JJ,JK)**(XCEXVT) * & ! McF&H
& MAX( 0.05E6,-0.15319E6-0.021454E6* &
& ALOG(PRHODREF(JI,JJ,JK)*ZQP(JI,JJ)) )**XEXCSEDI
ENDIF
ENDDO
DO JJ = IJB, IJE
DO JI = IIB, IIE
ZH=PDZZ(JI,JJ,JK)
ZP1 = MIN(1., ZWSEDW1(JI,JJ,JK) * PTSTEP / ZH )
IF (ZWSEDW2(JI,JJ,JK) /= 0.) THEN
ZP2 = MAX(0.,1 - ZH &
& / (PTSTEP*ZWSEDW2(JI,JJ,JK)) )
ELSE
ZP2 = 0.
ENDIF
ZWSED (JI,JJ,JK)=ZP1*PRHODREF(JI,JJ,JK)*&
&ZH*PRIS(JI,JJ,JK)&
&* ZINVTSTEP+ ZP2 * ZWSED (JI,JJ,JK+KKL)
ENDDO
ENDDO
DO JK = IKTB , IKTE
PRIS(:,:,JK) = PRIS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
ENDDO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,4)=ZWSED(:,:,JK)
ENDDO
ENDIF
PRIS(:,:,:) = PRIS(:,:,:) * ZINVTSTEP
!
!* 2.4 for aggregates/snow
!
PRSS(:,:,:) = PRSS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
ZWSEDW1(:,:,:) = 0.
ZWSEDW2(:,:,:) = 0.
! calculation of ZP1, ZP2 and sedimentation flux
1890
1891
1892
1893
1894
1895
1896
1897
1898
1899
1900
1901
1902
1903
1904
1905
1906
1907
1908
1909
1910
1911
1912
1913
DO JK = IKE , IKB, -1*KKL
!estimation of q' taking into account incomming ZWSED
ZQP(:,:)=ZWSED(:,:,JK+KKL)*ZW(:,:,JK)
JCOUNT=COUNTJV2((PRSS(:,:,JK) > ZRTMIN(5)) .OR. &
(ZQP(:,:) > ZRTMIN(5)),I1(:),I2(:))
DO JL=1, JCOUNT
JI=I1(JL)
JJ=I2(JL)
!calculation of w
IF (PRSS(JI,JJ,JK) > ZRTMIN(5) ) THEN
ZWSEDW1(JI,JJ,JK)=XFSEDS*(PRSS(JI,JJ,JK))**(XEXSEDS-1)*&
PRHODREF(JI,JJ,JK)**(XEXSEDS-XCEXVT-1)
ENDIF
IF ( ZQP(JI,JJ) > ZRTMIN(5) ) THEN
ZWSEDW2(JI,JJ,JK)=XFSEDS*(ZQP(JI,JJ))**(XEXSEDS-1)*&
PRHODREF(JI,JJ,JK)**(XEXSEDS-XCEXVT-1)
ENDIF
ENDDO
DO JJ = IJB, IJE
DO JI = IIB, IIE
ZH=PDZZ(JI,JJ,JK)
ZP1 = MIN(1., ZWSEDW1(JI,JJ,JK) * PTSTEP / ZH )
IF (ZWSEDW2(JI,JJ,JK) /= 0.) THEN
ZP2 = MAX(0.,1 - ZH&
/ (PTSTEP*ZWSEDW2(JI,JJ,JK)) )
ELSE
ZP2 = 0.
ENDIF
ZWSED (JI,JJ,JK)=ZP1*PRHODREF(JI,JJ,JK)*&
&ZH*PRSS(JI,JJ,JK)&
&* ZINVTSTEP+ ZP2 * ZWSED (JI,JJ,JK+KKL)
ENDDO
ENDDO
DO JK = IKTB , IKTE
PRSS(:,:,JK) = PRSS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
ENDDO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,5)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRS(:,:) = ZWSED(:,:,IKB)/XRHOLW ! in m/s
PRSS(:,:,:) = PRSS(:,:,:) * ZINVTSTEP
!
!* 2.5 for graupeln
!
PRGS(:,:,:) = PRGS(:,:,:) * PTSTEP
ZWSED(:,:,:) = 0.
ZWSEDW1(:,:,:) = 0.
ZWSEDW2(:,:,:) = 0.
! calculation of ZP1, ZP2 and sedimentation flux
1950
1951
1952
1953
1954
1955
1956
1957
1958
1959
1960
1961
1962
1963
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
DO JK = IKE, IKB, -1*KKL
!estimation of q' taking into account incomming ZWSED
ZQP(:,:)=ZWSED(:,:,JK+KKL)*ZW(:,:,JK)
JCOUNT=COUNTJV2((PRGS(:,:,JK) > ZRTMIN(6)) .OR. &
(ZQP(:,:) > ZRTMIN(6)),I1(:),I2(:))
DO JL=1, JCOUNT
JI=I1(JL)
JJ=I2(JL)
!calculation of w
IF ( PRGS(JI,JJ,JK) > ZRTMIN(6) ) THEN
ZWSEDW1 (JI,JJ,JK)= XFSEDG*(PRGS(JI,JJ,JK))**(XEXSEDG-1) * &
PRHODREF(JI,JJ,JK)**(XEXSEDG-XCEXVT-1)
ENDIF
IF ( ZQP(JI,JJ) > ZRTMIN(6) ) THEN
ZWSEDW2 (JI,JJ,JK)= XFSEDG*(ZQP(JI,JJ))**(XEXSEDG-1) * &
PRHODREF(JI,JJ,JK)**(XEXSEDG-XCEXVT-1)
ENDIF
ENDDO
DO JJ = IJB, IJE
DO JI = IIB, IIE
ZH=PDZZ(JI,JJ,JK)
ZP1 = MIN(1., ZWSEDW1(JI,JJ,JK) * PTSTEP / ZH )
IF (ZWSEDW2(JI,JJ,JK) /= 0.) THEN
ZP2 = MAX(0.,1 - ZH &
& / (PTSTEP*ZWSEDW2(JI,JJ,JK)) )
ELSE
ZP2 = 0.
ENDIF
ZWSED (JI,JJ,JK)=ZP1*PRHODREF(JI,JJ,JK)*&
&ZH*PRGS(JI,JJ,JK)&
&* ZINVTSTEP+ ZP2 * ZWSED (JI,JJ,JK+KKL)
ENDDO
ENDDO
DO JK = IKTB , IKTE
PRGS(:,:,JK) = PRGS(:,:,JK) + ZW(:,:,JK)*(ZWSED(:,:,JK+KKL)-ZWSED(:,:,JK))
ENDDO
IF (PRESENT(PFPR)) THEN
DO JK = IKTB , IKTE
PFPR(:,:,JK,6)=ZWSED(:,:,JK)
ENDDO
ENDIF
PINPRG(:,:) = ZWSED(:,:,IKB)/XRHOLW ! in m/s
PRGS(:,:,:) = PRGS(:,:,:) * ZINVTSTEP
!
!* 2.6 for hail