Skip to content
Snippets Groups Projects
Commit 9737946a authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

Philippe + Benoit Vié 11/06/2020: improve removal of negative scalar variables... 

Philippe + Benoit Vié 11/06/2020: improve removal of negative scalar variables + adapt the corresponding budgets
parent aba69365
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
src/MNH/advection_metsv.f90
+ 33
− 18
View file @ 9737946a
......@@ -139,6 +139,7 @@ END MODULE MODI_ADVECTION_METSV
!! 03/2020 (B.Vie) : LIMA negativity checks after turbulence, advection and
!! microphysics budgets
! P. Wautelet 11/06/2020: bugfix: correct PRSVS array indices
! P. Wautelet + Benoit Vié 11/06/2020: improve removal of negative scalar variables + adapt the corresponding budgets
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
......@@ -692,6 +693,20 @@ DO JSV=1,KSV
END DO
!
SELECT CASE ( HCLOUD )
CASE('KESS')
WHERE (PRRS(:,:,:,2) < 0.)
PRRS(:,:,:,1) = PRRS(:,:,:,1) + PRRS(:,:,:,2)
ZEXN(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
ZT(:,:,:)= PTHT(:,:,:)*ZEXN(:,:,:)
ZLV(:,:,:)=XLVTT +(XCPV-XCL) *(ZT(:,:,:)-XTT)
ZCPH(:,:,:)=XCPD +XCPV*PRT(:,:,:,1)
PRTHS(:,:,:) = PRTHS(:,:,:) - PRRS(:,:,:,2) * ZLV(:,:,:) / &
ZCPH(:,:,:) / ZEXN(:,:,:)
PRRS(:,:,:,2) = 0.0
END WHERE
!
PRRS(:,:,:,:) = MAX( 0.0,PRRS(:,:,:,:) )
!
CASE('ICE3','ICE4')
ZEXN(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
ZT(:,:,:)= PTHT(:,:,:)*ZEXN(:,:,:)
......@@ -731,6 +746,8 @@ SELECT CASE ( HCLOUD )
PRRS(:,:,:,4) = PRRS(:,:,:,4) -ZCOR(:,:,:)
END WHERE
END IF
!
PRRS(:,:,:,:) = MAX( 0.0,PRRS(:,:,:,:) )
!
CASE('C2R2','KHKO')
ZEXN(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
......@@ -750,8 +767,10 @@ SELECT CASE ( HCLOUD )
PRSVS(:,:,:,NSV_C2R2BEG-1+JSV) = 0.0
END WHERE
END DO
!
CASE('LIMA')
!
PRRS(:,:,:,:) = MAX( 0.0,PRRS(:,:,:,:) )
!
CASE('LIMA')
ZEXN(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
ZT(:,:,:)= PTHT(:,:,:)*ZEXN(:,:,:)
ZLV(:,:,:)=XLVTT +(XCPV-XCL) *(ZT(:,:,:)-XTT)
......@@ -793,7 +812,8 @@ SELECT CASE ( HCLOUD )
!
END SELECT
!
IF ((HCLOUD == 'ICE3') .OR. (HCLOUD == 'ICE4') .OR. (HCLOUD == 'KHKO') .OR. (HCLOUD == 'C2R2') .OR. (HCLOUD == 'LIMA') ) THEN
IF ( HCLOUD == 'KESS' .OR. HCLOUD == 'ICE3' .OR. HCLOUD == 'ICE4' .OR. &
HCLOUD == 'KHKO' .OR. HCLOUD == 'C2R2' .OR. HCLOUD == 'LIMA' ) THEN
IF (LBUDGET_TH) CALL BUDGET (PRTHS(:,:,:) , 4,'NEADV_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRRS(:,:,:,1), 6,'NEADV_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2), 7,'NEADV_BU_RRC')
......@@ -803,22 +823,17 @@ IF ((HCLOUD == 'ICE3') .OR. (HCLOUD == 'ICE4') .OR. (HCLOUD == 'KHKO') .OR. (HCL
IF (LBUDGET_RG) CALL BUDGET (PRRS(:,:,:,6),11,'NEADV_BU_RRG')
IF (LBUDGET_RH) CALL BUDGET (PRRS(:,:,:,7),12,'NEADV_BU_RRH')
END IF
IF (LBUDGET_SV .AND. (HCLOUD == 'LIMA')) THEN
IF (LWARM) CALL BUDGET (PRSVS(:,:,:,NSV_LIMA_NC),12+NSV_LIMA_NC,'NEADV_BU_RSV')
IF (LWARM.AND.LRAIN) CALL BUDGET (PRSVS(:,:,:,NSV_LIMA_NR),12+NSV_LIMA_NR,'NEADV_BU_RSV')
IF (LCOLD) CALL BUDGET (PRSVS(:,:,:,NSV_LIMA_NI),12+NSV_LIMA_NI,'NEADV_BU_RSV')
IF (NMOD_CCN.GE.1) THEN
DO JI=1, NMOD_CCN
CALL BUDGET ( PRSVS(:,:,:,NSV_LIMA_CCN_FREE+JI-1),12+NSV_LIMA_CCN_FREE+JI-1,'NEADV_BU_RSV')
END DO
END IF
IF (NMOD_IFN.GE.1) THEN
DO JI=1, NMOD_IFN
CALL BUDGET ( PRSVS(:,:,:,NSV_LIMA_IFN_FREE+JI-1),12+NSV_LIMA_IFN_FREE+JI-1,'NEADV_BU_RSV')
END DO
END IF
END IF
IF (LBUDGET_SV .AND. (HCLOUD == 'C2R2' .OR. HCLOUD == 'KHKO')) THEN
DO JI = 1, 3
CALL BUDGET ( PRSVS(:,:,:,NSV_C2R2BEG-1+JI),12+NSV_C2R2BEG-1+JI,'NEADV_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. HCLOUD == 'LIMA') THEN
DO JI = NSV_LIMA_BEG, NSV_LIMA_END
CALL BUDGET ( PRSVS(:,:,:,JI),12+JI,'NEADV_BU_RSV')
END DO
END IF
!-------------------------------------------------------------------------------
!
END SUBROUTINE ADVECTION_METSV
src/MNH/ini_budget.f90
+ 202
− 91
View file @ 9737946a
This diff is collapsed.
src/MNH/resolved_cloud.f90
+ 71
− 88
View file @ 9737946a
......@@ -272,6 +272,7 @@ END MODULE MODI_RESOLVED_CLOUD
!! B.Vié 03/03/2020 : use DTHRAD instead of dT/dt in Smax diagnostic computation
! P. Wautelet 11/06/2020: bugfix: correct ZSVS array indices
! P. Wautelet 11/06/2020: bugfix: add "Non local correction for precipitating species" for ICE4
! P. Wautelet + Benoit Vié 11/06/2020: improve removal of negative scalar variables + adapt the corresponding budgets
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
......@@ -635,25 +636,24 @@ SELECT CASE ( HCLOUD )
PRS(:,:,:,2) = 0.0
END WHERE
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
! CASE('C2R2','KHKO')
! CALL GET_HALO(PRS(:,:,:,2))
! CALL GET_HALO(ZSVS(:,:,:,2))
! WHERE (PRS(:,:,:,2) < 0. .OR. ZSVS(:,:,:,2) < 0.)
! ZSVS(:,:,:,1) = 0.0
! END WHERE
! DO JSV = 2, 3
! WHERE (PRS(:,:,:,JSV) < 0. .OR. ZSVS(:,:,:,JSV) < 0.)
! PRS(:,:,:,1) = PRS(:,:,:,1) + PRS(:,:,:,JSV)
! PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,JSV) * ZLV(:,:,:) / &
! ZCPH(:,:,:) / ZEXN(:,:,:)
! PRS(:,:,:,JSV) = 0.0
! ZSVS(:,:,:,JSV) = 0.0
! END WHERE
! ENDDO
! Commented 03/2013 O.Thouron
! (at least necessary to be commented for supersaturation variable)
! ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
!
CASE('C2R2','KHKO')
WHERE (PRS(:,:,:,2) < 0. .OR. ZSVS(:,:,:,2) < 0.)
ZSVS(:,:,:,1) = 0.0
END WHERE
DO JSV = 2, 3
WHERE (PRS(:,:,:,JSV) < 0. .OR. ZSVS(:,:,:,JSV) < 0.)
PRS(:,:,:,1) = PRS(:,:,:,1) + PRS(:,:,:,JSV)
PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,JSV) * ZLV(:,:,:) / &
ZCPH(:,:,:) / ZEXN(:,:,:)
PRS(:,:,:,JSV) = 0.0
ZSVS(:,:,:,JSV) = 0.0
END WHERE
ENDDO
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
!
CASE('ICE3','ICE4')
......@@ -690,6 +690,8 @@ SELECT CASE ( HCLOUD )
PRS(:,:,:,4) = PRS(:,:,:,4) -ZCOR(:,:,:)
END WHERE
END IF
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
CASE('C3R5')
WHERE (PRS(:,:,:,2) < 0. .OR. ZSVS(:,:,:,2) < 0.)
......@@ -704,14 +706,13 @@ SELECT CASE ( HCLOUD )
ZSVS(:,:,:,JSV) = 0.0
END WHERE
ENDDO
ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
! ice
WHERE (PRS(:,:,:,4) < 0.)
PRS(:,:,:,1) = PRS(:,:,:,1) + PRS(:,:,:,4)
PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,4) * ZLV(:,:,:) / &
ZCPH(:,:,:) / PEXNREF(:,:,:)
PRS(:,:,:,4) = 0.0
PSVS(:,:,:,4) = 0.0
ZSVS(:,:,:,4) = 0.0
END WHERE
! cloud
WHERE (PRS(:,:,:,2) < 0.)
......@@ -719,11 +720,13 @@ SELECT CASE ( HCLOUD )
PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,2) * ZLV(:,:,:) / &
ZCPH(:,:,:) / PEXNREF(:,:,:)
PRS(:,:,:,2) = 0.0
PSVS(:,:,:,2) = 0.0
ZSVS(:,:,:,2) = 0.0
END WHERE
PSVS(:,:,:,:) = MAX( 0.0,PSVS(:,:,:,:) )
!
CASE('LIMA')
ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
CASE('LIMA')
! Correction where rc<0 or Nc<0
IF (OWARM) THEN
WHERE (PRS(:,:,:,2) < YRTMIN(2)/PTSTEP .OR. ZSVS(:,:,:,NSV_LIMA_NC-NSV_LIMA_BEG+1) < YCTMIN(2)/PTSTEP)
......@@ -756,7 +759,7 @@ SELECT CASE ( HCLOUD )
END IF
!
ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
END SELECT
!
......@@ -764,35 +767,31 @@ END SELECT
!* 3.3 STORE THE BUDGET TERMS
! ----------------------
!
IF ((HCLOUD /= 'KHKO') .AND. (HCLOUD /= 'C2R2') ) THEN
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:) * PRHODJ(:,:,:), 4,'NEGA_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRS(:,:,:,1) * PRHODJ(:,:,:), 6,'NEGA_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRS(:,:,:,2) * PRHODJ(:,:,:), 7,'NEGA_BU_RRC')
END IF
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:) * PRHODJ(:,:,:), 4,'NEGA_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRS(:,:,:,1) * PRHODJ(:,:,:), 6,'NEGA_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRS(:,:,:,2) * PRHODJ(:,:,:), 7,'NEGA_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRS(:,:,:,3) * PRHODJ(:,:,:), 8,'NEGA_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRS(:,:,:,4) * PRHODJ(:,:,:) ,9,'NEGA_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRS(:,:,:,5) * PRHODJ(:,:,:),10,'NEGA_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRS(:,:,:,6) * PRHODJ(:,:,:),11,'NEGA_BU_RRG')
IF (LBUDGET_RH) CALL BUDGET (PRS(:,:,:,7) * PRHODJ(:,:,:),12,'NEGA_BU_RRH')
IF (LBUDGET_SV .AND. (HCLOUD == 'LIMA')) THEN
IF (OWARM) CALL BUDGET (ZSVS(:,:,:,NSV_LIMA_NC-NSV_LIMA_BEG+1) * PRHODJ(:,:,:),12+NSV_LIMA_NC,'NEGA_BU_RSV')
IF (OWARM.AND.ORAIN) CALL BUDGET (ZSVS(:,:,:,NSV_LIMA_NR-NSV_LIMA_BEG+1) * PRHODJ(:,:,:),12+NSV_LIMA_NR,'NEGA_BU_RSV')
IF (LCOLD) CALL BUDGET (ZSVS(:,:,:,NSV_LIMA_NI-NSV_LIMA_BEG+1) * PRHODJ(:,:,:),12+NSV_LIMA_NI,'NEGA_BU_RSV')
IF (NMOD_CCN.GE.1) THEN
DO JL=1, NMOD_CCN
CALL BUDGET ( ZSVS(:,:,:,NSV_LIMA_CCN_FREE+JL-NSV_LIMA_BEG)* &
PRHODJ(:,:,:),12+NSV_LIMA_CCN_FREE+JL-1,'NEGA_BU_RSV')
END DO
END IF
IF (NMOD_IFN.GE.1) THEN
DO JL=1, NMOD_IFN
CALL BUDGET ( ZSVS(:,:,:,NSV_LIMA_IFN_FREE+JL-NSV_LIMA_BEG)* &
PRHODJ(:,:,:),12+NSV_LIMA_IFN_FREE+JL-1,'NEGA_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. (HCLOUD == 'C2R2' .OR. HCLOUD == 'KHKO')) THEN
DO JSV = 1, 3
CALL BUDGET ( ZSVS(:,:,:,JSV),12+NSV_C2R2BEG-1+JSV,'NEGA_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. HCLOUD == 'C3R5') THEN
DO JSV = 1, 4
CALL BUDGET ( ZSVS(:,:,:,JSV),12+NSV_C2R2BEG-1+JSV,'NEGA_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. HCLOUD == 'LIMA') THEN
DO JSV = NSV_LIMA_BEG, NSV_LIMA_END
CALL BUDGET ( ZSVS(:,:,:,JSV-NSV_LIMA_BEG+1),12+JSV,'NEGA_BU_RSV')
END DO
END IF
!
!* 3.4 Limitations of Na and Nc to the CCN max number concentration
!
! Commented by O.Thouron 03/2013
......@@ -1148,26 +1147,7 @@ SELECT CASE ( HCLOUD )
PRS(:,:,:,2) = 0.0
END WHERE
!
!
! CASE('C2R2','KHKO')
! CALL GET_HALO(PRS(:,:,:,2))
! CALL GET_HALO(ZSVS(:,:,:,2))
! WHERE (PRS(:,:,:,2) < 0. .OR. ZSVS(:,:,:,2) < 0.)
! ZSVS(:,:,:,1) = 0.0
! END WHERE
! DO JSV = 2, 3
! WHERE (PRS(:,:,:,JSV) < 0. .OR. ZSVS(:,:,:,JSV) < 0.)
! PRS(:,:,:,1) = PRS(:,:,:,1) + PRS(:,:,:,JSV)
! PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,JSV) * ZLV(:,:,:) / &
! ZCPH(:,:,:) / ZEXN(:,:,:)
! PRS(:,:,:,JSV) = 0.0
! ZSVS(:,:,:,JSV) = 0.0
! END WHERE
! ENDDO
! Commented 03/2013 O.Thouron
! (at least necessary to be commented for supersaturation variable)
! ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
CASE('ICE3','ICE4')
WHERE (PRS(:,:,:,4) < 0.)
......@@ -1203,6 +1183,8 @@ SELECT CASE ( HCLOUD )
PRS(:,:,:,4) = PRS(:,:,:,4) -ZCOR(:,:,:)
END WHERE
END IF
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
CASE('C2R2','KHKO')
WHERE (PRS(:,:,:,2) < 0. .OR. ZSVS(:,:,:,2) < 0.)
......@@ -1217,7 +1199,9 @@ SELECT CASE ( HCLOUD )
ZSVS(:,:,:,JSV) = 0.0
END WHERE
ENDDO
!
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
CASE('C3R5')
WHERE (PRS(:,:,:,2) < 0. .OR. ZSVS(:,:,:,2) < 0.)
ZSVS(:,:,:,1) = 0.0
......@@ -1231,14 +1215,13 @@ SELECT CASE ( HCLOUD )
ZSVS(:,:,:,JSV) = 0.0
END WHERE
ENDDO
ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
! ice
WHERE (PRS(:,:,:,4) < 0.)
PRS(:,:,:,1) = PRS(:,:,:,1) + PRS(:,:,:,4)
PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,4) * ZLV(:,:,:) / &
ZCPH(:,:,:) / PEXNREF(:,:,:)
PRS(:,:,:,4) = 0.0
PSVS(:,:,:,4) = 0.0
ZSVS(:,:,:,4) = 0.0
END WHERE
! cloud
WHERE (PRS(:,:,:,2) < 0.)
......@@ -1246,9 +1229,10 @@ SELECT CASE ( HCLOUD )
PTHS(:,:,:) = PTHS(:,:,:) - PRS(:,:,:,2) * ZLV(:,:,:) / &
ZCPH(:,:,:) / PEXNREF(:,:,:)
PRS(:,:,:,2) = 0.0
PSVS(:,:,:,2) = 0.0
ZSVS(:,:,:,2) = 0.0
END WHERE
PSVS(:,:,:,:) = MAX( 0.0,PSVS(:,:,:,:) )
!
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
!
CASE('LIMA')
! Correction where rc<0 or Nc<0
......@@ -1283,7 +1267,7 @@ SELECT CASE ( HCLOUD )
END IF
!
ZSVS(:,:,:,:) = MAX( 0.0,ZSVS(:,:,:,:) )
PRS(:,:,:,:) = MAX( 0.0,PRS(:,:,:,:) )
PRS(:,:,:,:) = MAX( 0.0,PRS (:,:,:,:) )
!
END SELECT
!
......@@ -1299,22 +1283,21 @@ IF (LBUDGET_RI) CALL BUDGET (PRS(:,:,:,4) * PRHODJ(:,:,:) ,9,'NECON_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRS(:,:,:,5) * PRHODJ(:,:,:),10,'NECON_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRS(:,:,:,6) * PRHODJ(:,:,:),11,'NECON_BU_RRG')
IF (LBUDGET_RH) CALL BUDGET (PRS(:,:,:,7) * PRHODJ(:,:,:),12,'NECON_BU_RRH')
IF (LBUDGET_SV .AND. (HCLOUD == 'LIMA')) THEN
IF (OWARM) CALL BUDGET (ZSVS(:,:,:,NSV_LIMA_NC-NSV_LIMA_BEG+1) * PRHODJ(:,:,:),12+NSV_LIMA_NC,'NECON_BU_RSV')
IF (OWARM.AND.ORAIN) CALL BUDGET (ZSVS(:,:,:,NSV_LIMA_NR-NSV_LIMA_BEG+1) * PRHODJ(:,:,:),12+NSV_LIMA_NR,'NECON_BU_RSV')
IF (LCOLD) CALL BUDGET (ZSVS(:,:,:,NSV_LIMA_NI-NSV_LIMA_BEG+1) * PRHODJ(:,:,:),12+NSV_LIMA_NI,'NECON_BU_RSV')
IF (NMOD_CCN.GE.1) THEN
DO JL=1, NMOD_CCN
CALL BUDGET ( ZSVS(:,:,:,NSV_LIMA_CCN_FREE+JL-NSV_LIMA_BEG)* &
PRHODJ(:,:,:),12+NSV_LIMA_CCN_FREE+JL-1,'NECON_BU_RSV')
END DO
END IF
IF (NMOD_IFN.GE.1) THEN
DO JL=1, NMOD_IFN
CALL BUDGET ( ZSVS(:,:,:,NSV_LIMA_IFN_FREE+JL-NSV_LIMA_BEG)* &
PRHODJ(:,:,:),12+NSV_LIMA_IFN_FREE+JL-1,'NECON_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. (HCLOUD == 'C2R2' .OR. HCLOUD == 'KHKO')) THEN
DO JSV = 1, 3
CALL BUDGET ( ZSVS(:,:,:,JSV),12+NSV_C2R2BEG-1+JSV,'NECON_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. HCLOUD == 'C3R5') THEN
DO JSV = 1, 4
CALL BUDGET ( ZSVS(:,:,:,JSV),12+NSV_C2R2BEG-1+JSV,'NECON_BU_RSV')
END DO
END IF
IF (LBUDGET_SV .AND. HCLOUD == 'LIMA') THEN
DO JSV = NSV_LIMA_BEG, NSV_LIMA_END
CALL BUDGET ( ZSVS(:,:,:,JSV-NSV_LIMA_BEG+1),12+JSV,'NECON_BU_RSV')
END DO
END IF
!-------------------------------------------------------------------------------
!
......
src/MNH/turb.f90
+ 46
− 30
View file @ 9737946a
......@@ -342,6 +342,7 @@ END MODULE MODI_TURB
!! 01/2018 (Q.Rodier) Introduction of RM17
!! 03/2020 (B.Vie) : LIMA negativity checks after turbulence, advection and microphysics budgets
! P. Wautelet 11/06/2020: bugfix: correct PRSVS array indices
! P. Wautelet + Benoit Vié 11/06/2020: improve removal of negative scalar variables + adapt the corresponding budgets
!! --------------------------------------------------------------------------
!
!* 0. DECLARATIONS
......@@ -1120,6 +1121,22 @@ IF ( KRRL >= 1 ) THEN
END IF
!
SELECT CASE ( HCLOUD )
CASE('KESS')
WHERE (PRRS(:,:,:,2) < 0.)
ZEXNE(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
PRRS(:,:,:,1) = PRRS(:,:,:,1) + PRRS(:,:,:,2)
ZTT(:,:,:)= PTHLT(:,:,:)*ZEXNE(:,:,:)
ZLV(:,:,:)=XLVTT +(XCPV-XCL) *(ZTT(:,:,:)-XTT)
ZCPH(:,:,:)=XCPD +XCPV*PRT(:,:,:,1)
PRTHLS(:,:,:) = PRTHLS(:,:,:) - PRRS(:,:,:,2) * ZLV(:,:,:) / &
ZCPH(:,:,:) / ZEXNE(:,:,:)
PRRS(:,:,:,2) = 0.0
END WHERE
!
IF (LBUDGET_TH) CALL BUDGET (PRTHLS(:,:,:), 4,'NETUR_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRRS(:,:,:,1), 6,'NETUR_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2), 7,'NETUR_BU_RRC')
CASE('ICE3','ICE4')
ZEXNE(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
ZTT(:,:,:)= PTHLT(:,:,:)*ZEXNE(:,:,:)
......@@ -1159,6 +1176,11 @@ SELECT CASE ( HCLOUD )
PRRS(:,:,:,4) = PRRS(:,:,:,4) -ZCOR(:,:,:)
END WHERE
END IF
!
IF (LBUDGET_TH) CALL BUDGET (PRTHLS(:,:,:), 4,'NETUR_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRRS(:,:,:,1), 6,'NETUR_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2), 7,'NETUR_BU_RRC')
IF (LBUDGET_RI) CALL BUDGET (PRRS(:,:,:,4), 9,'NETUR_BU_RRI')
!
CASE('C2R2','KHKO')
ZEXNE(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
......@@ -1178,8 +1200,18 @@ SELECT CASE ( HCLOUD )
PRSVS(:,:,:,NSV_C2R2BEG-1+JSV) = 0.0
END WHERE
END DO
!
CASE('LIMA')
!
IF (LBUDGET_TH) CALL BUDGET (PRTHLS(:,:,:), 4,'NETUR_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRRS(:,:,:,1), 6,'NETUR_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2), 7,'NETUR_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:,3), 8,'NETUR_BU_RRR')
IF (LBUDGET_SV) THEN
DO JSV = 1, 3
CALL BUDGET (PRSVS(:,:,:,NSV_C2R2BEG-1+JSV),12+NSV_C2R2BEG-1+JSV,'NETUR_BU_RSV')
END DO
END IF
!
CASE('LIMA')
ZEXNE(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD)
ZTT(:,:,:)= PTHLT(:,:,:)*ZEXNE(:,:,:)
ZLV(:,:,:)=XLVTT +(XCPV-XCL) *(ZTT(:,:,:)-XTT)
......@@ -1217,36 +1249,20 @@ SELECT CASE ( HCLOUD )
END IF
!
PRSVS(:, :, :, NSV_LIMA_BEG:NSV_LIMA_END) = MAX( 0.0, PRSVS(:, :, :, NSV_LIMA_BEG:NSV_LIMA_END) )
PRRS(:,:,:,:) = MAX( 0.0,PRRS(:,:,:,:) )
!Commented out for B. Vie 15/05/2020 PRRS(:,:,:,:) = MAX( 0.0,PRRS(:,:,:,:) )
!
IF (LBUDGET_TH) CALL BUDGET (PRTHLS(:,:,:), 4,'NETUR_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRRS(:,:,:,1), 6,'NETUR_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2), 7,'NETUR_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:,3), 8,'NETUR_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRRS(:,:,:,4), 9,'NETUR_BU_RRI')
IF (LBUDGET_SV) THEN
DO JSV = NSV_LIMA_BEG, NSV_LIMA_END
CALL BUDGET (PRSVS(:,:,:,JSV),12+JSV,'NETUR_BU_RSV')
END DO
END IF
!
END SELECT
!
IF ((HCLOUD == 'ICE3') .OR. (HCLOUD == 'ICE4') .OR. (HCLOUD == 'KHKO') .OR. (HCLOUD == 'C2R2') .OR. (HCLOUD == 'LIMA') ) THEN
IF (LBUDGET_TH) CALL BUDGET (PRTHLS(:,:,:), 4,'NETUR_BU_RTH')
IF (LBUDGET_RV) CALL BUDGET (PRRS(:,:,:,1), 6,'NETUR_BU_RRV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2), 7,'NETUR_BU_RRC')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:,3), 8,'NETUR_BU_RRR')
IF (LBUDGET_RI) CALL BUDGET (PRRS(:,:,:,4), 9,'NETUR_BU_RRI')
IF (LBUDGET_RS) CALL BUDGET (PRRS(:,:,:,5),10,'NETUR_BU_RRS')
IF (LBUDGET_RG) CALL BUDGET (PRRS(:,:,:,6),11,'NETUR_BU_RRG')
IF (LBUDGET_RH) CALL BUDGET (PRRS(:,:,:,7),12,'NETUR_BU_RRH')
END IF
IF (LBUDGET_SV .AND. (HCLOUD == 'LIMA')) THEN
IF (LWARM) CALL BUDGET (PRSVS(:,:,:,NSV_LIMA_NC),12+NSV_LIMA_NC,'NETUR_BU_RSV')
IF (LWARM.AND.LRAIN) CALL BUDGET (PRSVS(:,:,:,NSV_LIMA_NR),12+NSV_LIMA_NR,'NETUR_BU_RSV')
IF (LCOLD) CALL BUDGET (PRSVS(:,:,:,NSV_LIMA_NI),12+NSV_LIMA_NI,'NETUR_BU_RSV')
IF (NMOD_CCN.GE.1) THEN
DO JI=1, NMOD_CCN
CALL BUDGET ( PRSVS(:,:,:,NSV_LIMA_CCN_FREE+JI-1),12+NSV_LIMA_CCN_FREE+JI-1,'NETUR_BU_RSV')
END DO
END IF
IF (NMOD_IFN.GE.1) THEN
DO JI=1, NMOD_IFN
CALL BUDGET ( PRSVS(:,:,:,NSV_LIMA_IFN_FREE+JI-1),12+NSV_LIMA_IFN_FREE+JI-1,'NETUR_BU_RSV')
END DO
END IF
END IF
!
!----------------------------------------------------------------------------
!
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme